Calcium Directly Regulates Phosphatidylinositol 4,5-Bisphosphate Headgroup Conformation and Recognition

Phosphoinositide switches in cell physiology - From molecular mechanisms to disease

Phosphoinositides are amphipathic lipid molecules derived from phosphatidylinositol that represent low abundance components of biological membranes. Rather than serving as mere structural elements of lipid bilayers, they represent molecular switches for a broad range of biological processes, including cell signaling, membrane dynamics and remodeling, and many other functions. Here, we focus on the molecular mechanisms that turn phosphoinositides into molecular switches and how the dysregulation of these processes can lead to disease.

In situ study of structural changes: Exploring the mechanism of protein corona transition from soft to hard

HYPOTHESIS

The process of protein corona changes has been widely believed to follow the Vroman effect, while protein structural change during the process is rarely reported, due to the lack of analytical methods. In-situ interpretation for protein structural change is critical to processes such as the recognition and transport of nanomaterials.

EXPERIMENTS

Molecular dynamics (MD) simulation was used to predict the deflection and twist of the protein tertiary structure. The structural changes of the surface protein corona during the interaction of nanoparticles (NPs) with lipid bilayer were probed in situ and real-time by sum frequency generation (SFG) spectroscopy.

FINDINGS

The ring tertiary structure of the protein corona is altered from vertical to horizontal on particle surface, a process of the soft-to-hard structural transition, which is contributed by the hydrogen bonding force between the protein and water molecules. The negatively charged protein corona can induce the redistribution of interfacial charge, leading to a more stable hydrogen bond network of the interfacial water. Our findings suggest that the structural change from flexible to rigid is a crucial process in the soft-to-hard transition of the protein corona, which will be a beneficial supplement to the Vroman effect of protein adsorption.

Deciphering Lipid Arrangement in Phosphatidylserine/Phosphatidylcholine Mixed Membranes: Simulations and Experiments

Phosphatidylserine (PS) exposure on the plasma membrane is crucial for many cellular processes including apoptotic cell recognition, blood clotting regulation, cellular signaling, and intercellular interactions. In this study, we investigated the arrangement of PS headgroups in mixed PS/phosphatidylcholine (PC) bilayers, serving as a simplified model of the outer leaflets of mammalian cell plasma membranes. Combining atomistic-scale molecular dynamics (MD) simulations with Langmuir monolayer experiments, we unraveled the mutual miscibility of POPC and POPS lipids and the intricate intermolecular interactions inherent to these membranes as well as the disparities in position and orientation of PC and PS headgroups. Our experiments revealed micrometer-scale miscibility at all mole fractions of POPC and POPS, marked by modest deviations from ideal mixing with no apparent microscale phase separation. The MD simulations, meanwhile, demonstrated that these deviations were due to strong electrostatic interactions between like-lipid pairs (POPC-POPC and POPS-POPS), culminating in lateral segregation and nanoscale clustering. Notably, PS headgroups profoundly affect the ordering of the lipid acyl chains, leading to lipid elongation and subtle PS protrusion above the zwitterionic membrane. In addition, PC headgroups are more tilted with respect to the membrane normal, while PS headgroups align at a smaller angle, making them more exposed to the surface of the mixed PC/PS membranes. These findings provide a detailed molecular-level account of the organization of mixed PC/PS membranes, corroborated by experimental data. The insights gained here extend our comprehension of the physiological role of PSs.

Solvation Behavior of Elastin-like Polypeptides in Divalent Metal Salt Solutions

The effects of CaCl2 and MgCl2 on the cloud point temperature of two different elastin-like polypeptides (ELPs) were studied using a combination of cloud point measurements, molecular dynamics simulations, and infrared spectroscopy. Changes in the cloud point for the ELPs in aqueous divalent metal cation solutions were primarily governed by two competing interactions: the cation-amide oxygen electrostatic interaction and the hydration of the cation. In particular, Ca2+ cations can more readily shed their hydration shells and directly contact two amide oxygens by the formation of ion bridges. By contrast, Mg2+ cations were more strongly hydrated and preferred to partition toward the amide oxygens along with their hydration shells. In fact, although hydrophilic ELP V5A2G3 was salted-out at low concentrations of MgCl2, it was salted-in at higher salt concentrations. By contrast, CaCl2 salted the ELP sharply out of solution at higher salt concentrations because of the bridging effect.

Efficient Simulations of Solvent Asymmetry Across Lipid Membranes Using Flat-Bottom Restraints

The routinely employed periodic boundary conditions complicate molecular simulations of physiologically relevant asymmetric lipid membranes together with their distinct solvent environments. Therefore, separating the extracellular fluid from its cytosolic counterpart has often been performed using a costly double-bilayer setup. Here, we demonstrate that the lipid membrane and solvent asymmetry can be efficiently modeled with a single lipid bilayer by applying an inverted flat-bottom potential to ions and other solute molecules, thereby restraining them to only interact with the relevant leaflet. We carefully optimized the parameters of the suggested method so that the results obtained using the flat-bottom and double-bilayer approaches become mutually indistinguishable. Then, we apply the flat-bottom approach to lipid bilayers with various compositions and solvent environments, covering ions and cationic peptides to validate the approach in a realistic use case. We also discuss the possible limitations of the method as well as its computational efficiency and provide a step-by-step guide on how to set up such simulations in a straightforward manner.

"Nonlinear" pursuit of understanding pollutant accumulation and chemistry at environmental and biological interfaces

Over the past few decades, the public recognition of the prevalence of certain classes of pollutants, such as perfluoroalkyl substances and nanoplastics, within the environment, has sparked growing concerns over their potential impact on environmental and human health. Within both environmental and biological systems, the adsorption and structural organization of pollutants at aqueous interfaces can greatly impact the chemical reactivity and transformation. Experimentally probing chemical behavior at interfaces can often pose a problem due to bulk solvated molecules convoluting molecular signatures from interfacial molecules. To solve this problem, there exist interface-specific nonlinear spectroscopy techniques that can directly probe both macroscopic planar interfaces and nanoplastic interfaces in aqueous environments. These techniques can provide essential information such as chemical adsorption, structure, and reactivity at interfaces. In this perspective, these techniques are presented with obvious advantages for studying the chemical properties of pollutants adsorbed to environmental and biological interfaces.

Boldine Alters Serum Lipidomic Signatures after Acute Spinal Cord Transection in Male Mice

Traumatic spinal cord injury (SCI) results in wide-ranging cellular and systemic dysfunction in the acute and chronic time frames after the injury. Chronic SCI has well-described secondary medical consequences while acute SCI has unique metabolic challenges as a result of physical trauma, in-patient recovery and other post-operative outcomes. Here, we used high resolution mass spectrometry approaches to describe the circulating lipidomic and metabolomic signatures using blood serum from mice 7 d after a complete SCI. Additionally, we probed whether the aporphine alkaloid, boldine, was able to prevent SCI-induced changes observed using these 'omics platforms'. We found that SCI resulted in large-scale changes to the circulating lipidome but minimal changes in the metabolome, with boldine able to reverse or attenuate SCI-induced changes in the abundance of 50 lipids. Multiomic integration using xMWAS demonstrated unique network structures and community memberships across the groups.

Amyloid-β Tetramers and Divalent Cations at the Membrane/Water Interface: Simple Models Support a Functional Role

Charge polarization at the membrane interface is a fundamental process in biology. Despite the lower concentration compared to the abundant monovalent ions, the relative abundance of divalent cations (Ca2+, Mg2+, Zn2+, Fe2+, Cu2+) in particular spaces, such as the neuron synapse, raised many questions on the possible effects of free multivalent ions and of the required protection of membranes by the eventual defects caused by the free forms of the cations. In this work, we first applied a recent realistic model of divalent cations to a well-investigated model of a polar lipid bilayer, di-myristoyl phosphatidyl choline (DMPC). The full atomistic model allows a fairly good description of changes in the hydration of charged and polar groups upon the association of cations to lipid atoms. The lipid-bound configurations were analyzed in detail. In parallel, amyloid-β 1-42 (Aβ42) peptides assembled into tetramers were modeled at the surface of the same bilayer. Two of the protein tetramers' models were loaded with four Cu2+ ions, the latter bound as in DMPC-free Aβ42 oligomers. The two Cu-bound models differ in the binding topology: one with each Cu ion binding each of the monomers in the tetramer; one with pairs of Cu ions linking two monomers into dimers, forming tetramers as dimers of dimers. The models here described provide hints on the possible role of Cu ions in synaptic plasticity and of Aβ42 oligomers in storing the same ions away from lipids. The release of structurally disordered peptides in the synapse can be a mechanism to recover ion homeostasis and lipid membranes from changes in the divalent cation concentration.

Phosphoinositides and intracellular calcium signaling: novel insights into phosphoinositides and calcium coupling as negative regulators of cellular signaling

Intracellular calcium (Ca2+) and phosphoinositides (PIPs) are crucial for regulating cellular activities such as metabolism and cell survival. Cells maintain precise intracellular Ca2+ and PIP levels via the actions of a complex system of Ca2+ channels, transporters, Ca2+ ATPases, and signaling effectors, including specific lipid kinases, phosphatases, and phospholipases. Recent research has shed light on the complex interplay between Ca2+ and PIP signaling, suggesting that elevated intracellular Ca2+ levels negatively regulate PIP signaling by inhibiting the membrane localization of PIP-binding proteins carrying specific domains, such as the pleckstrin homology (PH) and Ca2+-independent C2 domains. This dysregulation is often associated with cancer and metabolic diseases. PIPs recruit various proteins with PH domains to the plasma membrane in response to growth hormones, which activate signaling pathways regulating metabolism, cell survival, and growth. However, abnormal PIP signaling in cancer cells triggers consistent membrane localization and activation of PIP-binding proteins. In the context of obesity, an excessive intracellular Ca2+ level prevents the membrane localization of the PIP-binding proteins AKT, IRS1, and PLCδ via Ca2+-PIPs, contributing to insulin resistance and other metabolic diseases. Furthermore, an excessive intracellular Ca2+ level can cause functional defects in subcellular organelles such as the endoplasmic reticulum (ER), lysosomes, and mitochondria, causing metabolic diseases. This review explores how intracellular Ca2+ overload negatively regulates the membrane localization of PIP-binding proteins.

Engineering phosphatidylinositol-4,5-bisphosphate model membranes enriched in endocytic cargo: A neutron reflectometry, AFM and QCM-D structural study

The combination of in vitro models of biological membranes based on solid-supported lipid bilayers (SLBs) and of surface sensitive techniques, such as neutron reflectometry (NR), atomic force microscopy (AFM) and quartz crystal microbalance with dissipation monitoring (QCM-D), is well suited to provide quantitative information about molecular level interactions and lipid spatial distributions. In this work, cellular plasma membranes have been mimicked by designing complex SLB, containing phosphatidylinositol 4,5-bisphosphate (PtdIns4,5P2) lipids as well as incorporating synthetic lipo-peptides that simulate the cytoplasmic tails of transmembrane proteins. The QCM-D results revealed that the adsorption and fusion kinetics of PtdIns4,5P2 are highly dependent of Mg2+. Additionally, it was shown that increasing concentrations of PtdIns4,5P2 leads to the formation of SLBs with higher homogeneity. The presence of PtdIns4,5P2 clusters was visualized by AFM. NR provided important insights about the structural organization of the various components within the SLB, highlighting that the leaflet symmetry of these SLBs is broken by the presence of CD4-derived cargo peptides. Finally, we foresee our study to be a starting point for more sophisticated in vitro models of biological membranes with the incorporation of inositol phospholipids and synthetic endocytic motifs.

Spherocytosis-Related L1340P Mutation in Ankyrin Affects Its Interactions with Spectrin

Previously, we reported a new missense mutation in the ANK1 gene that correlated with the hereditary spherocytosis phenotype. This mutation, resulting in L1340P substitution (HGMD CM149731), likely leads to the changes in the conformation of the ankyrin ZZUD domain important for ankyrin binding to spectrin. Here, we report the molecular and physiological effects of this mutation. First, we assessed the binding activity of human β-spectrin to the mutated ZZUDL1340P domain of ankyrin using two different experimental approaches-the study of association and dissociation responses of the spectrin-ankyrin binding domain and a sedimentation assay. In addition, we documented the changes in morphology caused by the overexpressed ankyrin ZZUD domain in human cell models. Our results prove the key role of the L1340 aa residue for the correct alignment of the ZZUD domain of ankyrin, which results in binding the latter with spectrin within the erythrocyte membrane. Replacing L1340 with a proline residue disrupts the spectrin-binding activity of ankyrin.

Phosphatidylinositol 4,5-Bisphosphate and Cholesterol Regulators of the Calcium-Activated Chloride Channels TMEM16A and TMEM16B

Chloride fluxes through homo-dimeric calcium-activated channels TMEM16A and TMEM16B are critical to blood pressure, gastrointestinal motility, hormone, fluid and electrolyte secretion, pain sensation, sensory transduction, and neuronal and muscle excitability. Their gating depends on the voltage-dependent binding of two intracellular calcium ions to a high-affinity site formed by acidic residues from α-helices 6-8 in each monomer. Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2), a low-abundant lipid of the inner leaflet, supports TMEM16A function; it allows TMEM16A to evade the down-regulation induced by calcium, poly-L-lysine, or PI(4,5)P2 5-phosphatase. In stark contrast, adding or removing PI(4,5)P2 diminishes or increases TMEM16B function, respectively. PI(4,5)P2-binding sites on TMEM16A, and presumably on TMEM16B, are on the cytosolic side of α-helices 3-5, opposite the calcium-binding sites. This modular structure suggested that PI(4,5)P2 and calcium cooperate to maintain the conductive state in TMEM16A. Cholesterol, the second-largest constituent of the plasma membrane, also regulates TMEM16A though the mechanism, functional outcomes, binding site(s), and effects on TMEM16A and TMEM16B remain unknown.

SNARE Proteins in Synaptic Vesicle Fusion

Neurotransmitters are stored in small membrane-bound vesicles at synapses; a subset of synaptic vesicles is docked at release sites. Fusion of docked vesicles with the plasma membrane releases neurotransmitters. Membrane fusion at synapses, as well as all trafficking steps of the secretory pathway, is mediated by SNARE proteins. The SNAREs are the minimal fusion machinery. They zipper from N-termini to membrane-anchored C-termini to form a 4-helix bundle that forces the apposed membranes to fuse. At synapses, the SNAREs comprise a single helix from syntaxin and synaptobrevin; SNAP-25 contributes the other two helices to complete the bundle. Unc13 mediates synaptic vesicle docking and converts syntaxin into the permissive "open" configuration. The SM protein, Unc18, is required to initiate and proofread SNARE assembly. The SNAREs are then held in a half-zippered state by synaptotagmin and complexin. Calcium removes the synaptotagmin and complexin block, and the SNAREs drive vesicle fusion. After fusion, NSF and alpha-SNAP unwind the SNAREs and thereby recharge the system for further rounds of fusion. In this chapter, we will describe the discovery of the SNAREs, their relevant structural features, models for their function, and the central role of Unc18. In addition, we will touch upon the regulation of SNARE complex formation by Unc13, complexin, and synaptotagmin.

Unravelling the orientation of the inositol-biphosphate ring and its dependence on phosphatidylinositol 4,5-bisphosphate cluster formation in model membranes

HYPOTHESIS

Inositol phospholipids are well known to form clusters in the cytoplasmic leaflet of the plasma membrane that are responsible for the interaction and recruitment of proteins involved in key biological processes like endocytosis, ion channel activation and secondary messenger production. Although their phosphorylated inositol ring headgroup plays an important role in protein binding, its orientation with respect to the plane of the membrane and its lateral packing density has not been previously described experimentally.

EXPERIMENTS

Here, we study phosphatidylinositol 4,5-bisphosphate (PIP2) planar model membranes in the form of Langmuir monolayers by surface pressure-area isotherms, Brewster angle microscopy and neutron reflectometry to elucidate the relation between lateral (in-plane) and perpendicular (out-of-plane) molecular organization of PIP2.

FINDINGS

Different surface areas were explored through monolayer compression, allowing us to correlate the formation of transient PIP2 clusters with the change in orientation of the inositol-biphosphate headgroup, which was experimentally determined by neutron reflectometry.

Structural determinants of REMORIN nanodomain formation in anionic membranes

Remorins are a family of multigenic plasma membrane phosphoproteins involved in biotic and abiotic plant interaction mechanisms, partnering in molecular signaling cascades. Signaling activity of remorins depends on their phosphorylation states and subsequent clustering into nanosized membrane domains. The presence of a coiled-coil domain and a C-terminal domain is crucial to anchor remorins to negatively charged membrane domains; however, the exact role of the N-terminal intrinsically disordered domain (IDD) on protein clustering and lipid interactions is largely unknown. Here, we combine chemical biology and imaging approaches to study the partitioning of group 1 remorin into anionic model membranes mimicking the inner leaflet of the plant plasma membrane. Using reconstituted membranes containing a mix of saturated and unsaturated phosphatidylcholine, phosphatidylinositol phosphates, and sterol, we investigate the clustering of remorins to the membrane and monitor the formation of nanosized membrane domains. REM1.3 promoted membrane nanodomain organization on the exposed external leaflet of both spherical lipid vesicles and flat supported lipid bilayers. Our results reveal that REM1.3 drives a mechanism allowing lipid reorganization, leading to the formation of remorin-enriched nanodomains. Phosphorylation of the N-terminal IDD by the calcium protein kinase CPK3 influences this clustering and can lead to the formation of smaller and more disperse domains. Our work reveals the phosphate-dependent involvement of the N-terminal IDD in the remorin-membrane interaction process by driving structural rearrangements at lipid-water interfaces.

The PH Domain and C-Terminal polyD Motif of Phafin2 Exhibit a Unique Concurrence in Animals

Phafin2, a member of the Phafin family of proteins, contributes to a plethora of cellular activities including autophagy, endosomal cargo transportation, and macropinocytosis. The PH and FYVE domains of Phafin2 play key roles in membrane binding, whereas the C-terminal poly aspartic acid (polyD) motif specifically autoinhibits the PH domain binding to the membrane phosphatidylinositol 3-phosphate (PtdIns3P). Since the Phafin2 FYVE domain also binds PtdIns3P, the role of the polyD motif remains unclear. In this study, bioinformatics tools and resources were employed to determine the concurrence of the PH-FYVE module with the polyD motif among Phafin2 and PH-, FYVE-, or polyD-containing proteins from bacteria to humans. FYVE was found to be an ancient domain of Phafin2 and is related to proteins that are present in both prokaryotes and eukaryotes. Interestingly, the polyD motif only evolved in Phafin2 and PH- or both PH-FYVE-containing proteins in animals. PolyD motifs are absent in PH domain-free FYVE-containing proteins, which usually display cellular trafficking or autophagic functions. Moreover, the prediction of the Phafin2-interacting network indicates that Phafin2 primarily cross-talks with proteins involved in autophagy, protein trafficking, and neuronal function. Taken together, the concurrence of the polyD motif with the PH domain may be associated with complex cellular functions that evolved specifically in animals.

Impact of Ca2+-Induced PI(4,5)P2 Clusters on PH-YFP Organization and Protein-Protein Interactions

Despite its low abundance, phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) is a key modulator of membrane-associated signaling events in eukaryotic cells. Temporal and spatial regulation of PI(4,5)P2 concentration can achieve localized increases in the levels of this lipid, which are crucial for the activation or recruitment of peripheral proteins to the plasma membrane. The recent observation of the dramatic impact of physiological divalent cation concentrations on PI(4,5)P2 clustering, suggests that protein anchoring to the plasma membrane through PI(4,5)P2 is likely not defined solely by a simple (monomeric PI(4,5)P2)/(protein bound PI(4,5)P2) equilibrium, but instead depends on complex protein interactions with PI(4,5)P2 clusters. The insertion of PI(4,5)P2-binding proteins within these clusters can putatively modulate protein–protein interactions in the membrane, but the relevance of such effects is largely unknown. In this work, we characterized the impact of Ca2+ on the organization and protein–protein interactions of PI(4,5)P2-binding proteins. We show that, in giant unilamellar vesicles presenting PI(4,5)P2, the membrane diffusion properties of pleckstrin homology (PH) domains tagged with a yellow fluorescent protein (YFP) are affected by the presence of Ca2+, suggesting direct interactions between the protein and PI(4,5)P2 clusters. Importantly, PH-YFP is found to dimerize in the membrane in the absence of Ca2+. This oligomerization is inhibited in the presence of physiological concentrations of the divalent cation. These results confirm that cation-dependent PI(4,5)P2 clustering promotes interactions between PI(4,5)P2-binding proteins and has the potential to dramatically influence the organization and downstream interactions of PI(4,5)P2-binding proteins in the plasma membrane.

Phospholipid acyl tail affects lipid headgroup orientation and membrane hydration

Biomembrane hydration is crucial for understanding processes at biological interfaces. While the effect of the lipid headgroup has been studied extensively, the effect (if any) of the acyl chain chemical structure on lipid-bound interfacial water has remained elusive. We study model membranes composed of phosphatidylethanolamine (PE) and phosphatidylcholine (PC) lipids, the most abundant lipids in biomembranes. We explore the extent to which the lipid headgroup packing and associated water organization are affected by the lipid acyl tail unsaturation and chain length. To this end, we employ a combination of surface-sensitive techniques, including sum-frequency generation spectroscopy, surface pressure measurements, and Brewster angle microscopy imaging. Our results reveal that the acyl tail structure critically affects the headgroup phosphate orientational distribution and lipid-associated water molecules, for both PE and PC lipid monolayers at the air/water interface. These insights reveal the importance of acyl chain chemistry in determining not only membrane fluidity but also membrane hydration.

Contact Ion Pair Formation Is Not Necessarily Stronger than Solvent Shared Ion Pairing

Vibrational sum frequency spectroscopy (VSFS) and pressure-area Langmuir trough measurements were used to investigate the binding of alkali metal cations to eicosyl sulfate (ESO₄) surfactants in monolayers at the air/water interface. The number density of sulfate groups could be tuned by mixing the anionic surfactant with eicosanol. The equilibrium dissociation constant for K⁺ to the fatty sulfate interface showed 10 times greater affinity than for Li⁺ and approximately 3 times greater than for Na⁺. All three cations formed solvent shared ion pairs when the mole fraction of ESO₄ was 0.33 or lower. Above this threshold charge density, Li⁺ formed contact ion pairs with the sulfate headgroups, presumably via bridging structures. By contrast, K⁺ only bound to the sulfate moieties in solvent shared ion pairing configurations. The behavior for Na⁺ was intermediate. These results demonstrate that there is not necessarily a correlation between contact ion pair formation and stronger binding affinity.

Examining the Effect of Charged Lipids on Mitochondrial Outer Membrane Dynamics Using Atomistic Simulations

The outer mitochondrial membrane (OMM) is involved in multiple cellular functions such as apoptosis, inflammation and signaling via its membrane-associated and -embedded proteins. Despite the central role of the OMM in these vital phenomena, the structure and dynamics of the membrane have regularly been investigated in silico using simple two-component models. Accordingly, the aim was to generate the realistic multi-component model of the OMM and inspect its properties using atomistic molecular dynamics (MD) simulations. All major lipid components, phosphatidylinositol (PI), phosphatidylcholine (PC), phosphatidylethanolamine (PE), and phosphatidylserine (PS), were included in the probed OMM models. Because increased levels of anionic PS lipids have potential effects on schizophrenia and, more specifically, on monoamine oxidase B enzyme activity, the effect of varying the PS concentration was explored. The MD simulations indicate that the complex membrane lipid composition (MLC) behavior is notably different from the two-component PC-PE model. The MLC changes caused relatively minor effects on the membrane structural properties such as membrane thickness or area per lipid; however, notable effects could be seen with the dynamical parameters at the water-membrane interface. Increase of PS levels appears to slow down lateral diffusion of all lipids and, in general, the presence of anionic lipids reduced hydration and slowed down the PE headgroup rotation. In addition, sodium ions could neutralize the membrane surface, when PI was the main anionic component; however, a similar effect was not seen for high PS levels. Based on these results, it is advisable for future studies on the OMM and its protein or ligand partners, especially when wanting to replicate the correct properties on the water-membrane interface, to use models that are sufficiently complex, containing anionic lipid types, PI in particular.

Recapitulation of cell-like KRAS4b membrane dynamics on complex biomimetic membranes

Understanding the spatiotemporal distribution and dynamics of RAS on the plasma membrane (PM) is the key for elucidating the molecular mechanisms of the RAS signaling pathway. Single particle tracking (SPT) experiments show that in cells, KRAS diffuses in at least three interchanging states on the cellular PM; however, KRAS remains monomeric and always shows homogeneous diffusion on artificial membranes. Here, we show for the first time on a supported lipid bilayer composed of heterogeneous lipid components that we can recapitulate the three-state diffusion of KRAS seen in cells. The use of a biologically relevant eight-lipid system opens a new frontier in the biophysical studies of RAS and other membrane associated proteins on a biomimetic system that recapitulates the complexity of a cellular PM.

•

KRAS4b shows homogeneous diffusion on simple 2-lipids bilayer

•

KRAS4b shows a cell-like, three-state diffusion on a complex 8-lipid bilayer

•

Phase separation in lipids favors the multi-state diffusion of KRAS4b

•

The complex lipid composition favors RAS nanoclustering irrespective of nucleotide state

Ca2+ entry through mechanotransduction channels localizes BAIAP2L2 to stereocilia tips

Brain-specific angiogenesis inhibitor 1-associated protein 2-like protein 2 (BAIAP2L2), a membrane-binding protein required for the maintenance of mechanotransduction in hair cells, is selectively retained at the tips of transducing stereocilia. BAIAP2L2 trafficked to stereocilia tips in the absence of EPS8, but EPS8 increased the efficiency of localization. A tripartite complex of BAIAP2L2, EPS8, and MYO15A formed efficiently in vitro, and these three proteins robustly targeted to filopodia tips when coexpressed in cultured cells. Mice lacking functional transduction channels no longer concentrated BAIAP2L2 at row 2 stereocilia tips, a result that was phenocopied by blocking channels with tubocurarine in cochlear explants. Transduction channels permit Ca²⁺ entry into stereocilia, and we found that membrane localization of BAIAP2L2 was enhanced in the presence of Ca²⁺. Finally, reduction of intracellular Ca²⁺ in hair cells using BAPTA-AM led to a loss of BAIAP2L2 at stereocilia tips. Taken together, our results show that a MYO15A-EPS8 complex transports BAIAP2L2 to stereocilia tips, and Ca²⁺ entry through open channels at row 2 tips retains BAIAP2L2 there.

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.

Simulations of Kindlin-2 PIP binding domains reveal protonation-dependent membrane binding modes

Kindlin-2, a member of the Kindlin family of peripheral membrane proteins, is important for integrin activation and stabilization of epidermal growth factor receptor. It associates with the cytoplasmic face of the plasma membrane via dedicated phosphatidylinositol phosphate binding domains located in the N-terminal F0 and Pleckstrin Homology domains. These domains have binding affinity for phosphatidylinositol 4,5-bisphosphate and, to a greater degree, phosphatidylinositol 3,4,5-trisphosphate. The biological significance of the differential binding of these phosphatidylinositol phosphates to Kindlin-2 and the mechanism by which they activate Kindlin-2 are not well understood. Recently, ssNMR identified the predominant protonation states of phosphatidylinositol 4,5-bisphosphate and phosphatidylinositol 3,4,5-trisphosphate near physiological pH in the presence of anionic lipids. Here, we perform atomistic simulation of the bound state of the Pleckstrin Homology and F0 domains of Kindlin-2 at membranes containing phosphatidylinositol 4,5-bisphosphate/phosphatidylinositol 3,4,5-trisphosphate with differing protonation states. This computational approach demonstrates that these two phosphatidylinositol phosphates differently modulate Kindlin-2 subdomain binding in a protonation-state-dependent manner. We speculate these variations in binding mode provide a mechanism for intracellular pH and Ca²⁺ influx to control the membrane binding behavior and activity of Kindlin-2.

Surfactant-induced wettability reversal on oil-wet calcite surfaces: Experimentation and molecular dynamics simulations with scaled-charges

HYPOTHESIS

Surfactant flooding is the leading approach for reversing the wettability of oil-wet carbonate reservoirs, which is critical for the recovery of the remaining oil. Combination of molecular dynamics (MD) simulations with experiments on simplified model systems can uncover the molecular mechanisms of wettability reversal and identify key molecular properties for systematic design of new, effective chemical formulations for the enhanced oil recovery.

EXPERIMENTS/SIMULATIONS

Wettability reversal by a series of surfactant solutions was studied experimentally using contact angle measurements on aged calcite chips, and a novel MD simulation methodology with scaled-charges that provides superior description of the ionic interactions in aqueous solutions.

FINDINGS

The MD simulation results were in excellent agreement with the experiments. Cationic surfactants were the most effective in reversing the calcite wettability, resulting in complete detachment of the oil from the surface. Some nonionic surfactants also altered the wettability, but to a lesser degree, while the amphoteric and anionic surfactants had no effect. From the tested cationic surfactants, the double-tailed one was the least effective, but the experiments were inconclusive due to its poor solubility. Contributions of specific interactions to the wettability reversal process and implications for the design and optimization of surfactants for the enhanced oil recovery are discussed.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Regulation of Store-Operated Ca2+ Entry by SARAF

Calcium (Ca²⁺) signaling plays a dichotomous role in cellular biology, controlling cell survival and proliferation on the one hand and cellular toxicity and cell death on the other. Store-operated Ca²⁺ entry (SOCE) by CRAC channels represents a major pathway for Ca²⁺ entry in non-excitable cells. The CRAC channel has two key components, the endoplasmic reticulum Ca²⁺ sensor stromal interaction molecule (STIM) and the plasma-membrane Ca²⁺ channel Orai. Physical coupling between STIM and Orai opens the CRAC channel and the resulting Ca²⁺ flux is regulated by a negative feedback mechanism of slow Ca²⁺ dependent inactivation (SCDI). The identification of the SOCE-associated regulatory factor (SARAF) and investigations of its role in SCDI have led to new functional and molecular insights into how SOCE is controlled. In this review, we provide an overview of the functional and molecular mechanisms underlying SCDI and discuss how the interaction between SARAF, STIM1, and Orai1 shapes Ca²⁺ signaling in cells.

Binding of Ca2+-independent C2 domains to lipid membranes: A multi-scale molecular dynamics study

C2 domains facilitate protein interactions with lipid bilayers in either a Ca²⁺-dependent or -independent manner. We used molecular dynamics (MD) simulations to explore six Ca²⁺-independent C2 domains, from KIBRA, PI3KC2α, RIM2, PTEN, SHIP2, and Smurf2. In coarse-grained MD simulations these C2 domains formed transient interactions with zwitterionic bilayers, compared with longer-lived interactions with anionic bilayers containing phosphatidylinositol bisphosphate (PIP₂). Type I C2 domains bound non-canonically via the front, back, or side of the β sandwich, whereas type II C2 domains bound canonically, via the top loops. C2 domains interacted strongly with membranes containing PIP₂, causing bound anionic lipids to cluster around the protein. Binding modes were refined via atomistic simulations. For PTEN and SHIP2, CG simulations of their phosphatase plus C2 domains with PIP₂-containing bilayers were also performed, and the roles of the two domains in membrane localization compared. These studies establish a simulation protocol for membrane-recognition proteins.

•

Binding of Ca²⁺-independent C2 domains to membranes was explored by MD simulation

•

C2 domains from KIBRA, PI3KC2α, RIM2, PTEN, SHIP2, and Smurf2 were compared

•

C2 domains formed longer-lived interactions with lipid bilayers containing PIP₂

•

For PTEN and SHIP2, simulations of their phosphatase plus C2 domains were performed

Distinct EH domains of the endocytic TPLATE complex confer lipid and protein binding

Clathrin-mediated endocytosis (CME) is the gatekeeper of the plasma membrane. In contrast to animals and yeasts, CME in plants depends on the TPLATE complex (TPC), an evolutionary ancient adaptor complex. However, the mechanistic contribution of the individual TPC subunits to plant CME remains elusive. In this study, we used a multidisciplinary approach to elucidate the structural and functional roles of the evolutionary conserved N-terminal Eps15 homology (EH) domains of the TPC subunit AtEH1/Pan1. By integrating high-resolution structural information obtained by X-ray crystallography and NMR spectroscopy with all-atom molecular dynamics simulations, we provide structural insight into the function of both EH domains. Both domains bind phosphatidic acid with a different strength, and only the second domain binds phosphatidylinositol 4,5-bisphosphate. Unbiased peptidome profiling by mass-spectrometry revealed that the first EH domain preferentially interacts with the double N-terminal NPF motif of a previously unidentified TPC interactor, the integral membrane protein Secretory Carrier Membrane Protein 5 (SCAMP5). Furthermore, we show that AtEH/Pan1 proteins control the internalization of SCAMP5 via this double NPF peptide interaction motif. Collectively, our structural and functional studies reveal distinct but complementary roles of the EH domains of AtEH/Pan1 in plant CME and connect the internalization of SCAMP5 to the TPLATE complex.

Transmembrane potential of physiologically relevant model membranes: Effects of membrane asymmetry

Transmembrane potential difference (V_m) plays important roles in regulating various biological processes. At the macro level, V_m can be experimentally measured or calculated using the Nernst or Goldman-Hodgkin-Katz equation. However, the atomic details responsible for its generation and impact on protein and lipid dynamics still need to be further elucidated. In this work, we performed a series of all-atom molecular dynamics (MD) simulations of symmetric model membranes of various lipid compositions and cation contents to evaluate the relationship between membrane asymmetry and V_m. Specifically, we studied the impact of the asymmetric distribution of POPS (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine), PIP2 (phosphatidylinositol 4,5-bisphosphate), as well as Na⁺ and K⁺ on V_m using atomically detailed MD simulations of symmetric model membranes. The results suggest that, for an asymmetric POPC-POPC/POPS bilayer in the presence of NaCl, the presence of the monovalent anionic lipid POPS in the inner leaflet polarizes the membrane (ΔV_m < 0). Intriguingly, replacing a third of the POPS lipids by the polyvalent anionic signaling lipid PIP2 counteracts this effect, resulting in a smaller negative membrane potential. We also found that replacing Na⁺ ions in the inner region by K⁺ depolarizes the membrane (ΔV_m > 0). These divergent effects arise from variations in the strength of cation-lipid interactions and are correlated with changes in lipid chain order and head-group orientation.

α-Synuclein plasma membrane localization correlates with cellular phosphatidylinositol polyphosphate levels

The Parkinson's disease protein α-synuclein (αSyn) promotes membrane fusion and fission by interacting with various negatively charged phospholipids. Despite postulated roles in endocytosis and exocytosis, plasma membrane (PM) interactions of αSyn are poorly understood. Here, we show that phosphatidylinositol 4,5-bisphosphate (PIP₂) and phosphatidylinositol 3,4,5-trisphosphate (PIP₃), two highly acidic components of inner PM leaflets, mediate PM localization of endogenous pools of αSyn in A2780, HeLa, SK-MEL-2, and differentiated and undifferentiated neuronal SH-SY5Y cells. We demonstrate that αSyn binds to reconstituted PIP₂ membranes in a helical conformation in vitro and that PIP₂ synthesizing kinases and hydrolyzing phosphatases reversibly redistribute αSyn in cells. We further delineate that αSyn-PM targeting follows phosphoinositide-3 kinase (PI3K)-dependent changes of cellular PIP₂ and PIP₃ levels, which collectively suggests that phosphatidylinositol polyphosphates contribute to αSyn's function(s) at the plasma membrane.

Role of Water in CaCO3 Biomineralization

Biomineralization occurs in aqueous environments. Despite the ubiquity and relevance of CaCO₃ biomineralization, the role of water in the biomineralization process has remained elusive. Here, we demonstrate that water reorganization accompanies CaCO₃ biomineralization for sea urchin spine generation in a model system. Using surface-specific vibrational spectroscopy, we probe the water at the interface of the spine-associated protein during CaCO₃ mineralization. Our results show that, while the protein structure remains unchanged, the structure of interfacial water is perturbed differently in the presence of both Ca²⁺ and CO₃^2– compared to the addition of only Ca²⁺. This difference is attributed to the condensation of prenucleation mineral species. Our findings are consistent with a nonclassical mineralization pathway for sea urchin spine generation and highlight the importance of protein hydration in biomineralization.

PI(4,5)P2 Clustering and Its Impact on Biological Functions

Located at the inner leaflet of the plasma membrane (PM), phosphatidyl-inositol 4,5-bisphosphate [PI(4,5)P₂] composes only 1-2 mol% of total PM lipids. With its synthesis and turnover both spatially and temporally regulated, PI(4,5)P₂ recruits and interacts with hundreds of cellular proteins to support a broad spectrum of cellular functions. Several factors contribute to the versatile and dynamic distribution of PI(4,5)P₂ in membranes. Physiological multivalent cations such as Ca²⁺ and Mg²⁺ can bridge between PI(4,5)P₂ headgroups, forming nanoscopic PI(4,5)P₂-cation clusters. The distinct lipid environment surrounding PI(4,5)P₂ affects the degree of PI(4,5)P₂ clustering. In addition, diverse cellular proteins interacting with PI(4,5)P₂ can further regulate PI(4,5)P₂ lateral distribution and accessibility. This review summarizes the current understanding of PI(4,5)P₂ behavior in both cells and model membranes, with emphasis on both multivalent cation- and protein-induced PI(4,5)P₂ clustering. Understanding the nature of spatially separated pools of PI(4,5)P₂ is fundamental to cell biology.

Structural Effects of Cation Binding to DPPC Monolayers

Ions at the two sides of the plasma membrane maintain the transmembrane potential, participate in signaling, and affect the properties of the membrane itself. The extracellular leaflet is particularly enriched in phosphatidylcholine lipids and under the influence of Na⁺, Ca²⁺, and Cl^- ions. In this work, we combined molecular dynamics simulations performed using state-of-the-art models with vibrational sum frequency generation (VSFG) spectroscopy to study the effects of these key ions on the structure of dipalmitoylphosphatidylcholine. We used lipid monolayers as a proxy for membranes, as this approach enabled a direct comparison between simulation and experiment. We find that the effects of Na⁺ are minor. Ca²⁺, on the other hand, strongly affects the lipid headgroup conformations and induces a tighter packing of lipids, thus promoting the liquid condensed phase. It does so by binding to both the phosphate and carbonyl oxygens via direct and water-mediated binding modes, the ratios of which depend on the monolayer packing. Clustering analysis performed on simulation data revealed that changes in area per lipid or CaCl₂ concentration both affect the headgroup conformations, yet their effects are anticorrelated. Cations at the monolayer surface also attract Cl^-, which at large CaCl₂ concentrations penetrates deep to the monolayer. This phenomenon coincides with a radical change in the VSFG spectra of the phosphate group, thus indicating the emergence of a new binding mode.

Zn2+ Binds to Phosphatidylserine and Induces Membrane Blebbing

Herein, we show that Zn²⁺ binds to phosphatidylserine (PS) lipids in supported lipid bilayers (SLBs), forming a PS-Zn²⁺ complex with an equilibrium dissociation constant of ∼100 μM. Significantly, Zn²⁺ binding to SLBs containing more than 10 mol % PS induces extensive reordering of the bilayer. This reordering is manifest through bright spots of high fluorescence intensity that can be observed when the bilayer contains a dye-labeled lipid. Measurements using atomic force microscopy (AFM) reveal that these spots represent three-dimensional unilamellar blebs. Bleb formation is ion specific, inducible by exposing the bilayer to μM concentrations of Zn²⁺ but not Mg²⁺, Cu²⁺, Co²⁺, or Mn²⁺. Moreover, Ca²⁺ can induce some blebbing at mM concentrations but not nearly as effectively as Zn²⁺. The interactions of divalent metal cations with PS lipids were further investigated by a combination of vibrational sum frequency spectroscopy (VSFS) and surface pressure-area isotherm measurements. VSFS revealed that Zn²⁺ and Ca²⁺ were bound to the phosphate and carboxylate moieties on PS via contact ion pairing, dehydrating the lipid headgroup, whereas Mg²⁺ and Cu²⁺ were bound without perturbing the hydration of these functional groups. Additionally, Zn²⁺ was found to dramatically reduce the area per lipid in lipid monolayers, while Mg²⁺ and Cu²⁺ did not. Ca²⁺ could also reduce the area per lipid but only when significantly higher surface pressures were applied. These measurements suggest that Zn²⁺ caused lipid blebbing by decreasing the area per lipid on the side of the bilayer to which the salt was exposed. Such findings have implications for blebbing, fusion, oxidation, and related properties of PS-rich membranes in biological systems where Zn²⁺ concentrations are asymmetrically distributed.

Microscopic Picture of Calcium-Assisted Lipid Demixing and Membrane Remodeling Using Multiscale Simulations

The specificity of anionic phospholipids-calcium ion interaction and lipid demixing has been established as a key regulatory mechanism in several cellular signaling processes. The mechanism and implications of this calcium-assisted demixing have not been elucidated from a microscopic point of view. Here, we present an overview of atomic interactions between calcium and phospholipids that can drive nonideal mixing of lipid molecules in a model lipid bilayer composed of zwitterionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)) and anionic (1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS)) lipids with computer simulations at multiple resolutions. Lipid nanodomain formation and growth were driven by calcium-enabled lipid bridging of the charged phosphatidylserine (PS) headgroups, which were favored against inter-POPS dipole interactions. Consistent with several experimental studies of calcium-associated membrane sculpting, our analyses also suggest modifications in local membrane curvature and cross-leaflet couplings as a response to such induced lateral heterogeneity. In addition, reverse mapping to a complementary atomistic description revealed structural insights in the presence of anionic nanodomains, at timescales not accessed by previous computational studies. This work bridges information across multiple scales to reveal a mechanistic picture of calcium ion's impact on membrane biophysics.

Structure and Lateral Organization of Phosphatidylinositol 4,5-bisphosphate

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂) is a minor but ubiquitous component of the inner leaflet of the plasma membrane of eukaryotic cells. However, due to its particular complex biophysical properties, it stands out from its neighboring lipids as one of the most important regulators of membrane-associated signaling events. Despite its very low steady-state concentration, PI(4,5)P₂ is able to engage in a multitude of simultaneous cellular functions that are temporally and spatially regulated through the presence of localized transient pools of PI(4,5)P₂ in the membrane. These pools are crucial for the recruitment, activation, and organization of signaling proteins and consequent regulation of downstream signaling. The present review showcases some of the most important PI(4,5)P2 molecular and biophysical properties as well as their impact on its membrane dynamics, lateral organization, and interactions with other biochemical partners.

Mechanisms of PI(4,5)P2 Enrichment in HIV-1 Viral Membranes

Phosphatidylinositol 4,5-bisphosphate (PIP2) is critical for HIV-1 virus assembly. The viral membrane is enriched in PIP2, suggesting that the virus assembles at PIP2-rich microdomains. We showed previously that in model membranes PIP2 can form nanoscopic clusters bridged by multivalent cations. Here, using purified proteins we quantitated the binding of HIV-1 Gag-related proteins to giant unilamellar vesicles containing either clustered or free PIP2. Myristoylated MA strongly preferred binding to clustered PIP2. By contrast, unmyristoylated HIV-1 MA, RSV MA, and a PH domain all preferred to interact with free PIP2. We also found that HIV-1 Gag multimerization promotes PIP2 clustering. Truncated Gag proteins comprising the MA, CA, and SP domains (MACASP) or the MA and CA domains (MACA) induced self-quenching of acyl chain-labeled fluorescent PIP2 in liposomes, implying clustering. However, HIV-1 MA itself did not induce PIP2 clustering. A CA inter-hexamer dimer interface mutation led to a loss of induced PIP2 clustering in MACA, indicating the importance of protein multimerization. Cryo-electron tomography of liposomes with bound MACA showed an amorphous protein layer on the membrane surface. Thus, it appears that while protein–protein interactions are required for PIP2 clustering, formation of a regular lattice is not. Protein-induced PIP2 clustering and multivalent cation-induced PIP2 clustering are additive. Taken together, these results provide the first evidence that HIV-1 Gag can selectively target pre-existing PIP2-enriched domains of the plasma membrane for viral assembly, and that Gag multimerization can further enrich PIP2 at assembly sites. These effects could explain the observed PIP2 enrichment in HIV-1.

The Glycosphingolipid GM3 Modulates Conformational Dynamics of the Glucagon Receptor

The extracellular domain (ECD) of class B1 G-protein-coupled receptors (GPCRs) plays a central role in signal transduction and is uniquely positioned to sense both the extracellular and membrane environments. Although recent studies suggest a role for membrane lipids in the modulation of class A and class F GPCR signaling properties, little is known about the effect of lipids on class B1 receptors. In this study, we employed multiscale molecular dynamics simulations to access the dynamics of the glucagon receptor (GCGR) ECD in the presence of native-like membrane bilayers. Simulations showed that the ECD could move about a hinge region formed by residues Q122–E126 to adopt both closed and open conformations relative to the transmembrane domain. ECD movements were modulated by binding of the glycosphingolipid GM3. These large-scale fluctuations in ECD conformation may affect the ligand binding and receptor activation properties. We also identify a unique phosphatidylinositol (4,5)-bisphosphate (PIP₂) interaction profile near intracellular loop (ICL) 2/TM3 at the G-protein-coupling interface, suggesting a mechanism of engaging G-proteins that may have a distinct dependence on PIP₂ compared with class A GPCRs. Given the structural conservation of class B1 GPCRs, the modulatory effects of GM3 and PIP₂ on GCGR may be conserved across these receptors, offering new insights into potential therapeutic targeting.

New opportunities and insights into Papaver self-incompatibility by imaging engineered Arabidopsis pollen

Pollen tube growth is essential for plant reproduction. Their rapid extension using polarized tip growth provides an exciting system for studying this specialized type of growth. Self-incompatibility (SI) is a genetically controlled mechanism to prevent self-fertilization. Mechanistically, one of the best-studied SI systems is that of Papaver rhoeas (poppy). This utilizes two S-determinants: stigma-expressed PrsS and pollen-expressed PrpS. Interaction of cognate PrpS-PrsS triggers a signalling network, causing rapid growth arrest and programmed cell death (PCD) in incompatible pollen. We previously demonstrated that transgenic Arabidopsis thaliana pollen expressing PrpS-green fluorescent protein (GFP) can respond to Papaver PrsS with remarkably similar responses to those observed in incompatible Papaver pollen. Here we describe recent advances using these transgenic plants combined with genetically encoded fluorescent probes to monitor SI-induced cellular alterations, including cytosolic calcium, pH, the actin cytoskeleton, clathrin-mediated endocytosis (CME), and the vacuole. This approach has allowed us to study the SI response in depth, using multiparameter live-cell imaging approaches that were not possible in Papaver. This lays the foundations for new opportunities to elucidate key mechanisms involved in SI. Here we establish that CME is disrupted in self-incompatible pollen. Moreover, we reveal new detailed information about F-actin remodelling in pollen tubes after SI.

The energetics of protein-lipid interactions as viewed by molecular simulations

Membranes are formed from a bilayer containing diverse lipid species with which membrane proteins interact. Integral, membrane proteins are embedded in this bilayer, where they interact with lipids from their surroundings, whilst peripheral membrane proteins bind to lipids at the surface of membranes. Lipid interactions can influence the function of membrane proteins, either directly or allosterically. Both experimental (structural) and computational approaches can reveal lipid binding sites on membrane proteins. It is, therefore, important to understand the free energies of these interactions. This affords a more complete view of the engagement of a particular protein with the biological membrane surrounding it. Here, we describe many computational approaches currently in use for this purpose, including recent advances using both free energy and unbiased simulation methods. In particular, we focus on interactions of integral membrane proteins with cholesterol, and with anionic lipids such as phosphatidylinositol 4,5-bis-phosphate and cardiolipin. Peripheral membrane proteins are exemplified via interactions of PH domains with phosphoinositide-containing membranes. We summarise the current state of the field and provide an outlook on likely future directions of investigation.

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

Numerous cellular functions mediated by phosphatidylinositol (4,5)-bisphosphate (PI(4,5)P2; PIP2) involve clustering of the lipid as well as colocalization with other lipids. Although the cation-mediated electrostatic interaction is regarded as the primary clustering mechanism, the ion-specific nature of the intermolecular network formation makes it challenging to characterize the clusters. Here we use all-atom molecular dynamics (MD) simulations of PIP2 monolayers and graph-theoretic analysis to gain insight into the phenomenon. MD simulations reveal that the intermolecular interactions preferentially occur between specific cations and phosphate groups (P1, P4, and P5) of the inositol headgroup with better-matched kosmotropic/chaotropic characters consistent with the law of matching water affinities (LMWA). Ca2+ is strongly attracted to P4/P5, while K+ preferentially binds to P1; Na+ interacts with both P4/P5 and P1. These specific interactions lead to the characteristic clustering patterns. Specificially, the size distributions and structures of PIP2 clusters generated by kosmotropic cations Ca2+ and Na+ are bimodal, with a combination of small and large clusters, while there is little clustering in the presence of only chaotropic K+; the largest clusters are obtained in systems with all three cations. The small-world network (a model with both local and long-range connections) best characterizes the clusters, followed by the random and the scale-free networks. More generally, the present results interpreted within the LMWA are consistent with the relative eukaryotic intracellular concentrations Ca2+ ≪ Na+ < Mg2+ < K+; that is, concentrations of Ca2+ and Na+ must be low to prevent damaging aggregation of lipids, DNA, RNA and phosphate-containing proteins.

Partial Metal Ion Saturation of C2 Domains Primes Synaptotagmin 1-Membrane Interactions

Synaptotagmin 1 (Syt1) is an integral membrane protein whose phospholipid-binding tandem C2 domains, C2A and C2B, act as Ca²⁺ sensors of neurotransmitter release. Our objective was to understand the role of individual metal-ion binding sites of these domains in the membrane association process. We used Pb²⁺, a structural and functional surrogate of Ca²⁺, to generate the protein states with well-defined protein-metal ion stoichiometry. NMR experiments revealed that binding of one divalent metal ion per C2 domain results in loss of conformational plasticity of the loop regions, potentially pre-organizing them for additional metal-ion and membrane-binding events. In C2A, a divalent metal ion in site 1 is sufficient to drive its weak association with phosphatidylserine-containing membranes, whereas in C2B, it enhances the interactions with the signaling lipid phosphatidylinositol-4,5-bisphosphate. In full-length Syt1, both Pb²⁺-complexed C2 domains associate with phosphatidylserine-containing membranes. Electron paramagnetic resonance experiments show that the extent of membrane insertion correlates with the occupancy of the C2 metal ion sites. Together, our results indicate that upon partial metal ion saturation of the intra-loop region, Syt1 adopts a dynamic, partially membrane-bound state. The properties of this state, such as conformationally restricted loop regions and positioning of C2 domains in close proximity to anionic lipid headgroups, "prime" Syt1 for cooperative binding of a full complement of metal ions and deeper membrane insertion.

Multiple lipid binding sites determine the affinity of PH domains for phosphoinositide-containing membranes

Association of peripheral proteins with lipid bilayers regulates membrane signaling and dynamics. Pleckstrin homology (PH) domains bind to phosphatidylinositol phosphate (PIP) molecules in membranes. The effects of local PIP enrichment on the interaction of PH domains with membranes is unclear. Molecular dynamics simulations allow estimation of the binding energy of GRP1 PH domain to PIP3-containing membranes. The free energy of interaction of the PH domain with more than two PIP3 molecules is comparable to experimental values, suggesting that PH domain binding involves local clustering of PIP molecules within membranes. We describe a mechanism of PH binding proceeding via an encounter state to two bound states which differ in the orientation of the protein relative to the membrane, these orientations depending on the local PIP concentration. These results suggest that nanoscale clustering of PIP molecules can control the strength and orientation of PH domain interaction in a concentration-dependent manner.

Phospholipase C

Phospholipase C (PLC) family members constitute a family of diverse enzymes. Thirteen different family members have been cloned. These family members have unique structures that mediate various functions. Although PLC family members all appear to signal through the bi-products of cleaving phospholipids, it is clear that each family member, and at times each isoform, contributes to unique cellular functions. This chapter provides a review of the current literature on PLC. In addition, references have been provided for more in-depth information regarding areas that are not discussed including tyrosine kinase activation of PLC. Understanding the roles of the individual PLC enzymes, and their distinct cellular functions, will lead to a better understanding of the physiological roles of these enzymes in the development of diseases and the maintenance of homeostasis.