Characterization of Specific Ion Effects on PI(4,5)P2 Clustering: Molecular Dynamics Simulations and Graph-Theoretic Analysis

The impact of volatile anesthetics and propofol on phosphatidylinositol 4,5-bisphosphate signaling

Phosphatidylinositol 4,5-bisphosphate (PIP2), as well as other anionic phospholipids, play a pivotal role in various cellular processes, including ion channel regulation, receptor trafficking, and intracellular signaling pathways. The binding of volatile anesthetics and propofol to PIP2 leads to alterations in PIP2-mediated signaling causing modulation of ion channels such as ɣ-aminobutyric acid type A (GABAA) receptors, voltage-gated calcium channels, and potassium channels through various mechanisms. Additionally, the interaction between anionic phospholipids and G protein-coupled receptors plays a critical role in various anesthetic pathways, with these anesthetic-induced changes impacting PIP2 levels which cause cascading effects on receptor trafficking, including GABAA receptor internalization. This comprehensive review of various mechanisms of interaction provides insights into the intricate interplay between PIP2 signaling and anesthetic-induced changes, shedding light on the molecular mechanisms underlying anesthesia.

Biophysical Modeling of Synaptic Plasticity

Dendritic spines are small, bulbous compartments that function as postsynaptic sites and undergo intense biochemical and biophysical activity. The role of the myriad signaling pathways that are implicated in synaptic plasticity is well studied. A recent abundance of quantitative experimental data has made the events associated with synaptic plasticity amenable to quantitative biophysical modeling. Spines are also fascinating biophysical computational units because spine geometry, signal transduction, and mechanics work in a complex feedback loop to tune synaptic plasticity. In this sense, ideas from modeling cell motility can inspire us to develop multiscale approaches for predictive modeling of synaptic plasticity. In this article, we review the key steps in postsynaptic plasticity with a specific focus on the impact of spine geometry on signaling, cytoskeleton rearrangement, and membrane mechanics. We summarize the main experimental observations and highlight how theory and computation can aid our understanding of these complex processes.

PIP2 modulates TRPC3 activity via TRP helix and S4-S5 linker

The transient receptor potential canonical type 3 (TRPC3) channel plays a pivotal role in regulating neuronal excitability in the brain via its constitutive activity. The channel is intricately regulated by lipids and has previously been demonstrated to be positively modulated by PIP2. Using molecular dynamics simulations and patch clamp techniques, we reveal that PIP2 predominantly interacts with TRPC3 at the L3 lipid binding site, located at the intersection of pre-S1 and S1 helices. We demonstrate that PIP2 sensing involves a multistep mechanism that propagates from L3 to the pore domain via a salt bridge between the TRP helix and S4-S5 linker. Notably, we find that both stimulated and constitutive TRPC3 activity require PIP2. These structural insights into the function of TRPC3 are invaluable for understanding the role of the TRPC subfamily in health and disease, in particular for cardiovascular diseases, in which TRPC3 channels play a major role.

Ion-Induced PIP2 Clustering with Martini3: Modification of Phosphate-Ion Interactions and Comparison with CHARMM36

Phosphatidylinositol 4,5-bisphosphate (PIP2) is a critical lipid for cellular signaling. The specific phosphorylation of the inositol ring controls protein binding as well as clustering behavior. Two popular models to describe ion-mediated clustering of PIP2 are Martini3 (M3) and CHARMM36 (C36). Molecular dynamics simulations of PIP2-containing bilayers in solutions of potassium chloride, sodium chloride, and calcium chloride, and at two different resolutions are performed to understand the aggregation and the model parameters that drive it. The average M3 clusters of PIP2 in bilayers of 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine and PIP2 bilayers in the presence of K+, Na+, or Ca2+ contained 2.2, 2.6, and 6.4 times more PIP2 than C36 clusters, respectively. Indeed, the Ca2+-containing systems often formed a single large aggregate. Reparametrization of the M3 ion-phosphate Lennard-Jones interaction energies to reproduce experimental osmotic pressure of sodium dimethyl phosphate (DMP), K[DMP], and Ca[DMP]2 solutions, the same experimental target as C36, yielded comparably sized PIP2 clusters for the two models. Furthermore, C36 and the modified M3 predict similar saturation of the phosphate groups with increasing Ca2+, although the coarse-grained model does not capture the cooperativity between K+ and Ca2+. This characterization of the M3 behavior in the presence of monovalent and divalent ions lays a foundation to study cation/protein/PIP2 clustering.

FERM domains recruit ample PI(4,5)P2s to form extensive protein-membrane attachments

The four-point-one ezrin-radixin-moesin homology (FERM) protein domain is a multifunctional protein-lipid binding site, constituting an integral part of numerous membrane-associated proteins. Its interaction with the lipid phosphatidylinositol-4,5-bisphosphate (PIP2), located at the inner leaflet of eukaryotic plasma membranes, is important for localization, anchorage, and activation of FERM-containing proteins. FERM-PIP2 complexes structurally determined so far exclusively feature a 1:1 binding stoichiometry of protein and lipid, with a few basic FERM residues neutralizing the -4 charge of the bound PIP2. Whether this picture from static crystal structures also applies to the dynamic interaction of FERM domains on PIP2 membranes is unknown. We here quantified the stoichiometry of FERM-PIP2 binding in a lipid bilayer using atomistic molecular dynamics simulations and experiments on solid supported membranes for the FERM domains of focal adhesion kinase and ezrin. In contrast to the structural data, we find much higher average stoichiometries of FERM-PIP2 binding, amounting to 1:3 or 1:4 ratios, respectively. In simulations, the full set of basic residues at the membrane interface, 7 and 15 residues for focal adhesion kinase and ezrin, respectively, engages in PIP2 interactions. In addition, Na ions enter the FERM-membrane binding interface, compensating negative PIP2 charges in case of high charge surpluses from bound PIP2. We propose the multivalent binding of FERM domains to PIP2 in lipid bilayers to significantly enhance the stability of FERM-membrane binding and to render the FERM-membrane linkage highly adjustable.

Exploring CRD mobility during RAS/RAF engagement at the membrane

During the activation of mitogen-activated protein kinase (MAPK) signaling, the RAS-binding domain (RBD) and cysteine-rich domain (CRD) of RAF bind to active RAS at the plasma membrane. The orientation of RAS at the membrane may be critical for formation of the RAS-RBDCRD complex and subsequent signaling. To explore how RAS membrane orientation relates to the protein dynamics within the RAS-RBDCRD complex, we perform multiscale coarse-grained and all-atom molecular dynamics (MD) simulations of KRAS4b bound to the RBD and CRD domains of RAF-1, both in solution and anchored to a model plasma membrane. Solution MD simulations describe dynamic KRAS4b-CRD conformations, suggesting that the CRD has sufficient flexibility in this environment to substantially change its binding interface with KRAS4b. In contrast, when the ternary complex is anchored to the membrane, the mobility of the CRD relative to KRAS4b is restricted, resulting in fewer distinct KRAS4b-CRD conformations. These simulations implicate membrane orientations of the ternary complex that are consistent with NMR measurements. While a crystal structure-like conformation is observed in both solution and membrane simulations, a particular intermolecular rearrangement of the ternary complex is observed only when it is anchored to the membrane. This configuration emerges when the CRD hydrophobic loops are inserted into the membrane and helices α3-5 of KRAS4b are solvent exposed. This membrane-specific configuration is stabilized by KRAS4b-CRD contacts that are not observed in the crystal structure. These results suggest modulatory interplay between the CRD and plasma membrane that correlate with RAS/RAF complex structure and dynamics, and potentially influence subsequent steps in the activation of MAPK signaling.

Molecular mechanisms of spontaneous curvature and softening in complex lipid bilayer mixtures

Membrane reshaping is an essential biological process. The chemical composition of lipid membranes determines their mechanical properties and thus the energetics of their shape. Hundreds of distinct lipid species make up native bilayers, and this diversity complicates efforts to uncover what compositional factors drive membrane stability in cells. Simplifying assumptions, therefore, are used to generate quantitative predictions of bilayer dynamics based on lipid composition. One assumption commonly used is that "per lipid" mechanical properties are both additive and constant-that they are an intrinsic property of lipids independent of the surrounding composition. Related to this is the assumption that lipid bulkiness, or "shape," determines its curvature preference, independently of context. In this study, all-atom molecular dynamics simulations on three separate multilipid systems were used to explicitly test these assumptions, applying methodology recently developed to isolate properties of single lipids or nanometer-scale patches of lipids. The curvature preference experienced by populations of lipid conformations were inferred from their redistribution on a dynamically fluctuating bilayer. Representative populations were extracted by both structural similarity and semi-automated hidden Markov model analysis. The curvature preferences of lipid dimers were then determined and compared with an additive model that combines the monomer curvature preference of both the individual lipids. In all three systems, we identified conformational subpopulations of lipid dimers that showed non-additive curvature preference, in each case mediated by a special chemical interaction (e.g., hydrogen bonding). Our study highlights the importance of specific chemical interactions between lipids in multicomponent bilayers and the impact of interactions on bilayer stiffness. We identify two mechanisms of bilayer softening: diffusional softening, driven by the dynamic coupling between lipid distributions and membrane undulations, and conformational softening, driven by the inter-conversion between distinct dimeric conformations.

Design principles of PI(4,5)P2 clustering under protein-free conditions: Specific cation effects and calcium-potassium synergy

SignificanceClustering of phosphatidylinositol 4,5-bisphosphate (PIP₂) with proteins into what are known as "PIP₂ rafts" is a critical component of intracellular signaling, yet little is known about PIP₂ clusters at the atomic level. Using molecular dynamics simulations and network theory, this paper shows that Ca²⁺ generates large clusters by linking PIP₂ dimers already formed by doubly charged P4/P5 phosphates, while monovalent cations form smaller and less-stable clusters by adding PIP₂ monomers preferentially via weaker interactions with P4/P5 (for Na⁺) or with glycerol P1 (for K⁺). Synergy arises between K⁺ and Ca²⁺ because each ion forms linkages with different phosphates, thereby giving clusters more ways to grow. This explains why Ca²⁺ is pumped into cells by ion channels to form PIP₂ rafts.

Effects of Cholesterol and PIP2 on Interactions between Glycophorin A and Band 3 in Lipid Bilayers

In the erythrocyte membrane, the interactions between Glycophorin A (GPA) and Band 3 are associated strongly with the biological function of the membrane and several blood disorders. In this work, using coarse-grained molecular dynamics simulation, we systematically investigate the effects of cholesterol and phosphatidylinositol-4,5-bisphosphate (PIP2) on the interactions of GPA with Band 3 in the model erythrocyte membranes. We examine the dynamics of the interactions of GPA with Band 3 in different lipid bilayers on the microsecond time scale and calculate the binding free energy between GPA and Band 3. The results indicate that cholesterols thermodynamically favor the binding of GPA to Band 3 by increasing the thickness of the lipid bilayer and by producing an effective attraction between the proteins due to the depletion effect. Cholesterols also slow the kinetics of the binding of GPA to Band 3 by reducing the lateral mobility of the lipids and proteins and may influence the binding sites between the proteins. The anionic PIP2 lipids prefer binding to the surface of the proteins through electrostatic attraction between the PIP2 headgroup and the positively charged residues on the protein surface. Ions in the solvent facilitate the PIP2 aggregation which promotes the binding of GPA to Band 3.

Two-Dimensional Interfacial Exchange Diffusion Has the Potential to Augment Spatiotemporal Precision of Ca2+ Signaling

Nano-junctions between the endoplasmic reticulum and cytoplasmic surfaces of the plasma membrane and other organelles shape the spatiotemporal features of biological Ca²⁺ signals. Herein, we propose that 2D Ca²⁺ exchange diffusion on the negatively charged phospholipid surface lining nano-junctions participates in guiding Ca²⁺ from its source (channel or carrier) to its target (transport protein or enzyme). Evidence provided by in vitro Ca²⁺ flux experiments using an artificial phospholipid membrane is presented in support of the above proposed concept, and results from stochastic simulations of Ca²⁺ trajectories within nano-junctions are discussed in order to substantiate its possible requirements. Finally, we analyze recent literature on Ca²⁺ lipid interactions, which suggests that 2D interfacial Ca²⁺ diffusion may represent an important mechanism of signal transduction in biological systems characterized by high phospholipid surface to aqueous volume ratios.

Improved Parameterization of Phosphatidylinositide Lipid Headgroups for the Martini 3 Coarse-Grain Force Field

Phosphoinositides are a family of membrane phospholipids that play crucial roles in membrane regulatory events. As such, these lipids are often a key part of molecular dynamics simulation studies of biological membranes, in particular of those employing coarse-grain models because of the potential long times and sizes of the involved membrane processes. Version 3 of the widely used Martini coarse-grain force field has been recently published, greatly refining many aspects of biomolecular interactions. In order to properly use it for lipid membrane simulations with phosphoinositides, we put forth the Martini 3-specific parameterization of inositol, phosphatidylinositol, and seven physiologically relevant phosphorylated derivatives of phosphatidylinositol. Compared to parameterizations for earlier Martini versions, focus was put on a more accurate reproduction of the behavior seen in both atomistic simulations and experimental studies, including the signaling-relevant phosphoinositide interaction with divalent cations. The models that we develop improve upon the conformational dynamics of phosphoinositides in the Martini force field and provide stable topologies at typical Martini time steps. They are able to reproduce experimentally known protein-binding poses as well as phosphoinositide aggregation tendencies. The latter was tested both in the presence and absence of calcium and included correct behavior of PI(4,5)P₂ calcium-induced clusters, which can be of relevance for regulation.

Molecular and mesoscopic geometries in autophagosome generation. A review

Autophagy is an essential process in cell self-repair and survival. The centre of the autophagic event is the generation of the so-called autophagosome (AP), a vesicle surrounded by a double membrane (two bilayers). The AP delivers its cargo to a lysosome, for degradation and re-use of the hydrolysis products as new building blocks. AP formation is a very complex event, requiring dozens of specific proteins, and involving numerous instances of membrane biogenesis and architecture, including membrane fusion and fission. Many stages of AP generation can be rationalised in terms of curvature, both the molecular geometry of lipids interpreted in terms of 'intrinsic curvature', and the overall mesoscopic curvature of the whole membrane, as observed with microscopy techniques. The present contribution intends to bring together the worlds of biophysics and cell biology of autophagy, in the hope that the resulting cross-pollination will generate abundant fruit.

Membrane signalosome: where biophysics meets systems biology

We opine on the recent advances in experiments and modeling of modular signaling complexes assembled on mammalian cell membranes (membrane signalosomes) in the context of several applications including intracellular trafficking, cell migration, and immune response. Characterizing the individual components of the membrane assemblies at the nanoscale, ranging from protein-lipid and protein-protein interactions, to membrane morphology, and the energetics of emergent assemblies at the subcellular to cellular scales pose significant challenges. Overcoming these challenges through the iterative coupling of multiscale modeling and experiment can be transformative in terms of addressing the gaps between structural biology and super-resolution microscopy, as it holds the key to the discovery of fundamental mechanisms behind the emergence of function in the membrane signalosome.

The degree and position of phosphorylation determine the impact of toxic and trace metals on phosphoinositide containing model membranes

This work assessed effects of metal binding on membrane fluidity, liposome size, and lateral organization in biomimetic membranes composed of 1 mol% of selected phosphorylated phosphoinositides in each system. Representative examples of phosphoinositide phosphate, bisphosphate and triphosphate were investigated. These include phosphatidylinositol-(4,5)-bisphosphate, an important signaling lipid constituting a minor component in plasma membranes whereas phosphatidylinositol-(4,5)-bisphosphate clusters support the propagation of secondary messengers in numerous signaling pathways. The high negative charge of phosphoinositides facilitates electrostatic interactions with metals. Lipids are increasingly identified as toxicological targets for divalent metals, which potentially alter lipid packing and domain formation. Exposure to heavy metals, such as lead and cadmium or elevated levels of essential metals, like cobalt, nickel, and manganese, implicated with various toxic effects were investigated.  Phosphatidylinositol-(4)-phosphate and phosphatidylinositol-(3,4,5)-triphosphate containing membranes are rigidified by lead, cobalt, and manganese whilst cadmium and nickel enhanced fluidity of membranes containing phosphatidylinositol-(4,5)-bisphosphate. Only cobalt induced liposome aggregation. All metals enhanced lipid clustering in phosphatidylinositol-(3,4,5)-triphosphate systems, cobalt in phosphatidylinositol-(4,5)-bisphosphate systems, while all metals showed limited changes in lateral film organization in phosphatidylinositol-(4)-phosphate matrices. These observed changes are relevant from the biophysical perspective as interference with the spatiotemporal formation of intricate domains composed of important signaling lipids may contribute to metal toxicity.

PLD2-PI(4,5)P2 interactions in fluid phase membranes: Structural modeling and molecular dynamics simulations

Interaction of phospholipase D2 (PLD2) with phosphatidylinositol (4,5)-bisphosphate (PIP2) is regarded as the critical step of numerous physiological processes. Here we build a full-length model of human PLD2 (hPLD2) combining template-based and ab initio modeling techniques and use microsecond all-atom molecular dynamics (MD) simulations of the protein in contact with a complex membrane to determine hPLD2-PIP2 interactions. MD simulations reveal that the intermolecular interactions preferentially occur between specific PIP2 phosphate groups and hPLD2 residues; the most strongly interacting residues are arginine at the pbox consensus sequence (PX) and pleckstrin homology (PH) domain. Interaction networks indicate formation of clusters at the protein-membrane interface consisting of amino acids, PIP2, and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidic acid (POPA); the largest cluster was in the PH domain.

Divalent cations bind to phosphoinositides to induce ion and isomer specific propensities for nano-cluster initiation in bilayer membranes

We report all-atom molecular dynamics simulations of asymmetric bilayers containing phosphoinositides in the presence of monovalent and divalent cations. We have characterized the molecular mechanism by which these divalent cations interact with phosphoinositides. Ca²⁺ desolvates more readily, consistent with single-molecule calculations, and forms a network of ionic-like bonds that serve as a 'molecular glue' that allows a single ion to coordinate with up to three phosphatidylinositol-(4,5)-bisphosphate (PI(4, 5)P₂) lipids. The phosphatidylinositol-(3,5)-bisphosphate isomer shows no such effect and neither does PI(4, 5)P₂ in the presence of Mg²⁺. The resulting network of Ca²⁺-mediated lipid-lipid bonds grows to span the entire simulation space and therefore has implications for the lateral distribution of phosophoinositides in the bilayer. We observe context-specific differences in lipid diffusion rates, lipid surface densities and bilayer structure. The molecular-scale delineation of ion-lipid arrangements reported here provides insight into similar nanocluster formation induced by peripheral proteins to regulate the formation of functional signalling complexes on the membrane.