Improved binding of raf to Ras.GDP is correlated with biological activity

Novel RAF-directed approaches to overcome current clinical limits and block the RAS/RAF node

Mutations in the RAS-RAF-MEK-ERK pathway are frequent alterations in cancer and RASopathies, and while RAS oncogene activation alone affects 19% of all patients and accounts for approximately 3.4 million new cases every year, less frequent alterations in the cascade's downstream effectors are also involved in cancer etiology. RAS proteins initiate the signaling cascade by promoting the dimerization of RAF kinases, which can act as oncoproteins as well: BRAFV600E is the most common oncogenic driver, mutated in the 8% of all malignancies. Research in this field led to the development of drugs that target the BRAFV600-like mutations (Class I), which are now utilized in clinics, but cause paradoxical activation of the pathway and resistance development. Furthermore, they are ineffective against non-BRAFV600E malignancies that dimerize and could be either RTK/RAS independent or dependent (Class II and III, respectively), which are still lacking an effective treatment. This review discusses the recent advances in anti-RAF therapies, including paradox breakers, dimer-inhibitors, immunotherapies, and other novel approaches, critically evaluating their efficacy in overcoming the therapeutic limitations, and their putative role in blocking the RAS pathway.

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.

Kinetic and thermodynamic allostery in the Ras protein family

Allostery, the transfer of information between distant parts of a macromolecule, is a fundamental feature of protein function and regulation. However, allosteric mechanisms are usually not explained by protein structure, requiring information on correlated fluctuations uniquely accessible to molecular simulation. Existing work to extract allosteric pathways from molecular dynamics simulations has focused on thermodynamic correlations. Here, we show how kinetic correlations encode complementary information essential to explain observed variations in allosteric regulation. We applied kinetic and thermodynamic correlation analysis on atomistic simulations of H, K, and NRas isoforms in the apo, GTP, and GDP-bound states of Ras protein, with and without complexing to its downstream effector, Raf. We show that switch I and switch II are the primary components of thermodynamic and kinetic allosteric networks, consistent with the key roles of these two motifs. These networks connect the switches to an allosteric loop recently discovered from a crystal structure of HRas. This allosteric loop is inactive in KRas, but is coupled to the hydrolysis arm switch II in NRas and HRas. We find that the mechanism in the latter two isoforms are thermodynamic and kinetic, respectively. Binding of Raf-RBD further activates thermodynamic allostery in HRas and KRas but has limited effect on NRas. These results indicate that kinetic and thermodynamic correlations are both needed to explain protein function and allostery. These two distinct channels of allosteric regulation, and their combinatorial variability, may explain how subtle mutational differences can lead to diverse regulatory profiles among enzymatic proteins.

Multiplexed imaging for probing RAS-RAF interactions in living cells

GTP-bound RAS interacts with its protein effectors in response to extracellular stimuli, leading to chemical inputs for downstream pathways. Significant progress has been made in measuring these reversible protein-protein interactions (PPIs) in various cell-free environments. Yet, acquiring high sensitivity in heterogeneous solutions remains challenging. Here, using an intermolecular fluorescence resonance energy transfer (FRET) biosensing approach, we develop a method to visualize and localize HRAS-CRAF interactions in living cells. We demonstrate that the EGFR activation and the HRAS-CRAF complex formation can be concurrently probed in a single cell. This biosensing strategy discriminates EGF-stimulated HRAS-CRAF interactions at the cell and organelle membranes. In addition, we provide quantitative FRET measurements for assessing these transient PPIs in a cell-free environment. Finally, we prove the utility of this approach by showing that an EGFR-binding compound is a potent inhibitor of HRAS-CRAF interactions. The outcomes of this work form a fundamental basis for further explorations of the spatiotemporal dynamics of various signaling networks.

Targeting Ras with protein engineering

Ras proteins are small GTPases that regulate cell growth and division. Mutations in Ras genes are associated with many types of cancer, making them attractive targets for cancer therapy. Despite extensive efforts, targeting Ras proteins with small molecules has been extremely challenging due to Ras's mostly flat surface and lack of small molecule-binding cavities. These challenges were recently overcome by the development of the first covalent small-molecule anti-Ras drug, sotorasib, highlighting the efficacy of Ras inhibition as a therapeutic strategy. However, this drug exclusively inhibits the Ras G12C mutant, which is not a prevalent mutation in most cancer types. Unlike the G12C variant, other Ras oncogenic mutants lack reactive cysteines, rendering them unsuitable for targeting via the same strategy. Protein engineering has emerged as a promising method to target Ras, as engineered proteins have the ability to recognize various surfaces with high affinity and specificity. Over the past few years, scientists have engineered antibodies, natural Ras effectors, and novel binding domains to bind to Ras and counteract its carcinogenic activities via a variety of strategies. These include inhibiting Ras-effector interactions, disrupting Ras dimerization, interrupting Ras nucleotide exchange, stimulating Ras interaction with tumor suppressor genes, and promoting Ras degradation. In parallel, significant advancements have been made in intracellular protein delivery, enabling the delivery of the engineered anti-Ras agents into the cellular cytoplasm. These advances offer a promising path for targeting Ras proteins and other challenging drug targets, opening up new opportunities for drug discovery and development.

Classification of GTP-dependent K-Ras4B active and inactive conformational states

Classifying reliably active and inactive molecular conformations of wildtype (WT) and mutated oncogenic proteins is a key, ongoing challenge in molecular cancer studies. Here, we probe the GTP-bound K-Ras4B conformational dynamics using long-time atomistic molecular dynamics (MD) simulations. We extract and analyze the detailed underlying free energy landscape of WT K-Ras4B. We use two key reaction coordinates, labeled d₁ and d₂ (i.e., distances coordinating the P_β atom of the GTP ligand with two key residues, T35 and G60), shown to correlate closely with activities of WT and mutated K-Ras4B. However, our new K-Ras4B conformational kinetics study reveals a more complex network of equilibrium Markovian states. We show that a new reaction coordinate is required to account for the orientation of acidic K-Ras4B sidechains such as D38 with respect to the interface with binding effector RAF1 and rationalize the activation/inactivation propensities and the corresponding molecular binding mechanisms. We use this understanding to unveil how a relatively conservative mutation (i.e., D33E, in the switch I region) can lead to significantly different activation propensities compared with WT K-Ras4B. Our study sheds new light on the ability of residues near the K-Ras4B-RAF1 interface to modulate the network of salt bridges at the binding interface with the RAF1 downstream effector and, thus, to influence the underlying GTP-dependent activation/inactivation mechanism. Altogether, our hybrid MD-docking modeling approach enables the development of new in silico methods for quantitative assessment of activation propensity changes (e.g., due to mutations or local binding environment). It also unveils the underlying molecular mechanisms and facilitates the rational design of new cancer drugs.

A conformation-selective protein binder for a KRAS mutant inhibits the interaction between RAS and RAF

Small GTPases are key signaling nodes that regulate the cellular processes and subcellular events, and their abnormal activities and dysregulations are closely linked with diverse cancers. Here, we report the development of conformation-selective protein binders for a KRAS mutant. The conformation-specific protein binders were selected from a repebody scaffold composed of LRR (Leucine-rich repeat) modules through phage display and modular engineering against constitute active conformation of KRAS. Epitope of the selected binders was mapped to be located close to switch I of KRAS. The conformation-selective protein binders were shown to effectively block the interaction between active KRAS and RAS-binding domain of BRAF, suppressing the KRAS-mediated downstream signaling.

Lifetime of actin-dependent protein nanoclusters

Protein nanoclusters (PNCs) are dynamic collections of a few proteins that spatially organize in nanometer-length clusters. PNCs are one of the principal forms of spatial organization of membrane proteins, and they have been shown or hypothesized to be important in various cellular processes, including cell signaling. PNCs show remarkable diversity in size, shape, and lifetime. In particular, the lifetime of PNCs can vary over a wide range of timescales. The diversity in size and shape can be explained by the interaction of the clustering proteins with the actin cytoskeleton or the lipid membrane, but very little is known about the processes that determine the lifetime of the nanoclusters. In this paper, using mathematical modeling of the cluster dynamics, we model the biophysical processes that determine the lifetime of actin-dependent PNCs. In particular, we investigated the role of actin aster fragmentation, which had been suggested to be a key determinant of the PNC lifetime, and we found that it is important only for a small class of PNCs. A simple extension of our model allowed us to investigate the kinetics of protein-ligand interaction near PNCs. We found an anomalous increase in the lifetime of ligands near PNCs, which agrees remarkably well with experimental data on RAS-RAF kinetics. In particular, analysis of the RAS-RAF data through our model provides falsifiable predictions and novel hypotheses that will not only shed light on the role of RAS-RAF kinetics in various cancers, but also will be useful in studying membrane protein clustering in general.

Titration of apparent in-cellula affinities of protein-protein interactions

A genetic assay permits simultaneous quantification of two interacting proteins and their bound fraction at the single-cell level using flow cytometry. Apparent in-cellula affinities of protein-protein interactions can be extracted from the acquired data through a titration-like analysis. The applicability of this approach is demonstrated on a diverse set of interactions with proteins from different families and organisms and with in-vitro dissociation constants ranging from picomolar to micromolar.

Engineered variants of the Ras effector protein RASSF5 (NORE1A) promote anticancer activities in lung adenocarcinoma

Within the superfamily of small GTPases, Ras appears to be the master regulator of such processes as cell cycle progression, cell division, and apoptosis. Several oncogenic Ras mutations at amino acid positions 12, 13, and 61 have been identified that lose their ability to hydrolyze GTP, giving rise to constitutive signaling and eventually development of cancer. While disruption of the Ras/effector interface is an attractive strategy for drug design to prevent this constitutive activity, inhibition of this interaction using small molecules is impractical due to the absence of a cavity to which such molecules could bind. However, proteins and especially natural Ras effectors that bind to the Ras/effector interface with high affinity could disrupt Ras/effector interactions and abolish procancer pathways initiated by Ras oncogene. Using a combination of computational design and in vitro evolution, we engineered high-affinity Ras-binding proteins starting from a natural Ras effector, RASSF5 (NORE1A), which is encoded by a tumor suppressor gene. Unlike previously reported Ras oncogene inhibitors, the proteins we designed not only inhibit Ras-regulated procancer pathways, but also stimulate anticancer pathways initiated by RASSF5. We show that upon introduction into A549 lung carcinoma cells, the engineered RASSF5 mutants decreased cell viability and mobility to a significantly greater extent than WT RASSF5. In addition, these mutant proteins induce cellular senescence by increasing acetylation and decreasing phosphorylation of p53. In conclusion, engineered RASSF5 variants provide an attractive therapeutic strategy able to oppose cancer development by means of inhibiting of procancer pathways and stimulating anticancer processes.

Divergent Mechanisms Activating RAS and Small GTPases Through Post-translational Modification

RAS is a founding member of the RAS superfamily of GTPases. These small 21 kDa proteins function as molecular switches to initialize signaling cascades involved in various cellular processes, including gene expression, cell growth, and differentiation. RAS is activated by GTP loading and deactivated upon GTP hydrolysis to GDP. Guanine nucleotide exchange factors (GEFs) and GTPase-activating proteins (GAPs) accelerate GTP loading and hydrolysis, respectively. These accessory proteins play a fundamental role in regulating activities of RAS superfamily small GTPase via a conserved guanine binding (G)-domain, which consists of five G motifs. The Switch regions lie within or proximal to the G2 and G3 motifs, and undergo dynamic conformational changes between the GDP-bound "OFF" state and GTP-bound "ON" state. They play an important role in the recognition of regulatory factors (GEFs and GAPs) and effectors. The G4 and G5 motifs are the focus of the present work and lie outside Switch regions. These motifs are responsible for the recognition of the guanine moiety in GTP and GDP, and contain residues that undergo post-translational modifications that underlie new mechanisms of RAS regulation. Post-translational modification within the G4 and G5 motifs activates RAS by populating the GTP-bound "ON" state, either through enhancement of intrinsic guanine nucleotide exchange or impairing GAP-mediated down-regulation. Here, we provide a comprehensive review of post-translational modifications in the RAS G4 and G5 motifs, and describe the role of these modifications in RAS activation as well as potential applications for cancer therapy.

A comprehensive analysis of RAS-effector interactions reveals interaction hotspots and new binding partners

RAS effectors specifically interact with GTP-bound RAS proteins to link extracellular signals to downstream signaling pathways. These interactions rely on two types of domains, called RAS-binding (RB) and RAS association (RA) domains, which share common structural characteristics. Although the molecular nature of RAS-effector interactions is well-studied for some proteins, most of the RA/RB-domain-containing proteins remain largely uncharacterized. Here, we searched through human proteome databases, extracting 41 RA domains in 39 proteins and 16 RB domains in 14 proteins, each of which can specifically select at least one of the 25 members in the RAS family. We next comprehensively investigated the sequence–structure–function relationship between different representatives of the RAS family, including HRAS, RRAS, RALA, RAP1B, RAP2A, RHEB1, and RIT1, with all members of RA domain family proteins (RASSFs) and the RB-domain-containing CRAF. The binding affinity for RAS-effector interactions, determined using fluorescence polarization, broadly ranged between high (0.3 μM) and very low (500 μM) affinities, raising interesting questions about the consequence of these variable binding affinities in the regulation of signaling events. Sequence and structural alignments pointed to two interaction hotspots in the RA/RB domains, consisting of an average of 19 RAS-binding residues. Moreover, we found novel interactions between RRAS1, RIT1, and RALA and RASSF7, RASSF9, and RASSF1, respectively, which were systematically explored in sequence–structure–property relationship analysis, and validated by mutational analysis. These data provide a set of distinct functional properties and putative biological roles that should now be investigated in the cellular context.

The Ins and Outs of RAS Effector Complexes

RAS oncogenes are among the most commonly mutated proteins in human cancers. They regulate a wide range of effector pathways that control cell proliferation, survival, differentiation, migration and metabolic status. Including aberrations in these pathways, RAS-dependent signaling is altered in more than half of human cancers. Targeting mutant RAS proteins and their downstream oncogenic signaling pathways has been elusive. However, recent results comprising detailed molecular studies, large scale omics studies and computational modeling have painted a new and more comprehensive portrait of RAS signaling that helps us to understand the intricacies of RAS, how its physiological and pathophysiological functions are regulated, and how we can target them. Here, we review these efforts particularly trying to relate the detailed mechanistic studies with global functional studies. We highlight the importance of computational modeling and data integration to derive an actionable understanding of RAS signaling that will allow us to design new mechanism-based therapies for RAS mutated cancers.

Novel, provable algorithms for efficient ensemble-based computational protein design and their application to the redesign of the c-Raf-RBD:KRas protein-protein interface

The K* algorithm provably approximates partition functions for a set of states (e.g., protein, ligand, and protein-ligand complex) to a user-specified accuracy ε. Often, reaching an ε-approximation for a particular set of partition functions takes a prohibitive amount of time and space. To alleviate some of this cost, we introduce two new algorithms into the osprey suite for protein design: fries, a Fast Removal of Inadequately Energied Sequences, and EWAK*, an Energy Window Approximation to K*. fries pre-processes the sequence space to limit a design to only the most stable, energetically favorable sequence possibilities. EWAK* then takes this pruned sequence space as input and, using a user-specified energy window, calculates K* scores using the lowest energy conformations. We expect fries/EWAK* to be most useful in cases where there are many unstable sequences in the design sequence space and when users are satisfied with enumerating the low-energy ensemble of conformations. In combination, these algorithms provably retain calculational accuracy while limiting the input sequence space and the conformations included in each partition function calculation to only the most energetically favorable, effectively reducing runtime while still enriching for desirable sequences. This combined approach led to significant speed-ups compared to the previous state-of-the-art multi-sequence algorithm, BBK*, while maintaining its efficiency and accuracy, which we show across 40 different protein systems and a total of 2,826 protein design problems. Additionally, as a proof of concept, we used these new algorithms to redesign the protein-protein interface (PPI) of the c-Raf-RBD:KRas complex. The Ras-binding domain of the protein kinase c-Raf (c-Raf-RBD) is the tightest known binder of KRas, a protein implicated in difficult-to-treat cancers. fries/EWAK* accurately retrospectively predicted the effect of 41 different sets of mutations in the PPI of the c-Raf-RBD:KRas complex. Notably, these mutations include mutations whose effect had previously been incorrectly predicted using other computational methods. Next, we used fries/EWAK* for prospective design and discovered a novel point mutation that improves binding of c-Raf-RBD to KRas in its active, GTP-bound state (KRas^GTP). We combined this new mutation with two previously reported mutations (which were highly-ranked by osprey) to create a new variant of c-Raf-RBD, c-Raf-RBD(RKY). fries/EWAK* in osprey computationally predicted that this new variant binds even more tightly than the previous best-binding variant, c-Raf-RBD(RK). We measured the binding affinity of c-Raf-RBD(RKY) using a bio-layer interferometry (BLI) assay, and found that this new variant exhibits single-digit nanomolar affinity for KRas^GTP, confirming the computational predictions made with fries/EWAK*. This new variant binds roughly five times more tightly than the previous best known binder and roughly 36 times more tightly than the design starting point (wild-type c-Raf-RBD). This study steps through the advancement and development of computational protein design by presenting theory, new algorithms, accurate retrospective designs, new prospective designs, and biochemical validation.

Computational structure-based protein design is an innovative tool for redesigning proteins to introduce a particular or novel function. One such function is improving the binding of one protein to another, which can increase our understanding of important protein systems. Herein we introduce two novel, provable algorithms, fries and EWAK*, for more efficient computational structure-based protein design as well as their application to the redesign of the c-Raf-RBD:KRas protein-protein interface. These new algorithms speed-up computational structure-based protein design while maintaining accurate calculations, allowing for larger, previously infeasible protein designs. Additionally, using fries and EWAK* within the osprey suite, we designed the tightest known binder of KRas, a heavily studied cancer target that interacts with a number of different proteins. This previously undiscovered variant of a KRas-binding domain, c-Raf-RBD, has potential to serve as a tool to further probe the protein-protein interface of KRas with its effectors and its discovery alone emphasizes the potential for more successful applications of computational structure-based protein design.

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.

A Quantitative Tri-fluorescent Yeast Two-hybrid System: From Flow Cytometry to In cellula Affinities

We present a technological advancement for the estimation of the affinities of Protein-Protein Interactions (PPIs) in living cells. A novel set of vectors is introduced that enables a quantitative yeast two-hybrid system based on fluorescent fusion proteins. The vectors allow simultaneous quantification of the reaction partners (Bait and Prey) and the reporter at the single-cell level by flow cytometry. We validate the applicability of this system on a small but diverse set of PPIs (eleven protein families from six organisms) with different affinities; the dissociation constants range from 117 pm to 17 μm After only two hours of reaction, expression of the reporter can be detected even for the weakest PPI. Through a simple gating analysis, it is possible to select only cells with identical expression levels of the reaction partners. As a result of this standardization of expression levels, the mean reporter levels directly reflect the affinities of the studied PPIs. With a set of PPIs with known affinities, it is straightforward to construct an affinity ladder that permits rapid classification of PPIs with thus far unknown affinities. Conventional software can be used for this analysis. To permit automated analysis, we provide a graphical user interface for the Python-based FlowCytometryTools package.

Conformation-specific inhibitors of activated Ras GTPases reveal limited Ras dependency of patient-derived cancer organoids

The small GTPases H, K, and NRAS are molecular switches indispensable for proper regulation of cellular proliferation and growth. Several mutations in the genes encoding members of this protein family are associated with cancer and result in aberrant activation of signaling processes caused by a deregulated recruitment of downstream effector proteins. In this study, we engineered variants of the Ras-binding domain (RBD) of the C-Raf proto-oncogene, Ser/Thr kinase (CRAF). These variants bound with high affinity with the effector-binding site of Ras in an active conformation. Structural characterization disclosed how the newly identified RBD mutations cooperate and thereby enhance affinity with the effector-binding site in Ras compared with WT RBD. The engineered RBD variants closely mimicked the interaction mode of naturally occurring Ras effectors and acted as dominant-negative affinity reagents that block Ras signal transduction. Experiments with cancer cells showed that expression of these RBD variants inhibits Ras signaling, reducing cell growth and inducing apoptosis. Using these optimized RBD variants, we stratified patient-derived colorectal cancer organoids with known Ras mutational status according to their response to Ras inhibition. These results revealed that the presence of Ras mutations was insufficient to predict sensitivity to Ras inhibition, suggesting that not all of these tumors required Ras signaling for proliferation. In summary, by engineering the Ras/Raf interface of the CRAF-RBD, we identified potent and selective inhibitors of Ras in its active conformation that outcompete binding of Ras-signaling effectors.

Structural snapshots of RAF kinase interactions

RAF (rapidly accelerated fibrosarcoma) Ser/Thr kinases (ARAF, BRAF, and CRAF) link the RAS (rat sarcoma) protein family with the MAPK (mitogen-activated protein kinase) pathway and control cell growth, differentiation, development, aging, and tumorigenesis. Their activity is specifically modulated by protein-protein interactions, post-translational modifications, and conformational changes in specific spatiotemporal patterns via various upstream regulators, including the kinases, phosphatase, GTPases, and scaffold and modulator proteins. Dephosphorylation of Ser-259 (CRAF numbering) and dissociation of 14-3-3 release the RAF regulatory domains RAS-binding domain and cysteine-rich domain for interaction with RAS-GTP and membrane lipids. This, in turn, results in RAF phosphorylation at Ser-621 and 14-3-3 reassociation, followed by its dimerization and ultimately substrate binding and phosphorylation. This review focuses on structural understanding of how distinct binding partners trigger a cascade of molecular events that induces RAF kinase activation.

RAS ubiquitylation modulates effector interactions

RAS proteins function as molecular switches that regulate cellular growth by cycling between active GTP- and inactive GDP bound states. While RAS activity is modulated by factors (guanine nucleotide exchange and GTPase activating proteins) that control levels of active Ras-GTP, RAS proteins also undergo a number of post-translational modifications that regulate their function. One such modification is ubiquitylation. Monoubiquitylation of KRAS at lysine 147 (mUbRAS) enhances Ras activation and promotes signaling through the RAF and Phosphoinositide 3-Kinase (PI3K) signaling pathways. We have previously shown that mUbRAS leads to activation of RAS through a defect in GTPase activating protein (GAP) mediated downregulation, similar to the action of most oncogenic mutations. Consistent with these findings, we now show that mUbRASimpairsRAS binding to the p120 GAP catalytic domain. Mutations in activated G12V RAS that prevent ubiquitylaton at 147 show a decrease in tumorigenesis, suggesting that in addition to activating KRAS, monoubiquitylation at this site may promote downstream signaling and transformation. To investigate whether mUbRAS alters RAS effector interactions, we chemically ubiquitylated KRAS at residue 147 and characterized binding of mUbRAS to RAS binding domains (RBDs) from three distinct downstream effectors that play key roles in RAS-mediated transformation. Results from these studies show a decrease in binding of mUbRAS (7-10-fold) relative to the CRAF RAS Binding Domain (RBD), the catalytic subunit of Phosphoinositide 3-Kinase catalytic gamma (PI3Kcγ) and RALGDS RBD. Intriguingly, we find that mUbRAS shows greatly enhanced (> 40-fold) binding to the CRAF RBD when bound to GDP. These findings, taken together, suggest that mUbRASmay promoteactivation of RAS through a GAP defect, and facilitate RAF association and MAPK signaling in a nucleotide independent manner.

Antibody targeting intracellular oncogenic Ras mutants exerts anti-tumour effects after systemic administration

Oncogenic Ras mutants, frequently detected in human cancers, are high-priority anticancer drug targets. However, direct inhibition of oncogenic Ras mutants with small molecules has been extremely challenging. Here we report the development of a human IgG1 format antibody, RT11, which internalizes into the cytosol of living cells and selectively binds to the activated GTP-bound form of various oncogenic Ras mutants to block the interactions with effector proteins, thereby suppressing downstream signalling and exerting anti-proliferative effects in a variety of tumour cells harbouring oncogenic Ras mutants. When systemically administered, an RT11 variant with an additional tumour-associated integrin binding moiety for tumour tissue targeting significantly inhibits the in vivo growth of oncogenic Ras-mutated tumour xenografts in mice, but not wild-type Ras-harbouring tumours. Our results demonstrate the feasibility of developing therapeutic antibodies for direct targeting of cytosolic proteins that are inaccessible using current antibody technology.

Deconstruction of the Ras switching cycle through saturation mutagenesis

Ras proteins are highly conserved signaling molecules that exhibit regulated, nucleotide-dependent switching between active and inactive states. The high conservation of Ras requires mechanistic explanation, especially given the general mutational tolerance of proteins. Here, we use deep mutational scanning, biochemical analysis and molecular simulations to understand constraints on Ras sequence. Ras exhibits global sensitivity to mutation when regulated by a GTPase activating protein and a nucleotide exchange factor. Removing the regulators shifts the distribution of mutational effects to be largely neutral, and reveals hotspots of activating mutations in residues that restrain Ras dynamics and promote the inactive state. Evolutionary analysis, combined with structural and mutational data, argue that Ras has co-evolved with its regulators in the vertebrate lineage. Overall, our results show that sequence conservation in Ras depends strongly on the biochemical network in which it operates, providing a framework for understanding the origin of global selection pressures on proteins.

Secreted Glioblastoma Nanovesicles Contain Intracellular Signaling Proteins and Active Ras Incorporated in a Farnesylation-dependent Manner

Glioblastomas (GBMs) are malignant brain tumors with a median survival of less than 18 months. Redundancy of signaling pathways represented within GBMs contributes to their therapeutic resistance. Exosomes are extracellular nanovesicles released from cells and present in human biofluids that represent a possible biomarker of tumor signaling state that could aid in personalized treatment. Herein, we demonstrate that mouse GBM cell-derived extracellular nanovesicles resembling exosomes from an H-RasV12 myr-Akt mouse model for GBM are enriched for intracellular signaling cascade proteins (GO: 0007242) and Ras protein signal transduction (GO: 0007265), and contain active Ras. Active Ras isolated from human and mouse GBM extracellular nanovesicles lysates using the Ras-binding domain of Raf also coprecipitates with ESCRT (endosomal sorting complex required for transport)-associated exosome proteins Vps4a and Alix. Although we initially hypothesized a role for active Ras protein signaling in exosome biogenesis, we found that GTP binding of K-Ras was dispensable for its packaging within extracellular nanovesicles and for the release of Alix. By contrast, farnesylation of K-Ras was required for its packaging within extracellular nanovesicles, yet expressing a K-Ras farnesylation mutant did not decrease the number of nanovesicles or the amount of Alix protein released per cell. Overall, these results emphasize the primary importance of membrane association in packaging of extracellular nanovesicle factors and indicate that screening nanovesicles within human fluids could provide insight into tissue origin and the wiring of signaling proteins at membranes to predict onset and behavior of cancer and other diseases linked to deregulated membrane signaling states.

Observing Extremely Weak Protein-Protein Interactions with Conventional Single-Molecule Fluorescence Microscopy

Extremely weak protein-protein interactions (PPIs), signified by micromolar or even millimolar dissociation constants, are one of the keys to understanding the rapid responses of cellular systems. Although single-molecule methods are particularly useful in determining kinetics of biological processes, their application is largely limited to rather strong interactions because of the diffraction-limited observation volume. In this study, we report a single-molecule method that allows the characterization of PPIs using a prey concentration 4 orders of magnitude lower than the dissociation constant. Instead of increasing the concentration of diffusing molecules, which is inevitably limited by the optical diffraction limit, we employed an increased density of surface bait protein. The low occupancy of the surface baits permitted determination of the kinetics with single-molecule resolution. We used this approach to study a PPI network consisting of Ras and its downstream proteins including full-length Rafs and catalytic subunits of phosphoinositide 3-kinase.

Inferring the microscopic surface energy of protein-protein interfaces from mutation data

Mutations at protein-protein recognition sites alter binding strength by altering the chemical nature of the interacting surfaces. We present a simple surface energy model, parameterized with empirical ΔΔG values, yielding mean energies of -48 cal mol(-1) Å(-2) for interactions between hydrophobic surfaces, -51 to -80 cal mol(-1) Å(-2) for surfaces of complementary charge, and 66-83 cal mol(-1) Å(-2) for electrostatically repelling surfaces, relative to the aqueous phase. This places the mean energy of hydrophobic surface burial at -24 cal mol(-1) Å(-2) . Despite neglecting configurational entropy and intramolecular changes, the model correlates with empirical binding free energies of a functionally diverse set of rigid-body interactions (r = 0.66). When used to rerank docking poses, it can place near-native solutions in the top 10 for 37% of the complexes evaluated, and 82% in the top 100. The method shows that hydrophobic burial is the driving force for protein association, accounting for 50-95% of the cohesive energy. The model is available open-source from http://life.bsc.es/pid/web/surface_energy/ and via the CCharpPPI web server http://life.bsc.es/pid/ccharppi/.

Multiplex single-molecule interaction profiling of DNA-barcoded proteins

In contrast with advances in massively parallel DNA sequencing, high-throughput protein analyses are often limited by ensemble measurements, individual analyte purification and hence compromised quality and cost-effectiveness. Single-molecule protein detection using optical methods is limited by the number of spectrally non-overlapping chromophores. Here we introduce a single-molecular-interaction sequencing (SMI-seq) technology for parallel protein interaction profiling leveraging single-molecule advantages. DNA barcodes are attached to proteins collectively via ribosome display or individually via enzymatic conjugation. Barcoded proteins are assayed en masse in aqueous solution and subsequently immobilized in a polyacrylamide thin film to construct a random single-molecule array, where barcoding DNAs are amplified into in situ polymerase colonies (polonies) and analysed by DNA sequencing. This method allows precise quantification of various proteins with a theoretical maximum array density of over one million polonies per square millimetre. Furthermore, protein interactions can be measured on the basis of the statistics of colocalized polonies arising from barcoding DNAs of interacting proteins. Two demanding applications, G-protein coupled receptor and antibody-binding profiling, are demonstrated. SMI-seq enables 'library versus library' screening in a one-pot assay, simultaneously interrogating molecular binding affinity and specificity.

Quantitative in vivo fluorescence cross-correlation analyses highlight the importance of competitive effects in the regulation of protein-protein interactions

Computer-assisted simulation is a promising approach for clarifying complicated signaling networks. However, this approach is currently limited by a deficiency of kinetic parameters determined in living cells. To overcome this problem, we applied fluorescence cross-correlation spectrometry (FCCS) to measure dissociation constant (Kd) values of signaling molecule complexes in living cells (in vivo Kd). Among the pairs of fluorescent molecules tested, that of monomerized enhanced green fluorescent protein (mEGFP) and HaloTag-tetramethylrhodamine was most suitable for the measurement of in vivo Kd by FCCS. Using this pair, we determined 22 in vivo Kd values of signaling molecule complexes comprising the epidermal growth factor receptor (EGFR)-Ras-extracellular signal-regulated kinase (ERK) mitogen-activated protein (MAP) kinase pathway. With these parameters, we developed a kinetic simulation model of the EGFR-Ras-ERK MAP kinase pathway and uncovered a potential role played by stoichiometry in Shc binding to EGFR during the peak activations of Ras, MEK, and ERK. Intriguingly, most of the in vivo Kd values determined in this study were higher than the in vitro Kd values reported previously, suggesting the significance of competitive bindings inside cells. These in vivo Kd values will provide a sound basis for the quantitative understanding of signal transduction.

Flexibility and binding affinity in protein-ligand, protein-protein and multi-component protein interactions: limitations of current computational approaches

The recognition process between a protein and a partner represents a significant theoretical challenge. In silico structure-based drug design carried out with nothing more than the three-dimensional structure of the protein has led to the introduction of many compounds into clinical trials and numerous drug approvals. Central to guiding the discovery process is to recognize active among non-active compounds. While large-scale computer simulations of compounds taken from a library (virtual screening) or designed de novo are highly desirable in the post-genomic area, many technical problems remain to be adequately addressed. This article presents an overview and discusses the limits of current computational methods for predicting the correct binding pose and accurate binding affinity. It also presents the performances of the most popular algorithms for exploring binary and multi-body protein interactions.

Cu2+-cyclen as probe to identify conformational states in guanine nucleotide binding proteins

(31)P NMR spectroscopy is a suitable method for identifying conformational states in the active site of guanine nucleotide binding proteins detecting the nucleotide placed there. Because there is no labeling necessary, this method is gaining increasing interest. By (31)P NMR spectroscopy two major conformational states, namely state 1(T) and state 2(T), can be detected in active Ras protein characterized by different chemical shifts. Depending on the conformational state Ras shows clearly different physiological properties. Meanwhile analogous conformational equilibria could also be shown for other members of the Ras superfamily. It is often difficult to determine the conformational states of the proteins on the basis of chemical shift alone; therefore, direct detection would be a great advantage. With the use of Cu(2+)-cyclen which selectively interacts only with one of the major conformational states (state 1) one has a probe to distinguish between the two states, because only proteins existing in conformational state 1 interact with the Cu(2+)-cyclen at low millimolar concentrations. The suitability was proven using Ras(wt) and Ras mutants, Ras complexed with GTP, GppNHp, or GTPγS, as well as two further members of the Ras superfamily namely Arf1 and Ran.

Protein binding specificity versus promiscuity

Interactions between macromolecules in general, and between proteins in particular, are essential for any life process. Examples include transfer of information, inhibition or activation of function, molecular recognition as in the immune system, assembly of macromolecular structures and molecular machines, and more. Proteins interact with affinities ranging from millimolar to femtomolar and, because affinity determines the concentration required to obtain 50% binding, the amount of different complexes formed is very much related to local concentrations. Although the concentration of a specific binding partner is usually quite low in the cell (nanomolar to micromolar), the total concentration of other macromolecules is very high, allowing weak and non-specific interactions to play important roles. In this review we address the question of binding specificity, that is, how do some proteins maintain monogamous relations while others are clearly polygamous. We examine recent work that addresses the molecular and structural basis for specificity versus promiscuity. We show through examples how multiple solutions exist to achieve binding via similar interfaces and how protein specificity can be tuned using both positive and negative selection (specificity by demand). Binding of a protein to numerous partners can be promoted through variation in which residues are used for binding, conformational plasticity and/or post-translational modification. Natively unstructured regions represent the extreme case in which structure is obtained only upon binding. Many natively unstructured proteins serve as hubs in protein-protein interaction networks and such promiscuity can be of functional importance in biology.

Ras homolog enriched in brain (Rheb) enhances apoptotic signaling

Rheb is a homolog of Ras GTPase that regulates cell growth, proliferation, and regeneration via mammalian target of rapamycin (mTOR). Because of the well established potential of activated Ras to promote survival, we sought to investigate the ability of Rheb signaling to phenocopy Ras. We found that overexpression of lipid-anchored Rheb enhanced the apoptotic effects induced by UV light, TNFα, or tunicamycin in an mTOR complex 1 (mTORC1)-dependent manner. Knocking down endogenous Rheb or applying rapamycin led to partial protection, identifying Rheb as a mediator of cell death. Ras and c-Raf kinase opposed the apoptotic effects induced by UV light or TNFα but did not prevent Rheb-mediated apoptosis. To gain structural insight into the signaling mechanisms, we determined the structure of Rheb-GDP by NMR. The complex adopts the typical canonical fold of RasGTPases and displays the characteristic GDP-dependent picosecond to nanosecond backbone dynamics of the switch I and switch II regions. NMR revealed Ras effector-like binding of activated Rheb to the c-Raf-Ras-binding domain (RBD), but the affinity was 1000-fold lower than the Ras/RBD interaction, suggesting a lack of functional interaction. shRNA-mediated knockdown of apoptosis signal-regulating kinase 1 (ASK-1) strongly reduced UV or TNFα-induced apoptosis and suppressed enhancement by Rheb overexpression. In conclusion, Rheb-mTOR activation not only promotes normal cell growth but also enhances apoptosis in response to diverse toxic stimuli via an ASK-1-mediated mechanism. Pharmacological regulation of the Rheb/mTORC1 pathway using rapamycin should take the presence of cellular stress into consideration, as this may have clinical implications.

Zoledronic acid potentiates mTOR inhibition and abolishes the resistance of osteosarcoma cells to RAD001 (Everolimus): pivotal role of the prenylation process

Despite recent improvements in therapeutic management of osteosarcoma, ongoing challenges in improving the response to chemotherapy warrants new strategies still needed to improve overall patient survival. In this study, we investigated in vivo the effects of RAD001 (Everolimus), a new orally available mTOR inhibitor, on the growth of human and mouse osteosarcoma cells either alone or in combination with zoledronate (ZOL), an anti-osteoporotic drug used to treat bone metastases. RAD001 inhibited osteosarcoma cell proliferation in a dose- and time-dependent manner with no modification of cell-cycle distribution. Combination with ZOL augmented this inhibition of cell proliferation, decreasing PI3K/mTOR signaling compared with single treatments. Notably, in contrast to RAD001, ZOL downregulated isoprenylated membrane-bound Ras concomitantly with an increase of nonisoprenylated cytosolic Ras in sensitive and resistant osteosarcoma cell lines to both drugs. Moreover, ZOL and RAD001 synergized to decrease Ras isoprenylation and GTP-bound Ras levels. Further, the drug combination reduced tumor development in two murine models of osteoblastic or osteolytic osteosarcoma. We found that ZOL could reverse RAD001 resistance in osteosarcoma, limiting osteosarcoma cell growth in combination with RAD001. Our findings rationalize further study of the applications of mTOR and mevalonate pathway inhibitors that can limit protein prenylation pathways.

What makes Ras an efficient molecular switch: a computational, biophysical, and structural study of Ras-GDP interactions with mutants of Raf

Ras is a small GTP-binding protein that is an essential molecular switch for a wide variety of signaling pathways including the control of cell proliferation, cell cycle progression and apoptosis. In the GTP-bound state, Ras can interact with its effectors, triggering various signaling cascades in the cell. In the GDP-bound state, Ras looses its ability to bind to known effectors. The interaction of the GTP-bound Ras (Ras(GTP)) with its effectors has been studied intensively. However, very little is known about the much weaker interaction between the GDP-bound Ras (Ras(GDP)) and Ras effectors. We investigated the factors underlying the nucleotide-dependent differences in Ras interactions with one of its effectors, Raf kinase. Using computational protein design, we generated mutants of the Ras-binding domain of Raf kinase (Raf) that stabilize the complex with Ras(GDP). Most of our designed mutations narrow the gap between the affinity of Raf for Ras(GTP) and Ras(GDP), producing the desired shift in binding specificity towards Ras(GDP). A combination of our best designed mutation, N71R, with another mutation, A85K, yielded a Raf mutant with a 100-fold improvement in affinity towards Ras(GDP). The Raf A85K and Raf N71R/A85K mutants were used to obtain the first high-resolution structures of Ras(GDP) bound to its effector. Surprisingly, these structures reveal that the loop on Ras previously termed the switch I region in the Ras(GDP).Raf mutant complex is found in a conformation similar to that of Ras(GTP) and not Ras(GDP). Moreover, the structures indicate an increased mobility of the switch I region. This greater flexibility compared to the same loop in Ras(GTP) is likely to explain the natural low affinity of Raf and other Ras effectors to Ras(GDP). Our findings demonstrate that an accurate balance between a rigid, high-affinity conformation and conformational flexibility is required to create an efficient and stringent molecular switch.