Segregation of negatively charged phospholipids by the polycationic and farnesylated membrane anchor of Kras

From disorder comes function: Regulation of small GTPase function by intrinsically disordered lipidated membrane anchor

The intrinsically disordered, lipid-modified membrane anchor of small GTPases is emerging as a critical modulator of function through its ability to sort lipids in a conformation-dependent manner. We reviewed recent computational and experimental studies that have begun to shed light on the sequence-ensemble-function relationship in this unique class of lipidated intrinsically disordered regions (LIDRs).

Intrinsically Disordered Membrane Anchors of Rheb, RhoA, and DiRas3 Small GTPases: Molecular Dynamics, Membrane Organization, and Interactions

Protein structure has been well established to play a key role in determining function; however, intrinsically disordered proteins and regions (IDPs and IDRs) defy this paradigm. IDPs and IDRs exist as an ensemble of structures rather than a stable 3D structure yet play essential roles in many cell-signaling processes. Nearly all Ras superfamily GTPases are tethered to membranes by a lipid tail at the end of a flexible IDR. The sequence of the IDR is a key determinant of membrane localization, and interaction between the IDR and the membrane has been shown to affect signaling in RAS proteins through the modulation of dynamic membrane organization. Here, we utilized atomistic molecular dynamics simulations to study the membrane interaction, conformational dynamics, and lipid sorting of three IDRs from small GTPases Rheb, RhoA, and DiRas3 in model membranes representing their physiological target membranes. We found that complementarity between the lipidated IDR sequence and target membrane lipid composition is a determinant of conformational plasticity. We also show that electrostatic interactions between anionic lipids and basic residues on IDRs are correlated with sampling of semistable conformational substates, and lack of these interactions is associated with greater conformational diversity. Finally, we show that small GTPase IDRs with a polybasic domain alter local lipid composition by segregating anionic lipids and, in some cases, excluding other lipids from their immediate vicinity in favor of anionic lipids.

Lysophosphatidylcholine acyltransferase 1 suppresses nanoclustering and function of KRAS

KRAS is frequently mutated in cancer, contributing to 20% of all human cancer especially pancreatic, colorectal and lung cancer. Signaling of the constitutively active KRAS oncogenic mutants is mostly compartmentalized to proteolipid nanoclusters on the plasma membrane (PM). Signaling nanoclusters of many KRAS mutants selectively enrich phosphatidylserine (PS) lipids with unsaturated sn-2 acyl chains, but not the fully saturated PS species. Thus, remodeling PS acyl chains may suppress KRAS oncogenesis. Lysophosphatidylcholine acyltransferases (LPCATs) remodel sn-2 acyl chains of phospholipids, with LPCAT1 preferentially generating the fully saturated lipids. Here, we show that stable expression of LPCAT1 depletes major PS species with unsaturated sn-2 chains while decreasing minor phosphatidylcholine (PC) species with the corresponding acyl chains. LPCAT1 expression more effectively disrupts the nanoclustering of oncogenic GFP-KRASG12V, which is restored by acute addback of exogenous unsaturated PS. LPCAT1 expression compromises signaling and oncogenic activities of the KRAS-dependent pancreatic tumor lines. LPCAT1 expression sensitizes human pancreatic tumor MiaPaCa-2 cells to KRASG12C specific inhibitor, Sotorasib. Statistical analyses of patient data further reveal that pancreatic cancer patients with KRAS mutations express less LPCAT1. Higher LPCAT1 expression also improves survival probability of pancreatic and lung adenocarcinoma patients with KRAS mutations. Thus, PS acyl chain remodeling selectively suppresses KRAS oncogenesis.

Intrinsically disordered membrane anchors of Rheb, RhoA and DiRas3 small GTPases: Molecular dynamics, membrane organization, and interactions

Protein structure has been well established to play a key role in determining function; however, intrinsically disordered proteins and regions (IDPs and IDRs) defy this paradigm. IDPs and IDRs exist as an ensemble of structures rather than a stable 3D structure yet play essential roles in many cell signaling processes. Nearly all Ras Superfamily GTPases are tethered to membranes by a lipid tail at the end of a flexible IDR. The sequence of these IDRs are key determinants of membrane localization, and interactions between the IDR and the membrane have been shown to affect signaling in RAS proteins through modulation of dynamic membrane organization. Here we utilized atomistic molecular dynamics simulations to study the membrane interactions, conformational dynamics, and lipid sorting of three IDRs from small GTPases Rheb, RhoA and DiRas3 in model membranes representing their physiological target membranes. We found that complementarity between lipidated IDR sequence and target membrane lipid composition is a determinant of conformational plasticity. We also show that electrostatic interactions between anionic lipids and basic residues on IDRs generate semi-stable conformational sub-states, and a lack of these residues leads to greater conformational diversity. Finally, we show that small GTPase IDRs with a polybasic domain alter local lipid composition by segregating anionic membrane lipids, and, in some cases, excluding other lipids from their immediate proximity in favor of anionic lipids.

Oliveira AS, Mulholland AJ, Serapian SA, Colombo G. Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations

Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.

Conformational ensemble-dependent lipid recognition and segregation by prenylated intrinsically disordered regions in small GTPases

We studied diverse prenylated intrinsically disordered regions (PIDRs) of Ras and Rho family small GTPases using long timescale atomistic molecular dynamics simulations in an asymmetric model membrane of phosphatidylcholine (PC) and phosphatidylserine (PS) lipids. Here we show that conformational plasticity is a key determinant of lipid sorting by polybasic PIDRs and provide evidence for lipid sorting based on both headgroup and acyl chain structures. We further show that conformational ensemble-based lipid recognition is generalizable to all polybasic PIDRs, and that the sequence outside the polybasic domain (PBD) modulates the conformational plasticity, bilayer adsorption, and interactions of PIDRs with membrane lipids. Specifically, we find that palmitoylation, the ratio of basic to acidic residues, and the hydrophobic content of the sequence outside the PBD significantly impact the diversity of conformational substates and hence the extent of conformation-dependent lipid interactions. We thus propose that the PBD is required but not sufficient for the full realization of lipid sorting by prenylated PBD-containing membrane anchors, and that the membrane anchor is not only responsible for high affinity membrane binding but also directs the protein to the right target membrane where it participates in lipid sorting.

Conformational ensemble dependent lipid recognition and segregation by prenylated intrinsically disordered regions in small GTPases

We studied diverse prenylated intrinsically disordered regions (PIDRs) of Ras and Rho family small GTPases using long timescale atomistic molecular dynamics simulations in an asymmetric model membrane of phosphatidylcholine (PC) and phosphatidylserine (PS) lipids. We show that conformational plasticity is a key determinant of lipid sorting by polybasic PIDRs and provide evidence for lipid sorting based on both headgroup and acyl chain structures. We further show that conformational ensemble-based lipid recognition is generalizable to all polybasic PIDRs, and that the sequence outside the polybasic domain (PBD) modulates the conformational plasticity, bilayer adsorption, and interactions of PIDRs with membrane lipids. Specifically, we found that palmitoylation, the ratio of basic to acidic residues, and the hydrophobic content of the sequence outside the PBD significantly impact the diversity of conformational substates and hence the extent of conformation-dependent lipid interactions. We thus propose that the PBD is required but not sufficient for the full realization of lipid sorting by prenylated PBD-containing membrane anchors, and that the membrane anchor is not only responsible for high affinity membrane binding but also directs the protein to the right target membrane where it participates in lipid sorting.

RAS nanoclusters are cell surface transducers that convert extracellular stimuli to intracellular signalling

Mutations of rat sarcoma virus (RAS) oncogenes (HRAS, KRAS and NRAS) can contribute to the development of cancers and genetic disorders (RASopathies). The spatiotemporal organization of RAS is an important property that warrants further investigation. In order to function, wild-type or oncogenic mutants of RAS must be localized to the inner leaflet of the plasma membrane (PM), which is driven by interactions between their C-terminal membrane-anchoring domains and PM lipids. The isoform-specific RAS-lipid interactions promote the formation of nanoclusters on the PM. As main sites for effector recruitment, these nanoclusters are biologically important. Since the spatial distribution of lipids is sensitive to changing environments, such as mechanical and electrical perturbations, RAS nanoclusters act as transducers to convert external stimuli to intracellular mitogenic signalling. As such, effective inhibition of RAS oncogenesis requires consideration of the complex interplay between RAS nanoclusters and various cell surface and extracellular stimuli. In this review, we discuss in detail how, by sorting specific lipids in the PM, RAS nanoclusters act as transducers to convert external stimuli into intracellular signalling.

Remodeling of the Plasma Membrane by Surface-Bound Protein Monomers and Oligomers: The Critical Role of Intrinsically Disordered Regions

The plasma membrane (PM) of cells is a dynamic structure whose morphology and composition is in constant flux. PM morphologic changes are particularly relevant for the assembly and disassembly of signaling platforms involving surface-bound signaling proteins, as well as for many other mechanochemical processes that occur at the PM surface. Surface-bound membrane proteins (SBMP) require efficient association with the PM for their function, which is often achieved by the coordinated interactions of intrinsically disordered regions (IDRs) and globular domains with membrane lipids. This review focuses on the role of IDR-containing SBMPs in remodeling the composition and curvature of the PM. The ability of IDR-bearing SBMPs to remodel the Gaussian and mean curvature energies of the PM is intimately linked to their ability to sort subsets of phospholipids into nanoclusters. We therefore discuss how IDRs of many SBMPs encode lipid-binding specificity or facilitate cluster formation, both of which increase their membrane remodeling capacity, and how SBMP oligomers alter membrane shape by monolayer surface area expansion and molecular crowding.

Phosphatidylserine and Phosphatidylethanolamine Asymmetry Have a Negligible Effect on the Global Structure, Dynamics, and Interactions of the KRAS Lipid Anchor

The intrinsically disordered C-terminus of the prominent oncogenic protein KRAS-4B (KRAS) selectively interacts and clusters with phosphatidylserine (PS) lipids in the plasma membrane (PM). This 11-residue segment, called tK, contains a polybasic domain (PBD) of six contiguous lysine residues and a farnesylated cysteine. Previous molecular dynamics (MD) simulation studies of tK in phosphatidylcholine (PC)/PS bilayers have suggested that backbone conformational dynamics modulate tK-PS interactions. These simulations have been conducted in symmetric membranes whereas the PM is compositionally asymmetric, with the inner leaflet, where KRAS is localized, being enriched with PS and phosphatidylethanolamine (PE) lipids. To examine if bilayer asymmetry affects tK conformational dynamics and interaction with lipids, we conducted two 10 μs long MD simulations of tK bound to a PC/PS and a PC/PS/PE bilayer in which the PS and PE lipids are distributed in one leaflet. We found that, first, these compositional asymmetries caused differences in acyl chain dynamics between leaflets, but the equilibrium structural and dynamic properties of the two asymmetric bilayers are similar; second, in both systems tK is highly dynamic and samples at least two distinct conformational states; third, PS-tK hydrogen-bonding interactions vary with peptide backbone conformations, and lysine side chains in the PBD predominantly interact with the serine oxygens of PS. These results are in good agreement with previous observations of tK in symmetric membranes. The effects of POPS asymmetry or the presence of POPE on tK are limited to modulating the relative contribution of individual side chains to interactions with lipids and redistributing conformational substates. Additional observations include the larger flexibility of tK in the current simulations, which we attribute to the longer duration of the simulations and the use of the CHARMM36m force field, which more accurately models intrinsically disordered peptides such as tK.

Lipid Profiles of RAS Nanoclusters Regulate RAS Function

The lipid-anchored RAS (Rat sarcoma) small GTPases (guanosine triphosphate hydrolases) are highly prevalent in human cancer. Traditional strategies of targeting the enzymatic activities of RAS have been shown to be difficult. Alternatively, RAS function and pathology are mostly restricted to nanoclusters on the plasma membrane (PM). Lipids are important structural components of these signaling platforms on the PM. However, how RAS nanoclusters selectively enrich distinct lipids in the PM, how different lipids contribute to RAS signaling and oncogenesis and whether the selective lipid sorting of RAS nanoclusters can be targeted have not been well-understood. Latest advances in quantitative super-resolution imaging and molecular dynamic simulations have allowed detailed characterization RAS/lipid interactions. In this review, we discuss the latest findings on the select lipid composition (with headgroup and acyl chain specificities) within RAS nanoclusters, the specific mechanisms for the select lipid sorting of RAS nanoclusters on the PM and how perturbing lipid compositions within RAS nanoclusters impacts RAS function and pathology. We also describe different strategies of manipulating lipid composition within RAS nanoclusters on the PM.

RAS Nanoclusters Selectively Sort Distinct Lipid Headgroups and Acyl Chains

RAS proteins are lipid-anchored small GTPases that switch between the GTP-bound active and GDP-bound inactive states. RAS isoforms, including HRAS, NRAS and splice variants KRAS4A and KRAS4B, are some of the most frequently mutated proteins in cancer. In particular, constitutively active mutants of KRAS comprise ∼80% of all RAS oncogenic mutations and are found in 98% of pancreatic, 45% of colorectal and 31% of lung tumors. Plasma membrane (PM) is the primary location of RAS signaling in biology and pathology. Thus, a better understanding of how RAS proteins localize to and distribute on the PM is critical to better comprehend RAS biology and to develop new strategies to treat RAS pathology. In this review, we discuss recent findings on how RAS proteins sort lipids as they undergo macromolecular assembly on the PM. We also discuss how RAS/lipid nanoclusters serve as signaling platforms for the efficient recruitment of effectors and signal transduction, and how perturbing the PM biophysical properties affect the spatial distribution of RAS isoforms and their functions.

Unveiling the Dynamics of KRAS4b on Lipid Model Membranes

Small GTPase proteins are ubiquitous and responsible for regulating several processes related to cell growth and differentiation. Mutations that stabilize their active state can lead to uncontrolled cell proliferation and cancer. Although these proteins are well characterized at the cellular scale, the molecular mechanisms governing their functions are still poorly understood. In addition, there is limited information about the regulatory function of the cell membrane which supports their activity. Thus, we have studied the dynamics and conformations of the farnesylated KRAS4b in various membrane model systems, ranging from binary fluid mixtures to heterogeneous raft mimics. Our approach combines long time-scale coarse-grained (CG) simulations and Markov state models to dissect the membrane-supported dynamics of KRAS4b. Our simulations reveal that protein dynamics is mainly modulated by the presence of anionic lipids and to some extent by the nucleotide state (activation) of the protein. In addition, our results suggest that both the farnesyl and the polybasic hypervariable region (HVR) are responsible for its preferential partitioning within the liquid-disordered (Ld) domains in membranes, potentially enhancing the formation of membrane-driven signaling platforms.

The KRAS and other prenylated polybasic domain membrane anchors recognize phosphatidylserine acyl chain structure

KRAS interacts with the inner leaflet of the plasma membrane (PM) using a hybrid anchor that comprises a lysine-rich polybasic domain (PBD) and a C-terminal farnesyl chain. Electrostatic interactions have been envisaged as the primary determinant of interactions between KRAS and membranes. Here, we integrated molecular dynamics (MD) simulations and superresolution spatial analysis in mammalian cells and systematically compared four equally charged KRAS anchors: the wild-type farnesyl hexa-lysine and engineered mutants comprising farnesyl hexa-arginine, geranylgeranyl hexa-lysine, and geranylgeranyl hexa-arginine. MD simulations show that these equally charged KRAS mutant anchors exhibit distinct interactions and packing patterns with different phosphatidylserine (PtdSer) species, indicating that prenylated PBD-bilayer interactions extend beyond electrostatics. Similar observations were apparent in intact cells, where each anchor exhibited binding specificities for PtdSer species with distinct acyl chain compositions. Acyl chain composition determined responsiveness of the spatial organization of different PtdSer species to diverse PM perturbations, including transmembrane potential, cholesterol depletion, and PM curvature. In consequence, the spatial organization and PM binding of each KRAS anchor precisely reflected the behavior of its preferred PtdSer ligand to these same PM perturbations. Taken together these results show that small GTPase PBD-prenyl anchors, such as that of KRAS, have the capacity to encode binding specificity for specific acyl chains as well as lipid headgroups, which allow differential responses to biophysical perturbations that may have biological and signaling consequences for the anchored GTPase.

Long-lasting Salt Bridges Provide the Anchoring Mechanism of Oncogenic Kirsten Rat Sarcoma Proteins at Cell Membranes

RAS proteins work as GDP-GTP binary switches and regulate cytoplasmic signaling networks that are able to control several cellular processes, playing an essential role in signal transduction pathways involved in cell growth, differentiation, and survival, so that overacting RAS signaling can lead to cancer. One of the hardest challenges to face is the design of mutation-selective therapeutic strategies. In this work, a G12D-mutated farnesylated GTP-bound Kirsten RAt sarcoma (KRAS) protein has been simulated at the interface of a DOPC/DOPS/cholesterol model anionic cell membrane. A specific long-lasting salt bridge connection between farnesyl and the hypervariable region of the protein has been identified as the main mechanism responsible for the binding of oncogenic farnesylated KRAS-4B to the cell membrane. Free-energy landscapes allowed us to characterize local and global minima of KRAS-4B binding to the cell membrane, revealing the main pathways between anchored and released states.

Biology, pathology, and therapeutic targeting of RAS

RAS was identified as a human oncogene in the early 1980s and subsequently found to be mutated in nearly 30% of all human cancers. More importantly, RAS plays a central role in driving tumor development and maintenance. Despite decades of effort, there remain no FDA approved drugs that directly inhibit RAS. The prevalence of RAS mutations in cancer and the lack of effective anti-RAS therapies stem from RAS' core role in growth factor signaling, unique structural features, and biochemistry. However, recent advances have brought promising new drugs to clinical trials and shone a ray of hope in the field. Here, we will exposit the details of RAS biology that illustrate its key role in cell signaling and shed light on the difficulties in therapeutically targeting RAS. Furthermore, past and current efforts to develop RAS inhibitors will be discussed in depth.

Blue light-triggered photochemistry and cytotoxicity of retinal

The chemical- and photo- toxicity of chromophore retinal on cells have long been debated. Although we recently showed that retinal and blue light exposure interrupt cellular signaling, a comprehensive study examining molecular underpinnings of this perturbation and its consequences to cellular fate is lacking. Here, we report molecular evidence for blue light excited-retinal induced oxidative damage of polyunsaturated lipid anchors in membrane-interacting signaling molecules and DNA damage in cells using live-cell imaging and in vitro experimentation. The incurred molecular damage irreversibly disrupted subcellular localization of these molecules, a crucial criterion for their signaling. We further show retinal accumulation in lipid-bilayers of cell membranes could enhance the lifetime of retinal in cells. Comparative response-signatures suggest that retinal triggers reactions upon photoexcitation similar to photodynamic therapy agents and generate reactive oxygen species in cells. Additionally, data also shows that exposing retinal-containing cells to sunlight induces substantial cytotoxicity. Collectively, our results explain a likely in vivo mechanism and reaction conditions under which bio-available retinal in physiological light conditions damages cells.

Probing the Conformational and Energy Landscapes of KRAS Membrane Orientation

Membrane reorientation of oncogenic RAS proteins is emerging as an important modulator of their functions. Previous studies have shown that the most common orientations include those with either the three C-terminal α-helices (OS1) or N-terminal β-strands (OS2) of the catalytic domain facing the membrane. OS1 and OS2 differ by the degree to which the effector-interacting surface is occluded by the membrane. However, the relative stability of these states and the rates of transition between them remained undetermined. How mutations might modulate preferences for specific orientation states is also far from clear. The current work attempted to address these questions through a comprehensive analysis of two 20 μs-long atomistic molecular dynamics simulations. The simulations were conducted on the oncogenic G12D and Q61H KRAS mutants bound to an anionic lipid bilayer. G12D and Q61H are among the most prevalent cancer-causing mutations at the P-loop and switch 2 regions of KRAS, respectively. We found that both mutants fluctuate in a similar manner between OS1 and OS2 via an intermediate orientation OS0, and both favor the signaling competent OS1 and OS0 over the occluded OS2. However, they differ in the details, such as in the extent to which they sample OS1. Analysis of the orientation free-energy landscapes estimated from the simulations indicate that OS1 and OS2 are the most stable states. However, the overall free energy surface is rugged, indicating a large diversity of conformations including at least two substates in each orientation state that differ in stability only by about 0.5-1.0 kcal/mol. Reversible transitions between OS1 and OS2 occur via two well-defined pathways that traverse OS0. In the minimum energy path, helix 4 remains close to the membrane as the angle of the catalytic domain from the membrane plane changes, resulting in a barrier of ∼1 kcal/mol for OS1/OS2 interconversions. Estimation of the rates of the various transitions based on survival probabilities yielded two rate constants in the order of 10⁷ and 10⁶ s^-1, which we attribute to intrinsic protein conformational dynamics and transient protein-lipid interactions, respectively. The faster process dominates every transition, confirming a previous suggestion that RAS membrane reorientation is driven by conformational fluctuations rather than protein-lipid interactions.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Interplay between Membrane Curvature and Cholesterol: Role of Palmitoylated Caveolin-1

Caveolin-1 (cav-1) is an important player in cell signaling and endocytosis that has been shown to colocalize with cholesterol-rich membrane domains. Experimental studies with varying cav-1 constructs have suggested that it can induce both cholesterol clustering and membrane curvature. Here, we probe the molecular origin of membrane curvature and cholesterol clustering by cav-1 by using coarse-grain molecular dynamics simulations. We have performed a series of simulations of a functionally important cav-1 construct, comprising the membrane-interacting domains and a C-terminal palmitoyl tail. Our results suggest that cav-1 is able to induce cholesterol clustering in the membrane leaflet to which it is bound as well as the opposing leaflet. A positive membrane curvature is observed upon cav-1 binding in cholesterol-containing bilayers. Interestingly, we observe an interplay between cholesterol clustering and membrane curvature such that cav-1 is able to induce higher membrane curvature in cholesterol-rich membranes. The role of the cav-1 palmitoyl tail is less clear and appears to increase the membrane contacts. Further, we address the importance of the secondary structure of cav-1 domains and show that it could play an important role in membrane curvature and cholesterol clustering. Our work is an important step toward a molecular picture of caveolae and vesicular endocytosis.

Cooperative Nonbonded Forces Control Membrane Binding of the pH-Low Insertion Peptide pHLIP

Peptides with the ability to bind and insert into the cell membrane have immense potential in biomedical applications. pH (low) insertion peptide (pHLIP), a water-soluble polypeptide derived from helix C of bacteriorhodopsin, can insert into a membrane at acidic pH to form a stable transmembrane α-helix. The insertion process takes place in three stages: pHLIP is unstructured and soluble in water at neutral pH (state I), unstructured and bound to the surface of a membrane at neutral pH (state II), and inserted into the membrane as an α-helix at low pH (state III). Using molecular dynamics simulations, we have modeled state II of pHLIP and a fast-folding variant of pHLIP, in which each peptide is bound to a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer surface. Our results provide strong support for recently published spectroscopic studies, namely that pHLIP preferentially binds to the bilayer surface as a function of location of anionic amino acids and that backbone dehydration occurs upon binding. Unexpectedly, we also observed several instances of segments of pHLIP folding into a stable helical turn. Our results provide a molecular level of detail that is essential to providing new insights into pHLIP function and to facilitate design of variants with improved membrane-active capabilities.

Deciphering lipid codes: K-Ras as a paradigm

The cell plasma membrane (PM) is a highly dynamic and heterogeneous lipid environment, driven by complex hydrophobic and electrostatic interactions among the hundreds of types of lipid species. Although the biophysical processes governing lipid lateral segregation in the cell PM have been established in vitro, biological implications of lipid heterogeneity are poorly understood. Of particular interest is how membrane proteins potentially utilize transient spatial clustering of PM lipids to regulate function. This review focuses on a lipid-anchored small GTPase K-Ras as an example to explore how its C-terminal membrane-anchoring domain, consisting of a contiguous hexa-lysine polybasic domain and an adjacent farnesyl anchor, possesses a complex coding mechanism for highly selective lipid sorting on the PM. How this lipid specificity modulates K-Ras signal transmission will also be discussed.

K-Ras4B Remains Monomeric on Membranes over a Wide Range of Surface Densities and Lipid Compositions

Ras is a membrane-anchored signaling protein that serves as a hub for many signaling pathways and also plays a prominent role in cancer. The intrinsic behavior of Ras on the membrane has captivated the biophysics community in recent years, especially the possibility that it may form dimers. In this article, we describe results from a comprehensive series of experiments using fluorescence correlation spectroscopy and single-molecule tracking to probe the possible dimerization of natively expressed and fully processed K-Ras4B in supported lipid bilayer membranes. Key to these studies is the fact that K-Ras4B has its native membrane anchor, including both the farnesylation and methylation of the terminal cysteine, enabling detailed exploration of possible effects of cholesterol and lipid composition on K-Ras4B membrane organization. The results from all conditions studied indicate that full-length K-Ras4B lacks intrinsic dimerization capability. This suggests that any lateral organization of Ras in living cell membranes likely stems from interactions with other factors.

Direct inhibition of RAS: Quest for the Holy Grail?

RAS GTPases (H-, K-, and N-RAS) are the most frequently mutated oncoprotein family in human cancer. However, the relatively smooth surface architecture of RAS and its picomolar affinity for nucleotide have given rise to the assumption that RAS is an "undruggable" target. Recent advancements in drug screening, molecular modeling, and a greater understanding of RAS function have led to a resurgence in efforts to pharmacologically target this challenging foe. This review focuses on the state of the art of RAS inhibition, the approaches taken to achieve this goal, and the challenges of translating these discoveries into viable therapeutics.

Geodesic curvature driven surface microdomain formation

Lipid bilayer membranes are not uniform and clusters of lipids in a more ordered state exist within the generally disorder lipid milieu of the membrane. These clusters of ordered lipids microdomains are now referred to as lipid rafts. Recent reports attribute the formation of these microdomains to the geometrical and molecular mechanical mismatch of lipids of different species on the boundary. Here we introduce the geodesic curvature to characterize the geometry of the domain boundary, and develop a geodesic curvature energy model to describe the formation of these microdomains as a result of energy minimization. Our model accepts the intrinsic geodesic curvature of any binary lipid mixture as an input, and will produce microdomains of the given geodesic curvature as demonstrated by three sets of numerical simulations. Our results are in contrast to the surface phase separation predicted by the classical surface Cahn-Hilliard equation, which tends to generate large domains as a result of the minimizing line tension. Our model provides a direct and quantified description of the structure inhomogeneity of lipid bilayer membrane, and can be coupled to the investigations of biological processes on membranes for which such inhomogeneity plays essential roles.

Influence of isoform-specific Ras lipidation motifs on protein partitioning and dynamics in model membrane systems of various complexity

The partitioning of the lipidated signaling proteins N-Ras and K-Ras4B into various membrane systems, ranging from single-component fluid bilayers, binary fluid mixtures, heterogeneous raft model membranes up to complex native-like lipid mixtures (GPMVs) in the absence and presence of integral membrane proteins have been explored in the last decade in a combined chemical-biological and biophysical approach. These studies have revealed pronounced isoform-specific differences regarding the lateral distribution in membranes and formation of protein-rich membrane domains. In this context, we will also discuss the effects of lipid head group structure and charge density on the partitioning behavior of the lipoproteins. Moreover, the dynamic properties of N-Ras and K-Ras4B have been studied in different model membrane systems and native-like crowded milieus. Addition of crowding agents such as Ficoll and its monomeric unit, sucrose, gradually favors clustering of Ras proteins in forming small oligomers in the bulk; only at very high crowder concentrations association is disfavored.

Membrane localization and dynamics of geranylgeranylated Rab5 hypervariable region

The small GTPase Rab5 is a key regulator of endosomal trafficking processes and a marker for the early endosome. The C-terminal hypervariable region (HVR) of Rab5 is post-translationally modified at residues Cys²¹² and Cys²¹³ to accommodate two geranylgeranyl anchors (C20 carbon chain length) in order to associate Rab5 with the membrane. The structural role of the HVR regarding protein-early endosome membrane recruitment is not resolved due to its high degree of flexibility and lack of crystallographic information. Here, full-atomistic and coarse-grained molecular dynamics simulations of the truncated Rab5 HVR^206-215 in three model membranes of increasing complexity (pure phospholipid bilayer, ternary membrane with cholesterol, six-component early endosome) were performed. Specific electrostatic interactions between the HVR^206-215 Arg²⁰⁹ residue and the phosphate group of the inositol ring of PI(3)P were detected. This shows that PI(3)P acts as a first contact site of protein recruitment to the early endosome. The free energy change of HVR^206-215 extraction from the bilayer was largest for the physiological negatively charged membrane. 5μs coarse-grained simulations revealed an active recruitment of PI(3)P to the HVR^206-215 supporting the formation of Rab5- and PI(3)P enriched signaling platforms.

Modulating short tryptophan- and arginine-rich peptides activity by substitution with histidine

BACKGROUND

High antimicrobial efficacy of short tryptophan-and arginine-rich peptides makes them good candidates in the fight against pathogens. Substitution of tryptophan and arginine by histidine could be used to modulate the peptides efficacy by optimizing their structures.

METHODS

The peptide (RRWWRWWRR), reported to showed good antimicrobial efficacy, was used as template, seven new analogs being designed substituting tryptophan or arginine with histidine. The peptides' efficacy was tested against E. coli, B. subtilis and S. aureus. The cytotoxicity and hemolytic effect were evaluated and the therapeutic index was inferred for each peptide. Atomic force microscopy and molecular simulation were used to analyze the effects of peptides on bacterial membrane.

RESULTS

The substitution of tryptophan by histidine proved to strongly modulate the antimicrobial activity, mainly by changing the peptide-to-membrane binding energy. The substitution of arginine has low effect on the antimicrobial efficacy. The presence of histidine residue reduced the cytotoxic and hemolytic activity of the peptides in some cases maintaining the same efficacy against bacteria. The peptides' antimicrobial activity was correlated to the 3D-hydrophobic moment and to a simple structure-based packing parameter.

CONCLUSION

The results show that some of these peptides have the potential to become good candidates to fight against bacteria. The substitution by histidine proved to fine tune the therapeutic index allowing the optimization of the peptide structure mainly by changing its binding energy and 3D-hydrophobic moment.

GENERAL SIGNIFICANCE

The short tryptophan reach peptides therapeutic index can be maximized using the histidine substitution to optimize their structure.

Probing the druggability of membrane-bound Rab5 by molecular dynamics simulations

Rab5 is a small GTPase and a key regulator in early endosomal trafficking. Rab5 and its effectors are involved in a large number of infectious diseases and certain types of cancer. We performed µs atomistic molecular dynamics simulations of inactive and active full-length Rab5 anchored to a complex model bilayer with composition of the early endosome membrane. Direct interactions between the Rab5 G domain and the bilayer were observed. We found two dominant nucleotide-dependent orientations characterised by a different accessibility of the switch regions. The “buried switch” orientation was mainly associated with inactive Rab5 accompanied with a rather extended structure of the hypervariable C-terminal region. Active Rab5 preferred an orientation in which the switch regions are accessible to effector proteins. These structural differences may provide an opportunity to selectively target one Rab5 state and lead to new approaches in the development of Rab5-specific therapies.

Lipoprotein insertion into membranes of various complexity: lipid sorting, interfacial adsorption and protein clustering

In a combined chemical-biological and biophysical approach we explored the membrane partitioning of the lipidated signaling proteins N-Ras and K-Ras4B into membrane systems of different complexity, ranging from one-component lipid bilayers and anionic binary and ternary heterogeneous membrane systems even up to partitioning studies on protein-free and protein-containing giant plasma membrane vesicles (GPMVs). To yield a pictorial view of the localization process, imaging using confocal laser scanning and atomic force microscopy was performed. The results reveal pronounced isoform-specific differences regarding the lateral distribution and formation of protein-rich membrane domains. Line tension is one of the key parameters controlling not only the size and dynamic properties of segregated lipid domains but also the partitioning process of N-Ras that acts as a lineactant. The formation of N-Ras protein clusters is even recorded for almost vanishing hydrophobic mismatch. Conversely, for K-Ras4B, selective localization and clustering are electrostatically mediated by its polybasic farnesylated C-terminus. The formation of K-Ras4B clusters is also observed for the multi-component GPMV membrane, i.e., it seems to be a general phenomenon, largely independent of the details of the membrane composition, including the anionic charge density of lipid headgroups. Our data indicate that unspecific and entropy-driven membrane-mediated interactions play a major role in the partitioning behavior, thus relaxing the need for a multitude of fine-tuned interactions. Such a scenario seems also to be reasonable recalling the high dynamic nature of cellular membranes. Finally, we note that even relatively simple models of heterogeneous membranes are able to reproduce many of the properties of much more complex biological membranes.

Lipid-Sorting Specificity Encoded in K-Ras Membrane Anchor Regulates Signal Output

K-Ras is targeted to the plasma membrane by a C-terminal membrane anchor that comprises a farnesyl-cysteine-methyl-ester and a polybasic domain. We used quantitative spatial imaging and atomistic molecular dynamics simulations to examine molecular details of K-Ras plasma membrane binding. We found that the K-Ras anchor binds selected plasma membrane anionic lipids with defined head groups and lipid side chains. The precise amino acid sequence and prenyl group define a combinatorial code for lipid binding that extends beyond simple electrostatics; within this code lysine and arginine residues are non-equivalent and prenyl chain length modifies nascent polybasic domain lipid preferences. The code is realized by distinct dynamic tertiary structures of the anchor on the plasma membrane that govern amino acid side-chain-lipid interactions. An important consequence of this specificity is the ability of such anchors when aggregated to sort subsets of phospholipids into nanoclusters with defined lipid compositions that determine K-Ras signaling output.

Oncogenic K-Ras Binds to an Anionic Membrane in Two Distinct Orientations: A Molecular Dynamics Analysis

K-Ras is a membrane-associated GTPase that cycles between active and inactive conformational states to regulate a variety of cell signaling pathways. Somatic mutations in K-Ras are linked to 15-20% of all human tumors. K-Ras attaches to the inner leaflet of the plasma membrane via a farnesylated polybasic domain; however, the structural details of the complex remain poorly understood. Based on extensive (7.5 μs total) atomistic molecular dynamics simulations here we show that oncogenic mutant K-Ras interacts with a negatively charged lipid bilayer membrane in multiple orientations. Of these, two highly populated orientations account for ∼54% of the conformers whose catalytic domain directly interacts with the bilayer. In one of these orientation states, membrane binding involves helices 3 and 4 of the catalytic domain in addition to the farnesyl and polybasic motifs. In the other orientation, β-strands 1-3 and helix 2 on the opposite face of the catalytic domain contribute to membrane binding. Flexibility of the linker region was found to be important for the reorientation. The biological significance of these observations was evaluated by initial experiments in cells overexpressing mutant K-Ras as well as by an analysis of Ras-effector complex structures. The results suggest that only one of the two major orientation states is capable of effector binding. We propose that the different modes of membrane binding may be exploited in structure-based drug design efforts for cancer therapy.

AMPK and Endothelial Nitric Oxide Synthase Signaling Regulates K-Ras Plasma Membrane Interactions via Cyclic GMP-Dependent Protein Kinase 2

K-Ras must localize to the plasma membrane and be arrayed in nanoclusters for biological activity. We show here that K-Ras is a substrate for cyclic GMP-dependent protein kinases (PKGs). In intact cells, activated PKG2 selectively colocalizes with K-Ras on the plasma membrane and phosphorylates K-Ras at Ser181 in the C-terminal polybasic domain. K-Ras phosphorylation by PKG2 is triggered by activation of AMP-activated protein kinase (AMPK) and requires endothelial nitric oxide synthase and soluble guanylyl cyclase. Phosphorylated K-Ras reorganizes into distinct nanoclusters that retune the signal output. Phosphorylation acutely enhances K-Ras plasma membrane affinity, but phosphorylated K-Ras is progressively lost from the plasma membrane via endocytic recycling. Concordantly, chronic pharmacological activation of AMPK → PKG2 signaling with mitochondrial inhibitors, nitric oxide, or sildenafil inhibits proliferation of K-Ras-positive non-small cell lung cancer cells. The study shows that K-Ras is a target of a metabolic stress-signaling pathway that can be leveraged to inhibit oncogenic K-Ras function.

Ras Dimer Formation as a New Signaling Mechanism and Potential Cancer Therapeutic Target

The K-, N-, and HRas small GTPases are key regulators of cell physiology and are frequently mutated in human cancers. Despite intensive research, previous efforts to target hyperactive Ras based on known mechanisms of Ras signaling have been met with little success. Several studies have provided compelling evidence for the existence and biological relevance of Ras dimers, establishing a new mechanism for regulating Ras activity in cells additionally to GTP-loading and membrane localization. Existing data also start to reveal how Ras proteins dimerize on the membrane. We propose a dimer model to describe Ras-mediated effector activation, which contrasts existing models of Ras signaling as a monomer or as a 5-8 membered multimer. We also discuss potential implications of this model in both basic and translational Ras biology.

The higher level of complexity of K-Ras4B activation at the membrane

Is nucleotide exchange sufficient to activate K-Ras4B? To signal, oncogenic rat sarcoma (Ras) anchors in the membrane and recruits effectors by exposing its effector lobe. With the use of NMR and molecular dynamics (MD) simulations, we observed that in solution, farnesylated guanosine 5'-diphosphate (GDP)-bound K-Ras4B is predominantly autoinhibited by its hypervariable region (HVR), whereas the GTP-bound state favors an activated, HVR-released state. On the anionic membrane, the catalytic domain adopts multiple orientations, including parallel (∼180°) and perpendicular (∼90°) alignments of the allosteric helices, with respect to the membrane surface direction. In the autoinhibited state, the HVR is sandwiched between the effector lobe and the membrane; in the active state, with membrane-anchored farnesyl and unrestrained HVR, the catalytic domain fluctuates reinlessly, exposing its effector-binding site. Dimerization and clustering can reduce the fluctuations. This achieves preorganized, productive conformations. Notably, we also observe HVR-autoinhibited K-Ras4B-GTP states, with GDP-bound-like orientations of the helices. Thus, we propose that the GDP/GTP exchange may not be sufficient for activation; instead, our results suggest that the GDP/GTP exchange, HVR sequestration, farnesyl insertion, and orientation/localization of the catalytic domain at the membrane conjointly determine the active or inactive state of K-Ras4B. Importantly, K-Ras4B-GTP can exist in active and inactive states; on its own, GTP binding may not compel K-Ras4B activation.-Jang, H., Banerjee, A., Chavan, T. S, Lu, S., Zhang, J., Gaponenko, V., Nussinov, R. The higher level of complexity of K-Ras4B activation at the membrane.

Ras Conformational Ensembles, Allostery, and Signaling

Ras proteins are classical members of small GTPases that function as molecular switches by alternating between inactive GDP-bound and active GTP-bound states. Ras activation is regulated by guanine nucleotide exchange factors that catalyze the exchange of GDP by GTP, and inactivation is terminated by GTPase-activating proteins that accelerate the intrinsic GTP hydrolysis rate by orders of magnitude. In this review, we focus on data that have accumulated over the past few years pertaining to the conformational ensembles and the allosteric regulation of Ras proteins and their interpretation from our conformational landscape standpoint. The Ras ensemble embodies all states, including the ligand-bound conformations, the activated (or inactivated) allosteric modulated states, post-translationally modified states, mutational states, transition states, and nonfunctional states serving as a reservoir for emerging functions. The ensemble is shifted by distinct mutational events, cofactors, post-translational modifications, and different membrane compositions. A better understanding of Ras biology can contribute to therapeutic strategies.

pMD-Membrane: A Method for Ligand Binding Site Identification in Membrane-Bound Proteins

Probe-based or mixed solvent molecular dynamics simulation is a useful approach for the identification and characterization of druggable sites in drug targets. However, thus far the method has been applied only to soluble proteins. A major reason for this is the potential effect of the probe molecules on membrane structure. We have developed a technique to overcome this limitation that entails modification of force field parameters to reduce a few pairwise non-bonded interactions between selected atoms of the probe molecules and bilayer lipids. We used the resulting technique, termed pMD-membrane, to identify allosteric ligand binding sites on the G12D and G13D oncogenic mutants of the K-Ras protein bound to a negatively charged lipid bilayer. In addition, we show that differences in probe occupancy can be used to quantify changes in the accessibility of druggable sites due to conformational changes induced by membrane binding or mutation.

Ras nanoclusters: Versatile lipid-based signaling platforms

Ras proteins assemble into transient nanoclusters on the plasma membrane. Nanoclusters are the sites of Ras effector recruitment and activation and are therefore essential for signal transmission. The dynamics of nanocluster formation and disassembly result in interesting emergent properties including high-fidelity signal transmission. More recently the lipid structure of Ras nanoclusters has been reported and shown to contribute to isoform-specific Ras signaling. In addition specific lipids play critical roles in mediating the formation, stability and dynamics of Ras nanoclusters. In consequence the spatiotemporal organization of these lipids has emerged as important and novel regulators of Ras function. This article is part of a Special Issue entitled: Nanoscale membrane organisation and signalling.

Overview of simulation studies on the enzymatic activity and conformational dynamics of the GTPase Ras

Over the last 40 years, we have learnt a great deal about the Ras onco-proteins. These intracellular molecular switches are essential for the function of a variety of physiological processes, including signal transduction cascades responsible for cell growth and proliferation. Molecular simulations and free energy calculations have played an essential role in elucidating the conformational dynamics and energetics underlying the GTP hydrolysis reaction catalysed by Ras. Here we present an overview of the main lessons from molecular simulations on the GTPase reaction and conformational dynamics of this important anti-cancer drug target. In the first part, we summarise insights from quantum mechanical and combined quantum mechanical/molecular mechanical simulations as well as other free energy methods and highlight consensus viewpoints as well as remaining controversies. The second part provides a very brief overview of new insights emerging from large-scale molecular dynamics simulations. We conclude with a perspective regarding future studies of Ras where computational approaches will likely play an active role.

Deformation of a Two-domain Lipid Bilayer due to Asymmetric Insertion of Lipid-modified Ras Peptides

Ras proteins are attached to the inner leaflet of the plasma membrane via a lipid-modified anchor. Membrane-bound Ras proteins laterally segregate into nanoscale signaling platforms called nanoclusters. It has been shown that the membrane domain preference of Ras nanoclusters varies with the nature of lipidation but their effect on the membrane has not been well understood. To investigate the effect of Ras insertion on membrane structure, we carried out numerous coarse-grained molecular dynamics (CGMD) simulations on a two-domain DPPC/DLiPC/cholesterol lipid bilayer in which different number and type of H-ras peptides were attached on one side. We have shown previously that this lipid mixture forms co-existing liquid-ordered/liquid-disordered (L o /L d ) domains and that different H-ras peptides form clusters that variously accumulate at the L o /L d regions or the boundary between them. Here we show that asymmetric insertion of each of these peptides induces a vertical relative displacement of the domains and deforms the bilayer, with the domain boundary serving as the center of deformation. The extent of the deformation, however, varies with the type and number of lipid modification. This is because the number and type of the Ras lipid tails determines the degree to which the stress caused by asymmetric peptide insertion is relieved by inter-leaflet cholesterol transfer and lipid tilt. In addition, we have characterized the mechanism of bilayer deformation based on the collective effect of the Ras peptides on inter-leaflet surface area, pressure profile and line tension differences. This allowed us to elucidate how Ras lipid modification affects membrane geometry and how a two-domain bilayer adjusts its shape through boundary deformation. The result contributes to a better understanding of Ras signaling platforms and highlights some of the mechanisms by which a multi-domain membrane responds to external perturbation.

Lessons from computer simulations of Ras proteins in solution and in membrane

BACKGROUND

A great deal has been learned over the last several decades about the function of Ras proteins in solution and membrane environments. While much of this knowledge has been derived from a plethora of experimental techniques, computer simulations have also played a substantial role.

SCOPE OF REVIEW

Our goal here is to summarize the contribution of molecular simulations to our current understanding of normal and aberrant Ras function. We focus on lessons from molecular dynamics simulations in aqueous and membrane environments.

MAJOR CONCLUSIONS

The central message is that a close interaction between theory and simulation on the one hand and cell-biological, spectroscopic and other experimental approaches on the other has played, and will likely continue to play, a vital role in Ras research.

GENERAL SIGNIFICANCE

Atomistic insights emerging from detailed simulations of Ras in solution and in bilayers may be the key to unlock the secret that to date prevented development of selective anti-Ras inhibitors for cancer therapy.

Staurosporines disrupt phosphatidylserine trafficking and mislocalize Ras proteins

Oncogenic mutant Ras is frequently expressed in human cancers, but no anti-Ras drugs have been developed. Since membrane association is essential for Ras biological activity, we developed a high content assay for inhibitors of Ras plasma membrane localization. We discovered that staurosporine and analogs potently inhibit Ras plasma membrane binding by blocking endosomal recycling of phosphatidylserine, resulting in redistribution of phosphatidylserine from plasma membrane to endomembrane. Staurosporines are more active against K-Ras than H-Ras. K-Ras is displaced to endosomes and undergoes proteasomal-independent degradation, whereas H-Ras redistributes to the Golgi and is not degraded. K-Ras nanoclustering on the plasma membrane is also inhibited. Ras mislocalization does not correlate with protein kinase C inhibition or induction of apoptosis. Staurosporines selectively abrogate K-Ras signaling and proliferation of K-Ras-transformed cells. These results identify staurosporines as novel inhibitors of phosphatidylserine trafficking, yield new insights into the role of phosphatidylserine and electrostatics in Ras plasma membrane targeting, and validate a new target for anti-Ras therapeutics.

Phosphatidic acid regulation of PIPKI is critical for actin cytoskeletal reorganization

Type I phosphatidylinositol-4-phosphate 5-kinase (PIPKI) is the main enzyme generating the lipid second messenger phosphatidylinositol-4,5-bisphosphate [PI(4,5)P2], which has critical functions in many cellular processes, such as cytoskeletal reorganization, membrane trafficking, and signal transduction. All three members of the PIPKI family are activated by phosphatidic acid (PA). However, how PA regulates the activity and functions of PIPKI have not been fully elucidated. In this study, we identify a PA-binding site on PIPKIγ. Mutation of this site inhibited the PA-stimulated activity and membrane localization of PIPKIγ as well as the formation of actin comets and foci induced by PIPKIγ. We also demonstrate that phospholipase D (PLD) generates a pool of PA involved in PIPKIγ regulation by showing that PLD inhibitors blocked the membrane localization of PIPKIγ and its ability to induce actin cytoskeletal reorganization. Targeting the PIPKIγ PA-binding-deficient mutant to membranes by a membrane localization sequence failed to restore the actin reorganization activity of PIPKIγ, suggesting that PA binding is not only involved in recruiting PIPKIγ to membranes but also may induce a conformational change. Taken together, these results reveal a new molecular mechanism through which PA regulates PIPKI and provides direct evidence that PA is important for the localization and functions of PIPKI in intact cells.

What drives the clustering of membrane-bound Ras?

The dynamic assembly and lateral organization of Ras proteins on the plasma membrane has been the focus of much research in recent years. It has been shown that different isoforms of Ras proteins share a nearly identical catalytic domain, yet form distinct and non-overlapping nanoclusters. Though this difference in the clustering behavior of Ras proteins has been attributed largely to their different C terminal lipid modification, its precise physical basis was not determined. Recently, we used computer simulations to study the mechanism by which the triply lipid-modified membrane-anchor of H-ras, and its partially de-lipidated variants, form nanoclusters in a model lipid bilayer. We found that the specific nature of the lipid modification is less important for cluster formation, but plays a key role for the domain-specific distribution of the nanoclusters. Here we provide additional details on the interplay between bilayer structure perturbation and peptide-peptide association that provide the physical driving force for clustering. We present some thoughts about how enthalpic (i.e., interaction) and entropic effects might regulate nanocluster size and stability.