K-Ras G-domain binding with signaling lipid phosphatidylinositol (4,5)-phosphate (PIP2): membrane association, protein orientation, and function

RAS G-domains allosterically contribute to the recognition of lipid headgroups and acyl chains

Mutant RAS are major contributors to cancer and signal primarily from nanoclusters on the plasma membrane (PM). Their C-terminal membrane anchors are main features of membrane association. However, the same RAS isoform bound to different guanine nucleotides spatially segregate. Different RAS nanoclusters all enrich a phospholipid, phosphatidylserine (PS). These findings suggest more complex membrane interactions. Our electron microscopy-spatial analysis shows that wild-types, G12V mutants, and membrane anchors of isoforms HRAS, KRAS4A, and KRAS4B prefer distinct PS species. Mechanistically, reorientation of KRAS4B G-domain exposes distinct residues, such as Arg 135 in orientation state 1 (OS1) and Arg 73/Arg 102 in OS2, to the PM and differentially facilitates the recognition of PS acyl chains. Allele-specific oncogenic mutations of KRAS4B also shift G-domain reorientation equilibrium. Indeed, KRAS4BG12V, KRAS4BG12D, KRAS4BG12C, KRAS4BG13D, and KRAS4BQ61H associate with PM lipids with headgroup and acyl chain specificities. Distribution of these KRAS4B oncogenic mutants favors different nanoscale membrane topography. Thus, RAS G-domains allosterically facilitate membrane lateral distribution.

Recent advances in liposome development for studying protein-lipid interactions

Protein-lipid interactions are crucial for various cellular biological processes like intracellular signaling, membrane transport, and cytoskeletal dynamics. Therefore, studying these interactions is essential to understand and unravel their specific functions. Nevertheless, the interacting proteins of many lipids are poorly understood and still require systematic study. Liposomes are the most well-known and familiar biomimetic systems used to study protein-lipid interactions. Although liposomes have been widely used for studying protein-lipid interactions in classical methods such as the co-flotation assay (CFA), co-sedimentation assay (CSA), and flow cytometric assay (FCA), an overview of their current applications and developments in high-throughput methods is not yet available. Here, we summarize the liposome development in low and high-throughput methods to study protein-lipid interactions. Besides, a constructive comment for each platform is presented to stimulate the advancement of these technologies in the future.

Direct K-Ras Inhibitors to Treat Cancers: Progress, New Insights, and Approaches to Treat Resistance

Here we discuss approaches to K-Ras inhibition and drug resistance scenarios. A breakthrough offered a covalent drug against K-RasG12C. Subsequent innovations harnessed same-allele drug combinations, as well as cotargeting K-RasG12C with a companion drug to upstream regulators or downstream kinases. However, primary, adaptive, and acquired resistance inevitably emerge. The preexisting mutation load can explain how even exceedingly rare mutations with unobservable effects can promote drug resistance, seeding growth of insensitive cell clones, and proliferation. Statistics confirm the expectation that most resistance-related mutations are in cis, pointing to the high probability of cooperative, same-allele effects. In addition to targeted Ras inhibitors and drug combinations, bifunctional molecules and innovative tri-complex inhibitors to target Ras mutants are also under development. Since the identities and potential contributions of preexisting and evolving mutations are unknown, selecting a pharmacologic combination is taxing. Collectively, our broad review outlines considerations and provides new insights into pharmacology and resistance.

Revealing KRas4b topology on the membrane surface

KRas4b is a membrane-bound regulatory protein belonging to the family of small GTPases that function as a molecular switch, facilitating signal transduction from activated membrane receptors to intracellular pathways controlling cell growth and proliferation. Oncogenic mutations locking KRas4b in the active GTP state are responsible for nearly 85% of all Ras-driven cancers. Understanding the membrane-bound state of KRas4b is crucial for designing new therapeutic approaches targeting oncogenic KRas-driven signaling pathways. Extensive research demonstrates the significant involvement of the membrane bilayer in Ras-effector interactions, with anionic lipids playing a critical role in determining protein conformations The preferred topology of KRas4b for interacting with signaling partners has been a long-time question. Computational studies suggest a membrane-proximal conformation, while other biophysical methods like neutron reflectivity propose a membrane-distal conformation. To address these gaps, we employed FRET measurements to investigate the conformation of KRas4b. Using fully post-translationally modified KRas4b, we designed a Nanodisc based FRET assay to study KRas4b-membrane interactions. We suggest an extended conformation of KRas4b relative to the membrane surface. Measurement of FRET donor - acceptor distances reveal that a negatively charged membrane surface weakly favors closer association with the membrane surface. Our findings provide insights into the role of anionic lipids in determining the dynamic conformations of KRas4b and shed light on the predominant conformation of its topology on lipid headgroups.

Transmembrane dimers of type 1 receptors sample alternate configurations: MD simulations using coarse grain Martini 3 versus AlphaFold2 Multimer

Structures and dynamics of transmembrane (TM) receptor regions are key to understanding their signaling mechanism across membranes. Here we examine configurations of TM region dimers, assembled using the recent Martini 3 force field for coarse-grain (CG) molecular dynamics simulations. At first glance, our results show only a reasonable agreement with ab initio predictions using PREDDIMER and AlphaFold2 Multimer and with nuclear magnetic resonance (NMR)-derived structures. 5 of 11 CG TM structures are similar to the NMR structures (within <3.5 Å root-mean-square deviation [RMSD]) compared with 10 and 9 using PREDDIMER and AlphaFold2, respectively (with 8 structures of the later within 1.5 Å). Surprisingly, AlphaFold2 predictions are closer to NMR structures when the 2001 instead of 2020 database is used for training. The CG simulations reveal that alternative configurations of TM dimers readily interconvert with a predominant population. The implications for transmembrane signaling are discussed, including for the development of peptide-based pharmaceuticals.

Phosphatidylinositol 4,5-Bisphosphate Sensing Lipid Raft via Inter-Leaflet Coupling Regulated by Acyl Chain Length of Sphingomyelin

Phosphatidylinositol 4,5-bisphosphate (PIP₂) is an important molecule located at the inner leaflet of cell membrane, where it serves as anchoring sites for a cohort of membrane-associated molecules and as a broad-reaching signaling intermediate. The lipid raft is thought as the major platform recruiting proteins for signal transduction and also known to mediate PIP₂ accumulation across the membrane. While the significance of this cross-membrane coupling is increasingly appreciated, it remains unclear whether and how PIP₂ senses the dynamic change of the ordered lipid domains over the packed hydrophobic core of the bilayer. Herein, by means of molecular dynamic simulation, we reveal that inner PIP₂ molecules can sense the outer lipid domain via inter-leaflet coupling, and the coupling manner is dictated by the acyl chain length of sphingomyelin (SM) partitioned to the lipid domain. Shorter SM promotes membrane domain registration, whereby PIP₂ accumulates beneath the domain across the membrane. In contrast, the anti-registration is thermodynamically preferred if the lipid domain has longer SM due to the hydrophobic mismatch between the corresponding acyl chains in SM and PIP₂. In this case, PIP₂ is expelled by the domain with a higher diffusivity. These results provide molecular insights into the regulatory mechanism of correlation between the outer lipid domain and inner PIP₂, both of which are critical components for cell signal transduction.

Quantitative Paramagnetic NMR-Based Analysis of Protein Orientational Dynamics on Membranes: Dissecting the KRas4B-Membrane Interactions

Peripheral membrane proteins can adopt distinct orientations on the surfaces of lipid bilayers that are often short-lived and challenging to characterize by conventional experimental methods. Here we describe a robust approach for mapping protein orientational landscapes through quantitative interpretation of paramagnetic relaxation enhancement (PRE) data arising from membrane mimetics with spin-labeled lipids. Theoretical analysis, followed by experimental verification, reveals insights into the distinct properties of the PRE observables that are generally distorted in the case of stably membrane-anchored proteins. To suppress the artifacts, we demonstrate that undistorted Γ₂ values can be obtained via transient membrane anchoring, based on which a computational framework is established for deriving accurate orientational ensembles obeying Boltzmann statistics. Application of the approach to KRas4B, a classical peripheral membrane protein whose orientations are critical for its functions and drug design, reveals four distinct orientational states that are close but not identical to those reported previously. Similar orientations are also found for a truncated KRas4B without the hypervariable region (HVR) that can sample a broader range of orientations, suggesting a confinement role of the HVR geometrically prohibiting severe tilting. Comparison of the KRas4B Γ₂ rates measured using nanodiscs containing different types of anionic lipids reveals identical Γ₂ patterns for the G-domain but different ones for the HVR, indicating only the latter is able to selectively interact with anionic lipids.

PI(4,5)P2 and Cholesterol: Synthesis, Regulation, and Functions

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂) is the most abundant membrane phosphoinositide and cholesterol is an essential component of the plasma membrane (PM). Both lipids play key roles in a variety of cellular functions including as signaling molecules and major regulators of protein function. This chapter provides an overview of these two important lipids. Starting from a brief description of their structure, synthesis, and regulation, the chapter continues to describe the primary functions and signaling processes in which PI(4,5)P₂ and cholesterol are involved. While PI(4,5)P₂ and cholesterol can act independently, they often act in concert or affect each other's impact. The chapters in this volume on "Cholesterol and PI(4,5)P₂ in Vital Biological Functions: From Coexistence to Crosstalk" focus on the emerging relationship between cholesterol and PI(4,5)P₂ in a variety of biological systems and processes. In this chapter, the next section provides examples from the ion channel field demonstrating that PI(4,5)P₂ and cholesterol can act via common mechanisms. The chapter ends with a discussion of future directions.

Transmembrane Membrane Readers form a Novel Class of Proteins That Include Peripheral Phosphoinositide Recognition Domains and Viral Spikes

Membrane proteins are broadly classified as transmembrane (TM) or peripheral, with functions that pertain to only a single bilayer at a given time. Here, we explicate a class of proteins that contain both transmembrane and peripheral domains, which we dub transmembrane membrane readers (TMMRs). Their transmembrane and peripheral elements anchor them to one bilayer and reversibly attach them to another section of bilayer, respectively, positioning them to tether and fuse membranes while recognizing signals such as phosphoinositides (PIs) and modifying lipid chemistries in proximity to their transmembrane domains. Here, we analyze full-length models from AlphaFold2 and Rosetta, as well as structures from nuclear magnetic resonance (NMR) spectroscopy and X-ray crystallography, using the Membrane Optimal Docking Area (MODA) program to map their membrane-binding surfaces. Eukaryotic TMMRs include phospholipid-binding C1, C2, CRAL-TRIO, FYVE, GRAM, GTPase, MATH, PDZ, PH, PX, SMP, StART and WD domains within proteins including protrudin, sorting nexins and synaptotagmins. The spike proteins of SARS-CoV-2 as well as other viruses are also TMMRs, seeing as they are anchored into the viral membrane while mediating fusion with host cell membranes. As such, TMMRs have key roles in cell biology and membrane trafficking, and include drug targets for diseases such as COVID-19.

Exploring CRD mobility during RAS/RAF engagement at the membrane

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

Nanoscopic Spatial Association between Ras and Phosphatidylserine on the Cell Membrane Studied with Multicolor Super Resolution Microscopy

Recent work suggests that Ras small GTPases interact with the anionic lipid phosphatidylserine (PS) in an isoform-specific manner, with direct implications for their biological functions. Studies on PS-Ras associations in cells, however, have relied on immuno-EM imaging of membrane sheets. To study their spatial relationships in intact cells, we have combined the use of Lact-C2-GFP, a biosensor for PS, with multicolor super resolution imaging based on DNA-PAINT. At ~20 nm spatial resolution, the resulting super resolution images clearly show the nonuniform molecular distribution of PS on the cell membrane and its co-enrichment with caveolae, as well as with unidentified membrane structures. Two-color imaging followed by spatial analysis shows that KRas-G12D and HRas-G12V both co-enrich with PS in model U2OS cells, confirming previous observations, yet exhibit clear differences in their association patterns. Whereas HRas-G12V is almost always co-enriched with PS, KRas-G12D is strongly co-enriched with PS in about half of the cells, with the other half exhibiting a more moderate association. In addition, perturbations to the actin cytoskeleton differentially impact PS association with the two Ras isoforms. These results suggest that PS-Ras association is context-dependent and demonstrate the utility of multiplexed super resolution imaging in defining the complex interplay between Ras and the membrane.

Structural insights into Ras regulation by SIN1

Both the mTORC2 and Ras-ERK pathways respond to growth factor stimulation and play critical roles in cell growth and proliferation, disarray of these pathways leads to many diseases, especially cancer. These two signaling pathways crosstalk at many levels; recently it's become clear that the SIN1 component of mTORC2 could interact with Ras family small GTPases, but how these two proteins interact at the molecular level and the functional outcomes of this interaction remain to be addressed. In this work we determined the high-resolution structure of Ras-SIN1 complexes and revealed the detailed interaction mechanism. We also showed that Ras-SIN1 association inhibits insulin-induced ERK activation. Insights from this work could improve our understanding of the disease-causing mechanism of errant mTORC2 or Ras proteins.

Over the years it has been established that SIN1, a key component of mTORC2, could interact with Ras family small GTPases through its Ras-binding domain (RBD). The physical association of Ras and SIN1/mTORC2 could potentially affect both mTORC2 and Ras-ERK pathways. To decipher the precise molecular mechanism of this interaction, we determined the high-resolution structures of HRas/KRas-SIN1 RBD complexes, showing the detailed interaction interface. Mutation of critical interface residues abolished Ras-SIN1 interaction and in SIN1 knockout cells we demonstrated that Ras-SIN1 association promotes SGK1 activity but inhibits insulin-induced ERK activation. With structural comparison and competition fluorescence resonance energy transfer (FRET) assays we showed that HRas-SIN1 RBD association is much weaker than HRas-Raf1 RBD but is slightly stronger than HRas-PI3K RBD interaction, providing a possible explanation for the different outcome of insulin or EGF stimulation. We also found that SIN1 isoform lacking the PH domain binds stronger to Ras than other longer isoforms and the PH domain appears to have an inhibitory effect on Ras-SIN1 binding. In addition, we uncovered a Ras dimerization interface that could be critical for Ras oligomerization. Our results advance our understanding of Ras-SIN1 association and crosstalk between growth factor-stimulated pathways.

Membrane Composition and Raf[CRD]-Membrane Attachment Are Driving Forces for K-Ras4B Dimer Stability

Ras proteins are membrane-anchored GTPases that regulate key cellular signaling networks. It has been recently shown that different anionic lipid types can affect the properties of Ras in terms of dimerization/clustering on the cell membrane. To understand the effects of anionic lipids on key spatiotemporal properties of dimeric K-Ras4B, we perform all-atom molecular dynamics simulations of the dimer K-Ras4B in the presence and absence of Raf[RBD/CRD] effectors on two model anionic lipid membranes: one containing 78% mol DOPC, 20% mol DOPS, and 2% mol PIP2 and another one with enhanced concentration of anionic lipids containing 50% mol DOPC, 40% mol DOPS, and 10% mol PIP2. Analysis of our results unveils the orientational space of dimeric K-Ras4B and shows that the stability of the dimer is enhanced on the membrane containing a high concentration of anionic lipids in the absence of Raf effectors. This enhanced stability is also observed in the presence of Raf[RBD/CRD] effectors although it is not influenced by the concentration of anionic lipids in the membrane, but rather on the ability of Raf[CRD] to anchor to the membrane. We generate dominant K-Ras4B conformations by Markov state modeling and yield the population of states according to the K-Ras4B orientation on the membrane. For the membrane containing anionic lipids, we observe correlations between the diffusion of K-Ras4B and PIP2 and anchoring of anionic lipids to the Raf[CRD] domain. We conclude that the presence of effectors with the Raf[CRD] domain anchoring on the membrane as well as the membrane composition both influence the conformational stability of the K-Ras4B dimer, enabling the preservation of crucial interface interactions.

Ras Multimers on the Membrane: Many Ways for a Heart-to-Heart Conversation

Formation of Ras multimers, including dimers and nanoclusters, has emerged as an exciting, new front of research in the 'old' field of Ras biomedicine. With significant advances made in the past few years, we are beginning to understand the structure of Ras multimers and, albeit preliminary, mechanisms that regulate their formation in vitro and in cells. Here we aim to synthesize the knowledge accrued thus far on Ras multimers, particularly the presence of multiple globular (G-) domain interfaces, and discuss how membrane nanodomain composition and structure would influence Ras multimer formation. We end with some general thoughts on the potential implications of Ras multimers in basic and translational biology.

Machine learning-driven multiscale modeling reveals lipid-dependent dynamics of RAS signaling proteins

RAS is a signaling protein associated with the cell membrane that is mutated in up to 30% of human cancers. RAS signaling has been proposed to be regulated by dynamic heterogeneity of the cell membrane. Investigating such a mechanism requires near-atomistic detail at macroscopic temporal and spatial scales, which is not possible with conventional computational or experimental techniques. We demonstrate here a multiscale simulation infrastructure that uses machine learning to create a scale-bridging ensemble of over 100,000 simulations of active wild-type KRAS on a complex, asymmetric membrane. Initialized and validated with experimental data (including a new structure of active wild-type KRAS), these simulations represent a substantial advance in the ability to characterize RAS-membrane biology. We report distinctive patterns of local lipid composition that correlate with interfacially promiscuous RAS multimerization. These lipid fingerprints are coupled to RAS dynamics, predicted to influence effector binding, and therefore may be a mechanism for regulating cell signaling cascades.

Mechanism of activation and the rewired network: New drug design concepts

Precision oncology benefits from effective early phase drug discovery decisions. Recently, drugging inactive protein conformations has shown impressive successes, raising the cardinal questions of which targets can profit and what are the principles of the active/inactive protein pharmacology. Cancer driver mutations have been established to mimic the protein activation mechanism. We suggest that the decision whether to target an inactive (or active) conformation should largely rest on the protein mechanism of activation. We next discuss the recent identification of double (multiple) same-allele driver mutations and their impact on cell proliferation and suggest that like single driver mutations, double drivers also mimic the mechanism of activation. We further suggest that the structural perturbations of double (multiple) in cis mutations may reveal new surfaces/pockets for drug design. Finally, we underscore the preeminent role of the cellular network which is deregulated in cancer. Our structure-based review and outlook updates the traditional Mechanism of Action, informs decisions, and calls attention to the intrinsic activation mechanism of the target protein and the rewired tumor-specific network, ushering innovative considerations in precision medicine.

Oncogenic KRAS G12D mutation promotes dimerization through a second, phosphatidylserine-dependent interface: a model for KRAS oligomerization

KRAS forms transient dimers and higher-order multimers (nanoclusters) on the plasma membrane, which drive MAPK signaling and cell proliferation. KRAS is a frequently mutated oncogene, and while it is well known that the most prevalent mutation, G12D, impairs GTP hydrolysis, thereby increasing KRAS activation, G12D has also been shown to enhance nanoclustering. Elucidating structures of dynamic KRAS assemblies on a membrane has been challenging, thus we have refined our NMR approach that uses nanodiscs to study KRAS associated with membranes. We incorporated paramagnetic relaxation enhancement (PRE) titrations and interface mutagenesis, which revealed that, in addition to the symmetric ‘α–α’ dimerization interface shared with wild-type KRAS, the G12D mutant also self-associates through an asymmetric ‘α–β’ interface. The ‘α–β’ association is dependent on the presence of phosphatidylserine lipids, consistent with previous reports that this lipid promotes KRAS self-assembly on the plasma membrane in cells. Experiments using engineered mutants to spoil each interface, together with PRE probes attached to the membrane or free in solvent, suggest that dimerization through the primary ‘α–α’ interface releases β interfaces from the membrane promoting formation of the secondary ‘α–β’ interaction, potentially initiating nanoclustering. In addition, the small molecule BI-2852 binds at a β–β interface, stabilizing a new dimer configuration that outcompetes native dimerization and blocks the effector-binding site. Our data indicate that KRAS self-association involves a delicately balanced conformational equilibrium between transient states, which is sensitive to disease-associated mutation and small molecule inhibitors. The methods developed here are applicable to biologically important transient interactions involving other membrane-associated proteins.

Studies of membrane-dependent dimerization of KRAS on nanodiscs using paramagnetic NMR titrations and mutagenesis revealed a novel asymmetric ‘α–β’ interface that provides a potential mechanism for the enhanced assembly of KRAS–G12D nanoclusters.

Membrane-bound KRAS approximates an entropic ensemble of configurations

KRAS4B is a membrane-anchored signaling protein and a primary target in cancer research. Predictions from molecular dynamics simulations that have previously shaped our mechanistic understanding of KRAS signaling disagree with recent experimental results from neutron reflectometry, NMR, and thermodynamic binding studies. To gain insight into these discrepancies, we compare this body of biophysical data to back-calculated experimental results from a series of molecular simulations that implement different subsets of molecular interactions. Our results show that KRAS4B approximates an entropic ensemble of configurations at model membranes containing 30% phosphatidylserine lipids, which is not significantly shaped by interactions between the globular G-domain of KRAS4B and the lipid membrane. These findings revise our understanding of KRAS signaling and promote a model in which the protein samples the accessible conformational space in a near-uniform manner while being available to bind to effector proteins.

Intrinsically disordered proteins and membranes: a marriage of convenience for cell signalling?

The structure-function paradigm has guided investigations into the molecules involved in cellular signalling for decades. The peripheries of this paradigm, however, start to unravel when considering the co-operation between proteins and the membrane in signalling processes. Intrinsically disordered regions hold distinct advantages over folded domains in terms of their binding promiscuity, sensitivity to their particular environment and their ease of modulation through post-translational modifications. Low sequence complexity and bias towards charged residues are also favourable for the multivalent electrostatic interactions that occur at the surfaces of lipid bilayers. This review looks at the principles behind the successful marriage between protein disorder and membranes in addition to the role of this partnership in modifying and regulating signalling in cellular processes. The HVR (hypervariable region) of small GTPases is highlighted as a well-studied example of the nuanced role a short intrinsically disordered region can play in the fine-tuning of signalling pathways.

The phosphoinositide code is read by a plethora of protein domains

INTRODUCTION

The proteins that decipher nucleic acid- and protein-based information are well known, however, those that read membrane-encoded information remain understudied. Here we report 70 different human, microbial and viral protein folds that recognize phosphoinositides (PIs), comprising the readers of a vast membrane code.

AREAS COVERED

Membrane recognition is best understood for FYVE, PH and PX domains, which exemplify hundreds of PI code readers. Comparable lipid interaction mechanisms may be mediated by kinases, adjacent C1 and C2 domains, trafficking arrestin, GAT and VHS modules, membrane-perturbing annexin, BAR, CHMP, ENTH, HEAT, syntaxin and Tubby helical bundles, multipurpose FERM, EH, MATH, PHD, PDZ, PROPPIN, PTB and SH2 domains, as well as systems that regulate receptors, GTPases and actin filaments, transfer lipids and assembled bacterial and viral particles.

EXPERT OPINION

The elucidation of how membranes are recognized has extended the genetic code to the PI code. Novel discoveries include PIP-stop and MET-stop residues to which phosphates and metabolites are attached to block phosphatidylinositol phosphate (PIP) recognition, memteins as functional membrane protein apparatuses, and lipidons as lipid "codons" recognized by membrane readers. At least 5% of the human proteome senses such membrane signals and allows eukaryotic organelles and pathogens to operate and replicate.

Structural basis for the association of PLEKHA7 with membrane-embedded phosphatidylinositol lipids

PLEKHA7 (pleckstrin homology domain containing family A member 7) plays key roles in intracellular signaling, cytoskeletal organization, and cell adhesion, and is associated with multiple human cancers. The interactions of its pleckstrin homology (PH) domain with membrane phosphatidyl-inositol-phosphate (PIP) lipids are critical for proper cellular localization and function, but little is known about how PLEKHA7 and other PH domains interact with membrane-embedded PIPs. Here we describe the structural basis for recognition of membrane-bound PIPs by PLEHA7. Using X-ray crystallography, nuclear magnetic resonance, molecular dynamics simulations, and isothermal titration calorimetry, we show that the interaction of PLEKHA7 with PIPs is multivalent, distinct from a discrete one-to-one interaction, and induces PIP clustering. Our findings reveal a central role of the membrane assembly in mediating protein-PIP association and provide a roadmap for understanding how the PH domain contributes to the signaling, adhesion, and nanoclustering functions of PLEKHA7.

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.

RAS Nanoclusters: Dynamic Signaling Platforms Amenable to Therapeutic Intervention

RAS proteins are mutated in approximately 20% of all cancers and are generally associated with poor clinical outcomes. RAS proteins are localized to the plasma membrane and function as molecular switches, turned on by partners that receive extracellular mitogenic signals. In the on-state, they activate intracellular signal transduction cascades. Membrane-bound RAS molecules segregate into multimers, known as nanoclusters. These nanoclusters, held together through weak protein-protein and protein-lipid associations, are highly dynamic and respond to cellular input signals and fluctuations in the local lipid environment. Disruption of RAS nanoclusters results in downregulation of RAS-mediated mitogenic signaling. In this review, we discuss the propensity of RAS proteins to display clustering behavior and the interfaces that are associated with these assemblies. Strategies to therapeutically disrupt nanocluster formation or the stabilization of signaling incompetent RAS complexes are discussed.

Inhibition of Nonfunctional Ras

Intuitively, functional states should be targeted; not nonfunctional ones. So why could drugging the inactive K-Ras4BG12Cwork-but drugging the inactive kinase will likely not? The reason is the distinct oncogenic mechanisms. Kinase driver mutations work by stabilizing the active state and/or destabilizing the inactive state. Either way, oncogenic kinases are mostly in the active state. Ras driver mutations work by quelling its deactivation mechanisms, GTP hydrolysis, and nucleotide exchange. Covalent inhibitors that bind to the inactive GDP-bound K-Ras4BG12C conformation can thus work. By contrast, in kinases, allosteric inhibitors work by altering the active-site conformation to favor orthosteric drugs. From the translational standpoint this distinction is vital: it expedites effective pharmaceutical development and extends the drug classification based on the mechanism of action. Collectively, here we postulate that drug action relates to blocking the mechanism of activation, not to whether the protein is in the active or inactive state.

Influence of Cholesterol on the Orientation of the Farnesylated GTP-Bound KRas-4B Binding with Anionic Model Membranes

The Ras family of proteins is tethered to the inner leaflet of the cell membranes which plays an essential role in signal transduction pathways that promote cellular proliferation, survival, growth, and differentiation. KRas-4B, the most mutated Ras isoform in different cancers, has been under extensive study for more than two decades. Here we have focused our interest on the influence of cholesterol on the orientations that KRas-4B adopts with respect to the plane of the anionic model membranes. How cholesterol in the bilayer might modulate preferences for specific orientation states is far from clear. Herein, after analyzing data from in total 4000 ns-long molecular dynamics (MD) simulations for four KRas-4B systems, properties such as the area per lipid and thickness of the membrane as well as selected radial distribution functions, penetration of different moieties of KRas-4B, and internal conformational fluctuations of flexible moieties in KRas-4B have been calculated. It has been shown that high cholesterol content in the plasma membrane (PM) favors one orientation state (OS1), exposing the effector-binding loop for signal transduction in the cell from the atomic level. We confirm that high cholesterol in the PM helps KRas-4B mutant stay in its constitutively active state, which suggests that high cholesterol intake can increase mortality and may promote cancer progression for cancer patients. We propose that during the treatment of KRas-4B-related cancers, reducing the cholesterol level in the PM and sustaining cancer progression by controlling the plasma cholesterol intake might be taken into account in anti-cancer therapies.

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.

Targeting KRAS mutant cancers by preventing signaling transduction in the MAPK pathway

KRAS genes are the most commonly mutated oncogenes in cancer. Unfortunately, effective therapeutic strategies for targeting KRAS mutant cancers have proven to be difficult to obtain. A key reason for this setback is due to the lack of success direct KRAS mutant inhibitors have received. Researchers have turned their efforts away from targeting the KRAS nucleotide-binding site directly and towards targeting other areas of the MAPK signaling pathway to block KRAS function. Researchers found that inhibiting enzymes and protein-protein interactions involved in the MAPK signaling pathway inhibit the activation of KRAS mutant therefore can lead to a potential therapeutic for KRAS mutated cancers. Throughout the past two decades, various indirect inhibitors have been designed and tested. EGFR and MEK inhibitors have presented with less success; however, significant advances have been made when targeting the plasma membrane localization process and the allosteric site of KRAS mutant. Farnesyltransferase and allosteric inhibitors have both advanced to human clinical trials. This comprehensive review presents the most recent developments of direct and indirect KRAS mutant inhibitors. This review summarizes published data on the inhibitory and anti-cancer activity of compounds that target KRAS activation as well as highlights the most promising strategies for targeting KRAS mutant cancers.

Mechanisms of Ras Membrane Organization and Signaling: Ras Rocks Again

Ras is the most frequently mutated oncogene and recent drug development efforts have spurred significant new research interest. Here we review progress toward understanding how Ras functions in nanoscale, proteo-lipid signaling complexes on the plasma membrane, called nanoclusters. We discuss how G-domain reorientation is plausibly linked to Ras-nanoclustering and -dimerization. We then look at how these mechanistic features could cooperate in the engagement and activation of RAF by Ras. Moreover, we show how this structural information can be integrated with microscopy data that provide nanoscale resolution in cell biological experiments. Synthesizing the available data, we propose to distinguish between two types of Ras nanoclusters, an active, immobile RAF-dependent type and an inactive/neutral membrane anchor-dependent. We conclude that it is possible that Ras reorientation enables dynamic Ras dimerization while the whole Ras/RAF complex transits into an active state. These transient di/oligomer interfaces of Ras may be amenable to pharmacological intervention. We close by highlighting a number of open questions including whether all effectors form active nanoclusters and whether there is an isoform specific composition of Ras nanocluster.

Uncovering a membrane-distal conformation of KRAS available to recruit RAF to the plasma membrane

The small GTPase KRAS is localized at the plasma membrane where it functions as a molecular switch, coupling extracellular growth factor stimulation to intracellular signaling networks. In this process, KRAS recruits effectors, such as RAF kinase, to the plasma membrane where they are activated by a series of complex molecular steps. Defining the membrane-bound state of KRAS is fundamental to understanding the activation of RAF kinase and in evaluating novel therapeutic opportunities for the inhibition of oncogenic KRAS-mediated signaling. We combined multiple biophysical measurements and computational methodologies to generate a consensus model for authentically processed, membrane-anchored KRAS. In contrast to the two membrane-proximal conformations previously reported, we identify a third significantly populated state using a combination of neutron reflectivity, fast photochemical oxidation of proteins (FPOP), and NMR. In this highly populated state, which we refer to as "membrane-distal" and estimate to comprise ∼90% of the ensemble, the G-domain does not directly contact the membrane but is tethered via its C-terminal hypervariable region and carboxymethylated farnesyl moiety, as shown by FPOP. Subsequent interaction of the RAF1 RAS binding domain with KRAS does not significantly change G-domain configurations on the membrane but affects their relative populations. Overall, our results are consistent with a directional fly-casting mechanism for KRAS, in which the membrane-distal state of the G-domain can effectively recruit RAF kinase from the cytoplasm for activation at the membrane.

Trop2 inhibition of P16 expression and the cell cycle promotes intracellular calcium release in OSCC

Trop2 is an intracellular calcium signal transducer and a prognostic biomarker in many cancers. P16 is a cell cycle gene that negatively regulates cell proliferation and division in most human cancers. Oral squamous cell carcinoma (OSCC) is a common malignant tumor subgroup of head and neck squamous cell carcinoma worldwide. Both Ca²⁺-dependent and cell cycle signaling pathways play vital roles in OSCC, although the molecular mechanisms remain to be elucidated. We aimed to examine the function of Trop2 and P16 in regulating intracellular calcium ions and the cell cycle in OSCC cell lines. Furtherly, the mRNA and protein expression levels of Trop2 and P16 in OSCC tissue samples were assessed, and their function was evaluated as potential clinical prognostic biomarkers. Trop2 promoted intracellular calcium ion release in OSCC and induced S phase of the cell cycle. Moreover, Trop2-mediated Ca²⁺ inhibited P16 expression through the AMP-activated protein kinase pathway in OSCC. Interestingly, P16 overexpression could not reverse these phenomena in vitro. We also demonstrated that human OSCC tissues showed high Trop2 mRNA and protein expression, and Trop2+/P16- expression is an independent prognostic marker for OSCC patients. Our data suggest that Trop2+/P16- may be a valuable prognostic marker for OSCC and that Trop2 inhibits P16 expression and induces S phase by promoting intracellular calcium release in OSCC.

A regulatory role of membrane by direct modulation of the catalytic kinase domain

Cell membrane modulates the function and activity of specific proteins and acts more than just a non-specific scaffolding machinery. In this review, I focus on studies that highlight a direct membrane-mediated modulation of the catalytic kinase domain of a variety of kinases thereby regulating the kinase activity. It emerges that membrane provides a second level of regulation once kinase domain is relieved of its inactive auto-inhibitory state. For the first time a generalized regulatory role of membrane is proposed that governs the kinase activity by modulating the catalytic kinase domain. Striking similarities among a variety of multi-domain kinases as well as single-domain lipidated enzymes such as RAS proteins are presented.

How Anionic Lipids Affect Spatiotemporal Properties of KRAS4B on Model Membranes

RAS proteins are small membrane-anchored GTPases that regulate key cellular signaling networks. It has been recently shown that different anionic lipid types can affect the spatiotemporal properties of RAS through dimerization/clustering and signaling fidelity. To understand the effects of anionic lipids on key spatiotemporal properties of RAS, we dissected 1 ms of data from all-atom molecular dynamics simulations for KRAS4B on two model anionic lipid membranes that have 30% of POPS mixed with neutral POPC and 8% of PIP2 mixed with POPC. We unveiled the orientation space of KRAS4B, whose kinetics were slower and more distinguishable on the membrane containing PIP2 than the membrane containing POPS. Particularly, the PIP2-mixed membrane can differentiate a third kinetic orientation state from the other two known orientation states. We observed that each orientation state may yield different binding modes with an RAF kinase, which is required for activating the MAPK/ERK signaling pathway. However, an overall occluded probability, for which RAF kinases cannot bind KRAS4B, remains unchanged on the two different membranes. We identified rare fast diffusion modes of KRAS4B that appear coupled with orientations exposed to cytosolic RAF. Particularly, on the membrane having PIP2, we found nonlinear correlations between the orientation states and the conformations of the cationic farnesylated hypervariable region, which acts as an anchor in the membrane. Using diffusion coefficients estimated from the all-atom simulations, we quantified the effect of PIP2 and POPS on the KRAS4B dimerization via Green's function reaction dynamics simulations, in which the averaged dimerization rate is 12.5% slower on PIP2-mixed membranes.

Multivalent assembly of KRAS with the RAS-binding and cysteine-rich domains of CRAF on the membrane

Membrane anchoring of farnesylated KRAS is critical for activation of RAF kinases, yet our understanding of how these proteins interact on the membrane is limited to isolated domains. The RAS-binding domain (RBD) and cysteine-rich domain (CRD) of RAF engage KRAS and the plasma membrane, unleashing the kinase domain from autoinhibition. Due to experimental challenges, structural insight into this tripartite KRAS:RBD-CRD:membrane complex has relied on molecular dynamics simulations. Here, we report NMR studies of the KRAS:CRAF RBD-CRD complex. We found that the nucleotide-dependent KRAS-RBD interaction results in transient electrostatic interactions between KRAS and CRD, and we mapped the membrane interfaces of the CRD, RBD-CRD, and the KRAS:RBD-CRD complex. RBD-CRD exhibits dynamic interactions with the membrane through the canonical CRD lipid-binding site (CRD β7-8), as well as an alternative interface comprising β6 and the C terminus of CRD and β2 of RBD. Upon complex formation with KRAS, two distinct states were observed by NMR: State A was stabilized by membrane association of CRD β7-8 and KRAS α4-α5 while state B involved the C terminus of CRD, β3-5 of RBD, and part of KRAS α5. Notably, α4-α5, which has been proposed to mediate KRAS dimerization, is accessible only in state B. A cancer-associated mutation on the state B membrane interface of CRAF RBD (E125K) stabilized state B and enhanced kinase activity and cellular MAPK signaling. These studies revealed a dynamic picture of the assembly of the KRAS-CRAF complex via multivalent and dynamic interactions between KRAS, CRAF RBD-CRD, and the membrane.

ACLY is the novel signaling target of PIP2/PIP3 and Lyn in acute myeloid leukemia

A fundamental feature of tumor progression is reprogramming of metabolic pathways. ATP citrate lyase (ACLY) is a key metabolic enzyme that catalyzes the generation of Acetyl-CoA and is upregulated in cancer cells and required for their growth. The phosphoinositide 3-kinase (PI3K) and Src-family kinase (SFK) Lyn are constitutively activate in many cancers. We show here, for the first time, that both the substrate and product of PI3K, phosphatidylinositol-(4,5)-bisphosphate (PIP₂) and phosphatidylinositol-(3,4,5)-trisphosphate (PIP₃), respectively, bind to ACLY in Acute Myeloid Leukemia (AML) patient-derived, but not normal donor-derived cells. We demonstrate the binding of PIP₂ to the CoA-binding domain of ACLY and identify the six tyrosine residues of ACLY that are phosphorylated by Lyn. Three of them (Y682, Y252, Y227) can be also phosphorylated by Src and they are located in catalytic, citrate binding and ATP binding domains, respectively. PI3K and Lyn inhibitors reduce the ACLY enzyme activity, ACLY-mediated Acetyl-CoA synthesis, phospholipid synthesis, histone acetylation and cell growth. Thus, PIP₂/PIP₃ binding and Src tyrosine kinases-mediated stimulation of ACLY links oncogenic pathways to Acetyl-CoA-dependent pro-growth and survival metabolic pathways in cancer cells. These results indicate a novel function for Lyn, as a regulator of Acetyl-CoA-mediated metabolic pathways.

Comparative effects of oncogenic mutations G12C, G12V, G13D, and Q61H on local conformations and dynamics of K-Ras

K-Ras is the most frequently mutated protein in human cancers. However, until very recently, its oncogenic mutants were viewed as undruggable. To develop inhibitors that directly target oncogenic K-Ras mutants, we need to understand both their mutant-specific and pan-mutant dynamics and conformations. Recently, we have investigated how the most frequently observed K-Ras mutation in cancer patients, G12D, changes its local dynamics and conformations (Vatansever et al., 2019). Here, we extend our analysis to study and compare the local effects of other frequently observed oncogenic mutations, G12C, G12V, G13D and Q61H. For this purpose, we have performed Molecular Dynamics (MD) simulations of each mutant when active (GTP-bound) and inactive (GDP-bound), analyzed their trajectories, and compared how each mutant changes local residue conformations, inter-protein distance distributions, local flexibility and residue pair correlated motions. Our results reveal that in the four active oncogenic mutants we have studied, the α2 helix moves closer to the C-terminal of the α3 helix. However, P-loop mutations cause α3 helix to move away from Loop7, and only G12 mutations change the local conformational state populations of the protein. Furthermore, the motions of coupled residues are mutant-specific: G12 mutations lead to new negative correlations between residue motions, while Q61H destroys them. Overall, our findings on the local conformational states and protein dynamics of oncogenic K-Ras mutants can provide insights for both mutant-selective and pan-mutant targeted inhibition efforts.

The Plasma Membrane as a Competitive Inhibitor and Positive Allosteric Modulator of KRas4B Signaling

Mutant Ras proteins are important drivers of human cancers, yet no approved drugs act directly on this difficult target. Over the last decade, the idea has emerged that oncogenic signaling can be diminished by molecules that drive Ras into orientations in which effector-binding interfaces are occluded by the cell membrane. To support this approach to drug discovery, we characterize the orientational preferences of membrane-bound K-Ras4B in 1.45-ms aggregate time of atomistic molecular dynamics simulations. Individual simulations probe active or inactive states of Ras on membranes with or without anionic lipids. We find that the membrane orientation of Ras is relatively insensitive to its bound guanine nucleotide and activation state but depends strongly on interactions with anionic phosphatidylserine lipids. These lipids slow Ras' translational and orientational diffusion and promote a discrete population in which small changes in orientation control Ras' competence to bind multiple regulator and effector proteins. Our results suggest that compound-directed conversion of constitutively active mutant Ras into functionally inactive forms may be accessible via subtle perturbations of Ras' orientational preferences at the membrane surface.

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.

Membrane curvature sensing of the lipid-anchored K-Ras small GTPase

Plasma membrane (PM) curvature defines cell shape and intracellular organelle morphologies and is a fundamental cell property. Growth/proliferation is more stimulated in flatter cells than the same cells in elongated shapes. PM-anchored K-Ras small GTPase regulates cell growth/proliferation and plays key roles in cancer. The lipid-anchored K-Ras form nanoclusters selectively enriched with specific phospholipids, such as phosphatidylserine (PS), for efficient effector recruitment and activation. K-Ras function may, thus, be sensitive to changing lipid distribution at membranes with different curvatures. Here, we used complementary methods to manipulate membrane curvature of intact/live cells, native PM blebs, and synthetic liposomes. We show that the spatiotemporal organization and signaling of an oncogenic mutant K-Ras G12V favor flatter membranes with low curvature. Our findings are consistent with the more stimulated growth/proliferation in flatter cells. Depletion of endogenous PS abolishes K-Ras G12V PM curvature sensing. In cells and synthetic bilayers, only mixed-chain PS species, but not other PS species tested, mediate K-Ras G12V membrane curvature sensing. Thus, K-Ras nanoclusters act as relay stations to convert mechanical perturbations to mitogenic signaling.

Membrane Recognition and Binding by the Phosphatidylinositol Phosphate Kinase PIP5K1A: A Multiscale Simulation Study

Phosphatidylinositol phosphates (PIPs) are lipid signaling molecules that play key roles in many cellular processes. PIP5K1A kinase catalyzes phosphorylation of PI4P to form PIP₂, which in turn interacts with membrane and membrane-associated proteins. We explore the mechanism of membrane binding by the PIP5K1A kinase using a multiscale molecular dynamics approach. Coarse-grained simulations show binding of monomeric PIP5K1A to a model cell membrane containing PI4P. PIP5K1A did not bind to zwitterionic or anionic membranes lacking PIP molecules. Initial encounter of kinase and bilayer was followed by reorientation to enable productive binding to the PI4P-containing membrane. The simulations suggest that unstructured regions may be important for the preferred orientation for membrane binding. Atomistic simulations indicated that the dimeric kinase could not bind to the membrane via both active sites at the same time, suggesting a conformational change in the protein and/or bilayer distortion may be needed for dual-site binding to occur.

•

PIP5K1A kinase interacts with PIP-containing membranes via its activation loop

•

PIP5K1A does not bind to zwitterionic or anionic membranes lacking PIP molecules

•

Initial encounter of protein and bilayer is followed by reorientation and binding

•

Dimeric PIP5K1A binds with membrane contacts via only one catalytic site at a time