Skip to main content
  • Research article
  • Open access
  • Published:

Langevin dynamics simulations of charged model phosphatidylinositol lipids in the presence of diffusion barriers: toward an atomic level understanding of corralling of PIP2 by protein fences in biological membranes



The polyvalent acidic lipid phosphatidylinositol, 4,5-bisphosphate (PIP2) is important for many cellular functions. It has been suggested that different pools of PIP2 exist in the cytoplasmic leaflet of the plasma membrane, and that such pooling could play a role in the regulation of PIP2. The mechanism of fencing, however, is not understood.


This study presents the results of Langevin dynamics simulations of PIP2 to elucidate some of the molecular level considerations that must be applied to models for fencing. For each simulation, a pool of PIP2 (modeled as charged spheres) was placed in containments with boundaries modeled as a single row of rods (steric or electrostatic) or rigid protein filaments. It is shown that even a small gap (20 Å, which is 1.85 times larger than the diameter of a PIP2 sphere) leads to poor steric blocking, and that electrostatic blockage is only effective at very high charge density. Filaments of human septin, yeast septin, and actin also failed to provide adequate blockage when placed on the membrane surface. The two septins do provide high blockage consistent with experiment and with phenomenological considerations of permeability when they are buried 9 Å and 12 Å below the membrane surface, respectively. In contrast, burial does not improve blockage by the “arch-shaped” actin filaments. Free energy estimates using implicit membrane-solvent models indicate that burial of the septins to about 10 Å can be achieved without penetration of charged residues into the hydrophobic region of the membrane.


These results imply that a functioning fence assembled from protein filaments must either be buried well below the membrane surface, have more than a single row, or contain additional components that fill small gaps in the filaments.


Phosphatidylinositol 4,5-bisphosphate (PIP2) participates in numerous cellular processes such as generation of second messengers, activation of ion channels, endocytosis, and exocytosis [1],[2]. Although PIP2 occupies only ~2% of the phospholipids on the inner leaflet, there is experimental evidence that the local concentration is significantly enhanced in the regions where exo-, endo-, and phagocytosis occur [1],[2]; in the case of the forming phagosome, the local concentration is increased about 3-fold to ~6%. More generally, the pooling of PIP2 into specific compartments could allow it to participate in many different cellular functions.

Two different mechanisms of PIP2 pooling have been proposed: protein-fencing and protein-binding. In the so-called “protein-fence” hypothesis [3],[4], membrane-bound or transmembrane proteins form a boundary that greatly hinders escape of PIP2 from the corral. Protein fences have been implicated in restriction in diffusion of membrane proteins and lipids [5]. The alternative “protein-binding” hypothesis, or “reduced diffusion constant hypothesis”, is that peptides or proteins within the corral bind PIP2, which substantially reduces its diffusion constant, and thereby prevents its escape [6]. For example, PIP2 could be bound through simple electrostatic interactions to clusters of basic residues, which exist on many proteins involved in endo-, exo-, and phagocytosis, such as syntaxin and the MARCKS protein [7],[8]. Recent measurements of PIP2 diffusion in nascent phagosomes by Golebiewska et al. [4],[9] have provided unambiguous support for the protein-fence hypothesis in at least this system; diffusion within the corral is similar to other regions in the plasma membrane (and not very different from pure lipid liposomes), but diffusion out of the corral is reduced by a factor of at least 100 (1% of free diffusion). The molecular composition of the fence, however, was not determined.

What makes a good fence? An arrangement of actin filaments would seem to be a suitable candidate. Actin lies on the surface of the nascent phagosomes noted above and is highly negatively charged, as is PIP2. Hence, an unbroken row of actin filaments could in principle provide both steric and electrostatic barriers to PIP2 diffusion. However, fencing remained after actin was experimentally removed from nascent phagosomes, and computer simulations of PIP2-like charged spheres on a membrane surface indicated essentially unimpeded diffusion through an atomic-level model fence composed of actin [4].

This paper extends the aforementioned simulations to evaluate both rod-like and atomic models of protein fences. The PIP2 molecules are represented as charged spheres and simulated by Langevin dynamics (LD), which involves generating stochastic trajectories for individual particles consistent with the Langevin equation [10],[11]. LD is a close variant of Brownian dynamics [12],[13] and the results of the simulations (ratios of relaxation times) would be similar for the two methods [14]. The rod-like fences consist of a row of spheres each comparable to the size of PIP2. Steric interactions are isolated by placing the rod-like fences with different sizes of gaps (missing spheres) on the diffusion plane. Electrostatic interactions are probed by charging the spheres (either negatively or positively) and raising the row above the plane to eliminate steric interactions.

While the results of these simulations could be obtained using diffusion equations, a simulation-based solution is pedagogical and provides insight into more complex shapes and arrangements. In addition, the atomic-level particle simulations were performed with actin and two other fence candidates, human septin and yeast septin, at different levels of burial with respect to the diffusion plane. All-atom simulations suitable for studying these systems are not presently possible in terms of time and length scales. Hence, the protein fence is represented by a field on a grid and assumed to be rigid; water and all other membrane components are treated implicitly; i.e., only the PIP2-like spheres are simulated. For simplicity, only the single-row fence is considered and the possibility that peptide or other membrane components can bind to the protein filaments and fill in gaps is not explored.

By way of outline, the details of the modeling and simulation are described in the following section. The Results and Discussion section presents and analyzes the decay functions obtained from the simulations, and relates them with the experimental permeability.


Langevin dynamics of PIP2

The general simulation system is depicted in Figure 1a. PIP2 is modeled as a charged sphere (q = -4e, Lennard-Jones parameters ε LJ = 0.5 kcal/mol and r min/2 = 5.4 Å, and mass m = 1,043 amu) and confined to the diffusion plane (Z = 0, XY-plane) by a planar harmonic restraint potential with a force constant of 0.6 kcal/mol/Å2. The diffusion plane is located 5 Å above the membrane surface; this is the separation of the phosphates on the inositol ring of PIP2 and the phosphate plane of the bilayer obtained from molecular dynamics simulations [15],[16].

Figure 1
figure 1

Schematic representation of the Langevin dynamics system. (a) Schematic representation of the simulation systems. Each system is a rectangular box with hard-wall boundaries and a fence along the Y-axis centered at X = 0. PIP2 molecules are modeled as charged spheres and initially placed inside the pool (X < 0) at a concentration of 6%; the 329 PIP2 from one of the septin simulations are shown explicitly (see Figure S3 for snapshots at later times). (B-C) Schematic representations of two types of rod-like fences, with PIP2 and fence particles enlarged for clarity (b) A porous fence is generated from uncharged blocking atoms on the diffusion plane with a variation in the fence opening length (L open). (c) An electrostatic fence is generated from the single-lined charged atoms placed above the diffusion plane with a variation in the height (h) from the diffusion plane and charge per atom (q).

The systems are enclosed by hard wall boundaries modeled by a harmonic restraint potential with a force constant of 100 kcal/mol/Å2 with XY dimensions of 1,400 Å (-475 < X < 925 Å) × 1,000 Å for rod-like steric fence, 1,400 Å × 450 Å for rod-like electrostatic fence, and 1,400 Å × 750 Å for the protein fences. Fences were located parallel to the Y-axis at X = 0 (Figure 1a). The “corral” extends from -475 Å < X < 0 Å for all the systems. Simulations were initialized with 438, 197, and 329 PIP2-like spheres in the rod-like steric, rod-like electrostatic, and protein corrals (X < 0), respectively; these values correspond to the 6% concentrations of PIP2 if restricted inside the corral or 2% if distributed uniformly in the whole PIP2-accessible region. LD simulations were carried out with a position independent collision frequency γ = 1 ps-1, temperature T = 300 K, and a time step of 5 fs. The PIP2 spheres were equilibrated for 25 ns and constrained to remain in the corral, and then simulated for 10.0 μs without constraints. Different random seeds were used to generate 20 independent trajectories for the rod-like fences and 30 for the protein fences. The diffusion constant D = k B T/ = 2.4 ×10-5 cm2/s is over 10 times faster than the experimental value [4]. However, because only the ratios of hindered to free diffusion are considered, this large value of D does not alter the conclusions, and is computationally efficient. Note also that D is the same throughout the region, although the effective diffusion constant may be enhanced or retarded in the vicinity of the fence, and is larger in the beginning of the trajectory arising from repulsion of the changed particles in the corral.

Lastly, Brownian dynamics would have been equally acceptable for this study, but it is not available in CHARMM [17], the simulation package used.

Rod-like fences

As noted in the Background section, the rod-like models were simulated to quantify the general blockade characteristics for purely steric and electrostatic potentials. Steric fences (Figure 1b) of length L = 1,000 Å were developed by placing impenetrable spheres with no charge and gaps specified by the fence opening length (L open) on the diffusion plane. Simulations were performed with L open = 10, 20, 30, 40, 50, 100, 200, and 500 Å to model increasingly porous fences, and L open = 1,000 Å to model free diffusion. The electrostatic fences (Figure 1c) are single-lined charged spheres with charges (q = ±0.05e, ±0.10e, ±0.50e, and ±1.00e) and heights above the diffusion plane (h = 2, 5, 7, and 10 Å).

Protein fences

PIP2 pools in cells are several hundred nanometers on a side [4],[9] and appear to be bounded by long filamentous structures. We thereby assume that corral-building proteins have two physical properties. First, the fence proteins bind to or interact strongly with the anionic lipid membrane and/or PIP2, and second, the fence is built from filament formation of its subunits. Among peripheral membrane proteins satisfying these two conditions, actin and septins (human and yeast) were selected as protein fence candidates (Figure 2). Actin is one of the major components of the mammalian cytoskeleton. Septin has been reported to colocalize with actin filaments. All three proteins have been found near PIP2-abundant regions [18]-[20].

Figure 2
figure 2

Electrostatic potential on the molecular surfaces. (a) Actin, (b) Human septin, and (c) Yeast septin. Half of the biological units of the septins (Sept7-Sept6-Sept2 for human septin and Cdc11-Cdc12-Cdc3-Cdc10 for yeast septin) are shown here for simplicity. The blue regions shown in the bottom view of the septins are the polybasic regions, which are assumed to be the binding surface to the membrane in the simulation.

An actin fence was built from an X-ray fiber based model (PDB:2ZWH) of Oda et al. [21] and by adding additional monomers with 166.4° rotations and 27.6 Å translations along the fence axis to generate an actin filament. A human septin fence was built from the hexameric biological unit (Sept7-Sept6-Sept2-Sept2-Sept6-Sept7) in PDB:2QAG [22]. Missing coordinates of some residues in PDB:2QAG were reconstructed using I-TASSER [23], a protein structure prediction tool. The model structure was then aligned to the original PDB to determine the best-fit structure. Additional biological units were placed by 253.2 Å translations along the fence axis to generate a human septin filament.

There is no reported PDB structure for octameric yeast septin, but the full sequence and the octameric biological unit (Cdc11-Cdc12-Cdc3-Cdc10-Cdc10-Cdc3-Cdc12-Cdc11) of the protein is available. Atomic coordinates of the full sequence were generated using I-TASSER and the best-fit structure was selected from the alignment to the human septin model. The sequence identity between the subunits of yeast septin and human septin ranges from 32% to 38% (comparison in Protein BLAST [24]), which supports the similar structures to each other. Additional biological units were placed by 338.7 Å translations along the fence axis to generate a yeast septin filament.

The resulting septin structures form a long filamentous rod. The binding surfaces of the septins to the lipid membrane are assumed to be the one with the polybasic region that is conserved among septin family and its subunits (Additional file 1: Figure S1), thereby electrostatic attraction is expected to enhance the binding of the septins to the anionic lipid [18],[19]. The septins contain long C-terminal extensions that are believed to play a key role in binding with the other septins to make parallel filament bundle [25]. The C-terminal extensions were excluded here because they appear to contribute little to the single filament formation and binding of the septins to the membrane.

Continuum electrostatics based effective potentials

The interactions between PIP2-like charged spheres and the environments (water, membrane, and protein) were approximated by static electrostatic (ϕ PHIX) and repulsive core (U CORE) potentials [26]. Using the Poisson-Boltzmann PBEQ module [27],[28] in CHARMM [17], ϕ PHIX was calculated with an implicit membrane with a thickness of 36 Å and a dielectric constant of 2, the solution region with a dielectric constant of 80 and a salt concentration of 150 mM, and the proteins or artificial fence atoms with the optimized PB atomic radii [29] and a dielectric constant of 2. The electrostatic potential was first calculated with a coarse grid (grid-spacing of 1.0 Å), the result of which was then used to set the potential on a finer grid (grid-spacing of 0.5 Å). The PB radii augmented by 4.16 Å were used to setup the molecular surface by which a core repulsion potential map was built. U CORE was set to zero in all PIP2-accessible regions (outside proteins or artificial fence atoms) and 100 kcal/mol otherwise. The electrostatic and repulsive core potentials were stored on a grid (2,801 × 2,001 × 61 for the rod-like steric fence, 2,801 × 901 × 61 for the rod-like electrostatic fence, and 2,801 × 1,501 × 61 for the protein fence with a grid-spacing of 0.5 Å). Therefore, in this system representation, PIP2 molecules are the only explicitly simulated entities and the PIP2 energy is calculated by

U P I P 2 r i = j i ε i j r min / r i j 12 2 r min / r i j 6 + q i q j 80 r i j + U CORE r i + q i Φ PHIX r i

Note that ϕ PHIX and U CORE were computed once before the LD simulations, and their energies and forces were calculated using the 3rd-order B-spline interpolation during the simulations [30],[31]. The preparation of the simulation system (PDB manipulation), calculation of the potential maps (PBEQ module), and LD simulations were performed using CHARMM [17] and CHARMM-GUI ( [32].


The depletion of PIP2 population, i.e., the concentration decrease in the PIP2 pool from the initial 6%, characterizes the PIP2 diffusion through the fence. The number of particles in the pool was evaluated as a function of time and averaged over the independent runs. All decay curves were well fit by the stretched exponential function:

C t = C 0 - C exp - t / τ β + C

where C , τ, and β are iteratively determined from the simulation results with C 0 = 6%. The average relaxation time from the stretched exponential curve is given by

<τ>= τ β .Γ 1 / β

where Γ is the gamma function defined as Γ x = 0 t x 1 e - t dt. The retardation of diffusion is then quantified by the concentration relaxation time ratio ξ = < τ free >/< τ >, where < τ free > is the relaxation time in the absence of a barrier, i.e., free diffusion.

The potential of mean force (PMF) was calculated from position dependent concentrations of PIP2 after equilibrium had been reached [33]:

W 1 D x = k B T.ln C eq x / C r e f

where C ref is set to the bulk concentration under equilibrium (2%). The PMF represents the potential profile in the presence of PIP2 molecules in the equilibrated simulation system.

Results and Discussion

Rod-like fences

Additional file 1: Figure S2 shows the characteristics of the rod-like fences. The steric fences range from nearly continuous to 50% breached (Additional file 1: Figure S2A), and the electrostatic fences from strongly negative to strongly positive (Additional file 1: Figure S2B). The positive (attractive potential well) and negative charges (repulsive potential barrier) yield symmetric electrostatic potential profiles (Additional file 1: Figure S2C).

The retardation of PIP2 diffusion through the fence was estimated by the relative change of the relaxation time with respect to the free diffusion case (L open = 1,000 Å): ξ = < τ free >/< τ >. The time-series of the PIP2 concentration inside the pool (depletion) for the porous fences and its fitted curves are shown in Figure 3a and the resulting relaxation times with respect to free diffusion (ξ defined above) are shown in Figure 3b (see Additional file 1: Figure S3 for representative snapshots from the trajectory). There is little retardation of diffusion for the larger gaps, and even for L open = 20 Å (which is 1.85 times larger than the size of PIP2) the relaxation times with respect to free diffusion are only reduced by a factor of 6. Hence, the gap in the fence should be very close to the size of PIP2 to reproduce the experimentally observed factor of 100 reduction in diffusion [4]. This result highlights ineffectiveness of purely steric repulsion for blocking diffusion; once a diffusing particle is in the vicinity of a gap, it has multiple opportunities to escape through the gap due to its stochastic motions. (For ballistic trajectories, passage through a barrier is simply proportional to the fraction of the barrier that is not blocked.) Additionally, as deduced from a comparison with a continuum calculation for free diffusion, repulsion between the PIP2 particles (q = -4e) in the corral provides as the second driving force for dispersal.

Figure 3
figure 3

Results for rod - like steric fences. (a) Time-series of PIP2 concentration in the pool with different fence opening ratio (L open/L) with L = 1,000 Å. (b) Ratio of the relaxation time with respect to free diffusion (<τ free > for L open = L). The dotted line labeled as Continuum is obtained from the one-dimensional diffusion of continuum matter (unpublished).

Turning to the electrostatic models, the concentration relaxation time ratios ξ for assorted cases are plotted in Figure 4 (the depletion curves in Additional file 1: Figure S4). The general trends are not surprising. For both negatively and positively charged rods, blockage increases with increasing charge and decreasing height above the surface. However, the negative rods are more effective at the same charge magnitude and height. The difference arises because negatively charged rod directly repels PIP2 (there is an energy barrier at X = 0, directly under the rod), while the positively charged rod attracts PIP2 (there is an energy well at this position). In particular, from Figure 4, complete blockade of PIP2 diffusion is observed for q = -1.0e when h ≤ 7 Å, q ≤ -0.5e when h = 5 Å, and |q| ≥ 0.5e when h ≤ 2 Å. These conditions correspond to electrostatic energy barriers larger than 7.5 kcal/mol and energy wells deeper than 12.5 kcal/mol, as calculated for the interaction of an individual PIP2 with the rod.

Figure 4
figure 4

Concentration relaxation time ratios for rod - like electrostatic fences. Concentration relaxation time ratio for the electrostatic fence simulations, calculated from Figure 4: (a) q < 0 and (b) q > 0. < τ free > refers to the free diffusion case (q = 0).

The explanation for the difference requires consideration of the interactions between PIP2. The potential of mean force (PMF) provides the entry point for this analysis. Figure 5 plots the PMFs calculated from equilibrium simulations of initial uniform 2% PIP2 distribution using Eq. (4). Negative charges on the rod lead to a nearly impermeable fence when the PMF barrier is larger than 6.5 kcal/mol. The difference between the electrostatic potential barrier of 7.5 kcal/mol (at X = 0) and the PMF barrier of 6.5 kcal/mol is due to the lack of PIP2 sampling near X = 0. Positive charges result in accumulation of PIP2 in the potential well, which, in turn, can exert a potential barrier for the other PIP2 to enter and then move through the fence; i.e., there is both a well under the rod and a barrier on each side. A nearly complete blockade occurs when the trapped PIP2 in the potential well below the fence exerts an energy barrier of around 6.5 kcal/mol to the PIP2 outside of the fence (Figure 5, h = 2 Å).

Figure 5
figure 5

Potentials of mean force for rod - like electrostatic fences. Potentials of mean force (PMFs) for h = 5 Å and 2 Å, calculated from equilibrium simulations with the electrostatic fences: (a) q < 0 and (b) q > 0.

Protein fences

Figure 2 shows the electrostatic potentials from the protein charges projected on the molecular surfaces of (A) actin, (B) human septin, and (C) yeast septin. The polybasic amino acids of the septins (Additional file 1: Figure S1) are assumed to be the binding surfaces to the membrane (blue regions in “bottom view”). These electrostatic potentials are combined with the surrounding dielectrics and steric repulsive potentials from the proteins to develop the total potentials. Figure 6 maps the potentials for PIP2 on the diffusion plane at Zmin = 0, -5, -10, and -15 Å, where Zmin is the minimum location of the protein atoms to the membrane surface (Z = 0; the lipid phosphate plane). Figure 6d shows the different components of the potential at Zmin = -9 Å for human septin and Zmin = -12 Å for yeast septin where a complete blockade is observed from trajectory analysis (Figure 7). Importantly, when Zmin = 0 Å (at least one protein atom is in contact with the diffusion plane), the steric and electrostatic potentials yield a relatively porous projection on the surface.

Figure 6
figure 6

Potential maps of the rigid - protein fences. The potential maps of the rigid-protein fences at four depths (Zmin) of protein binding to the lipid membrane, Zmin = 0, -5, -10, and -15 Å (from left to right): (a) Actin, (b) Human septin, and (c) Yeast septin. (d) Comparison of the steric potential (U CORE) and the steric plus electrostatic potential (U CORE + qϕPHIX) for human septin with Zmin = -9 Å and yeast septin with Zmin = -12 Å.

Figure 7
figure 7

PIP 2 depletion for the rigid - protein fences. Comparison of the concentration relaxation time ratio for the rigid-protein fences at different burial depths (Zmin).

The concentration relaxation time ratio ξ is plotted as a function of Zmin in Figure 7. Consistent from the potential maps in Figure 6, fencing is ineffective when Zmin = 0 Å. This result is anticipated from the steric rod models, where even a small gap results in nearly free diffusion (Figure 3). Actin has a pronounced arch-shape (Figure 2a), which leaves unblocked fence even with the deepest burial (Figure 6a, Zmin = -15 v). Consequently, all the actin fence results show significant PIP2 escape out of corral (Figure 7) for all burial depths. Septins provide essentially complete blockage (ξ ≈ zero) when Zmin = -9 Å for the human septin fence and Zmin = -12 Å for the yeast septin fence. Specifically, the potential maps calculated at Zmin = -9 Å for the human septin fence and Zmin = -12 Å for the yeast septin fence (Figure 6d) show the “blanket coverage” required for blockage found in the rod-like fence studies. Interestingly, human and yeast septin block PIP2 diffusion in different ways as revealed by the comparison of the steric (U CORE) and electrostatic (ϕ PHIX) potentials. U CORE from the human septin fence at Zmin = -9 Å provides complete coverage. In contrast, narrow gaps in U CORE from yeast septin fence at Zmin = -12 Å are filled in by ϕ PHIX for complete blockade.

Given the preceding simulation results that the septin fences can block PIP2 diffusion when buried at certain depths from the membrane surface, the next question is the likelihood of having such burial depths. While it is challenging to quantify the penetration depth of peripheral membrane proteins in both experimental and computational studies, there are three qualitative arguments that such burial depths are possible. First, there is the polybasic region exposed on the putative membrane-binding protein surface, which enhances electrostatic binding to the anionic lipid membrane (Additional file 1: Figure S1). Second, there is abundance in hydrophobic residues (PHE, TYR, and TRP) on the membrane-binding protein surface, especially those with hydrophobic rings. As shown in Additional file 1: Figure S1, both human and yeast septins have hydrophobic residues in the loops near the polybasic region and the other basic residues. These hydrophobic rings can act as binding anchors into the hydrophobic core of the lipid membrane. It has been reported that anchoring of MARCKS effector domain to the lipid membrane is enhanced by the penetration of five PHE rings into the lipid hydrophobic core [8]. Third, there is no penetration of charged residues in the present rigid models into the hydrophobic region even at Zmin = -10 Å, i.e., 10 Å below the lipid phosphate plane (Figure 8). Thus, the charged residues can stay in the head region without invoking large free energy penalty due to their burial into the hydrophobic core. However, the remaining 2 Å penetration for yeast septin is not guaranteed in the current rigid-protein fence model. Further (computationally very demanding) studies would be required to explore the influence of the protein conformation change on membrane binding and to determine whether such charged residues can re-orient themselves to avoid the burial penalty.

Figure 8
figure 8

Number of charged residues inside the hydrophobic region. Number of charged residues inside the hydrophobic region as a function of the burial depths (Zmin) of the human and yeast septin fences. The number of charged residues is counted in 1 Å bin along the Z-axis.

Phenomenological treatment of permeability

Let us assume, for simplicity, a square corral of 1 μm per side. From Hilgemann's analysis [34], we can also assume that there are 100 kinases inside the corral and each kinase produces 100 PIP2 molecules per second, hence production rate is of 104 molecules/(μm2s). The permeability P through a fence region can be defined as P = k D/d, where k is the partition coefficient of PIP2 from the bulk membrane into the fence region, D is the diffusion coefficient, and d is the thickness of the fence region. Assuming d = 5 nm and D =1 μm2/s (in the bulk region [4]), P = 200 k μm/s. The flux per unit length F of PIP2 out of the fence can be simplified using Fick's law as F = P (C i - C o ), where C i (from our rigid-protein model, 105 μm-2) is the areal concentration of PIP2 inside the fence and C o (3 × 104 μm-2) that in the bulk. The steady-state PIP2 flux out of the corral can be estimated as F = 1.4 × 107k / (μms). As the condition for the corral to act as a barrier maintaining the concentrated pool of PIP2, the outward flux should be kept below of the production rate. If the outward flux is only through one side of the square, the partition coefficient k should be smaller than 7 × 10-3 (= 104 / (1.4 × 107)). In other words, if the fence is leaky to PIP2 above this value, the fence would not be an effective barrier to maintain the PIP2 pool. This impermeability calculation for an effective protein fence is in good agreement with our simulation results, which predicts virtually complete blockade from U CORE or U CORE + q PIP2ϕ PHIX barrier (Figures 6d and 7).


Langevin dynamics simulations of PIP2 on a membrane surface with rod-like steric and electrostatic as well as all-atom protein fence models were carried out to determine the general conditions for fencing PIP2 and the specific fencing ability of actin, human septin, and yeast septin. Simulations on the model systems indicate that even a small gap leads to ineffective blocking. Likewise, electrostatic blockage is only effective at very high charge density, and is unlikely to play a major role in blockage of PIP2 diffusion.

These observations place significant limitations on the abilities of individual proteins to form effective fences. In fact, single filaments of actin, human septin, and yeast septin provided little blockage when placed on the membrane surface. Even burial to 15 Å did not yield significant blocking by actin, as could be anticipated by its pronounced arch-like shape. However, the two septins did provide blockage consistent with experiment when the human septin is buried 9 Å and the yeast septin 12 Å below the membrane surface. Implicit membrane-solvent models indicate that burial to 10 Å can be achieved, though further penetration requires protein conformational changes. All-atom simulations with fully flexible proteins and explicit membrane and solvent will undoubtedly yield further insight to the mechanism of blockage.

Though not simulated here, a fence could be made more effective by adding more rows or more components. For example, if each of three rows provides 80% blockage and the rows are independent, the array leads to more than 99% blockage. Similarly, a small gap in a single row could be filled by a membrane-associated peptide. In closing, it was not the purpose of this study to determine the molecular level mechanism of corralling of PIP2 on the cell surface by protein fences. Rather, it was to elucidate some of the molecular level considerations that must be applied to models for fencing, and to stimulate more detailed simulation studies.

Authors’ contributions

KIL carried out the Langevin dynamics simulations and data analysis, and drafted the manuscript. WI and RWP conceived of the study, and participated in its design and coordination. WI and RWP also participated in the data analysis and preparation of the manuscript. The authors read and approved the final manuscript.

Additional file


  1. Di Paolo G, De Camilli P: Phosphoinositides in cell regulation and membrane dynamics. Nature. 2006, 443: 651-657. 10.1038/nature05185.

    Article  Google Scholar 

  2. McLaughlin S, Murray D: Plasma membrane phosphoinositide organization by protein electrostatics. Nature. 2005, 438: 605-611. 10.1038/nature04398.

    Article  Google Scholar 

  3. Fujiwara T, Ritchie K, Murakoshi H, Jacobson K, Kusumi A: Phospholipids undergo hop diffusion in compartmentalized cell membrane. J Cell Biol. 2002, 157: 1071-1081. 10.1083/jcb.200202050.

    Article  Google Scholar 

  4. Golebiewska U, Kay JG, Masters T, Grinstein S, Im W, Pastor RW, Scarlata S, McLaughlin S: Evidence for a fence that impedes the diffusion of phosphatidylinositol 4,5-bisphosphate out of the forming phagosomes of macrophages. Mol Biol Cell. 2011, 22: 3498-3507. 10.1091/mbc.E11-02-0114.

    Article  Google Scholar 

  5. Kusumi A, Shirai YM, Koyama-Honda I, Suzuki KGN, Fujiwara TK: Hierarchical organization of the plasma membrane: Investigations by single-molecule tracking vs. fluorescence correlation spectroscopy. FEBS Lett. 2010, 584: 1814-1823. 10.1016/j.febslet.2010.02.047.

    Article  Google Scholar 

  6. van den Bogaart G, Meyenberg K, Risselada HJ, Amin H, Willig KI, Hubrich BE, Dier M, Hell SW, Grubmuller H, Diederichsen U, Jahn R: Membrane protein sequestering by ionic protein-lipid interactions. Nature. 2011, 479: 552-555. 10.1038/nature10545.

    Article  Google Scholar 

  7. Arbuzova A, Wang L, Wang J, Hangyás-Mihályné G, Murray D, Honig B, McLaughlin S: Membrane binding of peptides containing both basic and aromatic residues. Experimental studies with peptides corresponding to the scaffolding region of caveolin and the effector region of MARCKS. Biochemistry. 2000, 39: 10330-10339. 10.1021/bi001039j.

    Article  Google Scholar 

  8. Ellena JF, Burnitz MC, Cafiso DS: Location of the myristoylated alanine-rich C-kinase substrate (MARCKS) effector domain in negatively charged phospholipid bicelles. Biophys J. 2003, 85: 2442-2448. 10.1016/S0006-3495(03)74667-0.

    Article  Google Scholar 

  9. Golebiewska UP, Kay J, Grinstein S, Im W, Pastor RW, Scarlata S, McLaughlin S: Don't Fence Me in: evidence for a ‘fence’ that impedes the diffusion of PIP 2into and Out of nascent phagosomes in macrophages. Biophys J 2010, 98:688a.,

  10. Chandrasekhar S: Stochastic problems in physics and astronomy. Rev Mod Phys. 1943, 15: 1-89. 10.1103/RevModPhys.15.1.

    Article  MATH  MathSciNet  Google Scholar 

  11. Zwanzig R: Langevin Theory of Polymer Dynamics and Dilute Solutions. Stochastic Processes in Chemical Physics. Volume 15. Edited by: Shuler KE. 1969, John Wiley and Sons, New York, 325-332.

    Google Scholar 

  12. Pastor RW, Brooks BR, Szabo A: An analysis of the accuracy of Langevin and molecular dynamics algorithms. Molecular Physics. 1988, 65: 1409-1419. 10.1080/00268978800101881.

    Article  Google Scholar 

  13. Pastor RW: Techniques and Applications of Langevin Dynamics Simulations. The Molecular Dynamics of Liquid Crystals. Volume 431. Edited by: Luckhurst GR, Veracini CA. 1994, NATO ASI Series], Springer Netherlands, 85-138. 10.1007/978-94-011-1168-3_5.

    Chapter  Google Scholar 

  14. Pastor RW, Karplus M: Inertial effects in butane stochastic dynamics. J Chem Phys. 1989, 91: 211-218. 10.1063/1.457508.

    Article  Google Scholar 

  15. Wu EL, Qi Y, Song KC, Klauda JB, Im W: Preferred orientations of phosphoinositides in bilayers and their implications in protein recognition mechanisms. J Phys Chem B. 2014, 118: 4315-4325. 10.1021/jp500610t.

    Article  Google Scholar 

  16. Li Z, Venable RM, Rogers LA, Murray D, Pastor RW: Molecular Dynamics simulations of PIP(2) and PIP(3) in lipid bilayers: determination of ring orientation, and the effects of surface roughness on a Poisson-Boltzmann description. Biophys J. 2009, 97: 155-163. 10.1016/j.bpj.2009.04.037.

    Article  Google Scholar 

  17. Brooks BR, Brooks CL, Mackerell AD, Nilsson L, Petrella RJ, Roux B, Won Y, Archontis G, Bartels C, Boresch S, Caflisch A, Caves L, Cui Q, Dinner AR, Feig M, Fischer S, Gao J, Hodoscek M, Im W, Kuczera K, Lazaridis T, Ma J, Ovchinnikov V, Paci E, Pastor RW, Post CB, Pu JZ, Schaefer M, Tidor B, Venable RM, et al: CHARMM: the biomolecular simulation program. J Comput Chem. 2009, 30: 1545-1614. 10.1002/jcc.21287.

    Article  Google Scholar 

  18. Zhang J, Kong C, Xie H, McPherson PS, Grinstein S, Trimble WS: Phosphatidylinositol polyphosphate binding to the mammalian septin H5 is modulated by GTP. Curr Biol. 1999, 9: 1458-1467. 10.1016/S0960-9822(00)80115-3.

    Article  Google Scholar 

  19. Casamayor A, Snyder M: Molecular dissection of a yeast septin: distinct domains are required for septin interaction, localization, and function. Mol Cell Biol. 2003, 23: 2762-2777. 10.1128/MCB.23.8.2762-2777.2003.

    Article  Google Scholar 

  20. Rodríguez-Escudero I, Roelants FM, Thorner J, Nombela C, Molina M, Cid VJ: Reconstitution of the mammalian PI3K/PTEN/Akt pathway in yeast. Biochem J. 2005, 390: 613-623. 10.1042/BJ20050574.

    Article  Google Scholar 

  21. Oda T, Iwasa M, Aihara T, Maeda Y, Narita A: The nature of the globular-to fibrous-actin transition. Nature. 2009, 457: 441-445. 10.1038/nature07685.

    Article  Google Scholar 

  22. Sirajuddin M, Farkasovsky M, Hauer F, Kuhlmann D, Macara IG, Weyand M, Stark H, Wittinghofer A: Structural insight into filament formation by mammalian septins. Nature. 2007, 449: 311-315. 10.1038/nature06052.

    Article  Google Scholar 

  23. Roy A, Kucukural A, Zhang Y: I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010, 5: 725-738. 10.1038/nprot.2010.5.

    Article  Google Scholar 

  24. Karlin S, Altschul SF: Applications and statistics for multiple high-scoring segments in molecular sequences. Proc Natl Acad Sci. 1993, 90: 5873-5877. 10.1073/pnas.90.12.5873.

    Article  Google Scholar 

  25. Oh Y, Bi E: Septin structure and function in yeast and beyond. Trends Cell Biol. 2011, 21: 141-148. 10.1016/j.tcb.2010.11.006.

    Article  Google Scholar 

  26. Rathinavelan T, Zhang L, Picking WL, Weis DD, De Guzman RN, Im W: A repulsive electrostatic mechanism for protein export through the type III secretion apparatus. Biophys J. 2010, 98: 452-461. 10.1016/j.bpj.2009.10.030.

    Article  Google Scholar 

  27. Im W, Beglov D, Roux B: Continuum solvation model: computation of electrostatic forces from numerical solutions to the Poisson-Boltzmann equation. Comput Phys Commun. 1998, 111: 59-75. 10.1016/S0010-4655(98)00016-2.

    Article  MATH  Google Scholar 

  28. Jo S, Vargyas M, Vasko-Szedlar J, Roux B, Im W: PBEQ-Solver for online visualization of electrostatic potential of biomolecules. Nucleic Acids Res. 2008, 36: W270-W275. 10.1093/nar/gkn314.

    Article  Google Scholar 

  29. Nina M, Im W, Roux B: Optimized atomic radii for protein continuum electrostatics solvation forces. Biophys Chem. 1999, 78: 89-96. 10.1016/S0301-4622(98)00236-1.

    Article  Google Scholar 

  30. Im W, Seefeld S, Roux B: A Grand Canonical Monte Carlo-Brownian dynamics algorithm for simulating ion channels. Biophys J. 2000, 79: 788-801. 10.1016/S0006-3495(00)76336-3.

    Article  Google Scholar 

  31. Im W, Roux BT: Ion permeation and selectivity of OmpF porin: a theoretical study based on molecular dynamics, Brownian dynamics, and continuum electrodiffusion theory. J Mol Biol. 2002, 322: 851-869. 10.1016/S0022-2836(02)00778-7.

    Article  Google Scholar 

  32. Jo S, Kim T, Iyer VG, Im W: CHARMM-GUI: a web-based graphical user interface for CHARMM. J Comput Chem. 2008, 29: 1859-1865. 10.1002/jcc.20945.

    Article  Google Scholar 

  33. Roux B: The calculation of the potential of mean force using computer-simulations. Comput Phys Commun. 1995, 91: 275-282. 10.1016/0010-4655(95)00053-I.

    Article  Google Scholar 

  34. Yaradanakul A, Feng S, Shen C, Lariccia V, Lin M-J, Yang J, Dong P, Yin HL, Albanesi JP, Hilgemann DW: Dual control of cardiac Na +-Ca 2+exchange by PIP 2: electrophysiological analysis of direct and indirect mechanisms. J Physiol 2007, 582:991-1010.,

Download references


We gratefully acknowledge Stuart McLaughlin for motivating the present study, and for helpful comments on the manuscript. This work was supported by the National Science Foundation (NSF MCB-1157677 to WI), XSEDE Resources (TG-MCB070009 to WI), and the Intramural Research Program of the NIH, National Heart, Lung, and Blood Institute (to RWP).

Author information

Authors and Affiliations


Corresponding authors

Correspondence to Wonpil Im or Richard W Pastor.

Additional information

Competing interests

The authors declare that they have no competing interests.

Electronic supplementary material


Additional file 1: Figure S1.: Basic and hydrophobic residues of the septins. Basic (blue: ARG and LYS) residues found in the polybasic regions and hydrophobic residues. with aromatic ring (magenta: PHE, green: TRP, and gray: TYR): (A) Human septin and (B). Yeast septin. Figure S2. Potential energy for rod-like fences. (A) Steric repulsive potential energy from the porous fences according to the fence opening. (L open) with the fence length L = 1,000 Å. (B-C) Potential energy from the electrostatic fences according to the charge per fence atom (q) and height (h) of the charged rod from the diffusion plane with L = 450 Å: (B) 2D potential map and (C) 1D profile along the X-axis (diffusion direction, Y = 0). Figure S3. Representative snapshots of the PIP2 depletion from the trajectory. Trajectory snapshots for human septin fence with Zmin = 8 Å taken at (A) initial configuration (t = 0), (B) t = 0.1 μs, (C) t = 1 μs, and (D) t = 5 μs. Figure S4. Results for rod-like electrostatic fences. Time-series of PIP2 concentration in the pool with different heights (h) of the charged bar from the diffusion plane and atomic charge (q): (A) h = 10 Å, (B) h = 7 Å, (C) h = 5 Å, and (D) h = 2 Å. The positive fence charge cases are shown in solid lines and negative fence charge cases in dotted lines. Black, red, green, and blue lines are for |q| = 0.05e, 0.10e, 0.50e, and 1.00e, respectively. (PDF 2 MB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lee, K.I., Im, W. & Pastor, R.W. Langevin dynamics simulations of charged model phosphatidylinositol lipids in the presence of diffusion barriers: toward an atomic level understanding of corralling of PIP2 by protein fences in biological membranes. BMC Biophys 7, 13 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: