Membrane-associated Ras dimers are isoform-specific: K-Ras dimers differ from H-Ras dimers

Conditional Cooperativity in RAS Assembly Pathways on Nanodiscs and Altered GTPase Cycling

Self-assemblies (i.e., nanoclusters) of the RAS GTPase on the membrane act as scaffolds that activate downstream RAF kinases and drive MAPK signaling for cell proliferation and tumorigenesis. However, the mechanistic details of nanoclustering remain largely unknown. Here, size-tunable nanodisc platforms and paramagnetic relaxation enhancement (PRE) analyses revealed the structural basis of the cooperative assembly processes of fully processed KRAS, mutated in a quarter of human cancers. The cooperativity is modulated by the mutation and nucleotide states of KRAS and the lipid composition of the membrane. Notably, the oncogenic mutants assemble in nonsequential pathways with two mutually cooperative 'α/α' and 'α/β' interfaces, while α/α dimerization of wild-type KRAS promotes the secondary α/β interaction sequentially. Mutation-based interface engineering was used to selectively trap the oligomeric intermediates of KRAS and probe their favorable interface interactions. Transiently exposed interfaces were available for the assembly. Real-time NMR demonstrated that higher-order oligomers retain higher numbers of active GTP-bound protomers in KRAS GTPase cycling. These data provide a deeper understanding of the nanocluster-enhanced signaling in response to the environment. Furthermore, our methodology is applicable to assemblies of many other membrane GTPases and lipid nanoparticle-based formulations of stable protein oligomers with enhanced cooperativity.

Membrane-Driven Dimerization of the Peripheral Membrane Protein KRAS: Implications for Downstream Signaling

Transient homo-dimerization of the RAS GTPase at the plasma membrane has been shown to promote the mitogen-activated protein kinase (MAPK) signaling pathway essential for cell proliferation and oncogenesis. To date, numerous crystallographic studies have focused on the well-defined GTPase domains of RAS isoforms, which lack the disordered C-terminal membrane anchor, thus providing limited structural insight into membrane-bound RAS molecules. Recently, lipid-bilayer nanodisc platforms and paramagnetic relaxation enhancement (PRE) analyses have revealed several distinct structures of the membrane-anchored homodimers of KRAS, an isoform that is most frequently mutated in human cancers. The KRAS dimerization interface is highly plastic and altered by biologically relevant conditions, including oncogenic mutations, the nucleotide states of the protein, and the lipid composition. Notably, PRE-derived structures of KRAS homodimers on the membrane substantially differ in terms of the relative orientation of the protomers at an "α-α" dimer interface comprising two α4-α5 regions. This interface plasticity along with the altered orientations of KRAS on the membrane impact the accessibility of KRAS to downstream effectors and regulatory proteins. Further, nanodisc platforms used to drive KRAS dimerization can be used to screen potential anticancer drugs that target membrane-bound RAS dimers and probe their structural mechanism of action.

Allosteric Activation of RhoA Complexed with p115-RhoGEF Deciphered by Conformational Dynamics

The Ras homologue family member A (RhoA) is a member of the Rho family, a subgroup of the Ras superfamily. RhoA interacts with the 115 kDa guanine nucleotide exchange factor (p115-RhoGEF), which assists in activation and binding with downstream effectors. Here, we use molecular dynamics (MD) simulations and essential dynamics analysis of the inactive RhoA-GDP and active RhoA-GTP, when bound to p115-RhoGEF to decipher the mechanism of RhoA activation at the structural level. We observe that inactive RhoA-GDP maintains its position near the catalytic site on the Dbl homology (DH) domain of p115-RhoGEF through the interaction of its Switch I region with the DH domain. We further show that the active RhoA-GTP is engaged in more interactions with the p115-RhoGEF membrane-bound Pleckstrin homology (PH) domain as compared to RhoA-GDP. We hypothesize that the role of the interactions between the active RhoA-GTP and the PH domain is to help release it from the DH domain upon activation. Our results support this premise, and our simulations uncover the beginning of this process and provide structural details. They also point to allosteric communication pathways that take part in RhoA activation to promote and strengthen the interaction between the active RhoA-GTP and the PH domain. Allosteric regulation also occurs among other members of the Rho superfamily. Collectively, we suggest that in the activation process, the role of the RhoA-GTP interaction with the PH domain is to release RhoA-GTP from the DH domain after activation, making it available to downstream effectors.

Structural insights into the complex of oncogenic KRas4BG12V and Rgl2, a RalA/B activator

About a quarter of total human cancers carry mutations in Ras isoforms. Accumulating evidence suggests that small GTPases, RalA, and RalB, and their activators, Ral guanine nucleotide exchange factors (RalGEFs), play an essential role in oncogenic Ras-induced signalling. We studied the interaction between human KRas4B and the Ras association (RA) domain of Rgl2 (Rgl2RA), one of the RA-containing RalGEFs. We show that the G12V oncogenic KRas4B mutation changes the interaction kinetics with Rgl2RA The crystal structure of the KRas4BG12V: Rgl2RA complex shows a 2:2 heterotetramer where the switch I and switch II regions of each KRasG12V interact with both Rgl2RA molecules. This structural arrangement is highly similar to the HRasE31K:RALGDSRA crystal structure and is distinct from the well-characterised Ras:Raf complex. Interestingly, the G12V mutation was found at the dimer interface of KRas4BG12V with its partner. Our study reveals a potentially distinct mode of Ras:effector complex formation by RalGEFs and offers a possible mechanistic explanation for how the oncogenic KRas4BG12V hyperactivates the RalA/B pathway.

Cancer and Autism: How PTEN Mutations Degrade Function at the Membrane and Isoform Expression in the Human Brain

Mutations causing loss of PTEN lipid phosphatase activity can promote cancer, benign tumors (PHTS), and neurodevelopmental disorders (NDDs). Exactly how they preferentially trigger distinct phenotypic outcomes has been puzzling. Here, we demonstrate that PTEN mutations differentially allosterically bias P loop dynamics and its connection to the catalytic site, affecting catalytic activity. NDD-related mutations are likely to sample conformations of the functional wild-type state, while sampled conformations for the strong, cancer-related driver mutation hotspots favor catalysis-primed conformations, suggesting that NDD mutations are likely to be weaker, and our large-scale simulations show why. Prenatal PTEN isoform expression data suggest exons 5 and 7, which harbor NDD mutations, as cancer-risk carriers. Since cancer requires more than a single mutation, our conformational and genomic analysis helps discover how same protein mutations can foster different clinical manifestations, articulates a role for co-occurring background latent driver mutations, and uncovers relationships of splicing isoform expression to life expectancy.

Crystal Packing Reveals a Potential Autoinhibited KRAS Dimer Interface and a Strategy for Small-Molecule Inhibition of RAS Signaling

KRAS GTPases harbor oncogenic mutations in more than 25% of human tumors. KRAS is considered to be largely undruggable due to the lack of a suitable small-molecule binding site. Here, we report a unique crystal structure of His-tagged KRASG12D that reveals a remarkable conformational change. The Switch I loop of one His-KRASG12D structure extends into the Switch I/II pocket of another His-KRASG12D in an adjacent unit cell to create an elaborate interface that is reminiscent of high-affinity protein-protein complexes. We explore the contributions of amino acids at this interface using alanine-scanning studies with alchemical free energy perturbation calculations based on explicit-solvent molecular dynamics simulations. Several interface amino acids were found to be hot spots as they contributed more than 1.5 kcal/mol to the protein-protein interaction. Computational analysis of the complex revealed the presence of two large binding pockets that possess physicochemical features typically found in pockets considered druggable. Small-molecule binding to these pockets may stabilize this autoinhibited structure of KRAS if it exists in cells to provide a new strategy to inhibit RAS signaling.

Protein conformational ensembles in function: roles and mechanisms

The sequence-structure-function paradigm has dominated twentieth century molecular biology. The paradigm tacitly stipulated that for each sequence there exists a single, well-organized protein structure. Yet, to sustain cell life, function requires (i) that there be more than a single structure, (ii) that there be switching between the structures, and (iii) that the structures be incompletely organized. These fundamental tenets called for an updated sequence-conformational ensemble-function paradigm. The powerful energy landscape idea, which is the foundation of modernized molecular biology, imported the conformational ensemble framework from physics and chemistry. This framework embraces the recognition that proteins are dynamic and are always interconverting between conformational states with varying energies. The more stable the conformation the more populated it is. The changes in the populations of the states are required for cell life. As an example, in vivo, under physiological conditions, wild type kinases commonly populate their more stable "closed", inactive, conformations. However, there are minor populations of the "open", ligand-free states. Upon their stabilization, e.g., by high affinity interactions or mutations, their ensembles shift to occupy the active states. Here we discuss the role of conformational propensities in function. We provide multiple examples of diverse systems, including protein kinases, lipid kinases, and Ras GTPases, discuss diverse conformational mechanisms, and provide a broad outlook on protein ensembles in the cell. We propose that the number of molecules in the active state (inactive for repressors), determine protein function, and that the dynamic, relative conformational propensities, rather than the rigid structures, are the hallmark of cell life.

Direct Observation of Membrane-Associated H-Ras in the Native Cellular Environment by In-Cell 19F-NMR Spectroscopy

Ras acts as a molecular switch to control intracellular signaling on the plasma membrane (PM). Elucidating how Ras associates with PM in the native cellular environment is crucial for understanding its control mechanism. Here, we used in-cell nuclear magnetic resonance (NMR) spectroscopy combined with site-specific ¹⁹F-labeling to explore the membrane-associated states of H-Ras in living cells. The site-specific incorporation of p-trifluoromethoxyphenylalanine (OCF₃Phe) at three different sites of H-Ras, i.e., Tyr32 in switch I, Tyr96 interacting with switch II, and Tyr157 on helix α5, allowed the characterization of their conformational states depending on the nucleotide-bound states and an oncogenic mutational state. Exogenously delivered ¹⁹F-labeled H-Ras protein containing a C-terminal hypervariable region was assimilated via endogenous membrane-trafficking, enabling proper association with the cell membrane compartments. Despite poor sensitivity of the in-cell NMR spectra of membrane-associated H-Ras, the Bayesian spectral deconvolution identified distinct signal components on three ¹⁹F-labeled sites, thus offering the conformational multiplicity of H-Ras on the PM. Our study may be helpful in elucidating the atomic-scale picture of membrane-associated proteins in living cells.

The Self-Association of the KRAS4b Protein is Altered by Lipid-Bilayer Composition and Electrostatics

KRAS is a peripheral membrane protein that regulates multiple signaling pathways, and is mutated in ≈30 % of cancers. Transient self-association of KRAS is essential for activation of the downstream effector RAF and oncogenicity. The presence of anionic phosphatidylserine (PS) lipids in the membrane was shown to promote KRAS self-assembly, however, the structural mechanisms remain elusive. Here, we employed nanodisc bilayers of defined lipid compositions, and probed the impact of PS concentration on KRAS self-association. Paramagnetic NMR experiments demonstrated the existence of two transient dimer conformations involving alternate electrostatic contacts between R135 and either D153 or E168 on the "α4/5-α4/5" interface, and revealed that lipid composition and salt modulate their dynamic equilibrium. These dimer interfaces were validated by charge-reversal mutants. This plasticity demonstrates how the dynamic KRAS dimerization interface responds to the environment, and likely extends to the assembly of other signaling complexes on the membrane.

How PTEN mutations degrade function at the membrane and life expectancy of carriers of mutations in the human brain

PTEN dysfunction, caused by loss of lipid phosphatase activity or deletion, promotes pathologies, cancer, benign tumors, and neurodevelopmental disorders (NDDs). Despite efforts, exactly how the mutations trigger distinct phenotypic outcomes, cancer or NDD, has been puzzling. It has also been unclear how to distinguish between mutations harbored by isoforms, are they cancer or NDDs-related. Here we address both. We demonstrate that PTEN mutations differentially allosterically bias P-loop dynamics and its connection to the catalytic site, affecting catalytic activity. NDD-related mutations are likely to sample conformations present in the wild-type, while sampled conformations sheltering cancer-related hotspots favor catalysis-prone conformations, suggesting that NDD mutations are weaker. Analysis of isoform expression data indicates that if the transcript has NDD-related mutations, alone or in combination with cancer hotspots, there is high prenatal expression. If no mutations within the measured days, low expression levels. Cancer mutations promote stronger signaling and cell proliferation; NDDs' are weaker, influencing brain cell differentiation. Further, exon 5 is impacted by NDD or non-NDD mutations, while exon 7 is exclusively impacted by NDD mutations. Our comprehensive conformational and genomic analysis helps discover how same allele mutations can foster different clinical manifestations and uncovers correlations of splicing isoform expression to life expectancy.

Higher-order interactions of Bcr-Abl can broaden chronic myeloid leukemia (CML) drug repertoire

Bcr-Abl, a nonreceptor tyrosine kinase, is associated with leukemias, especially chronic myeloid leukemia (CML). Deletion of Abl's N-terminal region, to which myristoyl is linked, renders the Bcr-Abl fusion oncoprotein constitutively active. The substitution of Abl's N-terminal region by Bcr enables Bcr-Abl oligomerization. Oligomerization is critical: it promotes clustering on the membrane, which is essential for potent MAPK signaling and cell proliferation. Here we decipher the Bcr-Abl specific, step-by-step oligomerization process, identify a specific packing surface, determine exactly how the process is structured and identify its key elements. Bcr's coiled coil (CC) domain at the N-terminal controls Bcr-Abl oligomerization. Crystallography validated oligomerization via Bcr-Abl dimerization between two Bcr CC domains, with tetramerization via tight packing between two binary assemblies. However, the structural principles guiding Bcr CC domain oligomerization are unknown, hindering mechanistic understanding and drugs exploiting it. Using molecular dynamics (MD) simulations, we determine that the binary complex of the Bcr CC domain serves as a basic unit in the quaternary complex providing a specific surface for dimer-dimer packing and higher-order oligomerization. We discover that the small α1-helix is the key. In the binary assembly, the helix forms interchain aromatic dimeric packing, and in the quaternary assembly, it contributes to the specific dimer-dimer packing. Our mechanism is supported by the experimental literature. It offers the key elements controlling this process which can expand the drug discovery strategy, including by Bcr CC-derived peptides, and candidate residues for small covalent drugs, toward quenching oligomerization, supplementing competitive and allosteric tyrosine kinase inhibitors.

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.

The structural basis of BCR-ABL recruitment of GRB2 in chronic myelogenous leukemia

BCR-ABL drives chronic myeloid leukemia (CML). BCR binding to GRB2 transduces signaling via the Ras/MAPK pathway. Despite considerable data confirming the binding, molecular-level understanding of exactly how the two proteins interact, and especially, what are the determinants of the specificity of the SH2^GRB2 domain-pBCR recognition are still open questions. Yet, this is vastly important for understanding binding selectivity, and for predicting the phosphorylated receptors, or peptides, that are likely to bind. Here, we uncover these determinants and ascertain to what extent they relate to the affinity of the interaction. Toward this end, we modeled the complexes of the phosphorylated BCR (pBCR) and SH2^GRB2 and other pY/Y-peptide-SH2 complexes and compared their specificity and affinity. We observed that pBCR's ₁₇₆FpYVNV₁₈₀ motif is favorable and specific to SH2^GRB2, similar to pEGFR, but not other complexes. SH2^GRB2 contains two binding pockets: pY-binding recognition pocket triggers binding, and the specificity pocket whose interaction is governed by N179 in pBCR and W121 in SH2^GRB2. Our proposed motif with optimal affinity to SH2^GRB2 is E/D-pY-E/V-N-I/L. Collectively, we provide the structural basis of BCR-ABL recruitment of GRB2, outline its specificity hallmarks, and delineate a blueprint for prediction of BCR-binding scaffolds and for therapeutic peptide design.

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.

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.

KARATE: PKA-induced KRAS4B-RHOA-mTORC2 supercomplex phosphorylates AKT in insulin signaling and glucose homeostasis

AKT is a serine/threonine kinase that plays an important role in metabolism, cell growth, and cytoskeletal dynamics. AKT is activated by two kinases, PDK1 and mTORC2. Although the regulation of PDK1 is well understood, the mechanism that controls mTORC2 is unknown. Here, by investigating insulin receptor signaling in human cells and biochemical reconstitution, we found that insulin induces the activation of mTORC2 toward AKT by assembling a supercomplex with KRAS4B and RHOA GTPases, termed KARATE (KRAS4B-RHOA-mTORC2 Ensemble). Insulin-induced KARATE assembly is controlled via phosphorylation of GTP-bound KRAS4B at S181 and GDP-bound RHOA at S188 by protein kinase A. By developing a KARATE inhibitor, we demonstrate that KRAS4B-RHOA interaction drives KARATE formation. In adipocytes, KARATE controls insulin-dependent translocation of the glucose transporter GLUT4 to the plasma membrane for glucose uptake. Thus, our work reveals a fundamental mechanism that activates mTORC2 toward AKT in insulin-regulated glucose homeostasis.

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.

A structural model of a Ras-Raf signalosome

The protein K-Ras functions as a molecular switch in signaling pathways regulating cell growth. In the human mitogen-activated protein kinase (MAPK) pathway, which is implicated in many cancers, multiple K-Ras proteins are thought to assemble at the cell membrane with Ras effector proteins from the Raf family. Here we propose an atomistic structural model for such an assembly. Our starting point was an asymmetric guanosine triphosphate-mediated K-Ras dimer model, which we generated using unbiased molecular dynamics simulations and verified with mutagenesis experiments. Adding further K-Ras monomers in a head-to-tail fashion led to a compact helical assembly, a model we validated using electron microscopy and cell-based experiments. This assembly stabilizes K-Ras in its active state and presents composite interfaces to facilitate Raf binding. Guided by existing experimental data, we then positioned C-Raf, the downstream kinase MEK1 and accessory proteins (Galectin-3 and 14-3-3σ) on and around the helical assembly. The resulting Ras-Raf signalosome model offers an explanation for a large body of data on MAPK signaling.

RAS Dimers: The Novice Couple at the RAS-ERK Pathway Ball

Signals conveyed through the RAS-ERK pathway constitute a pivotal regulatory element in cancer-related cellular processes. Recently, RAS dimerization has been proposed as a key step in the relay of RAS signals, critically contributing to RAF activation. RAS clustering at plasma membrane microdomains and endomembranes facilitates RAS dimerization in response to stimulation, promoting RAF dimerization and subsequent activation. Remarkably, inhibiting RAS dimerization forestalls tumorigenesis in cellular and animal models. Thus, the pharmacological disruption of RAS dimers has emerged as an additional target for cancer researchers in the quest for a means to curtail aberrant RAS activity.

Ras isoform-specific expression, chromatin accessibility, and signaling

The anchorage of Ras isoforms in the membrane and their nanocluster formations have been studied extensively, including their detailed interactions, sizes, preferred membrane environments, chemistry, and geometry. However, the staggering challenge of their epigenetics and chromatin accessibility in distinct cell states and types, which we propose is a major factor determining their specific expression, still awaits unraveling. Ras isoforms are distinguished by their C-terminal hypervariable region (HVR) which acts in intracellular transport, regulation, and membrane anchorage. Here, we review some isoform-specific activities at the plasma membrane from a structural dynamic standpoint. Inspired by physics and chemistry, we recognize that understanding functional specificity requires insight into how biomolecules can organize themselves in different cellular environments. Within this framework, we suggest that isoform-specific expression may largely be controlled by the chromatin density and physical compaction, which allow (or curb) access to "chromatinized DNA." Genes are preferentially expressed in tissues: proteins expressed in pancreatic cells may not be equally expressed in lung cells. It is the rule-not an exception, and it can be at least partly understood in terms of chromatin organization and accessibility state. Genes are expressed when they can be sufficiently exposed to the transcription machinery, and they are less so when they are persistently buried in dense chromatin. Notably, chromatin accessibility can similarly determine expression of drug resistance genes.

Novel MAPK/AKT-impairing germline NRAS variant identified in a melanoma-prone family

While several high-penetrance melanoma risk genes are known, variation in these genes fail to explain melanoma susceptibility in a large proportion of high-risk families. As part of a melanoma family sequencing study, including 435 families from Mediterranean populations we identified a novel NRAS variant (c.170A > C, p.D57A) in an Italian melanoma-prone family. This variant is absent in exomes in gnomAD, ESP, UKBiobank, and the 1000 Genomes Project, as well as in 11,273 Mediterranean individuals and 109 melanoma-prone families from the US and Australia. This variant occurs in the GTP-binding pocket of NRAS. Differently from other RAS activating alterations, NRAS D57A expression is unable to activate MAPK-pathway both constitutively and after stimulation but enhances EGF-induced PI3K-pathway signaling in serum starved conditions in vitro. Consistent with in vitro data demonstrating that NRAS D57A does not enrich GTP binding, molecular modeling suggests that the D57A substitution would be expected to impair Mg2 + binding and decrease nucleotide-binding and GTPase activity of NRAS. While we cannot firmly establish NRAS c.170A > C (p.D57A) as a melanoma susceptibility variant, further investigation of NRAS as a familial melanoma gene is warranted.

The structural basis of Akt PH domain interaction with calmodulin

Akt plays a key role in the Ras/PI3K/Akt/mTOR signaling pathway. In breast cancer, Akt translocation to the plasma membrane is enabled by the interaction of its pleckstrin homology domain (PHD) with calmodulin (CaM). At the membrane, the conformational change promoted by PIP3 releases CaM and facilitates Thr308 and Ser473 phosphorylation and activation. Here, using modeling and molecular dynamics simulations, we aim to figure out how CaM interacts with Akt's PHD at the atomic level. Our simulations show that CaM-PHD interaction is thermodynamically stable and involves a β-strand rather than an α-helix, in agreement with NMR data, and that electrostatic and hydrophobic interactions are critical. The PHD interacts with CaM lobes; however, multiple modes are possible. IP4, the polar head of PIP3, weakens the CaM-PHD interaction, implicating the release mechanism at the plasma membrane. Recently, we unraveled the mechanism of PI3Kα activation at the atomistic level and the structural basis for Ras role in the activation. Here, our atomistic structural data clarify the mechanism of how CaM interacts, delivers, and releases Akt-the next node in the Ras/PI3K pathway-at the plasma membrane.

Mechanistic Differences of Activation of Rac1P29S and Rac1A159V

Rac1 is a small GTPase that plays key roles in actin reorganization, cell motility, and cell survival/growth as well as in various cancer types and neurodegenerative diseases. Similar to other Ras superfamily GTPases, Rac1 switches between active GTP-bound and inactive GDP-bound states. Switch I and II regions open and close during GDP/GTP exchange. P29S and A159V (paralogous to K-RasA146) mutations are the two most common somatic mutations of Rac1. Rac1P29S is a known hotspot for melanoma, whereas Rac1A159V most commonly occurs in head and neck cancer. To investigate how these substitutions induce the Rac1 dynamics, we used atomistic molecular dynamics simulations on the wild-type Rac1 and two mutant systems (P29S and A159V) in the GTP bound state, and on the wild-type Rac1 and P29S mutated system in the GDP bound state. Here, we show that P29S and A159V mutations activate Rac1 with different mechanisms. In Rac1P29S-GTP, the substitution increases the flexibility of Switch I based on RMSF and dihedral angle calculations and leads to an open conformation. We propose that the open Switch I conformation is one of the underlying reasons for rapid GDP/GTP exchange of Rac1P29S. On the other hand, in Rac1A159V-GTP, some of the contacts of the guanosine ring of GTP with Rac1 are temporarily lost, enabling the guanosine ring to move toward Switch I and subsequently close the switch. Rac1A159V-GTP adopts a Ras state 2 like conformation, where both switch regions are in closed conformation and Thr35 forms a hydrogen bond with the nucleotide.

The mechanism of full activation of tumor suppressor PTEN at the phosphoinositide-enriched membrane

Tumor suppressor PTEN, the second most highly mutated protein in cancer, dephosphorylates signaling lipid PIP₃ produced by PI3Ks. Excess PIP₃ promotes cell proliferation. The mechanism at the membrane of this pivotal phosphatase is unknown hindering drug discovery. Exploiting explicit solvent simulations, we tracked full-length PTEN trafficking from the cytosol to the membrane. We observed its interaction with membranes composed of zwitterionic phosphatidylcholine, anionic phosphatidylserine, and phosphoinositides, including signaling lipids PIP₂ and PIP₃. We tracked its moving away from the zwitterionic and getting absorbed onto anionic membrane that harbors PIP₃. We followed it localizing on microdomains enriched in signaling lipids, as PI3K does, and observed PIP₃ allosterically unfolding the N-terminal PIP₂ binding domain, positioning it favorably for the polybasic motif interaction with PIP₂. Finally, we determined PTEN catalytic action at the membrane, all in line with experimental observations, deciphering the mechanisms of how PTEN anchors to the membrane and restrains cancer.

•

PTEN localizes on membrane microdomains enriched in phosphoinositides, as PI3K does

•

Full PTEN activation requires both signaling lipids, PIP₂ and PIP₃

•

Strong salt bridge interactions sustain stable PTEN membrane localization

•

Substrate-induced P loop conformational change implicates PTEN catalytic activity

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.

Active and Inactive Cdc42 Differ in Their Insert Region Conformational Dynamics

Cell division control protein 42 homolog (Cdc42) protein, a Ras superfamily GTPase, regulates cellular activities, including cancer progression. Using all-atom molecular dynamics (MD) simulations and essential dynamic analysis, we investigated the structure and dynamics of the catalytic domains of GDP-bound (inactive) and GTP-bound (active) Cdc42 in solution. We discovered substantial differences in the dynamics of the inactive and active forms, particularly in the "insert region" (residues 122-135), which plays a role in Cdc42 activation and binding to effectors. The insert region has larger conformational flexibility in the GDP-bound Cdc42 than in the GTP-bound Cdc42. The G2 loop and switch I at the effector lobe of the catalytic domain exhibit large conformational changes in both the GDP- and the GTP-bound systems, but in the GTP-bound Cdc42, the switch I interactions with GTP are retained. Oncogenic mutations were identified in the Ras superfamily. In Cdc42, the G12V and Q61L mutations decrease the GTPase activity. We simulated these mutations in both GDP- and GTP-bound Cdc42. Although the overall structural organization is quite similar between the wild type and the mutants, there are small differences in the conformational dynamics, especially in the two switch regions. Taken together, the G12V and Q61L mutations may play a role similar to their K-Ras counterparts in nucleotide binding and activation. The conformational differences, which are mainly in the insert region and, to a lesser extent, in the switch regions flanking the nucleotide binding site, can shed light on binding and activation. We propose that the differences are due to a network of hydrogen bonds that gets disrupted when Cdc42 is bound to GDP, a disruption that does not exist in other Rho GTPases. The differences in the dynamics between the two Cdc42 states suggest that the inactive conformation has reduced ability to bind to effectors.

Hetero-oligomerization of Rho and Ras GTPases Connects GPCR Activation to mTORC2-AKT Signaling

The activation of G-protein-coupled receptors (GPCRs) leads to the activation of mTORC2 in cell migration and metabolism. However, the mechanism that links GPCRs to mTORC2 remains unknown. Here, using Dictyostelium cells, we show that GPCR-mediated chemotactic stimulation induces hetero-oligomerization of phosphorylated GDP-bound Rho GTPase and GTP-bound Ras GTPase in directed cell migration. The Rho-Ras hetero-oligomers directly and specifically stimulate mTORC2 activity toward AKT in cells and after biochemical reconstitution using purified proteins in vitro. The Rho-Ras hetero-oligomers do not activate ERK/MAPK, another kinase that functions downstream of GPCRs and Ras. Human KRas4B functionally replace Dictyostelium Ras in mTORC2 activation. In contrast to GDP-Rho, GTP-Rho antagonizes mTORC2-AKT signaling by inhibiting the oligomerization of GDP-Rho with GTP-Ras. These data reveal that GPCR-stimulated hetero-oligomerization of Rho and Ras provides a critical regulatory step that controls mTORC2-AKT signaling.

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.

SOS1 interacts with Grb2 through regions that induce closed nSH3 conformations

Grb2 is an adaptor protein connecting the epidermal growth factor receptor and the downstream Son of sevenless 1 (SOS1), a Ras-specific guanine nucleotide exchange factor (RasGEF), which exchanges GDP by GTP. Grb2 contains three SH domains: N-terminal SH3 (nSH3), SH2, and C-terminal SH3 (cSH3). The C-terminal proline-rich (PR) domain of SOS1 regulates nSH3 open/closed conformations. Earlier, several nSH3 binding motifs were identified in the PR domain. More recently, we characterized by nuclear magnetic resonance and replica exchange simulations possible cSH3 binding regions. Among them, we discovered a cSH3-specific binding region. However, how PR binding at these sites regulates the nSH3/cSH3 conformation has been unclear. Here, we explore the nSH3/cSH3 interaction with linked and truncated PR segments using molecular dynamics simulations. Our 248 μs simulations include 620 distinct trajectories, each 400 ns. We construct the effective free energy landscape to validate the nSH3/cSH3 binding sites. The nSH3/cSH3-SOS1 peptide complex models indicate that strong peptide binders attract the flexible nSH3 n-Src loop, inducing a closed conformation of nSH3; by contrast, the cSH3 conformation remains unchanged. Inhibitors that disrupt the Ras-SOS1 interaction have been designed; the conformational details uncovered here may assist in the design of polypeptides inhibiting Grb2-SOS1 interaction, thus SOS1 recruitment to the membrane where Ras resides.

Nucleotide-Specific Autoinhibition of Full-Length K-Ras4B Identified by Extensive Conformational Sampling

K-Ras is one of the most frequently mutated oncogenes in human tumor cells. It consists of a well-conserved globular catalytic domain and a flexible tail-like hypervariable region (HVR) at its C-terminal end. It plays a key role in signaling networks in proliferation, differentiation, and survival, undergoing a conformational switch between the active and inactive states. It is regulated through the GDP-GTP cycle of the inactive GDP-bound and active GTP-bound states. Here, without imposing any prior constraints, we mapped the interaction pattern between the catalytic domain and the HVR using Molecular Dynamics with excited Normal Modes (MDeNM) starting from an initially extended HVR conformation for both states. Our sampling captured similar interaction patterns in both GDP- and GTP-bound states with shifted populations depending on the bound nucleotide. In the GDP-bound state, the conformations where the HVR interacts with the effector lobe are more populated than in the GTP-bound state, forming a buried thus autoinhibited catalytic site; in the GTP-bound state conformations where the HVR interacts with the allosteric lobe are more populated, overlapping the α3/α4 dimerization interface. The interaction of the GTP with Switch I and Switch II is stronger than that of the GDP in line with a decrease in the fluctuation upon GTP binding.

Two Distinct Structures of Membrane-Associated Homodimers of GTP- and GDP-Bound KRAS4B Revealed by Paramagnetic Relaxation Enhancement

KRAS homo-dimerization has been implicated in the activation of RAF kinases, however, the mechanism and structural basis remain elusive. We developed a system to study KRAS dimerization on nanodiscs using paramagnetic relaxation enhancement (PRE) NMR spectroscopy, and determined distinct structures of membrane-anchored KRAS dimers in the active GTP- and inactive GDP-loaded states. Both dimerize through an α4-α5 interface, but the relative orientation of the protomers and their contacts differ substantially. Dimerization of KRAS-GTP, stabilized by electrostatic interactions between R135 and E168, favors an orientation on the membrane that promotes accessibility of the effector-binding site. Remarkably, "cross"-dimerization between GTP- and GDP-bound KRAS molecules is unfavorable. These models provide a platform to elucidate the structural basis of RAF activation by RAS and to develop inhibitors that can disrupt the KRAS dimerization. The methodology is applicable to many other farnesylated small GTPases.

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.

The Mystery of Rap1 Suppression of Oncogenic Ras

Decades ago, Rap1, a small GTPase very similar to Ras, was observed to suppress oncogenic Ras phenotype, reverting its transformation. The proposed reason, persisting since, has been competition between Ras and Rap1 for a common target. Yet, none was found. There was also Rap1's puzzling suppression of Raf-1 versus activation of BRAF. Reemerging interest in Rap1 envisages capturing its Ras suppression action by inhibitors. Here, we review the literature and resolve the enigma. In vivo oncogenic Ras exists in isoform-distinct nanoclusters. The presence of Rap1 within the nanoclusters reduces the number of the clustered oncogenic Ras molecules, thus suppressing Raf-1 activation and mitogen-activated protein kinase (MAPK) signaling. Nanoclustering suggests that Rap1 suppression is Ras isoform dependent. Altogether, a potent Rap1-like inhibitor appears unlikely.

The quaternary assembly of KRas4B with Raf-1 at the membrane

Proximally located in the membrane, oncogenic Ras dimers (or nanoclusters) can recruit and promote Raf dimerization and MAPK (Raf/MEK/ERK) signaling. Among Ras isoforms, KRas4B is the most frequently mutated. Recent data on the binary KRas4B–Raf-1 complex suggested that Raf-1 CRD not only executes membrane anchorage, but also supports the high-affinity interaction of Raf-1 RBD with KRas4B catalytic domain. For a detailed mechanistic picture of Raf activation at the membrane, we employ explicit MD simulations of the quaternary KRas4B–Raf-1 complex. The complex contains two active GTP-bound KRas4B proteins forming a dimer through the allosteric lobe interface and two tandem RBD-CRD segments of Raf-1 interacting with the effector lobes at both ends of the KRas4B dimer. We show that Raf-1 RBD-CRD supports stable KRas4B dimer at preferred interface and orientation at the membrane, thereby cooperatively enhancing the affinity of the KRas4B–Raf-1 interaction. We propose that a Ras dimer at the membrane can increase the population of proximal Raf kinase domains, promoting kinase domain dimerization in the cytoplasm. Collectively, the dynamic Ras–Raf assembly promotes Raf activation not by allostery; instead, Ras activates Raf by shifting its ensemble toward kinase domain-accessible states through enhanced affinity at the membrane.

Presence or Absence of Ras Dimerization Shows Distinct Kinetic Signature in Ras-Raf Interaction

Initiations of cell signaling pathways often occur through the formation of multiprotein complexes that form through protein-protein interactions. Therefore, detecting their presence is central to understanding the function of a cell signaling pathway, aberration of which often leads to fatal diseases, including cancers. However, the multiprotein complexes are often difficult to detect using microscopes due to their small sizes. Therefore, currently, their presence can be only detected through indirect means. In this article, we propose to investigate the presence or absence of protein complexes through some easily measurable kinetic parameters, such as activation rates. As a proof of concept, we investigate the Ras-Raf system, a well-characterized cell signaling system. It has been hypothesized that Ras dimerization is necessary to create activated Raf dimers. Although there are circumstantial evidences supporting the Ras dimerization hypothesis, direct proof of Ras dimerization is still inconclusive. In the absence of conclusive direct experimental proof, this hypothesis can only be examined through indirect evidences of Ras dimerization. In this article, using a multiscale simulation technique, we provide multiple criteria that distinguishes an activation mechanism involving Ras dimerization from another mechanism that does not involve Ras dimerization. The provided criteria will be useful in the investigation of not only Ras-Raf interaction but also other two-protein interactions.

Oncogenic K-Ras4B Dimerization Enhances Downstream Mitogen-activated Protein Kinase Signaling

Ras recruits and activates effectors that transmit receptor-initiated signals. Monomeric Ras can bind Raf; however, Raf's activation requires dimerization, which can be facilitated by Ras dimerization. Previously, we showed that active K-Ras4B dimerizes in silico and in vitro through two major interfaces: (i) β-interface, mapped to Switch I and effector-binding regions, (ii) α-interface at the allosteric lobe. Here, we chose constitutively active K-Ras4B as our control and two double mutants (K101D and R102E; and R41E and K42D) in the α- and β-interfaces. Two of the mutations are from The Cancer Genome Atlas (TCGA) and the Catalogue Of Somatic Mutations In Cancer (COSMIC) data sets. R41 and R102 are found in several adenocarcinomas in Ras isoforms. We performed site-directed mutagenesis, cellular localization experiments, and molecular dynamics (MD) simulations to assess the impact of the mutations on K-Ras4B dimerization and function. α-interface K101D/R102E double mutations reduced dimerization but only slightly reduced downstream phosphorylated extracellular signal-regulated kinase (ERK) (pERK) levels. While β-interface R41E/K42D double mutations did not interfere with dimerization, they almost completely blocked K-Ras4B-mediated ERK phosphorylation. Both double mutations increased downstream phosphorylated Akt (pAkt) levels in cells. Changes in pERK and pAkt levels altered ERK- and Akt-regulated gene expressions, such as EGR1, JUN, and BCL2L11. These results underscore the role of the α-interface in K-Ras4B homodimerization and the β-surface in effector binding. MD simulations highlight that the membrane and hypervariable region (HVR) interact with both α- and β-interfaces of K-Ras4B mutants, respectively, inhibiting homodimerization and probably effector binding. Mutations at both interfaces interfered with mitogen-activated protein kinase (MAPK) and phosphoinositide 3-kinase signaling but in different forms and extents. We conclude that dimerization is not necessary but enhances downstream MAPK signaling.

KRAS: A Promising Therapeutic Target for Cancer Treatment

Kirsten rat sarcoma 2 viral oncogene homolog (KRAS) is the most commonly mutated oncogene in human cancer. The developments of many cancers depend on sustained expression and signaling of KRAS, which makes KRAS a high-priority therapeutic target. Scientists have not successfully developed drugs that target KRAS, although efforts have been made last three decades. In this review, we highlight the emerging experimental strategies of impairing KRAS membrane localization and the direct targeting of KRAS. We also conclude the combinatorial therapies and RNA interference technology for the treatment of KRAS mutant cancers. Moreover, the virtual screening approach to discover novel KRAS inhibitors and synthetic lethality interactors of KRAS are discussed in detail.

The Structural Basis of the Farnesylated and Methylated KRas4B Interaction with Calmodulin

Ca²⁺-calmodulin (CaM) extracts KRas4B from the plasma membrane, suggesting that KRas4B/CaM interaction plays a role in regulating Ras signaling. To gain mechanistic insight, we provide a computational model, supported by experimental structural data, of farnesylated/methylated KRas4B_1-185 interacting with CaM in solution and at anionic membranes including signaling lipids. Due to multiple interaction modes, we observe diverse conformational ensembles of the KRas4B-CaM complex. A highly populated conformation reveals the catalytic domain interacting with the N-lobe and the hypervariable region (HVR) wrapping around the linker with the farnesyl docking to the extended CaM's C-lobe pocket. Alternatively, KRas4B can interact with collapsed CaM with the farnesyl penetrating CaM's center. At anionic membranes, CaM interacts with the catalytic domain with large fluctuations, drawing the HVR. Signaling lipids establishing strong salt bridges with CaM prevent membrane departure. Membrane-interacting KRas4B-CaM complex can productively recruit phosphatidylinositol 3-kinase α (PI3Kα) to the plasma membrane, serving as a coagent in activating PI3Kα/Akt signaling.

Phosphorylated Rho-GDP directly activates mTORC2 kinase towards AKT through dimerization with Ras-GTP to regulate cell migration

mTORC2 plays critical roles in metabolism, cell survival and actin cytoskeletal dynamics through the phosphorylation of AKT. Despite its importance to biology and medicine, it is unclear how mTORC2-mediated AKT phosphorylation is controlled. Here, we identify an unforeseen principle by which a GDP-bound form of the conserved small G protein Rho GTPase directly activates mTORC2 in AKT phosphorylation in social amoebae (Dictyostelium discoideum) cells. Using biochemical reconstitution with purified proteins, we demonstrate that Rho-GDP promotes AKT phosphorylation by assembling a supercomplex with Ras-GTP and mTORC2. This supercomplex formation is controlled by the chemoattractant-induced phosphorylation of Rho-GDP at S192 by GSK-3. Furthermore, Rho-GDP rescues defects in both mTORC2-mediated AKT phosphorylation and directed cell migration in Rho-null cells in a manner dependent on phosphorylation of S192. Thus, in contrast to the prevailing view that the GDP-bound forms of G proteins are inactive, our study reveals that mTORC2-AKT signalling is activated by Rho-GDP.

KRAS-specific inhibition using a DARPin binding to a site in the allosteric lobe

Inhibiting the RAS oncogenic protein has largely been through targeting the switch regions that interact with signalling effector proteins. Here, we report designed ankyrin repeat proteins (DARPins) macromolecules that specifically inhibit the KRAS isoform by binding to an allosteric site encompassing the region around KRAS-specific residue histidine 95 at the helix α3/loop 7/helix α4 interface. We show that these DARPins specifically inhibit KRAS/effector interactions and the dependent downstream signalling pathways in cancer cells. Binding by the DARPins at that region influences KRAS/effector interactions in different ways, including KRAS nucleotide exchange and inhibiting KRAS dimerization at the plasma membrane. These results highlight the importance of targeting the α3/loop 7/α4 interface, a previously untargeted site in RAS, for specifically inhibiting KRAS function.

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

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

Review: Precision medicine and driver mutations: Computational methods, functional assays and conformational principles for interpreting cancer drivers

At the root of the so-called precision medicine or precision oncology, which is our focus here, is the hypothesis that cancer treatment would be considerably better if therapies were guided by a tumor’s genomic alterations. This hypothesis has sparked major initiatives focusing on whole-genome and/or exome sequencing, creation of large databases, and developing tools for their statistical analyses—all aspiring to identify actionable alterations, and thus molecular targets, in a patient. At the center of the massive amount of collected sequence data is their interpretations that largely rest on statistical analysis and phenotypic observations. Statistics is vital, because it guides identification of cancer-driving alterations. However, statistics of mutations do not identify a change in protein conformation; therefore, it may not define sufficiently accurate actionable mutations, neglecting those that are rare. Among the many thematic overviews of precision oncology, this review innovates by further comprehensively including precision pharmacology, and within this framework, articulating its protein structural landscape and consequences to cellular signaling pathways. It provides the underlying physicochemical basis, thereby also opening the door to a broader community.

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

Ras genes potently drive human cancers, with mutated proto-oncogene GTPase KRAS4B (K-Ras4B) being the most abundant isoform. Targeted inhibition of oncogenic gene products is considered the "holy grail" of present-day cancer therapy, and recent discoveries of small-molecule KRas4B inhibitors were made thanks to a deeper understanding of the structure and dynamics of this GTPase. Because interactions with biological membranes are key for Ras function, Ras-lipid interactions have become a major focus, especially because such interactions evidently involve both the Ras C terminus for lipid anchoring and its G-protein domain. Here, using NMR spectroscopy and molecular dynamics simulations complemented by biophysical- and cell-biology assays, we investigated the interaction between K-Ras4B with the signaling lipid phosphatidylinositol (4,5)-phosphate (PIP2). We discovered that the β2 and β3 strands as well as helices 4 and 5 of the GTPase G-domain bind to PIP2 and identified the specific residues in these structural elements employed in these interactions, likely occurring in two K-Ras4B orientation states relative to the membrane. Importantly, we found that some of these residues known to be oncogenic when mutated (D47K, D92N, K104M, and D126N) are critical for K-Ras-mediated transformation of fibroblast cells, but do not substantially affect basal and assisted nucleotide hydrolysis and exchange. Moreover, the K104M substitution abolished localization of K-Ras to the plasma membrane. The findings suggest that specific G-domain residues can critically regulate Ras function by mediating interactions with membrane-associated PIP2 lipids; these insights that may inform the future design of therapeutic reagents targeting Ras activity.