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.
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.
Display omitted
•The K-Ras anchor encodes previously unsuspected anionic lipid binding specificity•Lipid selectivity is realized by defined conformational structures of the anchor•Subtle changes to anchor sequence or prenylation profoundly alter lipid specificity•Lipid sorting capacity of the anchor is a key determinant of K-Ras signal output
Selective interactions with lipids shape Ras signaling output through formation of nanoclusters.
Cellular membranes are laterally organized into domains of distinct structures and compositions by the differential interaction affinities between various membrane lipids and proteins. A prominent ...example of such structures are lipid rafts, which are ordered, tightly packed domains that have been widely implicated in cellular processes. The functionality of raft domains is driven by their selective recruitment of specific membrane proteins to regulate their interactions and functions; however, there have been few general insights into the factors that determine the partitioning of membrane proteins between coexisting liquid domains. In this work, we used extensive coarse-grained and atomistic molecular dynamics simulations, potential of mean force calculations, and conceptual models to describe the partitioning dynamics and energetics of a model transmembrane domain from the linker of activation of T cells. We find that partitioning between domains is determined by an interplay between protein-lipid interactions and differential lipid packing between raft and nonraft domains. Specifically, we show that partitioning into ordered domains is promoted by preferential interactions between peptides and ordered lipids, mediated in large part by modification of the peptides by saturated fatty acids (i.e., palmitoylation). Ordered phase affinity is also promoted by elastic effects, specifically hydrophobic matching between the membrane and the peptide. Conversely, ordered domain partitioning is disfavored by the tight molecular packing of the lipids therein. The balance of these dominant drivers determines partitioning. In the case of the wild-type linker of activation of T cells transmembrane domain, these factors combine to yield enrichment of the peptide at Lo/Ld interfaces. These results define some of the general principles governing protein partitioning between coexisting membrane domains and potentially explain previous disparities among experiments and simulations across model systems.
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 107 and 106 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.
Eukaryotic plasma membranes are compartmentalized into functional lateral domains, including lipid-driven membrane rafts. Rafts are involved in most plasma membrane functions by selective recruitment ...and retention of specific proteins. However, the structural determinants of transmembrane protein partitioning to raft domains are not fully understood. Hypothesizing that protein transmembrane domains (TMDs) determine raft association, here we directly quantify raft affinity for dozens of TMDs. We identify three physical features that independently affect raft partitioning, namely TMD surface area, length, and palmitoylation. We rationalize these findings into a mechanistic, physical model that predicts raft affinity from the protein sequence. Application of these concepts to the human proteome reveals that plasma membrane proteins have higher raft affinity than those of intracellular membranes, consistent with raft-mediated plasma membrane sorting. Overall, our experimental observations and physical model establish general rules for raft partitioning of TMDs and support the central role of rafts in membrane traffic.
The plasma membrane (PM) serves as the functional interface between a cell and its environment, hosting extracellular signal transduction and nutrient transport among a variety of other processes. To ...support this extensive functionality, PMs are organized into lateral domains, including ordered, lipid-driven assemblies termed lipid rafts. Although the general requirements for ordered domain formation are well established, how these domains are regulated by cell-endogenous mechanisms or exogenous perturbations has not been widely addressed. In this context, an intriguing possibility is that dietary fats can incorporate into membrane lipids to regulate the properties and physiology of raft domains. Here, we investigate the effects of polyunsaturated fats on the organization of membrane domains across a spectrum of membrane models, including computer simulations, synthetic lipid membranes, and intact PMs isolated from mammalian cells. We observe that the ω-3 polyunsaturated fatty acid docosahexaenoic acid is robustly incorporated into membrane lipids, and this incorporation leads to significant remodeling of the PM lipidome. Across model systems, docosahexaenoic acid-containing lipids enhance the stability of ordered raft domains by increasing the order difference between them and coexisting nonraft domains. The relationship between interdomain order disparity and the stability of phase separation holds for a spectrum of different perturbations, including manipulation of cholesterol levels and high concentrations of exogenous amphiphiles, suggesting it as a general feature of the organization of biological membranes. These results demonstrate that polyunsaturated fats affect the composition and organization of biological membranes, suggesting a potential mechanism for the extensive effects of dietary fat on health and disease.
Recent experiments have shown that membrane-bound Ras proteins form transient, nanoscale signaling platforms that play a crucial role in high-fidelity signal transmission. However, a detailed ...characterization of these dynamic proteolipid substructures by high-resolution experimental techniques remains elusive. Here we use extensive semiatomic simulations to reveal the molecular basis for the formation and domain-specific distribution of Ras nanoclusters. As model systems, we chose the triply lipidated membrane targeting motif of H-ras (tH) and a large bilayer made up of di16:0-PC (DPPC), di18:2-PC (DLiPC), and cholesterol. We found that 4–10 tH molecules assemble into clusters that undergo molecular exchange in the sub-μs to μs time scale, depending on the simulation temperature and hence the stability of lipid domains. Driven by the opposite preference of tH palmitoyls and farnesyl for ordered and disordered membrane domains, clustered tH molecules segregate to the boundary of lipid domains. Additionally, a systematic analysis of depalmitoylated and defarnesylated tH variants allowed us to decipher the role of individual lipid modifications in domain-specific nanocluster localization and thereby explain why homologous Ras isoforms form nonoverlapping nanoclusters. Moreover, the localization of tH nanoclusters at domain boundaries resulted in a significantly lower line tension and increased membrane curvature. Taken together, these results provide a unique mechanistic insight into how protein assembly promoted by lipid-modification modulates bilayer shape to generate functional signaling platforms.
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.
Lipid-anchored Ras oncoproteins assemble into transient, nano-sized substructures on the plasma membrane. These substructures, called nanoclusters, were proposed to be crucial for high-fidelity ...signal transmission in cells. However, the molecular basis of Ras nanoclustering is poorly understood. In this work, we used coarse-grained (CG) molecular dynamics simulations to investigate the molecular mechanism by which full-length H-ras proteins form nanoclusters in a model membrane. We chose two different conformations of H-ras that were proposed to represent the active and inactive state of the protein, and a domain-forming model bilayer made up of di16:0-PC (DPPC), di18:2-PC (DLiPC) and cholesterol. We found that, irrespective of the initial conformation, Ras molecules assembled into a single large aggregate. However, the two binding modes, which are characterized by the different orientation of the G-domain with respect to the membrane, differ in dynamics and organization during and after aggregation. Some of these differences involve regions of Ras that are important for effector/modulator binding, which may partly explain observed differences in the ability of active and inactive H-ras nanoclusters to recruit effectors. The simulations also revealed some limitations in the CG force field to study protein assembly in solution, which we discuss in the context of proposed potential avenues of improvement.