Global conformational dynamics in ras

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.

Excited-state observation of active K-Ras reveals differential structural dynamics of wild-type versus oncogenic G12D and G12C mutants

Despite the prominent role of the K-Ras protein in many different types of human cancer, major gaps in atomic-level information severely limit our understanding of its functions in health and disease. Here, we report the quantitative backbone structural dynamics of K-Ras by solution nuclear magnetic resonance spectroscopy of the active state of wild-type K-Ras bound to guanosine triphosphate (GTP) nucleotide and two of its oncogenic P-loop mutants, G12D and G12C, using a new nanoparticle-assisted spin relaxation method, relaxation dispersion and chemical exchange saturation transfer experiments covering the entire range of timescales from picoseconds to milliseconds. Our combined experiments allow detection and analysis of the functionally critical Switch I and Switch II regions, which have previously remained largely unobservable by X-ray crystallography and nuclear magnetic resonance spectroscopy. Our data reveal cooperative transitions of K-Ras·GTP to a highly dynamic excited state that closely resembles the partially disordered K-Ras·GDP state. These results advance our understanding of differential GTPase activities and signaling properties of the wild type versus mutants and may thus guide new strategies for the development of therapeutics.

Reduced dynamic complexity allows structure elucidation of an excited state of KRASG13D

Localized dynamics of RAS, including regions distal to the nucleotide-binding site, is of high interest for elucidating the mechanisms by which RAS proteins interact with effectors and regulators and for designing inhibitors. Among several oncogenic mutants, methyl relaxation dispersion experiments reveal highly synchronized conformational dynamics in the active (GMPPNP-bound) KRAS^G13D, which suggests an exchange between two conformational states in solution. Methyl and ³¹P NMR spectra of active KRAS^G13D in solution confirm a two-state ensemble interconverting on the millisecond timescale, with a major P^γ atom peak corresponding to the dominant State 1 conformation and a secondary peak indicating an intermediate state different from the known State 2 conformation recognized by RAS effectors. High-resolution crystal structures of active KRAS^G13D and KRAS^G13D-RAF1 RBD complex provide snapshots of the State 1 and 2 conformations, respectively. We use residual dipolar couplings to solve and cross-validate the structure of the intermediate state of active KRAS^G13D, showing a conformation distinct from those of States 1 and 2 outside the known flexible switch regions. The dynamic coupling between the conformational exchange in the effector lobe and the breathing motion in the allosteric lobe is further validated by a secondary mutation in the allosteric lobe, which affects the conformational population equilibrium.

The Importance of Mg2+ -Free State in Nucleotide Exchange of Oncogenic K-Ras Mutants

For efficient targeting of oncogenic K-Ras interaction sites, a mechanistic picture of the Ras-cycle is necessary. Herein, we used NMR relaxation techniques and molecular dynamics simulations to decipher the role of slow dynamics in wild-type and three oncogenic P-loop mutants of K-Ras. Our measurements reveal a dominant two-state conformational exchange on the ms timescale in both GDP- and GTP-bound K-Ras. The identified low-populated higher energy state in GDP-loaded K-Ras has a conformation reminiscent of a nucleotide-bound/Mg²⁺ -free state characterized by shortened β2/β3-strands and a partially released switch-I region preparing K-Ras for the interaction with the incoming nucleotide exchange factor and subsequent reactivation. By providing insight into mutation-specific differences in K-Ras structural dynamics, our systematic analysis improves our understanding of prolonged K-Ras signaling and may aid the development of allosteric inhibitors targeting nucleotide exchange in K-Ras.

Insights into the Cross Talk between Effector and Allosteric Lobes of KRAS from Methyl Conformational Dynamics

KRAS is the most frequently mutated RAS protein in cancer patients, and it is estimated that about 20% of the cancer patients in the United States carried mutant RAS proteins. To accelerate therapeutic development, structures and dynamics of RAS proteins had been extensively studied by various biophysical techniques for decades. Although 31P NMR studies revealed population equilibrium of the two major states in the active GMPPNP-bound form, more complex conformational dynamics in RAS proteins and oncogenic mutants subtly modulate the interactions with their downstream effectors. We established a set of customized NMR relaxation dispersion techniques to efficiently and systematically examine the ms-μs conformational dynamics of RAS proteins. This method allowed us to observe varying synchronized motions that connect the effector and allosteric lobes in KRAS. We demonstrated the role of conformational dynamics of KRAS in controlling its interaction with the Ras-binding domain of the downstream effector RAF1, the first kinase in the MAPK pathway. This allows one to explain, as well as to predict, the altered binding affinities of various KRAS mutants, which was neither previously reported nor apparent from the structural perspective.

Intrinsically disordered proteins: Ensembles at the limits of Anfinsen's dogma

Intrinsically disordered proteins (IDPs) are proteins that lack rigid 3D structure. Hence, they are often misconceived to present a challenge to Anfinsen's dogma. However, IDPs exist as ensembles that sample a quasi-continuum of rapidly interconverting conformations and, as such, may represent proteins at the extreme limit of the Anfinsen postulate. IDPs play important biological roles and are key components of the cellular protein interaction network (PIN). Many IDPs can interconvert between disordered and ordered states as they bind to appropriate partners. Conformational dynamics of IDPs contribute to conformational noise in the cell. Thus, the dysregulation of IDPs contributes to increased noise and "promiscuous" interactions. This leads to PIN rewiring to output an appropriate response underscoring the critical role of IDPs in cellular decision making. Nonetheless, IDPs are not easily tractable experimentally. Furthermore, in the absence of a reference conformation, discerning the energy landscape representation of the weakly funneled IDPs in terms of reaction coordinates is challenging. To understand conformational dynamics in real time and decipher how IDPs recognize multiple binding partners with high specificity, several sophisticated knowledge-based and physics-based in silico sampling techniques have been developed. Here, using specific examples, we highlight recent advances in energy landscape visualization and molecular dynamics simulations to discern conformational dynamics and discuss how the conformational preferences of IDPs modulate their function, especially in phenotypic switching. Finally, we discuss recent progress in identifying small molecules targeting IDPs underscoring the potential therapeutic value of IDPs. Understanding structure and function of IDPs can not only provide new insight on cellular decision making but may also help to refine and extend Anfinsen's structure/function paradigm.

Novel approaches for the development of direct KRAS inhibitors: structural insights and drug design

INTRODUCTION

Hyperactivated RAS signaling is reported in 13% of all human cancers, in which ~80% resulted due to KRAS mutations alone. Direct inhibition of KRAS is an important aspect in treating KRAS-related tumors. Despite the efforts of more than four decades, not many KRAS inhibitors have been successful in obtaining clinical approval, except the very recent FDA approval for sotorasib. In recent years, the understanding of structural insights and allosteric pocket identification at catalytic sites of KRAS are likely to provide an excellent opportunity for the development of much more effective clinical candidates.

AREA COVERED

The presented review article mainly summarizes the developments of small molecule KRAS inhibitors as drug candidates and rational approaches that are being utilized for the selective targeting of KRAS signaling in the mutant cancer cells.

EXPERT OPINION

After the initial success in targeting the mutant KRAS G12C variants, the search has been shifted to address the challenges concerning the resistance and efficacy of small molecule KRAS inhibitors. However, the contribution of other KRAS mutations at G12V, G13C, and G13D variants causing cancers is much higher than the mutations at G12C. In view of this aspect, specific attention is required to target all other mutations as well. Accordingly, for the development of KRAS targeted therapies, the design of small molecule inhibitors that can inhibit KRAS signaling and as well as target inhibition of other signaling pathways like RAS-SOS and RAS-PI3K has to be explored extensively.

Unveiling the "invisible" druggable conformations of GDP-bound inactive Ras

The prevalent view on whether Ras is druggable has gradually changed in the recent decade with the discovery of effective inhibitors binding to cryptic sites unseen in the native structures. Despite the promising advances, therapeutics development toward higher potency and specificity is challenged by the elusive nature of these binding pockets. Here we derive a conformational ensemble of guanosine diphosphate (GDP)-bound inactive Ras by integrating spin relaxation-validated atomistic simulation with NMR chemical shifts and residual dipolar couplings, which provides a quantitative delineation of the intrinsic dynamics up to the microsecond timescale. The experimentally informed ensemble unequivocally demonstrates the preformation of both surface-exposed and buried cryptic sites in Ras•GDP, advocating design of inhibition by targeting the transient druggable conformers that are invisible to conventional experimental methods. The viability of the ensemble-based rational design has been established by retrospective testing of the ability of the Ras•GDP ensemble to identify known ligands from decoys in virtual screening.

Ras Isoforms from Lab Benches to Lives-What Are We Missing and How Far Are We?

The central protein in the oncogenic circuitry is the Ras GTPase that has been under intense scrutiny for the last four decades. From its discovery as a viral oncogene and its non-oncogenic contribution to crucial cellular functioning, an elaborate genetic, structural, and functional map of Ras is being created for its therapeutic targeting. Despite decades of research, there still exist lacunae in our understanding of Ras. The complexity of the Ras functioning is further exemplified by the fact that the three canonical Ras genes encode for four protein isoforms (H-Ras, K-Ras4A, K-Ras4B, and N-Ras). Contrary to the initial assessment that the H-, K-, and N-Ras isoforms are functionally similar, emerging data are uncovering crucial differences between them. These Ras isoforms exhibit not only cell-type and context-dependent functions but also activator and effector specificities on activation by the same receptor. Preferential localization of H-, K-, and N-Ras in different microdomains of the plasma membrane and cellular organelles like Golgi, endoplasmic reticulum, mitochondria, and endosome adds a new dimension to isoform-specific signaling and diverse functions. Herein, we review isoform-specific properties of Ras GTPase and highlight the importance of considering these towards generating effective isoform-specific therapies in the future.

Engineering subtilisin proteases that specifically degrade active RAS

We describe the design, kinetic properties, and structures of engineered subtilisin proteases that degrade the active form of RAS by cleaving a conserved sequence in switch 2. RAS is a signaling protein that, when mutated, drives a third of human cancers. To generate high specificity for the RAS target sequence, the active site was modified to be dependent on a cofactor (imidazole or nitrite) and protease sub-sites were engineered to create a linkage between substrate and cofactor binding. Selective proteolysis of active RAS arises from a 2-step process wherein sub-site interactions promote productive binding of the cofactor, enabling cleavage. Proteases engineered in this way specifically cleave active RAS in vitro, deplete the level of RAS in a bacterial reporter system, and also degrade RAS in human cell culture. Although these proteases target active RAS, the underlying design principles are fundamental and will be adaptable to many target proteins.

Dynamically encoded reactivity of Ras enzymes: opening new frontiers for drug discovery

Decoding molecular flexibility in order to understand and predict biological processes-applying the principles of dynamic-structure-activity relationships (DSAR)-becomes a necessity when attempting to design selective and specific inhibitors of a protein that has overlapping interaction surfaces with its upstream and downstream partners along its signaling cascade. Ras proteins are molecular switches that meet this definition perfectly. The close-lying P-loop and the highly flexible switch I and switch II regions are the site of nucleotide-, assisting-, and effector-protein binding. Oncogenic mutations that also appear in this region do not cause easily characterized overall structural changes, due partly to the inherent conformational heterogeneity and pliability of these segments. In this review, we present an overview of the results obtained using approaches targeting Ras dynamics, such as nuclear magnetic resonance (NMR) measurements and experiment-based modeling calculations (mostly molecular dynamics (MD) simulations). These methodologies were successfully used to decipher the mutant- and isoform-specific nature of certain transient states, far-lying allosteric sites, and the internal interaction networks, as well as the interconnectivity of the catalytic and membrane-binding regions. This opens new therapeutic potential: the discovered interaction hotspots present hitherto not targeted, selective sites for drug design efforts in diverse locations of the protein matrix.

Millisecond Allosteric Dynamics of Activated Ras Reproduced with a Slowly Hydrolyzable GTP Analogue

The millisecond timescale dynamics of activated Ras transiently sample a low-populated conformational state that has distinct surface property from the major state and represents a promising target for binding of small-molecule compounds. To avoid the complications of hydrolysis, dynamics and other properties of active Ras have so far been routinely investigated by using non-hydrolyzable GTP analogues, which, however, were previously reported to alter both the kinetics and distribution of the conformational exchange. In this study, we quantitatively measured and validated the internal dynamics of Ras complexed with a slowly hydrolyzable GTP analogue, GTPγS, which increases the lifetime of active Ras by 23 times relative to that of native GTP. It was found that GTPγS, in addition to its better mimicking of the exchange kinetics than the commonly used non-hydrolyzable analogues GppNHp and GppCH₂ p, can rigorously reproduce the natural dynamics network in active Ras, thus indicating its fitness for use in the development of allosteric inhibitors.

Inhibition of RAS: proven and potential vulnerabilities

RAS is a membrane localized small GTPase frequently mutated in human cancer. As such, RAS has been a focal target for developing cancer therapeutics since its discovery nearly four decades ago. However, efforts to directly target RAS have been challenging due to the apparent lack of readily discernable deep pockets for binding small molecule inhibitors leading many to consider RAS as undruggable. An important milestone in direct RAS inhibition was achieved recently with the groundbreaking discovery of covalent inhibitors that target the mutant Cys residue in KRAS(G12C). Surprisingly, these G12C-reactive compounds only target mutant RAS in the GDP-bound state thereby locking it in the inactive conformation and blocking its ability to couple with downstream effector pathways. Building on this success, several groups have developed similar compounds that selectively target KRAS(G12C), with AMG510 and MRTX849 the first to advance to clinical trials. Both have shown early promising results. Though the success with these compounds has reignited the possibility of direct pharmacological inhibition of RAS, these covalent inhibitors are limited to treating KRAS(G12C) tumors which account for <15% of all RAS mutants in human tumors. Thus, there remains an unmet need to identify more broadly efficacious RAS inhibitors. Here, we will discuss the current state of RAS(G12C) inhibitors and the potential for inhibiting additional RAS mutants through targeting RAS dimerization which has emerged as an important step in the allosteric regulation of RAS function.

In-depth NMR characterization of Rab4a structure, nucleotide exchange and hydrolysis kinetics reveals an atypical GTPase profile

Rab4a is a small GTPase associated with endocytic compartments and a key regulator of early endosomes recycling. Gathering evidence indicates that its expression and activation are required for the development of metastases. Rab4a-intrinsic GTPase properties that control its activity, i.e. nucleotide exchange and hydrolysis rates, have not yet been thoroughly studied. The determination of these properties is of the utmost importance to understand its functions and contributions to tumorigenesis. Here, we used the constitutively active (Rab4aQ67L) and dominant negative (Rab4aS22N) mutants to characterize the thermodynamical and structural determinants of the interaction between Rab4a and GTP (GTPγS) as well as GDP. We report the first ¹H, ¹³C, ¹⁵N backbone NMR assignments of a Rab GTPase family member with Rab4a in complex with GDP and GTPγS. We also provide a qualitative description of the extent of structural and dynamical changes caused by the Q67L and S22N mutations. Using a real-time NMR approach and the two aforementioned mutants as controls, we evaluated Rab4a intrinsic nucleotide exchange and hydrolysis rates. Compared to most small GTPases such as Ras, a rapid GTP exchange rate along with slow hydrolysis rate were observed. This suggests that, in a cellular context, Rab4a can self-activate and persist in an activated state in absence of regulatory mechanisms. This peculiar profile is uncommon among the Ras superfamily members, making Rab4a an atypical fast-cycling GTPase and may explain, at least in part, how it contributes to metastases.

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.

Biology, pathology, and therapeutic targeting of RAS

RAS was identified as a human oncogene in the early 1980s and subsequently found to be mutated in nearly 30% of all human cancers. More importantly, RAS plays a central role in driving tumor development and maintenance. Despite decades of effort, there remain no FDA approved drugs that directly inhibit RAS. The prevalence of RAS mutations in cancer and the lack of effective anti-RAS therapies stem from RAS' core role in growth factor signaling, unique structural features, and biochemistry. However, recent advances have brought promising new drugs to clinical trials and shone a ray of hope in the field. Here, we will exposit the details of RAS biology that illustrate its key role in cell signaling and shed light on the difficulties in therapeutically targeting RAS. Furthermore, past and current efforts to develop RAS inhibitors will be discussed in depth.

NMR-Derived Conformational Ensemble of State 1 of Activated Ras Reveals Insights into a Druggable Pocket

The lack of apparent pockets in the ground conformation of Ras has long challenged the rational design of inhibitors against this oncogenic protein. The sparsely populated, transiently formed state 1 of activated Ras, on the other hand, shows appreciable surface roughness and is increasingly recognized as a potential target for drug discovery. State 1, however, is extremely flexible, and a static structure cannot fully unveil its conformational space that can be exploited for drug design. Here, we present a conformational ensemble of state 1 that was derived using chemical shift-based modeling. The ensemble reveals the intrinsic plasticity of a druggable pocket in state 1 and demonstrates the mechanism of conformational selection for inhibitor recognition. The large set of structural templates in the ensemble, providing a comprehensive description of thermally accessible pocket conformations, is expected to significantly aid the rational design of anti-Ras drugs.

Therapeutic targeting of RAS: New hope for drugging the "undruggable"

RAS is the most frequently mutated oncogene in cancer and a critical driver of oncogenesis. Therapeutic targeting of RAS has been a goal of cancer research for more than 30 years due to its essential role in tumor formation and maintenance. Yet the quest to inhibit this challenging foe has been elusive. Although once considered "undruggable", the struggle to directly inhibit RAS has seen recent success with the development of pharmacological agents that specifically target the KRAS(G12C) mutant protein, which include the first direct RAS inhibitor to gain entry to clinical trials. However, the limited applicability of these inhibitors to G12C-mutant tumors demands further efforts to identify more broadly efficacious RAS inhibitors. Understanding allosteric influences on RAS may open new avenues to inhibit RAS. Here, we provide a brief overview of RAS biology and biochemistry, discuss the allosteric regulation of RAS, and summarize the various approaches to develop RAS inhibitors.

Oncogenic G12D mutation alters local conformations and dynamics of K-Ras

K-Ras is the most frequently mutated oncoprotein in human cancers, and G12D is its most prevalent mutation. To understand how G12D mutation impacts K-Ras function, we need to understand how it alters the regulation of its dynamics. Here, we present local changes in K-Ras structure, conformation and dynamics upon G12D mutation, from long-timescale Molecular Dynamics simulations of active (GTP-bound) and inactive (GDP-bound) forms of wild-type and mutant K-Ras, with an integrated investigation of atomistic-level changes, local conformational shifts and correlated residue motions. Our results reveal that the local changes in K-Ras are specific to bound nucleotide (GTP or GDP), and we provide a structural basis for this. Specifically, we show that G12D mutation causes a shift in the population of local conformational states of K-Ras, especially in Switch-II (SII) and α3-helix regions, in favor of a conformation that is associated with a catalytically impaired state through structural changes; it also causes SII motions to anti-correlate with other regions. This detailed picture of G12D mutation effects on the local dynamic characteristics of both active and inactive protein helps enhance our understanding of local K-Ras dynamics, and can inform studies on the development of direct inhibitors towards the treatment of K-Ras^G12D-driven cancers.

Extending the Lifetime of Native GTP-Bound Ras for Site-Resolved NMR Measurements: Quantifying the Allosteric Dynamics

Characterization of native GTP-bound Ras is important for an appreciation of its cellular signaling and for the design of inhibitors, which however has been depressed by its intrinsic instability. Herein, an effective approach for extending the lifetime of Ras⋅GTP samples by exploiting the active role of Son of Sevenless (Sos) is demonstrated that sustains the activated state of Ras. This approach, combined with a postprocessing method that suppresses residual Ras⋅GDP signals, is applied to the site-resolved NMR measurement of the allosteric dynamics of Ras⋅GTP. The observed network of concerted motions well covers the recently identified allosteric inhibitor-binding pockets, but the motions are more confined than those of Ras⋅GppNHp, advocating the use of native GTP for development of allosteric inhibitors. The Sos-based approach is anticipated to generally facilitate experiments on active Ras when native GTP is preferred.

Pharmacological targeting of RAS: Recent success with direct inhibitors

RAS has long been viewed as undruggable due to its lack of deep pockets for binding of small molecule inhibitors. However, recent successes in the development of direct RAS inhibitors suggest that the goal of pharmacological inhibition of RAS in patients may soon be realized. This review will discuss the role of RAS in cancer, the approaches used to develop direct RAS inhibitors, and highlight recent successes in the development of novel RAS inhibitory compounds that target different aspects of RAS biochemistry. In particular, this review will discuss the different properties of RAS that have been targeted by various inhibitors including membrane localization, the different activation states of RAS, effector binding, and nucleotide exchange. In addition, this review will highlight the recent success with mutation-specific inhibitors that exploit the unique biochemistry of the RAS(G12C) mutant. Although this mutation in KRAS accounts for 11% of all KRAS mutations in cancer, it is the most prominent KRAS mutant in lung cancer suggesting that G12C-specific inhibitors may provide a new approach for treating the subset of lung cancer patients harboring this mutant allele. Finally, this review will discuss the involvement of dimerization in RAS function and highlight new approaches to inhibit RAS by specifically interfering with RAS:RAS interaction.

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.

K-Ras Populates Conformational States Differently from Its Isoform H-Ras and Oncogenic Mutant K-RasG12D

Structures of wild-type K-Ras from crystals obtained in the presence of guanosine triphosphate (GTP) or its analogs have remained elusive. Of the K-Ras mutants, only K-RasG12D and K-RasQ61H are available in the PDB representing the activated form of the GTPase not in complex with other proteins. We present the crystal structure of wild-type K-Ras bound to the GTP analog GppCH₂p, with K-Ras in the state 1 conformation. Signatures of conformational states obtained by one-dimensional proton NMR confirm that K-Ras has a more substantial population of state 1 in solution than H-Ras, which predominantly favors state 2. The oncogenic mutant K-RasG12D favors state 2, changing the balance of conformational states in favor of interactions with effector proteins. Differences in the population of conformational states between K-Ras and H-Ras, as well as between K-Ras and its mutants, can provide a structural basis for focused targeting of the K-Ras isoform in cancer-specific strategies.

Study of Förster Resonance Energy Transfer to Lipid Domain Markers Ascertains Partitioning of Semisynthetic Lipidated N-Ras in Lipid Raft Nanodomains

Cellular membranes are heterogeneous planar lipid bilayers displaying lateral phase separation with the nanometer-scale liquid-ordered phase (also known as "lipid rafts") surrounded by the liquid-disordered phase. Many membrane-associated proteins were found to permanently integrate into the lipid rafts, which is critical for their biological function. Isoforms H and N of Ras GTPase possess a unique ability to switch their lipid domain preference depending on the type of bound guanine nucleotide (GDP or GTP). This behavior, however, has never been demonstrated in vitro in model bilayers with recombinant proteins and therefore has been attributed to the action of binding of Ras to other proteins at the membrane surface. In this paper, we report the observation of the nucleotide-dependent switch of lipid domain preferences of the semisynthetic lipidated N-Ras in lipid raft vesicles in the absence of additional proteins. To detect segregation of Ras molecules in raft and disordered lipid domains, we measured Förster resonance energy transfer between the donor fluorophore, mant, attached to the protein-bound guanine nucleotides, and the acceptor, rhodamine-conjugated lipid, localized into the liquid-disordered domains. Herein, we established that N-Ras preferentially populated raft domains when bound to mant-GDP, while losing its preference for rafts when it was associated with a GTP mimic, mant-GppNHp. At the same time, the isolated lipidated C-terminal peptide of N-Ras was found to be localized outside of the liquid-ordered rafts, most likely in the bulk-disordered lipid. Substitution of the N-terminal G domain of N-Ras with a homologous G domain of H-Ras disrupted the nucleotide-dependent lipid domain switch.

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.

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.

Dynamic studies of H-Ras•GTPγS interactions with nucleotide exchange factor Sos reveal a transient ternary complex formation in solution

The cycling between GDP- and GTP- bound forms of the Ras protein is partly regulated by the binding of Sos. The structural/dynamic behavior of the complex formed between activated Sos and Ras at the point of the functional cycle where the nucleotide exchange is completed has not been described to date. Here we show that solution NMR spectra of H-Ras∙GTPγS mixed with a functional fragment of Sos (Sos^Cat) at a 2:1 ratio are consistent with the formation of a rather dynamic assembly. H-Ras∙GTPγS binding was in fast exchange on the NMR timescale and retained a significant degree of molecular tumbling independent of Sos^Cat, while Sos^Cat also tumbled largely independently of H-Ras. Estimates of apparent molecular weight from both NMR data and SEC-MALS revealed that, at most, only one H-Ras∙GTPγS molecule appears stably bound to Sos. The weak transient interaction between Sos and the second H-Ras∙GTPγS may provide a necessary mechanism for complex dissociation upon the completion of the native GDP → GTP exchange reaction, but also explains measurable GTP → GTP exchange activity of Sos routinely observed in in vitro assays that use fluorescently-labelled analogs of GTP. Overall, the data presents the first dynamic snapshot of Ras functional cycle as controlled by Sos.

The Ras G Domain Lacks the Intrinsic Propensity to Form Dimers

Ras GTPase is a molecular switch controlling a number of cellular pathways including growth, proliferation, differentiation, and apoptosis. Recent reports indicated that Ras undergoes dimerization at the membrane surface through protein-protein interactions. If firmly established this property of Ras would require profound reassessment of a large amount of published data and modification of the Ras signaling paradigm. One proposed mechanism of dimerization involves formation of salt bridges between the two GTPase domains (G domains) leading to formation of a compact dimer as observed in Ras crystal structures. In this work, we interrogated the intrinsic ability of Ras to self-associate in solution by creating conditions of high local concentration through irreversibly tethering the two G domains together at their unstructured C-terminal tails. We evaluated possible self-association in this inverted tandem conjugate via analysis of the time-domain fluorescence anisotropy and NMR chemical shift perturbations. We did not observe the increased rotational correlation time expected for the G domain dimer. Variation of the ionic strength (to modulate stability of the salt bridges) did not affect the rotational correlation time in the tandem further supporting independent rotational diffusion of two G domains. In a parallel line of experiments to detect and map weak self-association of the G domains, we analyzed NMR chemical shifts perturbations at a number of sites near the crystallographic dimer interface. The nearly complete lack of chemical shift perturbations in the tandem construct supported a simple model with the independent G domains repelled from each other by their overall negative charge. These results lead us to the conclusion that self-association of the G domains cannot be responsible for homodimerization of Ras reported in the literature.

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."

Computational allosteric ligand binding site identification on Ras proteins

A number of computational techniques have been proposed to expedite the process of allosteric ligand binding site identification in inherently flexible and hence challenging drug targets. Some of these techniques have been instrumental in the discovery of allosteric ligand binding sites on Ras proteins, a group of elusive anticancer drug targets. This review provides an overview of these techniques and their application to Ras proteins. A summary of molecular docking and binding site identification is provided first, followed by a more detailed discussion of two specific techniques for binding site identification in ensembles of Ras conformations generated by molecular simulations.

Mapping the Conformation Space of Wildtype and Mutant H-Ras with a Memetic, Cellular, and Multiscale Evolutionary Algorithm

An important goal in molecular biology is to understand functional changes upon single-point mutations in proteins. Doing so through a detailed characterization of structure spaces and underlying energy landscapes is desirable but continues to challenge methods based on Molecular Dynamics. In this paper we propose a novel algorithm, SIfTER, which is based instead on stochastic optimization to circumvent the computational challenge of exploring the breadth of a protein's structure space. SIfTER is a data-driven evolutionary algorithm, leveraging experimentally-available structures of wildtype and variant sequences of a protein to define a reduced search space from where to efficiently draw samples corresponding to novel structures not directly observed in the wet laboratory. The main advantage of SIfTER is its ability to rapidly generate conformational ensembles, thus allowing mapping and juxtaposing landscapes of variant sequences and relating observed differences to functional changes. We apply SIfTER to variant sequences of the H-Ras catalytic domain, due to the prominent role of the Ras protein in signaling pathways that control cell proliferation, its well-studied conformational switching, and abundance of documented mutations in several human tumors. Many Ras mutations are oncogenic, but detailed energy landscapes have not been reported until now. Analysis of SIfTER-computed energy landscapes for the wildtype and two oncogenic variants, G12V and Q61L, suggests that these mutations cause constitutive activation through two different mechanisms. G12V directly affects binding specificity while leaving the energy landscape largely unchanged, whereas Q61L has pronounced, starker effects on the landscape. An implementation of SIfTER is made available at http://www.cs.gmu.edu/~ashehu/?q=OurTools. We believe SIfTER is useful to the community to answer the question of how sequence mutations affect the function of a protein, when there is an abundance of experimental structures that can be exploited to reconstruct an energy landscape that would be computationally impractical to do via Molecular Dynamics.

Some (dis)assembly required: partial unfolding in the Par-6 allosteric switch

Allostery is commonly described as a functional connection between two distant sites in a protein, where a binding event at one site alters affinity at the other. Here we review the conformational dynamics that encode an allosteric switch in the PDZ domain of Par-6. Par-6 is a scaffold protein that organizes other proteins into a complex required to initiate and maintain cell polarity. NMR measurements revealed that the PDZ domain samples an evolutionarily conserved unfolding intermediate allowing rearrangement of two adjacent loop residues that control ligand binding affinity. Cdc42 binding to Par-6 creates a novel interface between the PDZ domain and the adjoining CRIB motif that stabilizes the high-affinity PDZ conformation. Thermodynamic and kinetic studies suggest that partial PDZ unfolding is an integral part of the Par-6 switching mechanism. The Par-6 CRIB-PDZ module illustrates two important structural aspects of protein evolution: the interface between adjacent domains in the same protein can give rise to allosteric regulation, and thermodynamic stability may be sacrificed to increase the sampling frequency of an unfolding intermediate required for conformational switching.

Thermodynamic mapping of effector protein interfaces with RalA and RalB

RalA and RalB are members of the Ras family of small G proteins and are activated downstream of Ras via RalGEFs. The RalGEF-Ral axis represents one of the major effector pathways controlled by Ras and as such is an important pharmacological target. RalA and RalB are approximately 80% identical at the amino acid level; despite this, they have distinct roles both in normal cells and in the disease state. We have used our structure of RalB-RLIP76 to guide an analysis of Ral-effector interaction interfaces, creating panels of mutant proteins to probe the energetics of these interactions. The data provide a physical mechanism that underpins the effector selective mutations commonly employed to dissect Ral G protein function. Comparing the energetic landscape of the RalB-RLIP76 and RalB-Sec5 complexes reveals mutations in RalB that lead to differential binding of the two effector proteins. A panel of RLIP76 mutants was used to probe the interaction between RLIP76 and RalA and -B. Despite 100% sequence identity in the RalA and -B contact residues with RLIP76, differences still exist in the energetic profiles of the two complexes. Therefore, we have revealed properties that may account for some of the functional separation observed with RalA and RalB at the cellular level. Our mutations, in both the Ral isoforms and RLIP76, provide new tools that can be employed to parse the complex biology of Ral G protein signaling networks. The combination of these thermodynamic and structural data can also guide efforts to ablate RalA and -B activity with small molecules and peptides.

The Ras-Membrane Interface: Isoform-specific Differences in The Catalytic Domain

The small GTPase Ras is mutated in about 20% of human cancers, primarily at active site amino acid residues G12, G13, and Q61. Thus, structural biology research has focused on the active site, impairment of GTP hydrolysis by oncogenic mutants, and characterization of protein-protein interactions in the effector lobe half of the protein. The C-terminal hypervariable region has increasingly gained attention due to its importance in H-Ras, N-Ras, and K-Ras differences in membrane association. A high-resolution molecular view of the Ras-membrane interaction involving the allosteric lobe of the catalytic domain has lagged behind, although evidence suggests that it contributes to isoform specificity. The allosteric lobe has recently gained interest for harboring potential sites for more selective targeting of this elusive "undruggable" protein. The present review reveals critical insight that isoform-specific differences appear prominently at these potentially targetable sites and integrates these differences with knowledge of Ras plasma membrane localization, with the intent to better understand the structure-function relationships needed to design isoform-specific Ras inhibitors.

High Pressure NMR Methods for Characterizing Functional Substates of Proteins

Proteins usually exist in multiple conformational states in solution. High pressure NMR spectroscopy is a well-suited method to identify these states. In addition, these states can be characterized by their thermodynamic parameters, the free enthalpies at ambient pressure, the partial molar volumes, and the partial molar compressibility that can be obtained from the analysis of the high pressure NMR data. Two main types of states of proteins exist, functional states and folding states. There is a strong link between these two types, the functional states represent essential folding states (intermediates), other folding states may have no functional meaning (optional folding states). In this chapter, this concept is tested on the Ras protein, an important proto-oncogen in humans where all substates required by theory can be identified experimentally by high pressure NMR spectroscopy. Finally, we show how these data can be used to develop allosteric inhibitors of proteins.

Fluorescence2D: software for accelerated acquisition and analysis of two-dimensional fluorescence spectra

The Fluorescence2D is free software that allows analysis of two-dimensional fluorescence spectra obtained using the accelerated “triangular” acquisition schemes. The software is a combination of Python and MATLAB-based programs that perform conversion of the triangular data, display of the two-dimensional spectra, extraction of 1D slices at different wavelengths, and output in various graphic formats.

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.

Lessons from computer simulations of Ras proteins in solution and in membrane

BACKGROUND

A great deal has been learned over the last several decades about the function of Ras proteins in solution and membrane environments. While much of this knowledge has been derived from a plethora of experimental techniques, computer simulations have also played a substantial role.

SCOPE OF REVIEW

Our goal here is to summarize the contribution of molecular simulations to our current understanding of normal and aberrant Ras function. We focus on lessons from molecular dynamics simulations in aqueous and membrane environments.

MAJOR CONCLUSIONS

The central message is that a close interaction between theory and simulation on the one hand and cell-biological, spectroscopic and other experimental approaches on the other has played, and will likely continue to play, a vital role in Ras research.

GENERAL SIGNIFICANCE

Atomistic insights emerging from detailed simulations of Ras in solution and in bilayers may be the key to unlock the secret that to date prevented development of selective anti-Ras inhibitors for cancer therapy.

Allosteric activation of the Par-6 PDZ via a partial unfolding transition

Proteins exist in a delicate balance between the native and unfolded states, where thermodynamic stability may be sacrificed to attain the flexibility required for efficient catalysis, binding, or allosteric control. Partition-defective 6 (Par-6) regulates the Par polarity complex by transmitting a GTPase signal through the Cdc42/Rac interaction binding PSD-95/Dlg/ZO-1 (CRIB-PDZ) module that alters PDZ ligand binding. Allosteric activation of the PDZ is achieved by local rearrangement of the L164 and K165 side chains to stabilize the interdomain CRIB:PDZ interface and reposition a conserved element of the ligand binding pocket. However, microsecond to millisecond dynamics measurements revealed that L164/K165 exchange requires a larger rearrangement than expected. The margin of thermodynamic stability for the PDZ domain is modest (∼3 kcal/mol) and further reduced by transient interactions with the disordered CRIB domain. Measurements of local structural stability revealed that tertiary contacts within the PDZ are disrupted by a partial unfolding transition that enables interconversion of the L/K switch. The unexpected participation of partial PDZ unfolding in the allosteric mechanism of Par-6 suggests that native-state unfolding may be essential for the function of other marginally stable proteins.

Characterization of the second ion-binding site in the G domain of H-Ras

Ras is a small monomeric GTPase acting as molecular switch in multiple cellular processes. The N-terminal G domain of Ras binds GTP or GDP accompanied by a magnesium ion, which is strictly required for GTPase activity and performs a structural role. Another ion-binding site on the opposite face of the G domain has been recently observed to specifically associate with calcium acetate in the crystal [Buhrman, G., et al. (2010) Proc. Natl. Aacd. Sci. U.S.A. 107, 4931-4936]. In this article, we report thermodynamic measurements of the affinity and specificity of the remote ion-binding site in H-Ras as observed in solution. Using (15)N-(1)H nuclear magnetic resonance spectroscopy, we determined that, in contrast to the crystalline state, the remote site in solution is specific for a divalent cation, binding both calcium and magnesium with anions playing a minimal role. The affinity of the remote site for divalent cations is in the low millimolar range and remarkably different for GDP- and GppNHp-bound forms of the G domain, indicating that the GTP-binding pocket and the remote site are allosterically coupled through the distance of more than 25 Å. Considering that the remote site is oriented toward the membrane surface in vivo, we hypothesize that its cognate biological ligand might be a positively charged group extending from a lipid or an integral membrane protein.

Ras inhibition via direct Ras binding--is there a path forward?

Three decades after identification of the Ras oncogene, no effective treatments for Ras mutant tumors are available despite intensive drug discovery efforts. Here we critically review the attempts to inhibit Ras function via direct binding of small molecules at the Ras surface with the aim to disrupt its interaction with other proteins. Multiple binders at different binding sites have been discovered, and recent efforts afforded crystal structures of Ras-binder complexes. Albeit with low affinities, many of the binders were shown to impart inhibitory activities, and inhibition of nucleotide exchange as a consequence of disrupting the Ras-SOS interaction has been the most commonly identified mode of action. We see two key challenges in the development of these early starting points: Enhancing binding affinities and achieving selectivity, both against other GTPases and for mutant Ras over the wildtype form. In light of the large unmet medical need, we encourage the continued search for functionally active Ras binders, and we believe that integrated use of biophysical and biochemical tools will provide the highest chances for success. Given the failures experienced in the past and the significant hurdles ahead, we propose that this challenge be tackled through alliances between industry and academia.