Structural impact of GTP binding on downstream KRAS signaling

NMR 1H, 13C, 15N backbone resonance assignments of wild-type human K-Ras and its oncogenic mutants G12D and G12C bound to GTP

Human K-Ras protein, which is a member of the GTPase Ras family, hydrolyzes GTP to GDP and concomitantly converts from its active to its inactive state. It is a key oncoprotein, because several mutations, particularly those at residue position 12, occur with a high frequency in a wide range of human cancers. The K-Ras protein is therefore an important target for developing therapeutic anti-cancer agents. In this work we report the almost complete sequence-specific resonance assignments of wild-type and the oncogenic G12C and G12D mutants in the GTP-complexed active forms, including the functionally important Switch I and Switch II regions. These assignments serve as the basis for a comprehensive functional dynamics study of wild-type K-Ras and its G12 mutants.

Conserved allosteric perturbation of the GTPase domains by region 1 of Ras hypervariable regions

Ras proteins are important intracellular signaling hubs that can interact with numerous downstream effectors and upstream regulators through their GTPase domains (G-domains) anchored to plasma membranes by the C-terminal hypervariable regions (HVRs). The biological functions of Ras were proposed to be regulated at multiple levels including the intramolecular G-domain-HVR interactions, of which the exact mechanism and specificity are still controversial. Here, we demonstrate that the HVRs, instead of having direct contacts, can weakly perturb the G-domains via an allosteric interaction that is restricted to a ∼20 Å range and highly conserved in the tested Ras isoforms (HRas and KRas4B) and nucleotide-bound states. The origin of this allosteric perturbation has been localized to a short segment (residues 167-171) coinciding with region 1 of HVRs, which exhibits moderate to weak α-helical propensities. A charge-reversal mutation (E168K) of KRas4B in region 1, previously described in the Catalog of Somatic Mutations in Cancer database, was found to induce similar chemical shift perturbations as truncation of the HVR does. Further membrane paramagnetic relaxation enhancement (mPRE) data show that this region 1 mutation alters the membrane orientations of KRas4B and moderately increases the relative population of the signaling-compatible state.

Early pharmacological profiling of isatin derivatives as potent and selective cytotoxic agents

Isatin derivatives have attracted a lot of interest for their potential in the development of new anticancer drugs. A library of 38 isatin derivatives, created through an Ugi four-component reaction, underwent an initial screening in a panel of six human solid tumor cell lines. The four most active derivatives were then selected for further testing. These compounds showed selectivity towards the non-small cell lung cancer (NSCLC) cell line SW1573, whilst NSCLC A549 cells were barely affected. The combination of phenotypic assays, including wound healing, clonogenic and continuous live cell imaging provided a deeper understanding of the compounds' mode of action. In particular, the latter demonstrated that isatin derivatives were able to induce necroptosis in SW1573 cells. The kinetics of cell death showed that necroptosis appeared after 2.5 h of exposure, which could be delayed to 7 h when co-treated with necrostatin-1. Interaction between the isatin derivatives and the KRAS G12C protein variant was discarded after in silico studies. Further studies are warranted to identify the cellular target responsible for the observed selectivity among cell lines.

Revealing the mechanism of action of a first-in-class covalent inhibitor of KRASG12C (ON) and other functional properties of oncogenic KRAS by 31P NMR

Individual oncogenic KRAS mutants confer distinct differences in biochemical properties and signaling for reasons that are not well understood. KRAS activity is closely coupled to protein dynamics and is regulated through two interconverting conformations: state 1 (inactive, effector binding deficient) and state 2 (active, effector binding enabled). Here, we use 31P NMR to delineate the differences in state 1 and state 2 populations present in WT and common KRAS oncogenic mutants (G12C, G12D, G12V, G13D, and Q61L) bound to its natural substrate GTP or a commonly used nonhydrolyzable analog GppNHp (guanosine-5'-[(β,γ)-imido] triphosphate). Our results show that GppNHp-bound proteins exhibit significant state 1 population, whereas GTP-bound KRAS is primarily (90% or more) in state 2 conformation. This observation suggests that the predominance of state 1 shown here and in other studies is related to GppNHp and is most likely nonexistent in cells. We characterize the impact of this differential conformational equilibrium of oncogenic KRAS on RAF1 kinase effector RAS-binding domain and intrinsic hydrolysis. Through a KRAS G12C drug discovery, we have identified a novel small-molecule inhibitor, BBO-8956, which is effective against both GDP- and GTP-bound KRAS G12C. We show that binding of this inhibitor significantly perturbs state 1-state 2 equilibrium and induces an inactive state 1 conformation in GTP-bound KRAS G12C. In the presence of BBO-8956, RAF1-RAS-binding domain is unable to induce a signaling competent state 2 conformation within the ternary complex, demonstrating the mechanism of action for this novel and active-conformation inhibitor.

Far-reaching effects of tyrosine64 phosphorylation on Ras revealed with BeF3- complexes

Tyrosine phosphorylation on Ras by Src kinase is known to uncouple Ras from upstream regulation and downstream communication. However, the mechanisms by which phosphorylation modulates these interactions have not been detailed. Here, the major mono-phosphorylation level on tyrosine64 is quantified by 31P NMR and mutagenesis. Crystal structures of unphosphorylated and tyrosine64-phosphorylated Ras in complex with a BeF3- ground state analogue reveal "closed" Ras conformations very different from those of the "open" conformations previously observed for non-hydrolysable GTP analogue structures of Ras. They deliver new mechanistic and conformational insights into intrinsic GTP hydrolysis. Phosphorylation of tyrosine64 delivers conformational changes distant from the active site, showing why phosphorylated Ras has reduced affinity to its downstream effector Raf. 19F NMR provides evidence for changes in the intrinsic GTPase and nucleotide exchange rate and identifies the concurrent presence of a major "closed" conformation alongside a minor yet functionally important "open" conformation at the ground state of Ras. This study expands the application of metal fluoride complexes in revealing major and minor conformational changes of dynamic and modified Ras proteins.

Bioinformatics and Experimental Validation for Identifying Biomarkers Associated with AMG510 (Sotorasib) Resistance in KRASG12C-Mutated Lung Adenocarcinoma

The Kirsten rat sarcoma viral oncogene homolog (KRAS)G12C mutation is prevalent in lung adenocarcinoma (LUAD), driving tumor progression and indicating a poor prognosis. While the FDA-approved AMG510 (Sotorasib) initially demonstrated efficacy in treating KRASG12C-mutated LUAD, resistance emerged within months. Data from AMG510 treatment-resistant LUAD (GSE204753) and single-cell datasets (GSE149655) were analyzed. Gene set variation analysis (GSVA) and gene set enrichment analysis (GSEA) were used to explore enriched signaling pathways, nomogram models were constructed, and transcription factors predicting resistance biomarkers were predicted. CIBERSORT identified immune cell subpopulations, and their association with resistance biomarkers was assessed through single-cell analysis. AMG510-resistant LUAD cells (H358-AR) were constructed, and proliferative changes were evaluated using a CCK-8 assay. Key molecules for AMG510 resistance, including SLC2A1, TLE1, FAM83A, HMGA2, FBXO44, and MTRNR2L12, were recognized. These molecules impacted multiple signaling pathways and the tumor microenvironment and were co-regulated by various transcription factors. Single-cell analysis revealed a dampening effect on immune cell function, with associations with programmed cell death ligand 1 (PDL1) expression, cytokine factors, and failure factors. The findings indicate that these newly identified biomarkers are linked to the abnormal expression of PDL1 and have the potential to induce resistance through immunosuppression. These results highlight the need for further research and therapeutic intervention to address this issue effectively.

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

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

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.

Mechanism-Based Redesign of GAP to Activate Oncogenic Ras

Ras GTPases play a crucial role in cell signaling pathways. Mutations of the Ras gene occur in about one third of cancerous cell lines and are often associated with detrimental clinical prognosis. Hot spot residues Gly12, Gly13, and Gln61 cover 97% of oncogenic mutations, which impair the enzymatic activity in Ras. Using QM/MM free energy calculations, we present a two-step mechanism for the GTP hydrolysis catalyzed by the wild-type Ras.GAP complex. We found that the deprotonation of the catalytic water takes place via the Gln61 as a transient Brønsted base. We also determined the reaction profiles for key oncogenic Ras mutants G12D and G12C using QM/MM minimizations, matching the experimentally observed loss of catalytic activity, thereby validating our reaction mechanism. Using the optimized reaction paths, we devised a fast and accurate procedure to design GAP mutants that activate G12D Ras. We replaced GAP residues near the active site and determined the activation barrier for 190 single mutants. We furthermore built a machine learning for ultrafast screening, by fast prediction of the barrier heights, tested both on the single and double mutations. This work demonstrates that fast and accurate screening can be accomplished via QM/MM reaction path optimizations to design protein sequences with increased catalytic activity. Several GAP mutations are predicted to re-enable catalysis in oncogenic G12D, offering a promising avenue to overcome aberrant Ras-driven signal transduction by activating enzymatic activity instead of inhibition. The outlined computational screening protocol is readily applicable for designing ligands and cofactors analogously.

NMR-Chemical-Shift-Driven Protocol Reveals the Cofactor-Bound, Complete Structure of Dynamic Intermediates of the Catalytic Cycle of Oncogenic KRAS G12C Protein and the Significance of the Mg2+ Ion

In this work, catalytically significant states of the oncogenic G12C variant of KRAS, those of Mg2+-free and Mg2+-bound GDP-loaded forms, have been determined using CS-Rosetta software and NMR-data-driven molecular dynamics simulations. There are several Mg2+-bound G12C KRAS/GDP structures deposited in the Protein Data Bank (PDB), so this system was used as a reference, while the structure of the Mg2+-free but GDP-bound state of the RAS cycle has not been determined previously. Due to the high flexibility of the Switch-I and Switch-II regions, which also happen to be the catalytically most significant segments, only chemical shift information could be collected for the most important regions of both systems. CS-Rosetta was used to derive an "NMR ensemble" based on the measured chemical shifts, which, however, did not contain the nonprotein components of the complex. We developed a torsional restraint set for backbone torsions based on the CS-Rosetta ensembles for MD simulations, overriding the force-field-based parametrization in the presence of the reinserted cofactors. This protocol (csdMD) resulted in complete models for both systems that also retained the structural features and heterogeneity defined by the measured chemical shifts and allowed a detailed comparison of the Mg2+-bound and Mg2+-free states of G12C KRAS/GDP.

Quantitative Paramagnetic NMR-Based Analysis of Protein Orientational Dynamics on Membranes: Dissecting the KRas4B-Membrane Interactions

Peripheral membrane proteins can adopt distinct orientations on the surfaces of lipid bilayers that are often short-lived and challenging to characterize by conventional experimental methods. Here we describe a robust approach for mapping protein orientational landscapes through quantitative interpretation of paramagnetic relaxation enhancement (PRE) data arising from membrane mimetics with spin-labeled lipids. Theoretical analysis, followed by experimental verification, reveals insights into the distinct properties of the PRE observables that are generally distorted in the case of stably membrane-anchored proteins. To suppress the artifacts, we demonstrate that undistorted Γ₂ values can be obtained via transient membrane anchoring, based on which a computational framework is established for deriving accurate orientational ensembles obeying Boltzmann statistics. Application of the approach to KRas4B, a classical peripheral membrane protein whose orientations are critical for its functions and drug design, reveals four distinct orientational states that are close but not identical to those reported previously. Similar orientations are also found for a truncated KRas4B without the hypervariable region (HVR) that can sample a broader range of orientations, suggesting a confinement role of the HVR geometrically prohibiting severe tilting. Comparison of the KRas4B Γ₂ rates measured using nanodiscs containing different types of anionic lipids reveals identical Γ₂ patterns for the G-domain but different ones for the HVR, indicating only the latter is able to selectively interact with anionic lipids.

Identification of Dietary Bioflavonoids as Potential Inhibitors against KRAS G12D Mutant—Novel Insights from Computer-Aided Drug Discovery

The KRAS G12D mutation is very frequent in many cancers, such as pancreatic, colon and lung, and has remained undruggable for the past three decades, due to its smooth surface and lack of suitable pockets. Recent small pieces of evidence suggest that targeting the switch I/II of KRAS G12D mutant could be an efficient strategy. Therefore, in the present study, we targeted the switch I (residues 25–40) and switch II (residues 57–76) regions of KRAS G12D with dietary bioflavonoids in comparison with the reference KRAS SI/II inhibitor BI-2852. Initially, we screened 925 bioflavonoids based on drug-likeness properties, and ADME properties and selected 514 bioflavonoids for further studies. Molecular docking resulted in four lead bioflavonoids, namely 5-Dehydroxyparatocarpin K (L1), Carpachromene (L2), Sanggenone H (L3), and Kuwanol C (L4) with binding affinities of 8.8 Kcal/mol, 8.64 Kcal/mol, 8.62 Kcal/mol, and 8.58 Kcal/mol, respectively, in comparison with BI-2852 (−8.59 Kcal/mol). Further steered-molecular dynamics, molecular-dynamics simulation, toxicity, and in silico cancer-cell-line cytotoxicity predictions significantly support these four lead bioflavonoids as potential inhibitors of KRAS G12D SI/SII inhibitors. We finally conclude that these four bioflavonoids have potential inhibitory activity against the KRAS G12D mutant, and are further to be studied in vitro and in vivo, to evaluate their therapeutic potential and the utility of these compounds against KRAS G12D mutated cancers.

Ran GTPase and Its Importance in Cellular Signaling and Malignant Phenotype

Ran is a member of the Ras superfamily of proteins, which primarily regulates nucleocytoplasmic trafficking and mediates mitosis by regulating spindle formation and nuclear envelope (NE) reassembly. Therefore, Ran is an integral cell fate determinant. It has been demonstrated that aberrant Ran expression in cancer is a result of upstream dysregulation of the expression of various factors, such as osteopontin (OPN), and aberrant activation of various signaling pathways, including the extracellular-regulated kinase/mitogen-activated protein kinase (ERK/MEK) and phosphatidylinositol 3-kinase/Protein kinase B (PI3K/Akt) pathways. In vitro, Ran overexpression has severe effects on the cell phenotype, altering proliferation, adhesion, colony density, and invasion. Therefore, Ran overexpression has been identified in numerous types of cancer and has been shown to correlate with tumor grade and the degree of metastasis present in various cancers. The increased malignancy and invasiveness have been attributed to multiple mechanisms. Increased dependence on Ran for spindle formation and mitosis is a consequence of the upregulation of these pathways and the ensuing overexpression of Ran, which increases cellular dependence on Ran for survival. This increases the sensitivity of cells to changes in Ran concentration, with ablation being associated with aneuploidy, cell cycle arrest, and ultimately, cell death. It has also been demonstrated that Ran dysregulation influences nucleocytoplasmic transport, leading to transcription factor misallocation. Consequently, patients with tumors that overexpress Ran have been shown to have a higher malignancy rate and a shorter survival time compared to their counterparts.

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.

Thermal Shift Assay for Small GTPase Stability Screening: Evaluation and Suitability

Thermal unfolding methods are commonly used as a predictive technique by tracking the protein's physical properties. Inherent protein thermal stability and unfolding profiles of biotherapeutics can help to screen or study potential drugs and to find stabilizing or destabilizing conditions. Differential scanning calorimetry (DSC) is a 'Gold Standard' for thermal stability assays (TSA), but there are also a multitude of other methodologies, such as differential scanning fluorimetry (DSF). The use of an external probe increases the assay throughput, making it more suitable for screening studies, but the current methodologies suffer from relatively low sensitivity. While DSF is an effective tool for screening, interpretation and comparison of the results is often complicated. To overcome these challenges, we compared three thermal stability probes in small GTPase stability studies: SYPRO Orange, 8-anilino-1-naphthalenesulfonic acid (ANS), and the Protein-Probe. We studied mainly KRAS, as a proof of principle to obtain biochemical knowledge through TSA profiles. We showed that the Protein-Probe can work at lower concentration than the other dyes, and its sensitivity enables effective studies with non-covalent and covalent drugs at the nanomolar level. Using examples, we describe the parameters, which must be taken into account when characterizing the effect of drug candidates, of both small molecules and Designed Ankyrin Repeat Proteins.

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.

Identification of New KRAS G12D Inhibitors through Computer-Aided Drug Discovery Methods

Owing to several mutations, the oncogene Kirsten rat sarcoma 2 viral oncogene homolog (KRAS) is activated in the majority of cancers, and targeting it has been pharmacologically challenging. In this study, using an in silico approach comprised of pharmacophore modeling, molecular docking, and molecular dynamics simulations, potential KRAS G12D inhibitors were investigated. A ligand-based common feature pharmacophore model was generated to identify the framework necessary for effective KRAS inhibition. The chemical features in the selected pharmacophore model comprised two hydrogen bond donors, one hydrogen bond acceptor, two aromatic rings and one hydrophobic feature. This model was used for screening in excess of 214,000 compounds from InterBioScreen (IBS) and ZINC databases. Eighteen compounds from the IBS and ten from the ZINC database mapped onto the pharmacophore model and were subjected to molecular docking. Molecular docking results highlighted a higher affinity of four hit compounds towards KRAS G12D in comparison to the reference inhibitor, BI-2852. Sequential molecular dynamics (MD) simulation studies revealed all four hit compounds them possess higher KRAS G12D binding free energy and demonstrate stable polar interaction with key residues. Further, Principal Component Analysis (PCA) analysis of the hit compounds in complex with KRAS G12D also indicated stability. Overall, the research undertaken provides strong support for further in vitro testing of these newly identified KRAS G12D inhibitors, particularly Hit1 and Hit2.

Reciprocally-Coupled Gating: Strange Loops in Bioenergetics, Genetics, and Catalysis

Bioenergetics, genetic coding, and catalysis are all difficult to imagine emerging without pre-existing historical context. That context is often posed as a "Chicken and Egg" problem; its resolution is concisely described by de Grasse Tyson: "The egg was laid by a bird that was not a chicken". The concision and generality of that answer furnish no details-only an appropriate framework from which to examine detailed paradigms that might illuminate paradoxes underlying these three life-defining biomolecular processes. We examine experimental aspects here of five examples that all conform to the same paradigm. In each example, a paradox is resolved by coupling "if, and only if" conditions for reciprocal transitions between levels, such that the consequent of the first test is the antecedent for the second. Each condition thus restricts fluxes through, or "gates" the other. Reciprocally-coupled gating, in which two gated processes constrain one another, is self-referential, hence maps onto the formal structure of "strange loops". That mapping uncovers two different kinds of forces that may help unite the axioms underlying three phenomena that distinguish biology from chemistry. As a physical analog for Gödel's logic, biomolecular strange-loops provide a natural metaphor around which to organize a large body of experimental data, linking biology to information, free energy, and the second law of thermodynamics.

Lithium-7 qNMR as a method to quantify lithium content in brines using benchtop NMR

A novel ⁷Li quantitative NMR (qNMR) method to analyze lithium was developed to determine the lithium content in real brine samples using benchtop NMR instruments. The method was validated, and limits of detection and quantification of 40 and 100 ppm, respectively, were determined. Linearity, precision, and bias were also experimentally determined, and the results are presented herein. The results were compared to those obtained using atomic absorption (AA) spectroscopy, currently one of the few validated methods for the quantification of lithium. The method provides both accurate and precise results, as well as excellent correlation with AA. The absence of matrix effects, combined with no need for sample preparation or deuterated solvents, shows potential applicability in the mining industry.

Small molecule inhibitors of RAS proteins with oncogenic mutations

RAS proteins control a number of essential cellular processes as molecular switches in the human body. Presumably due to their important signalling role, RAS proteins are among the most frequently mutated oncogenes in human cancers. Hence, numerous efforts were done to develop appropriate therapies for RAS-mutant cancers in the last three decades. This review aimed to collect all of the reported small molecules that affect RAS signalling. These molecules can be divided in four main branches. First, we address approaches blocking RAS membrane association. Second, we focus on the stabilization efforts of non-productive RAS complexes. Third, we examine the approach to block RAS downstream signalling through disturbance of RAS-effector complex formation. Finally, we discuss direct inhibition; particularly the most recently reported covalent inhibitors, which are already advanced to human clinical trials.

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.

Comparative MD simulations and advanced analytics based studies on wild-type and hot-spot mutant A59G HRas

The Ras family of proteins is known to play an important role in cellular signal transduction. The oncoprotein Ras is also found to be mutated in ~90% of the pancreatic cancers, of which G12V, G13V, A59G and Q61L are the known hot-spot mutants. These ubiquitous proteins fall in the family of G-proteins, and hence switches between active GTP bound and inactive GDP bound states, which is hindered in most of its oncogenic mutant counterparts. Moreover, Ras being a GTPase has an intrinsic property to hydrolyze GTP to GDP, which is obstructed due to mutations and lends the mutants stuck in constitutively active state leading to oncogenic behavior. In this regard, the present study aims to understand the dynamics involved in the hot-spot mutant A59G-Ras using long 10μs classical MD simulations (5μs for each of the wild-type and mutant systems) and comparing the same with its wild-type counterpart. Advanced analytics using Markov State Model (MSM) based approach has been deployed to comparatively understand the transition path for the wild-type and mutant systems. Roles of crucial residues like Tyr32, Gln61 and Tyr64 have also been established using multivariate PCA analyses. Furthermore, this multivariate PCA analysis also provides crucial features which may be used as reaction coordinates for biased simulations for further studies. The absence of formation of pre-hydrolysis network is also reported for the mutant conformation, using the distance-based analyses (between crucial residues) of the conserved regions. The implications of this study strengthen the hypothesis that the disruption of the pre-hydrolysis network in the mutant A59G ensemble might lead to permanently active oncogenic conformation in the mutant conformers.