Shift in the equilibrium between on and off states of the allosteric switch in Ras-GppNHp affected by small molecules and bulk solvent composition

Allosteric site variants affect GTP hydrolysis on Ras

RAS GTPases are proto-oncoproteins that regulate cell growth, proliferation, and differentiation in response to extracellular signals. The signaling functions of RAS, and other small GTPases, are dependent on their ability to cycle between GDP-bound and GTP-bound states. Structural analyses suggest that GTP hydrolysis catalyzed by HRAS can be regulated by an allosteric site located between helices 3, 4, and loop 7. Here we explore the relationship between intrinsic GTP hydrolysis on HRAS and the position of helix 3 and loop 7 through manipulation of the allosteric site, showing that the two sites are functionally connected. We generated several hydrophobic mutations in the allosteric site of HRAS to promote shifts in helix 3 relative to helix 4. By combining crystallography and enzymology to study these mutants, we show that closure of the allosteric site correlates with increased hydrolysis of GTP on HRAS in solution. Interestingly, binding to the RAS binding domain of RAF kinase (RAF-RBD) inhibits GTP hydrolysis in the mutants. This behavior may be representative of a cluster of mutations found in human tumors, which potentially cooperate with RAF complex formation to stabilize the GTP-bound state of RAS.

Delineating The RAS Conformational Landscape

Mutations in RAS isoforms (KRAS, NRAS, and HRAS) are among the most frequent oncogenic alterations in many cancers, making these proteins high priority therapeutic targets. Effectively targeting RAS isoforms requires an exact understanding of their active, inactive, and druggable conformations. However, there is no structural catalog of RAS conformations to guide therapeutic targeting or examining the structural impact of RAS mutations. Here we present an expanded classification of RAS conformations based on analyses of the catalytic switch 1 (SW1) and switch 2 (SW2) loops. From 721 human KRAS, NRAS, and HRAS structures available in the Protein Data Bank (206 RAS-protein co-complexes, 190 inhibitor-bound, and 325 unbound, including 204 WT and 517 mutated structures), we created a broad conformational classification based on the spatial positions of Y32 in SW1 and Y71 in SW2. Clustering all well-modeled SW1 and SW2 loops using a density-based machine learning algorithm defined additional conformational subsets, some previously undescribed. Three SW1 conformations and nine SW2 conformations were identified, each associated with different nucleotide states (GTP-bound, nucleotide-free, and GDP-bound) and specific bound proteins or inhibitor sites. The GTP-bound SW1 conformation could be further subdivided based on the hydrogen bond type made between Y32 and the GTP γ-phosphate. Further analysis clarified the catalytic impact of G12D and G12V mutations and the inhibitor chemistries that bind to each druggable RAS conformation. Overall, this study has expanded our understanding of RAS structural biology, which could facilitate future RAS drug discovery.

Identification of functional substates of KRas during GTP hydrolysis with enhanced sampling simulations

As the hub of major signaling pathways, Ras proteins are implicated in 19% of tumor-caused cancers due to perturbations in their conformational and/or catalytic properties. Despite numerous studies, the functions of the conformational substates for the most important isoform, KRas, remain elusive. In this work, we perform an extensive simulation analysis on the conformational landscape of KRas in its various chemical states during the GTP hydrolysis cycle: the reactant state KRasGTP·Mg²⁺, the intermediate state KRasGDP·Pi·Mg²⁺ and the product state KRasGDP·Mg²⁺. The results from enhanced sampling simulations reveal that State 1 of KRasGTP·Mg²⁺ has multiple stable substates in solution, one of which might account for interacting with GEFs. State 2 of KRasGTP·Mg²⁺ features two substates "Tyr32_in" and "Tyr32_out", which are poised to interact with effectors and GAPs, respectively. For the intermediate state KRasGDP·Pi·Mg²⁺, Gln61 and Pi are found to assume a broad set of conformations, which might account for the weak oncogenic effect of Gln61 mutations in KRas in contrast to the situation in HRas and NRas. Finally, the product state KRasGDP·Mg²⁺ has more than two stable substates in solution, pointing to a conformation-selection mechanism for complexation with GEFs. Based on these results, some specific inhibition strategies for targeting the binding sites of the high-energy substates of KRas during GTP hydrolysis are discussed.

Equilibria between conformational states of the Ras oncogene protein revealed by high pressure crystallography

In this work, we experimentally investigate the allosteric transitions between conformational states on the Ras oncogene protein using high pressure crystallography. Ras protein is a small GTPase involved in central regulatory processes occurring in multiple conformational states. Ras acts as a molecular switch between active GTP-bound, and inactive GDP-bound states, controlling essential signal transduction pathways. An allosteric network of interactions between the effector binding regions and the membrane interacting regions is involved in Ras cycling. The conformational states which coexist simultaneously in solution possess higher Gibbs free energy than the ground state. Equilibria between these states can be shifted by applying pressure favouring conformations with lower partial molar volume, and has been previously analyzed by high-pressure NMR spectroscopy. High-pressure macromolecular crystallography (HPMX) is a powerful tool perfectly complementary to high-pressure NMR, allowing characterization at the molecular level with a high resolution the different allosteric states involved in the Ras cycling. We observe a transition above 300 MPa in the crystal leading to more stable conformers. Thus, we compare the crystallographic structures of Ras(wt)·Mg²⁺·GppNHp and Ras(D33K)·Mg²⁺·GppNHp at various high hydrostatic pressures. This gives insight into per-residue descriptions of the structural plasticity involved in allosteric equilibria between conformers. We have mapped out at atomic resolution the different segments of Ras protein which remain in the ground-state conformation or undergo structural changes, adopting excited-energy conformations corresponding to transient intermediate states. Such in crystallo phase transitions induced by pressure open the possibility to finely explore the structural determinants related to switching between Ras allosteric sub-states without any mutations nor exogenous partners.

The equilibria between structural states induced by pressure within the crystal structure of Ras are illustrated with different colors corresponding to different Ras substates.

Raf promotes dimerization of the Ras G-domain with increased allosteric connections

Ras dimerization is critical for Raf activation. Here we show that the Ras binding domain of Raf (Raf-RBD) induces robust Ras dimerization at low surface densities on supported lipid bilayers and, to a lesser extent, in solution as observed by size exclusion chromatography and confirmed by SAXS. Community network analysis based on molecular dynamics simulations shows robust allosteric connections linking the two Raf-RBD D113 residues located in the Galectin scaffold protein binding site of each Raf-RBD molecule and 85 Å apart on opposite ends of the dimer complex. Our results suggest that Raf-RBD binding and Ras dimerization are concerted events that lead to a high-affinity signaling complex at the membrane that we propose is an essential unit in the macromolecular assembly of higher order Ras/Raf/Galectin complexes important for signaling through the Ras/Raf/MEK/ERK pathway.

Isoform-Specific Destabilization of the Active Site Reveals a Molecular Mechanism of Intrinsic Activation of KRas G13D

Ras GTPases are mutated at codons 12, 13, and 61, with different frequencies in KRas, HRas, and NRas and in a cancer-specific manner. The G13D mutant appears in 25% of KRas-driven colorectal cancers, while observed only rarely in HRas or NRas. Structures of Ras G13D in the three isoforms show an open active site, with adjustments to the D13 backbone torsion angles and with disconnected switch regions. KRas G13D has unique features that destabilize the nucleotide-binding pocket. In KRas G13D bound to GDP, A59 is placed in the Mg²⁺ binding site, as in the HRas-SOS complex. Structure and biochemistry are consistent with an intermediate level of KRas G13D bound to GTP, relative to wild-type and KRas G12D, observed in genetically engineered mouse models. The results explain in part the elevated frequency of the G13D mutant in KRas over the other isoforms of Ras.

New Conformations of Acylation Adducts of Inhibitors of β-Lactamase from Mycobacterium tuberculosis

Mycobacterium tuberculosis (Mtb), the main causative agent of tuberculosis (TB), is naturally resistant to β-lactam antibiotics due to the production of the extended spectrum β-lactamase BlaC. β-Lactam/β-lactamase inhibitor combination therapies can circumvent the BlaC-mediated resistance of Mtb and are promising treatment options against TB. However, still little is known of the exact mechanism of BlaC inhibition by the β-lactamase inhibitors currently approved for clinical use, clavulanic acid, sulbactam, tazobactam, and avibactam. Here, we present the X-ray diffraction crystal structures of the acyl-enzyme adducts of wild-type BlaC with the four inhibitors. The +70 Da adduct derived from clavulanate and the trans-enamine acylation adducts of sulbactam and tazobactam are reported. BlaC in complex with avibactam revealed two inhibitor conformations. Preacylation binding could not be observed because inhibitor binding was not detected in BlaC variants carrying a substitution of the active site serine 70 to either alanine or cysteine, by crystallography, ITC or NMR. These results suggest that the catalytic serine 70 is necessary not only for enzyme acylation but also for increasing BlaC affinity for inhibitors in the preacylation state. The structure of BlaC with the serine to cysteine mutation showed a covalent linkage of the cysteine 70 Sγ atom to the nearby amino group of lysine 73. The differences of adduct conformations between BlaC and other β-lactamases are discussed.

Allostery and dynamics in small G proteins

The Ras family of small guanine nucleotide-binding proteins behave as molecular switches: they are switched off and inactive when bound to GDP but can be activated by GTP binding in response to signal transduction pathways. Early structural analysis showed that two regions of the protein, which change conformation depending on the nucleotide present, mediate this switch. A large number of X-ray, NMR and simulation studies have shown that this is an over-simplification. The switch regions themselves are highly dynamic and can exist in distinct sub-states in the GTP-bound form that have different affinities for other proteins. Furthermore, regions outside the switches have been found to be sensitive to the nucleotide state of the protein, indicating that allosteric change is more widespread than previously thought. Taken together, the accrued knowledge about small G protein structures, allostery and dynamics will be essential for the design and testing of the next generation of inhibitors, both orthosteric and allosteric, as well as for understanding their mode of action.

Specific Substates of Ras To Interact with GAPs and Effectors: Revealed by Theoretical Simulations and FTIR Experiments

The oncogenic Ras protein adopts various specific conformational states to execute its function in signal transduction. The large number of Ras structures obtained from X-ray and NMR experiments illustrates the diverse conformations that Ras adopts. It is difficult, however, to connect specific structural features with Ras functions. We report the free-energy landscape of Ras·GTP based on extensive explicit solvent simulations. The free-energy map clearly shows that the functional state 2 of Ras·GTP in fact has two distinct substates, denoted here as "Tyr32in" and "Tyr32out". Unbiased MD simulations show that the two substrates interconvert on the submicrosecond scale in solution, pointing to a novel mechanism for Ras·GTP to selectively interact with GAPs and effectors. This proposal is further supported by time-resolved FTIR experiments, which demonstrate that Tyr32 destabilizes the Ras·GAP complex and facilitates an efficient termination of Ras signaling.

The small GTPases K-Ras, N-Ras, and H-Ras have distinct biochemical properties determined by allosteric effects

H-Ras, K-Ras, and N-Ras are small GTPases that are important in the control of cell proliferation, differentiation, and survival, and their mutants occur frequently in human cancers. The G-domain, which catalyzes GTP hydrolysis and mediates downstream signaling, is 95% conserved between the Ras isoforms. Because of their very high sequence identity, biochemical studies done on H-Ras have been considered representative of all three Ras proteins. We show here that this is not a valid assumption. Using enzyme kinetic assays under identical conditions, we observed clear differences between the three isoforms in intrinsic catalysis of GTP by Ras in the absence and presence of the Ras-binding domain (RBD) of the c-Raf kinase protein (Raf-RBD). Given their identical active sites, isoform G-domain differences must be allosteric in origin, due to remote isoform-specific residues that affect conformational states. We present the crystal structure of N-Ras bound to a GTP analogue and interpret the kinetic data in terms of structural features specific for H-, K-, and N-Ras.

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.

Drugging Ras GTPase: a comprehensive mechanistic and signaling structural view

Ras proteins are small GTPases, cycling between inactive GDP-bound and active GTP-bound states. Through these switches they regulate signaling that controls cell growth and proliferation. Activating Ras mutations are associated with approximately 30% of human cancers, which are frequently resistant to standard therapies. Over the past few years, structural biology and in silico drug design, coupled with improved screening technology, led to a handful of promising inhibitors, raising the possibility of drugging Ras proteins. At the same time, the invariable emergence of drug resistance argues for the critical importance of additionally honing in on signaling pathways which are likely to be involved. Here we overview current advances in Ras structural knowledge, including the conformational dynamic of full-length Ras in solution and at the membrane, therapeutic inhibition of Ras activity by targeting its active site, allosteric sites, and Ras-effector protein-protein interfaces, Ras dimers, the K-Ras4B/calmodulin/PI3Kα trimer, and targeting Ras with siRNA. To mitigate drug resistance, we propose signaling pathways that can be co-targeted along with Ras and explain why. These include pathways leading to the expression (or activation) of YAP1 and c-Myc. We postulate that these and Ras signaling pathways, MAPK/ERK and PI3K/Akt/mTOR, act independently and in corresponding ways in cell cycle control. The structural data are instrumental in the discovery and development of Ras inhibitors for treating RAS-driven cancers. Together with the signaling blueprints through which drug resistance can evolve, this review provides a comprehensive and innovative master plan for tackling mutant Ras proteins.

Invited review: Small GTPases and their GAPs

Widespread utilization of small GTPases as major regulatory hubs in many different biological systems derives from a conserved conformational switch mechanism that facilitates cycling between GTP-bound active and GDP-bound inactive states under control of guanine nucleotide exchange factors (GEFs) and GTPase activating proteins (GAPs), which accelerate slow intrinsic rates of activation by nucleotide exchange and deactivation by GTP hydrolysis, respectively. Here we review developments leading to current understanding of intrinsic and GAP catalyzed GTP hydrolytic reactions in small GTPases from structural, molecular and chemical mechanistic perspectives. Despite the apparent simplicity of the GTPase cycle, the structural bases underlying the hallmark hydrolytic reaction and catalytic acceleration by GAPs are considerably more diverse than originally anticipated. Even the most fundamental aspects of the reaction mechanism have been challenging to decipher. Through a combination of experimental and in silico approaches, the outlines of a consensus view have begun to emerge for the best studied paradigms. Nevertheless, recent observations indicate that there is still much to be learned. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 431-448, 2016.

Molecular Mechanism for Conformational Dynamics of Ras·GTP Elucidated from In-Situ Structural Transition in Crystal

Ras•GTP adopts two interconverting conformational states, state 1 and state 2, corresponding to inactive and active forms, respectively. However, analysis of the mechanism for state transition was hampered by the lack of the structural information on wild-type Ras state 1 despite its fundamental nature conserved in the Ras superfamily. Here we solve two new crystal structures of wild-type H-Ras, corresponding to state 1 and state 2. The state 2 structure seems to represent an intermediate of state transition and, intriguingly, the state 1 crystal is successfully derived from this state 2 crystal by regulating the surrounding humidity. Structural comparison enables us to infer the molecular mechanism for state transition, during which a wide range of hydrogen-bonding networks across Switch I, Switch II and the α3-helix interdependently undergo gross rearrangements, where fluctuation of Tyr32, translocation of Gln61, loss of the functional water molecules and positional shift of GTP play major roles. The NMR-based hydrogen/deuterium exchange experiments also support this transition mechanism. Moreover, the unveiled structural features together with the results of the biochemical study provide a new insight into the physiological role of state 1 as a stable pool of Ras•GTP in the GDP/GTP cycle of Ras.

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.

Direct Attack on RAS: Intramolecular Communication and Mutation-Specific Effects

The crystal structure of RAS was first solved 25 years ago. In spite of tremendous and sustained efforts, there are still no drugs in the clinic that directly target this major driver of human cancers. Recent success in the discovery of compounds that bind RAS and inhibit signaling has fueled renewed enthusiasm, and in-depth understanding of the structure and function of RAS has opened new avenues for direct targeting. To succeed, we must focus on the molecular details of the RAS structure and understand at a high-resolution level how the oncogenic mutants impair function. Structural networks of intramolecular communication between the RAS active site and membrane-interacting regions on the G-domain are disrupted in oncogenic mutants. Although conserved across the isoforms, these networks are near hot spots of protein-ligand interactions with amino acid composition that varies among RAS proteins. These differences could have an effect on stabilization of conformational states of interest in attenuating signaling through RAS. The development of strategies to target these novel sites will add a fresh direction in the quest to conquer RAS-driven cancers. Clin Cancer Res; 21(8); 1810-8. ©2015 AACR. See all articles in this CCR Focus section, "Targeting RAS-Driven Cancers."

DT-Web: a web-based application for drug-target interaction and drug combination prediction through domain-tuned network-based inference

BACKGROUND

The identification of drug-target interactions (DTI) is a costly and time-consuming step in drug discovery and design. Computational methods capable of predicting reliable DTI play an important role in the field. Algorithms may aim to design new therapies based on a single approved drug or a combination of them. Recently, recommendation methods relying on network-based inference in connection with knowledge coming from the specific domain have been proposed.

DESCRIPTION

Here we propose a web-based interface to the DT-Hybrid algorithm, which applies a recommendation technique based on bipartite network projection implementing resources transfer within the network. This technique combined with domain-specific knowledge expressing drugs and targets similarity is used to compute recommendations for each drug. Our web interface allows the users: (i) to browse all the predictions inferred by the algorithm; (ii) to upload their custom data on which they wish to obtain a prediction through a DT-Hybrid based pipeline; (iii) to help in the early stages of drug combinations, repositioning, substitution, or resistance studies by finding drugs that can act simultaneously on multiple targets in a multi-pathway environment. Our system is periodically synchronized with DrugBank and updated accordingly. The website is free, open to all users, and available at http://alpha.dmi.unict.it/dtweb/.

CONCLUSIONS

Our web interface allows users to search and visualize information on drugs and targets eventually providing their own data to compute a list of predictions. The user can visualize information about the characteristics of each drug, a list of predicted and validated targets, associated enzymes and transporters. A table containing key information and GO classification allows the users to perform their own analysis on our data. A special interface for data submission allows the execution of a pipeline, based on DT-Hybrid, predicting new targets with the corresponding p-values expressing the reliability of each group of predictions. Finally, It is also possible to specify a list of genes tracking down all the drugs that may have an indirect influence on them based on a multi-drug, multi-target, multi-pathway analysis, which aims to discover drugs for future follow-up studies.

Enrichment of druggable conformations from apo protein structures using cosolvent-accelerated molecular dynamics

Here we describe the development of an improved workflow for utilizing experimental and simulated protein conformations in the structure-based design of inhibitors for anti-apoptotic Bcl-2 family proteins. Traditional structure-based approaches on similar targets are often constrained by the sparsity of available structures and difficulties in finding lead compounds that dock against flat, flexible protein-protein interaction surfaces. By employing computational docking of known small molecule inhibitors, we have demonstrated that structural ensembles derived from either accelerated MD (aMD) or MD in the presence of an organic cosolvent generally give better scores than those assessed from analogous conventional MD. Furthermore, conformations obtained from combined cosolvent aMD simulations started with the apo-Bcl-xL structure yielded better average and minimum docking scores for known binders than an ensemble of 72 experimental apo- and ligand-bound Bcl-xL structures. A detailed analysis of the simulated conformations indicates that the aMD effectively enhanced conformational sampling of the flexible helices flanking the main Bcl-xL binding groove, permitting the cosolvent acting as small ligands to penetrate more deeply into the binding pocket and shape ligand-bound conformations not evident in conventional simulations. We believe this approach could be useful for identifying inhibitors against other protein-protein interaction systems involving highly flexible binding sites, particularly for targets with less accumulated structural data.

Allosteric effects of the oncogenic RasQ61L mutant on Raf-RBD

The Ras/Raf/MEK/ERK signal transduction pathway is a major regulator of cell proliferation activated by Ras-guanosine triphosphate (GTP). The oncogenic mutant RasQ61L is not able to hydrolyze GTP in the presence of Raf and thus is a constitutive activator of this mitogenic pathway. The Ras/Raf interaction is essential for the activation of the Raf kinase domain through a currently unknown mechanism. We present the crystal structures of the Ras-GppNHp/Raf-RBD and RasQ61L-GppNHp/Raf-RBD complexes, which, in combination with MD simulations, reveal differences in allosteric interactions leading from the Ras/Raf interface to the Ras calcium-binding site and to the remote Raf-RBD loop L4. In the presence of Raf, the RasQ61L mutant has a rigid switch II relative to the wild-type and increased flexibility at the interface with switch I, which propagates across Raf-RBD. We show that in addition to local perturbations on Ras, RasQ61L has substantial long-range effects on the Ras allosteric lobe and on Raf-RBD.

Mechanism of the exchange reaction in HRAS from multiscale modeling

HRAS regulates cell growth promoting signaling processes by cycling between active (GTP-bound) and inactive (GDP-bound) states. Understanding the transition mechanism is central for the design of small molecules to inhibit the formation of RAS-driven tumors. Using a multiscale approach involving coarse-grained (CG) simulations, all-atom classical molecular dynamics (CMD; total of 3.02 µs), and steered molecular dynamics (SMD) in combination with Principal Component Analysis (PCA), we identified the structural features that determine the nucleotide (GDP) exchange reaction. We show that weakening the coupling between the SwitchI (residues 25–40) and SwitchII (residues 59–75) accelerates the opening of SwitchI; however, an open conformation of SwitchI is unstable in the absence of guanine nucleotide exchange factors (GEFs) and rises up towards the bound nucleotide to close the nucleotide pocket. Both I21 and Y32, play a crucial role in SwitchI transition. We show that an open SwitchI conformation is not necessary for GDP destabilization but is required for GDP/Mg escape from the HRAS. Further, we present the first simulation study showing displacement of GDP/Mg away from the nucleotide pocket. Both SwitchI and SwitchII, delays the escape of displaced GDP/Mg in the absence of GEF. Based on these results, a model for the mechanism of GEF in accelerating the exchange process is hypothesized.

Common hydrogen bond interactions in diverse phosphoryl transfer active sites

Phosphoryl transfer reactions figure prominently in energy metabolism, signaling, transport and motility. Prior detailed studies of selected systems have highlighted mechanistic features that distinguish different phosphoryl transfer enzymes. Here, a top-down approach is developed for comparing statistically the active site configurations between populations of diverse structures in the Protein Data Bank, and it reveals patterns of hydrogen bonding that transcend enzyme families. Through analysis of large samples of structures, insights are drawn at a level of detail exceeding the experimental precision of an individual structure. In phosphagen kinases, for example, hydrogen bonds with the O3β of the nucleotide substrate are revealed as analogous to those in unrelated G proteins. In G proteins and other enzymes, interactions with O3β have been understood in terms of electrostatic favoring of the transition state. Ground state quantum mechanical calculations on model compounds show that the active site interactions highlighted in our database analysis can affect substrate phosphate charge and bond length, in ways that are consistent with prior experimental observations, by modulating hyperconjugative orbital interactions that weaken the scissile bond. Testing experimentally the inference about the importance of O3β interactions in phosphagen kinases, mutation of arginine kinase Arg280 decreases kcat, as predicted, with little impact upon KM.

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.

'Pathway drug cocktail': targeting Ras signaling based on structural pathways

Tumors bearing Ras mutations are notoriously difficult to treat. Drug combinations targeting the Ras protein or its pathway have also not met with success. 'Pathway drug cocktails', which are combinations aiming at parallel pathways, appear more promising; however, to be usefully exploited, a repertoire of classified pathway combinations is desirable. This challenge would be facilitated by the availability of the structural network of signaling pathways. When integrated with functional and systems level clinical data, they can be powerful in advancing novel therapeutic platforms. Based on structural knowledge, drug cocktails may tear into multiple cellular processes that drive tumorigenesis, and help in deciphering the interrelationship between Ras mutations and the rewired Ras network. The pathway drug cocktail paradigm can be applied to other signaling protein targets.

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.