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

GDP-bound Rab27a regulates clathrin disassembly through HSPA8 after insulin secretion

Clathrin-dependent endocytosis is a key process for secretory cells, in which molecules on the plasma membrane are both degraded and recycled in a stimulus-dependent manner. There are many reports showing that disruption of endocytosis is involved in the onset of various diseases. Recently, it has been reported that such disruption in pancreatic β-cells causes impaired insulin secretion and might be associated with the pathology of diabetes mellitus. Compared with exocytosis, there are few reports on the molecular mechanism of endocytosis in pancreatic β-cells. We previously reported that GDP-bound Rab27a regulates endocytosis through its GDP-dependent effectors after insulin secretion. In this study, we identified heat shock protein family A member 8 (HSPA8) as a novel interacting protein for GDP-bound Rab27a. HSPA8 directly bound GDP-bound Rab27a via the β2 region of its substrate binding domain (SBD). The β2 fragment was capable of inhibiting the interaction between HSPA8 and GDP-bound Rab27a, and suppressed glucose-induced clathrin-dependent endocytosis in pancreatic β-cells. The region also affected clathrin dynamics on purified clathrin-coated vesicles (CCVs). These results suggest that the interaction between GDP-bound Rab27a and HSPA8 regulates clathrin disassembly from CCVs and subsequent vesicle transport. The regulatory stages in endocytosis by HSPA8 differ from those for other GDP-bound Rab27a effectors. This study shows that GDP-bound Rab27a dominantly regulates each stage in glucose-induced endocytosis through its specific effectors in pancreatic β-cells.

Targeting KRAS in Pancreatic Ductal Adenocarcinoma: The Long Road to Cure

Pancreatic ductal adenocarcinoma (PDAC) remains an important cause of cancer-related mortality, and it is expected to play an even bigger part in cancer burden in the years to come. Despite concerted efforts from scientists and physicians, patients have experienced little improvement in survival over the past decades, possibly because of the non-specific nature of the tested treatment modalities. Recently, the discovery of potentially targetable molecular alterations has paved the way for the personalized treatment of PDAC. Indeed, the central piece in the molecular framework of PDAC is starting to be unveiled. KRAS mutations are seen in 90% of PDACs, and multiple studies have demonstrated their pivotal role in pancreatic carcinogenesis. Recent investigations have shed light on the differences in prognosis as well as therapeutic implications of the different KRAS mutations and disentangled the relationship between KRAS and effectors of downstream and parallel signaling pathways. Additionally, the recognition of other mechanisms involving KRAS-mediated pathogenesis, such as KRAS dosing and allelic imbalance, has contributed to broadening the current knowledge regarding this molecular alteration. Finally, KRAS G12C inhibitors have been recently tested in patients with pancreatic cancer with relative success, and inhibitors of KRAS harboring other mutations are under clinical development. These drugs currently represent a true hope for a meaningful leap forward in this dreadful disease.

Defining bone fide effectors of RAS GTPases

RAS GTPases play essential roles in normal development and are direct drivers of human cancers. Three decades of study have failed to wholly characterize pathways stimulated by activated RAS, driven by engagement with 'effector' proteins that have RAS binding domains (RBDs). Bone fide effectors must bind directly to RAS GTPases in a nucleotide-dependent manner, and this interaction must impart a clear change in effector activity. Despite this, for most proteins currently deemed effectors there is little mechanistic understanding of how binding to the GTPase alters protein function. There has also been limited effort to comprehensively resolve the specificity of effector binding to the full array of RAS superfamily GTPase proteins. This review will summarize what is known about RAS-driven activation for an array of potential effector proteins, focusing on structural and mechanistic effects and highlighting how little is still known regarding this key paradigm of cellular signal transduction.

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.

RAS-targeted cancer therapy: Advances in drugging specific mutations

Rat sarcoma (RAS), as a frequently mutated oncogene, has been studied as an attractive target for treating RAS-driven cancers for over four decades. However, it is until the recent success of kirsten-RAS (KRAS)^G12C inhibitor that RAS gets rid of the title "undruggable". It is worth noting that the therapeutic effect of KRAS^G12C inhibitors on different RAS allelic mutations or even different cancers with KRAS^G12C varies significantly. Thus, deep understanding of the characteristics of each allelic RAS mutation will be a prerequisite for developing new RAS inhibitors. In this review, the structural and biochemical features of different RAS mutations are summarized and compared. Besides, the pathological characteristics and treatment responses of different cancers carrying RAS mutations are listed based on clinical reports. In addition, the development of RAS inhibitors, either direct or indirect, that target the downstream components in RAS pathway is summarized as well. Hopefully, this review will broaden our knowledge on RAS-targeting strategies and trigger more intensive studies on exploiting new RAS allele-specific inhibitors.

Des3PI: a fragment-based approach to design cyclic peptides targeting protein-protein interactions

Protein-protein interactions (PPIs) play crucial roles in many cellular processes and their deregulation often leads to cellular dysfunctions. One promising way to modulate PPIs is to use peptide derivatives that bind their protein target with high affinity and high specificity. Peptide modulators are often designed using secondary structure mimics. However, fragment-based design is an alternative emergent approach in the PPI field. Most of the reported computational fragment-based libraries targeting PPIs are composed of small molecules or already approved drugs, but, according to our knowledge, no amino acid based library has been reported yet. In this context, we developed a novel fragment-based approach called Des3PI (design of peptides targeting protein-protein interactions) with a library composed of natural amino acids. All the amino acids are docked into the target surface using Autodock Vina. The resulting binding modes are geometrically clustered, and, in each cluster, the most recurrent amino acids are identified and form the hotspots that will compose the designed peptide. This approach was applied on Ras and Mcl-1 proteins, as well as on A[Formula: see text] protofibril. For each target, at least five peptides generated by Des3PI were tested in silico: the peptides were first blindly docked on their target, and then, the stability of the successfully docked complexes was verified using 200 ns MD simulations. Des3PI shows very encouraging results by yielding at least 3 peptides for each protein target that succeeded in passing the two-step assessment.

Conformational Dynamics Allows Sampling of an "Active-like" State by Oncogenic K-Ras-GDP

Mutations in K-Ras GTPase replacing Gly12 with either Asp or Val are common in cancer. These mutations decelerate intrinsic and catalyzed GTP hydrolysis, leading to accumulation of K-Ras-GTP in cells. Signaling cascades initiated by K-Ras-GTP promote cell proliferation, survival, and invasion. Despite functional differences between the most frequent G12D mutation and the most aggressive and chemotherapy resistant G12V mutation, their long-suspected distinct structural features remain elusive. Using NMR, X-ray structures, and computational methods, we found that oncogenic mutants of K-Ras4B, the predominant splice variant of K-Ras, exhibit distinct conformational dynamics when GDP-bound, visiting the "active-like" conformational state similar to the one observed in GTP-bound K-Ras. This behavior distinguishes G12V from wild type and G12D K-Ras4B-GDP. The likely reason is hydrophobic interactions between the aliphatic sidechain of V12 and the Switch II region of K-Ras4B^G12V-GDP, which are distinct in K-Ras4B^G12D-GDP. In the X-ray structures, crystal contacts reduce the dynamics of the sidechain at position 12 by stabilizing the Switch I region of the protein. This explains why structural differences between G12V and G12D K-Ras have yet not been reported. Together, our results suggest a previously unknown mechanism of K-Ras activation. This mechanism relies on conformational dynamics caused by specific oncogenic mutations in the GDP-bound state. Our findings also imply that the therapeutic strategies decreasing the level of K-Ras-GTP by interfering with nucleotide exchange or by expediting GTP hydrolysis may work differently in different oncogenic mutants.

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.

Unveiling the Dynamics of KRAS4b on Lipid Model Membranes

Small GTPase proteins are ubiquitous and responsible for regulating several processes related to cell growth and differentiation. Mutations that stabilize their active state can lead to uncontrolled cell proliferation and cancer. Although these proteins are well characterized at the cellular scale, the molecular mechanisms governing their functions are still poorly understood. In addition, there is limited information about the regulatory function of the cell membrane which supports their activity. Thus, we have studied the dynamics and conformations of the farnesylated KRAS4b in various membrane model systems, ranging from binary fluid mixtures to heterogeneous raft mimics. Our approach combines long time-scale coarse-grained (CG) simulations and Markov state models to dissect the membrane-supported dynamics of KRAS4b. Our simulations reveal that protein dynamics is mainly modulated by the presence of anionic lipids and to some extent by the nucleotide state (activation) of the protein. In addition, our results suggest that both the farnesyl and the polybasic hypervariable region (HVR) are responsible for its preferential partitioning within the liquid-disordered (Ld) domains in membranes, potentially enhancing the formation of membrane-driven signaling platforms.

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.

Engineering of Biological Pathways: Complex Formation and Signal Transduction

The rational in silico design of interface mutations within protein complexes is a synthetic biology tool that enables-when introduced into biological systems-the artificial rewiring of biological pathways. Here we describe the three-dimensional structure-based design of "rewiring" mutations using the FoldX force field. Specifically, we provide the protocol for the design and selection of interface mutations in three Ras-effector complex structures (PDB entries 3KUD, 4K81, and 6AMB). Ras mutations that impair binding to some but not all interacting partners are selected.

The Ras switch in structural and historical perspective

Since its discovery as an oncogene more than 40 years ago, Ras has been and still is in the focus of many academic and pharmaceutical labs around the world. A huge amount of work has accumulated on its biology. However, many questions about the role of the different Ras isoforms in health and disease still exist and a full understanding will require more intensive work in the future. Here we try to survey some of the structural findings in a historical perspective and how it has influenced our understanding of structure-function and mechanistic relationships of Ras and its interactions. The structures show that Ras is a stable molecular machine that uses the dynamics of its switch regions for the interaction with all regulators and effectors. This conformational flexibility has been used to create small molecule drug candidates against this important oncoprotein.

Structural impact of GTP binding on downstream KRAS signaling

Oncogenic RAS proteins, involved in ∼30% of human tumors, are molecular switches of various signal transduction pathways. Here we apply a new protocol for the NMR study of KRAS in its (inactive) GDP- and (activated) GTP-bound form, allowing a comprehensive analysis of the backbone dynamics of its WT-, G12C- and G12D variants. We found that Tyr32 shows opposite mobility with respect to the backbone of its surroundings: it is more flexible in the GDP-bound form while more rigid in GTP-complexes (especially in WT- and G12D-GTP). Using the G12C/Y32F double mutant, we showed that the presence of the hydroxyl group of Tyr32 has a marked effect on the G12C-KRAS-GTP system as well. Molecular dynamics simulations indicate that Tyr32 is linked to the γ-phosphate of GTP in the activated states – an arrangement shown, using QM/MM calculations, to support catalysis. Anchoring Tyr32 to the γ-phosphate contributes to the capture of the catalytic waters participating in the intrinsic hydrolysis of GTP and supports a simultaneous triple proton transfer step (catalytic water → assisting water → Tyr32 → O1G of the γ-phosphate) leading to straightforward product formation. The coupled flip of negatively charged residues of switch I toward the inside of the effector binding pocket potentiates ligand recognition, while positioning of Thr35 to enter the coordination sphere of the Mg²⁺ widens the pocket. Position 12 mutations do not disturb the capture of Tyr32 by the γ-phosphate, but (partially) displace Gln61, which opens up the catalytic pocket and destabilizes catalytic water molecules thus impairing intrinsic hydrolysis.

Nucleotide exchange to the physiological, activated, GTP-bound form of KRAS results in the anchoring of Tyr32 within the active site.

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.

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.

Cysteine-based regulation of redox-sensitive Ras small GTPases

Reactive oxygen and nitrogen species (ROS and RNS, respectively) activate the redox-sensitive Ras small GTPases. The three canonical genes (HRAS, NRAS, and KRAS) are archetypes of the superfamily of small GTPases and are the most common oncogenes in human cancer. Oncogenic Ras is intimately linked to redox biology, mainly in the context of tumorigenesis. The Ras protein structure is highly conserved, especially in effector-binding regions. Ras small GTPases are redox-sensitive proteins thanks to the presence of the NKCD motif (Asn116-Lys 117-Cys118-Asp119). Notably, the ROS- and RNS-based oxidation of Cys118 affects protein stability, activity, and localization, and protein-protein interactions. Cys residues at positions 80, 181, 184, and 186 may also help modulate these actions. Moreover, oncogenic mutations of Gly12Cys and Gly13Cys may introduce additional oxidative centres and represent actionable drug targets. Here, the pathophysiological involvement of Cys-redox regulation of Ras proteins is reviewed in the context of cancer and heart and brain diseases.

In-cell single-molecule FRET measurements reveal three conformational state changes in RAF protein

BACKGROUND

The structures of proteins are intimately related to their functions. Significant efforts have been dedicated to the structural investigation of proteins, mainly those of purified proteins in in vitro environments. Proteins function in living cells and thus protein structures must be regulated by interactions with various molecules, some of which participate in reaction networks, depending on the states, conditions, or actions of the cell. Therefore, it is very important to understand the structural behavior of proteins in living cells.

METHODS

Single-molecule Förster resonance energy transfer (smFRET) measurements were conducted using the alternative laser excitation (ALEX) technique. smFRET distributions of cytosolic Rapidly Accelerated Fibrosarcoma (RAF) proteins in living HeLa cells were obtained with exclusion of the negative effects of photobleached fluorophores and incompletely labeled proteins on smFRET.

RESULTS

smFRET histograms of wildtype (wt) RAF in live cells exhibited two major peaks, whereas that of the S621A mutant, which has been thought to have an expanded structure, was almost single-peaked. A population shift involving the peaks for wt RAF was detected upon epidermal growth factor stimulation. Spontaneous transitions between the conformational states corresponding to the two peaks were also detected using the FRET-two-channel kernel-based density distribution estimator method in comparison to static double-stranded DNA samples.

CONCLUSIONS

Cytosolic CRAF has at least three conformational states; in addition to the closed and open forms, the fully-open form was distinctly specified. Based on the results, we propose a speculative structural model for CRAF.

GENERAL SIGNIFICANCE

Structural distribution and changes to proteins in live cells as a result of intracellular interactions were successfully identified. smFRET using ALEX is applicable to any other cytosolic proteins.

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.

PI3K: A Crucial Piece in the RAS Signaling Puzzle

RAS proteins are key signaling switches essential for control of proliferation, differentiation, and survival of eukaryotic cells. RAS proteins are mutated in 30% of human cancers. In addition, mutations in upstream or downstream signaling components also contribute to oncogenic activation of the pathway. RAS proteins exert their functions through activation of several signaling pathways and dissecting the contributions of these effectors in normal cells and in cancer is an ongoing challenge. In this review, we summarize our current knowledge about how RAS regulates type I phosphatidylinositol 3-kinase (PI3K), one of the main RAS effectors. RAS signaling through PI3K is necessary for normal lymphatic vasculature development and for RAS-induced transformation in vitro and in vivo, especially in lung cancer, where it is essential for tumor initiation and necessary for tumor maintenance.

Antibody Affinity Maturation by Computational Design

The immune systems protect our bodies from foreign molecules or antigens, where antibodies play important roles. Antibodies evolve over time upon antigen encounter by somatically mutating their genome sequences. The end result is a series of antibodies that display higher affinities and specificities to specific antigens. This process is called affinity maturation. Recent improvements in computer hardware and modeling algorithms now enable the rational design of protein structures and functions, and several works on computer-aided antibody design have been published. In this chapter, we briefly describe computational methods for antibody affinity maturation, focusing on methods for sampling antibody conformations and for scoring designed antibody variants. We also discuss lessons learned from the successful computer-aided design of antibodies.

Structural insight into the rearrangement of the switch I region in GTP-bound G12A K-Ras

K-Ras, a molecular switch that regulates cell growth, apoptosis and metabolism, is activated when it undergoes a conformation change upon binding GTP and is deactivated following the hydrolysis of GTP to GDP. Hydrolysis of GTP in water is accelerated by coordination to K-Ras, where GTP adopts a high-energy conformation approaching the transition state. The G12A mutation reduces intrinsic K-Ras GTP hydrolysis by an unexplained mechanism. Here, crystal structures of G12A K-Ras in complex with GDP, GTP, GTPγS and GppNHp, and of Q61A K-Ras in complex with GDP, are reported. In the G12A K-Ras-GTP complex, the switch I region undergoes a significant reorganization such that the Tyr32 side chain points towards the GTP-binding pocket and forms a hydrogen bond to the GTP γ-phosphate, effectively stabilizing GTP in its precatalytic state, increasing the activation energy required to reach the transition state and contributing to the reduced intrinsic GTPase activity of G12A K-Ras mutants.

Interface analysis of small GTP binding protein complexes suggests preferred membrane orientations

Crystal structures of small GTP binding protein complexes with their effectors and regulators reveal that one particularly flat side of the G domain that contains helix α4 and the C-terminal helix α5 is practically devoid of contacts. Although this observation seems trivial as the main binding targets are the switch I and II regions opposite of this side, the fact that all interacting proteins, even the largest ones, seem to avoid occupying this area (except for Ran, that does not localize to membranes) is very striking. An orientation with this 'flat' side parallel to the membrane was proposed before and would allow simultaneous interaction of the lipidated C-terminus and positive charges in the α4 helix with the membrane while being bound to effector or regulator molecules. Furthermore, this 'flat' side might be involved in regulatory mechanisms: a Ras dimer that is found in different crystal forms interacts exactly at this side. Additional interface analysis of GTPase complexes nicely confirms the effect of different flexibilities of the GTP and GDP forms. Besides Ran proteins, guanine nucleotide exchange factors (GEFs) bury the largest surface areas to provide the binding energy to open up the switch regions for nucleotide exchange.

RasIns: Genetically Encoded Intrabodies of Activated Ras Proteins

K- and H-Ras are the most commonly mutated genes in human tumors and are critical for conferring and maintaining the oncogenic phenotype in tumors with poor prognoses. Here, we design genetically encoded antibody-like ligands (intrabodies) that recognize active, GTP-bound K- and H-Ras. These ligands, which use the 10th domain of human fibronectin as their scaffold, are stable inside the cells and when fused with a fluorescent protein label, the constitutively active G12V mutant H-Ras. Primary selection of ligands against Ras with mRNA display resulted in an intrabody (termed RasIn1) that binds with a K_D of 2.1μM to H-Ras(G12V) (GTP), excellent state selectivity, and remarkable specificity for K- and H-Ras. RasIn1 recognizes residues in the Switch I region of Ras, similar to Raf-RBD, and competes with Raf-RBD for binding. An affinity maturation selection based on RasIn1 resulted in RasIn2, which binds with a K_D of 120nM and also retains excellent state selectivity. Both of these intrabodies colocalize with H-Ras, K-Ras, and G12V mutants inside the cells, providing new potential tools to monitor and modulate Ras-mediated signaling. Finally, RasIn1 and Rasin2 both display selectivity for the G12V mutants as compared with wild-type Ras providing a potential route for mutant selective recognition of Ras.

Protein Engineering by Combined Computational and In Vitro Evolution Approaches

Two alternative strategies are commonly used to study protein-protein interactions (PPIs) and to engineer protein-based inhibitors. In one approach, binders are selected experimentally from combinatorial libraries of protein mutants that are displayed on a cell surface. In the other approach, computational modeling is used to explore an astronomically large number of protein sequences to select a small number of sequences for experimental testing. While both approaches have some limitations, their combination produces superior results in various protein engineering applications. Such applications include the design of novel binders and inhibitors, the enhancement of affinity and specificity, and the mapping of binding epitopes. The combination of these approaches also aids in the understanding of the specificity profiles of various PPIs.

Ras Conformational Ensembles, Allostery, and Signaling

Ras proteins are classical members of small GTPases that function as molecular switches by alternating between inactive GDP-bound and active GTP-bound states. Ras activation is regulated by guanine nucleotide exchange factors that catalyze the exchange of GDP by GTP, and inactivation is terminated by GTPase-activating proteins that accelerate the intrinsic GTP hydrolysis rate by orders of magnitude. In this review, we focus on data that have accumulated over the past few years pertaining to the conformational ensembles and the allosteric regulation of Ras proteins and their interpretation from our conformational landscape standpoint. The Ras ensemble embodies all states, including the ligand-bound conformations, the activated (or inactivated) allosteric modulated states, post-translationally modified states, mutational states, transition states, and nonfunctional states serving as a reservoir for emerging functions. The ensemble is shifted by distinct mutational events, cofactors, post-translational modifications, and different membrane compositions. A better understanding of Ras biology can contribute to therapeutic strategies.

Combinatorial and Computational Approaches to Identify Interactions of Macrophage Colony-stimulating Factor (M-CSF) and Its Receptor c-FMS

The molecular interactions between macrophage colony-stimulating factor (M-CSF) and the tyrosine kinase receptor c-FMS play a key role in the immune response, bone metabolism, and the development of some cancers. Because no x-ray structure is available for the human M-CSF · c-FMS complex, the binding epitope for this complex is largely unknown. Our goal was to identify the residues that are essential for binding of the human M-CSF to c-FMS. For this purpose, we used a yeast surface display (YSD) approach. We expressed a combinatorial library of monomeric M-CSF (M-CSFM) single mutants and screened this library to isolate variants with reduced affinity for c-FMS using FACS. Sequencing yielded a number of single M-CSFM variants with mutations both in the direct binding interface and distant from the binding site. In addition, we used computational modeling to map the identified mutations onto the M-CSFM structure and to classify the mutations into three groups as follows: those that significantly decrease protein stability; those that destroy favorable intermolecular interactions; and those that decrease affinity through allosteric effects. To validate the YSD and computational data, M-CSFM and three variants were produced as soluble proteins; their affinity and structure were analyzed; and very good correlations with both YSD data and computational predictions were obtained. By identifying the M-CSFM residues critical for M-CSF · c-FMS interactions, we have laid down the basis for a deeper understanding of the M-CSF · c-FMS signaling mechanism and for the development of target-specific therapeutic agents with the ability to sterically occlude the M-CSF·c-FMS binding interface.

The RAS-Binding Domain of Human BRAF Protein Serine/Threonine Kinase Exhibits Allosteric Conformational Changes upon Binding HRAS

RAS binding is a critical step in the activation of BRAF protein serine/threonine kinase and stimulation of the mitogen-activated protein kinase signaling pathway. Mutations in both RAS and BRAF are associated with many human cancers. Here, we report the solution nuclear magnetic resonance (NMR) and X-ray crystal structures of the RAS-binding domain (RBD) from human BRAF. We further studied the complex between BRAF RBD and the GppNHp bound form of HRAS in solution. Backbone, side-chain, and (19)F NMR chemical shift perturbations reveal unexpected changes distal to the RAS-binding face that extend through the core of the RBD structure. Moreover, backbone amide hydrogen/deuterium exchange NMR data demonstrate conformational ensemble changes in the RBD core structure upon complex formation. These changes in BRAF RBD reveal a basis for allosteric regulation of BRAF structure and function, and suggest a mechanism by which RAS binding can signal the drastic domain rearrangements required for activation of BRAF kinase.

Biochemical and Structural Analysis of Common Cancer-Associated KRAS Mutations

KRAS mutations are the most common genetic abnormalities in cancer, but the distribution of specific mutations across cancers and the differential responses of patients with specific KRAS mutations in therapeutic clinical trials suggest that different KRAS mutations have unique biochemical behaviors. To further explain these high-level clinical differences and to explore potential therapeutic strategies for specific KRAS isoforms, we characterized the most common KRAS mutants biochemically for substrate binding kinetics, intrinsic and GTPase-activating protein (GAP)-stimulated GTPase activities, and interactions with the RAS effector, RAF kinase. Of note, KRAS G13D shows rapid nucleotide exchange kinetics compared with other mutants analyzed. This property can be explained by changes in the electrostatic charge distribution of the active site induced by the G13D mutation as shown by X-ray crystallography. High-resolution X-ray structures are also provided for the GDP-bound forms of KRAS G12V, G12R, and Q61L and reveal additional insight. Overall, the structural data and measurements, obtained herein, indicate that measurable biochemical properties provide clues for identifying KRAS-driven tumors that preferentially signal through RAF.

IMPLICATIONS

Biochemical profiling and subclassification of KRAS-driven cancers will enable the rational selection of therapies targeting specific KRAS isoforms or specific RAS effectors.

The Function of Embryonic Stem Cell-expressed RAS (E-RAS), a Unique RAS Family Member, Correlates with Its Additional Motifs and Its Structural Properties

E-RAS is a member of the RAS family specifically expressed in embryonic stem cells, gastric tumors, and hepatic stellate cells. Unlike classical RAS isoforms (H-, N-, and K-RAS4B), E-RAS has, in addition to striking and remarkable sequence deviations, an extended 38-amino acid-long unique N-terminal region with still unknown functions. We investigated the molecular mechanism of E-RAS regulation and function with respect to its sequence and structural features. We found that N-terminal extension of E-RAS is important for E-RAS signaling activity. E-RAS protein most remarkably revealed a different mode of effector interaction as compared with H-RAS, which correlates with deviations in the effector-binding site of E-RAS. Of all these residues, tryptophan 79 (arginine 41 in H-RAS), in the interswitch region, modulates the effector selectivity of RAS proteins from H-RAS to E-RAS features.

Knowledge-based modeling of peptides at protein interfaces: PiPreD

MOTIVATION

Protein-protein interactions (PPIs) underpin virtually all cellular processes both in health and disease. Modulating the interaction between proteins by means of small (chemical) agents is therefore a promising route for future novel therapeutic interventions. In this context, peptides are gaining momentum as emerging agents for the modulation of PPIs.

RESULTS

We reported a novel computational, structure and knowledge-based approach to model orthosteric peptides to target PPIs: PiPreD. PiPreD relies on a precompiled and bespoken library of structural motifs, iMotifs, extracted from protein complexes and a fast structural modeling algorithm driven by the location of native chemical groups on the interface of the protein target named anchor residues. PiPreD comprehensive and systematically samples the entire interface deriving peptide conformations best suited for the given region on the protein interface. PiPreD complements the existing technologies and provides new solutions for the disruption of selected interactions.

AVAILABILITY AND IMPLEMENTATION

Database and accessory scripts and programs are available upon request to the authors or at http://www.bioinsilico.org/PIPRED.

CONTACT

narcis.fernandez@gmail.com.

Src promotes GTPase activity of Ras via tyrosine 32 phosphorylation

Mutations in Ras GTPase and various other components of the Ras signaling pathways are among the most common genetic alterations in human cancers and also have been identified in several familial developmental syndromes. Over the past few decades it has become clear that the activity or the oncogenic potential of Ras is dependent on the nonreceptor tyrosine kinase Src to promote the Ras/Raf/MAPK pathway essential for proliferation, differentiation, and survival of eukaryotic cells. However, no direct relationship between Ras and Src has been established. We show here that Src binds to and phosphorylates GTP-, but not GDP-, loaded Ras on a conserved Y32 residue within the switch I region in vitro and that in vivo, Ras-Y32 phosphorylation markedly reduces the binding to effector Raf and concomitantly increases binding to GTPase-activating proteins and the rate of GTP hydrolysis. These results suggest that, in the context of predetermined crystallographic structures, Ras-Y32 serves as an Src-dependent keystone regulatory residue that modulates Ras GTPase activity and ensures unidirectionality to the Ras GTPase cycle.

DRoP: a water analysis program identifies Ras-GTP-specific pathway of communication between membrane-interacting regions and the active site

Ras GTPase mediates several cellular signal transduction pathways and is found mutated in a large number of cancers. It is active in the GTP-bound state, where it interacts with effector proteins, and at rest in the GDP-bound state. The catalytic domain is tethered to the membrane, with which it interacts in a nucleotide-dependent manner. Here we present the program Detection of Related Solvent Positions (DRoP) for crystallographic water analysis on protein surfaces and use it to study Ras. DRoP reads and superimposes multiple Protein Data Bank coordinates, transfers symmetry-related water molecules to the position closest to the protein surface, and ranks the waters according to how well conserved and tightly clustered they are in the set of structures. Coloring according to this rank allows visualization of the results. The effector-binding region of Ras is hydrated with highly conserved water molecules at the interface between the P-loop, switch I, and switch II, as well as at the Raf-RBD binding pocket. Furthermore, we discovered a new conserved water-mediated H-bonding network present in Ras-GTP, but not in Ras-GDP, that links the nucleotide sensor residues R161 and R164 on helix 5 to the active site. The double mutant RasN85A/N86A, where the final link between helix 5 and the nucleotide is not possible, is a severely impaired enzyme, while the single mutant RasN86A, with partial connection to the active site, has a wild-type hydrolysis rate. DRoP was instrumental in determining the water-mediated connectivity networks that link two lobes of the catalytic domain in Ras.

Structural basis for the interaction of the adaptor protein grb14 with activated ras

Grb14, a member of the Grb7-10-14 family of cytoplasmic adaptor proteins, is a tissue-specific negative regulator of insulin signaling. Grb7-10-14 contain several signaling modules, including a Ras-associating (RA) domain, a pleckstrin-homology (PH) domain, a family-specific BPS (between PH and SH2) region, and a C-terminal Src-homology-2 (SH2) domain. We showed previously that the RA and PH domains, along with the BPS region and SH2 domain, are necessary for downregulation of insulin signaling. Here, we report the crystal structure at 2.4-Å resolution of the Grb14 RA and PH domains in complex with GTP-loaded H-Ras (G12V). The structure reveals that the Grb14 RA and PH domains form an integrated structural unit capable of binding simultaneously to small GTPases and phosphoinositide lipids. The overall mode of binding of the Grb14 RA domain to activated H-Ras is similar to that of the RA domains of RalGDS and Raf1 but with important distinctions. The integrated RA-PH structural unit in Grb7-10-14 is also found in a second adaptor family that includes Rap1-interacting adaptor molecule (RIAM) and lamellipodin, proteins involved in actin-cytoskeleton rearrangement. The structure of Grb14 RA-PH in complex with H-Ras represents the first detailed molecular characterization of tandem RA-PH domains bound to a small GTPase and provides insights into the molecular basis for specificity.

A two-hybrid approach to identify inhibitors of the RAS-RAF interaction

MCP compounds were developed with the idea to inhibit RAS/RAF interaction. They were identified by carrying out high-throughput screens of chemical compounds for their ability to inhibit RAS/RAF interaction in the yeast two-hybrid assay. A number of compounds including MCP1, MCP53, and MCP110 were identified as active compounds. Their inhibition of the RAS signaling was demonstrated by examining RAF and MEK activities, phosphorylation of ERK as well as characterizing their effects on events downstream of RAF. Direct evidence for the inhibition of RAS/RAF interaction was obtained by carrying out co-IP experiments. MCP compounds inhibit proliferation of a wide range of human cancer cell lines. Combination studies with other drugs showed that MCP compounds synergize with MAPK pathway inhibitors as well as with microtubule-targeting chemotherapeutics. In particular, a strong synergy with paclitaxel was observed. Efficacy to inhibit tumor formation was demonstrated using mouse xenograft models. Combination of MCP110 and paclitaxel was particularly effective in inhibiting tumor growth in a mouse xenograft model of colorectal carcinoma.

NMR-based functional profiling of RASopathies and oncogenic RAS mutations

Defects in the RAS small G protein or its associated network of regulatory proteins that disrupt GTPase cycling are a major cause of cancer and developmental RASopathy disorders. Lack of robust functional assays has been a major hurdle in RAS pathway-targeted drug development. We used NMR to obtain detailed mechanistic data on RAS cycling defects conferred by oncogenic mutations, or full-length RASopathy-derived regulatory proteins. By monitoring the conformation of wild-type and oncogenic RAS in real-time, we show that opposing properties integrate with regulators to hyperactivate oncogenic RAS mutants. Q61L and G13D exhibited rapid nucleotide exchange and an unexpected susceptibility to GAP-mediated hydrolysis, in direct contrast with G12V, indicating different approaches must be taken to inhibit these oncoproteins. An NMR methodology was established to directly monitor RAS cycling by intact, multidomain proteins encoded by RASopathy genes in mammalian cell extracts. By measuring GAP activity from tumor cells, we demonstrate how loss of neurofibromatosis type 1 (NF1) increases RAS-GTP levels in NF1-derived cells. We further applied this methodology to profile Noonan Syndrome (NS)-derived SOS1 mutants. Combining NMR with cell-based assays allowed us to differentiate defects in catalysis, allosteric regulation, and membrane targeting of individual mutants, while revealing a membrane-dependent compensatory effect that attenuates dramatic increases in RAS activation shown by Y337C, L550P, and I252T. Our NMR method presents a precise and robust measure of RAS activity, providing mechanistic insights that facilitate discovery of therapeutics targeted against the RAS signaling network.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.

A comparison of successful and failed protein interface designs highlights the challenges of designing buried hydrogen bonds

The accurate design of new protein-protein interactions is a longstanding goal of computational protein design. However, most computationally designed interfaces fail to form experimentally. This investigation compares five previously described successful de novo interface designs with 158 failures. Both sets of proteins were designed with the molecular modeling program Rosetta. Designs were considered a success if a high-resolution crystal structure of the complex closely matched the design model and the equilibrium dissociation constant for binding was less than 10 μM. The successes and failures represent a wide variety of interface types and design goals including heterodimers, homodimers, peptide-protein interactions, one-sided designs (i.e., where only one of the proteins was mutated) and two-sided designs. The most striking feature of the successful designs is that they have fewer polar atoms at their interfaces than many of the failed designs. Designs that attempted to create extensive sets of interface-spanning hydrogen bonds resulted in no detectable binding. In contrast, polar atoms make up more than 40% of the interface area of many natural dimers, and native interfaces often contain extensive hydrogen bonding networks. These results suggest that Rosetta may not be accurately balancing hydrogen bonding and electrostatic energies against desolvation penalties and that design processes may not include sufficient sampling to identify side chains in preordered conformations that can fully satisfy the hydrogen bonding potential of the interface.

Specific Conformational States of Ras GTPase upon Effector Binding

To uncover the structural and dynamical determinants involved in the highly specific binding of Ras GTPase to its effectors, the conformational states of Ras in uncomplexed form and complexed to the downstream effectors Byr2, PI3Kγ, PLCε, and RalGDS were investigated using molecular dynamics and cross-comparison of the trajectories. The subtle changes in the dynamics and conformations of Ras upon effector binding require an analysis that targets local changes independent of global motions. Using a structural alphabet, a computational procedure is proposed to quantify local conformational changes. Positions detected by this approach were characterized as either specific for a particular effector, specific for an effector domain type, or as effector unspecific. A set of nine structurally connected residues (Ras residues 5–8, 32–35, 39–42, 55–59, 73–78, and 161–165), which link the effector binding site to the distant C-terminus, changed dynamics upon effector binding, indicating a potential effector-unspecific signaling route within the Ras structure. Additional conformational changes were detected along the N-terminus of the central β-sheet. Besides the Ras residues at the effector interface (e.g., D33, E37, D38, and Y40), which adopt effector-specific local conformations, the binding signal propagates from the interface to distant hot-spot residues, in particular to Y5 and D57. The results of this study reveal possible conformational mechanisms for the stabilization of the active state of Ras upon downstream effector binding and for the structural determinants responsible for effector specificity.

Structural basis for small G protein effector interaction of Ras-related protein 1 (Rap1) and adaptor protein Krev interaction trapped 1 (KRIT1)

Cerebral cavernous malformations (CCMs) affect 0.1-0.5% of the population resulting in leaky vasculature and severe neurological defects. KRIT1 (Krev interaction trapped-1) mutations associate with ∼40% of familial CCMs. KRIT1 is an effector of Ras-related protein 1 (Rap1) GTPase. Rap1 relocalizes KRIT1 from microtubules to cell membranes to impact integrin activation, potentially important for CCM pathology. We report the 1.95 Å co-crystal structure of KRIT1 FERM domain in complex with Rap1. Rap1-KRIT1 interaction encompasses an extended surface, including Rap1 Switch I and II and KRIT1 FERM F1 and F2 lobes. Rap1 binds KRIT1-F1 lobe using a GTPase-ubiquitin-like fold interaction but binds KRIT1-F2 lobe by a novel interaction. Point mutagenesis confirms the interaction. High similarity between KRIT1-F2/F3 and talin is revealed. Additionally, the mechanism for FERM domains acting as GTPase effectors is suggested. Finally, structure-based alignment of each lobe suggests classification of FERM domains as ERM-like and TMFK-like (talin-myosin-FAK-KRIT-like) and that FERM lobes resemble domain "modules."

Binding properties and stability of the Ras-association domain of Rap1-GTP interacting adapter molecule (RIAM)

The Rap1-GTP interacting adapter protein (RIAM) is an important protein in Rap1-mediated integrin activation. By binding to both Rap1 GTPase and talin, RIAM recruits talin to the cell membrane, thus facilitating talin-dependent integrin activation. In this article, we studied the role of the RIAM Ras-association (RA) and pleckstrin-homology (PH) domains in the interaction with Rap1. We found that the RA domain was sufficient for GTP-dependent interaction with Rap1B, and the addition of the PH domain did not change the binding affinity. We also detected GTP-independent interaction of Rap1B with the N-terminus of RIAM. In addition, we found that the PH domain stabilized the RA domain both in vitro and in cells.

RAS Mutations and Oncogenesis: Not all RAS Mutations are Created Equally

Mutation in RAS proteins is one of the most common genetic alterations observed in human and experimentally induced rodent cancers. In vivo, oncogenic mutations have been shown to occur at exons 12, 13, and 61, resulting in any 1 of 19 possible point mutations in a given tumor for a specific RAS isoform. While some studies have suggested a possible role of different mutant alleles in determining tumor severity and phenotype, no general consensus has emerged on the oncogenicity of different mutant alleles in tumor formation and progression. Part of this may be due to a lack of a single, signature pathway that shows significant alterations between different mutations. Rather, it is likely that subtle differences in the activation, or lack thereof, of downstream effectors by different RAS mutant alleles may determine the eventual outcome in terms of tumor phenotype. This paper reviews our current understanding of the potential role of different RAS mutations on tumorigenesis, highlights studies in model cell culture and in vivo systems, and discusses the potential of expression array and computational network modeling to dissect out differences in activated RAS genes in conferring a transforming phenotype.

Triathlon for energy functions: who is the winner for design of protein-protein interactions?

Computational prediction of stabilizing mutations into monomeric proteins has become an almost ordinary task. Yet, computational stabilization of protein–protein complexes remains a challenge. Design of protein–protein interactions (PPIs) is impeded by the absence of an energy function that could reliably reproduce all favorable interactions between the binding partners. In this work, we present three energy functions: one function that was trained on monomeric proteins, while the other two were optimized by different techniques to predict side-chain conformations in a dataset of PPIs. The performances of these energy functions are evaluated in three different tasks related to design of PPIs: predicting side-chain conformations in PPIs, recovering native binding-interface sequences, and predicting changes in free energy of binding due to mutations. Our findings show that both functions optimized on side-chain repacking in PPIs are more suitable for PPI design compared to the function trained on monomeric proteins. Yet, no function performs best at all three tasks. Comparison of the three energy functions and their performances revealed that (1) burial of polar atoms should not be penalized significantly in PPI design as in single-protein design and (2) contribution of electrostatic interactions should be increased several-fold when switching from single-protein to PPI design. In addition, the use of a softer van der Waals potential is beneficial in cases when backbone flexibility is important. All things considered, we define an energy function that captures most of the nuances of the binding energetics and hence, should be used in future for design of PPIs.

Allostery and population shift in drug discovery

Proteins can exist in a large number of conformations around their native states that can be characterized by an energy landscape. The landscape illustrates individual valleys, which are the conformational substates. From the functional standpoint, there are two key points: first, all functionally relevant substates pre-exist; and second, the landscape is dynamic and the relative populations of the substates will change following allosteric events. Allosteric events perturb the structure, and the energetic strain propagates and shifts the population. This can lead to changes in the shapes and properties of target binding sites. Here we present an overview of dynamic conformational ensembles focusing on allosteric events in signaling. We propose that combining equilibrium fluctuation concepts with genomic screens could help drug discovery.

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.