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.

Biophysical and Biochemical Characterization of Structurally Diverse Small Molecule Hits for KRAS Inhibition.

We describe six compounds as early hits for the development of direct inhibitors of KRAS, an important anticancer drug target. We show that these compounds bind to KRAS with affinities in the low micromolar range and exert different effects on its interactions with binding partners. Some of the compounds exhibit selective binding to the activated form of KRAS and inhibit signal transduction through both the MAPK or the phosphatidylinositide 3-kinase PI3K-protein kinase B (AKT) pathway in cells expressing mutant KRAS. Most inhibit intrinsic and/or SOS-mediated KRAS activation while others inhibit RAS-effector interaction. We propose these compounds as starting points for the development of non-covalent allosteric KRAS inhibitors.

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.

Intracellular Ca2+ regulation of H+/Ca2+ antiporter YfkE mediated by a Ca2+ mini-sensor.

The H+/Ca2+ (calcium ion) antiporter (CAX) plays an important role in maintaining cellular Ca2+ homeostasis in bacteria, yeast, and plants by promoting Ca2+ efflux across the cell membranes. However, how CAX facilitates Ca2+ balance in response to dynamic cytosolic Ca2+ perturbations is unknown. Here, we identified a type of Ca2+ "mini-sensor" in YfkE, a bacterial CAX homolog from Bacillus subtilis. The mini-sensor is formed by six tandem carboxylate residues within the transmembrane (TM)5-6 loop on the intracellular membrane surface. Ca2+ binding to the mini-sensor triggers the transition of the transport mode of YfkE from a high-affinity to a low-affinity state. Molecular dynamics simulation and fluorescence resonance energy transfer analysis suggest that Ca2+ binding to the mini-sensor causes an adjacent segment, namely, the exchanger inhibitory peptide (XIP), to move toward the Ca2+ translocation pathway to interact with TM2a in an inward-open cavity. The specific interaction was demonstrated with a synthetic peptide of the XIP, which inhibits YfkE transport and interrupts conformational changes mediated by the mini-sensor. By comparing the apo and Ca2+-bound CAX structures, we propose the following Ca2+ transport regulatory mechanism of YfkE: Ca2+ binding to the mini-sensor induces allosteric conformational changes in the Ca2+ translocation pathway via the XIP, resulting in a rearrangement of the Ca2+-binding transport site in the midmembrane. Since the Ca2+ mini-sensor and XIP sequences are also identified in other CAX homologs and/or Ca2+ transporters, including the mammalian Na+/Ca2+ exchanger (NCX), our study provides a regulatory mechanism for the Ca2+/cation transporter superfamily.

Conformational spread and dynamics in allostery of NMDA receptors.

Allostery can be manifested as a combination of repression and activation in multidomain proteins allowing for fine tuning of regulatory mechanisms. Here we have used single molecule fluorescence resonance energy transfer (smFRET) and molecular dynamics simulations to study the mechanism of allostery underlying negative cooperativity between the two agonists glutamate and glycine in the NMDA receptor. These data show that binding of one agonist leads to conformational flexibility and an increase in conformational spread at the second agonist site. Mutational and cross-linking studies show that the dimer-dimer interface at the agonist-binding domain mediates the allostery underlying the negative cooperativity. smFRET on the transmembrane segments shows that they are tightly coupled in the unliganded and single agonist-bound form and only upon binding both agonists the transmembrane domain explores looser packing which would facilitate activation.

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.

Antipsychotic phenothiazine drugs bind to KRAS in vitro.

We used NMR to show that the antipsychotic phenothiazine drugs promazine and promethazine bind to GDP-KRAS. Promazine also binds to oncogenic GDP-KRAS(G12D), and to wild type GppNHp-KRAS. A panel of additional phenothiazines bind to GDP-KRAS but with lower affinity than promazine or promethazine. Binding is most dependent on substitutions at C-2 of the tricyclic phenothiazine ring. Promazine was used to generate an NMR-driven HADDOCK model of the drug/GDP-KRAS complex. The structural model shows the tricyclic phenothiazine ring of promazine associates with the hydrophobic pocket p1 that is bordered by the central ß sheet and Switch II in KRAS. Binding appears to stabilize helix 2 in a conformation that is similar to that seen in KRAS bound to other small molecules. Association of phenothiazines with KRAS may affect normal KRAS signaling that could contribute to multiple biological activities of these antipsychotic drugs. Moreover, the phenothiazine ring represents a new core scaffold on which to design modulators of KRAS activity.

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

Transmembrane potential difference (Vm) plays important roles in regulating various biological processes. At the macro level, Vm 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 Vm. 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 Vm 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 (ΔVm < 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 (ΔVm > 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.

Dynamics of Membrane-Bound G12V-KRAS from Simulations and Single-Molecule FRET in Native Nanodiscs.

Recent studies have shown that the small GTPase KRAS adopts multiple orientations with respect to the plane of anionic model membranes, whereby either the three C-terminal helices or the three N-terminal ß-strands of the catalytic domain face the membrane. This has functional implications because, in the latter, the membrane occludes the effector-interacting surface. However, it remained unclear how membrane reorientation occurs and, critically, whether it occurs in the cell in which KRAS operates as a molecular switch in signaling pathways. Herein, using data from a 20 µs-long atomistic molecular dynamics simulation of the oncogenic G12V-KRAS mutant in a phosphatidylcholine/phosphatidylserine bilayer, we first show that internal conformational fluctuations of flexible regions in KRAS result in three distinct membrane orientations. We then show, using single-molecule fluorescence resonance energy transfer measurements in native lipid nanodiscs derived from baby hamster kidney cells, that G12V-KRAS samples three conformational states that correspond to the predicted orientations. The combined results suggest that relatively small energy barriers separate orientation states and that signaling-competent conformations dominate the overall population.

Protein Partitioning into Ordered Membrane Domains: Insights from Simulations.

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 potential and dynamics in a ternary lipid mixture: insights from molecular dynamics simulations.

Transmembrane potential (Vm) plays critical roles in cell signaling and other functions. However, the impact of Vm on the structure and dynamics of membrane lipids and proteins, which are critical for the regulation of signaling, is still an open question. All-atom molecular dynamics (MD) simulation is emerging as a useful technique to address this issue. Previous atomistic MD simulations of pure or binary model membranes indicated that both ion imbalance and electric field can be used to generate Vm, but both approaches failed to yield structural changes in lipids with statistical significance. We hypothesized that a possible reason for this could be oversimplified membrane composition or limited sampling. In this work, we tested if and how Vm modulates the structure and dynamics of lipids in a physiologically relevant model membrane. Using a detailed side-by-side comparison, we first show that while both ion imbalance and electric field generate Vm in our complex membranes, only the latter could produce physiologically relevant Vm. We further show that double bonds in lipid acyl chains have a relatively large sensitivity to Vm. A single-bilayer model with an electric field showed the highest sensitivity in simulations under the isothermal-isobaric (NPT) ensemble, reproducing expected responses of head-group dipoles to Vm and suggesting that this approach may be more suitable for studying the structural effects of Vm. Our findings also shed light on the relationship between the macroscopic Vm and its atomic-level underpinnings.

Distinct dynamics and interaction patterns in H- and K-Ras oncogenic P-loop mutants.

Despite years of study, the structural or dynamical basis for the differential reactivity and oncogenicity of Ras isoforms and mutants remains unclear. In this study, we investigated the effects of amino acid variations on the structure and dynamics of wild type and oncogenic mutants G12D, G12V, and G13D of H- and K-Ras proteins. Based on data from µs-scale molecular dynamics simulations, we show that the overall structure of the proteins remains similar but there are important differences in dynamics and interaction networks. We identified differences in residue interaction patterns around the canonical switch and distal loop regions, and persistent sodium ion binding near the GTP particularly in the G13D mutants. Our results also suggest that different Ras variants have distinct local structural features and interactions with the GTP, variations that have the potential to affect GTP release and hydrolysis. Furthermore, we found that H-Ras proteins and particularly the G12V and G13D variants are significantly more flexible than their K-Ras counterparts. Finally, while most of the simulated proteins sampled the effector-interacting state 2 conformational state, G12V and G13D H-Ras adopted an open switch state 1 conformation that is defective in effector interaction. These differences have implications for Ras GTPase activity, effector or exchange factor binding, dimerization and membrane interaction. Proteins 2017; 85:1618-1632. © 2017 Wiley Periodicals, Inc.

Spatiotemporal Analysis of K-Ras Plasma Membrane Interactions Reveals Multiple High Order Homo-oligomeric Complexes.

Self-assembly of plasma membrane-associated Ras GTPases has major implications to the regulation of cell signaling. However, the structural basis of homo-oligomerization and the fractional distribution of oligomeric states remained undetermined. We have addressed these issues by deciphering the distribution of dimers and higher-order oligomers of K-Ras4B, the most frequently mutated Ras isoform in human cancers. We focused on the constitutively active G12V K-Ras and two of its variants, K101E and K101C/E107C, which respectively destabilize and stabilize oligomers. Using raster image correlation spectroscopy and number and brightness analysis combined with fluorescence recovery after photobleaching, fluorescence correlation spectroscopy and electron microscopy in live cells, we show that G12V K-Ras exists as a mixture of monomers, dimers and larger oligomers, while the K101E mutant is predominantly monomeric and K101C/E107C is dominated by oligomers. This observation demonstrates the ability of K-Ras to exist in multiple oligomeric states whose population can be altered by interfacial mutations. Using molecular modeling and simulations we further show that K-Ras uses two partially overlapping interfaces to form compositionally and topologically diverse oligomers. Our results thus provide the first detailed insight into the multiplicity, structure, and membrane organization of K-Ras homomers.

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.

Polyunsaturated Lipids Regulate Membrane Domain Stability by Tuning Membrane Order.

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.

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.

Computational allosteric ligand binding site identification on Ras proteins.

A number of computational techniques have been proposed to expedite the process of allosteric ligand binding site identification in inherently flexible and hence challenging drug targets. Some of these techniques have been instrumental in the discovery of allosteric ligand binding sites on Ras proteins, a group of elusive anticancer drug targets. This review provides an overview of these techniques and their application to Ras proteins. A summary of molecular docking and binding site identification is provided first, followed by a more detailed discussion of two specific techniques for binding site identification in ensembles of Ras conformations generated by molecular simulations.

Andrographolide derivatives inhibit guanine nucleotide exchange and abrogate oncogenic Ras function.

Aberrant signaling by oncogenic mutant rat sarcoma (Ras) proteins occurs in ∼15% of all human tumors, yet direct inhibition of Ras by small molecules has remained elusive. Recently, several small-molecule ligands have been discovered that directly bind Ras and inhibit its function by interfering with exchange factor binding. However, it is unclear whether, or how, these ligands could lead to drugs that act against constitutively active oncogenic mutant Ras. Using a dynamics-based pocket identification scheme, ensemble docking, and innovative cell-based assays, here we show that andrographolide (AGP)--a bicyclic diterpenoid lactone isolated from Andrographis paniculata--and its benzylidene derivatives bind to transient pockets on Kirsten-Ras (K-Ras) and inhibit GDP-GTP exchange. As expected for inhibitors of exchange factor binding, AGP derivatives reduced GTP loading of wild-type K-Ras in response to acute EGF stimulation with a concomitant reduction in MAPK activation. Remarkably, however, prolonged treatment with AGP derivatives also reduced GTP loading of, and signal transmission by, oncogenic mutant K-RasG12V. In sum, the combined analysis of our computational and cell biology results show that AGP derivatives directly bind Ras, block GDP-GTP exchange, and inhibit both wild-type and oncogenic K-Ras signaling. Importantly, our findings not only show that nucleotide exchange factors are required for oncogenic Ras signaling but also demonstrate that inhibiting nucleotide exchange is a valid approach to abrogating the function of oncogenic mutant Ras.