Mapping the nucleotide and isoform-dependent structural and dynamical features of Ras proteins

Progress in RAS-targeted therapeutic strategies: From small molecule inhibitors to proteolysis targeting chimeras

As a widely considerable target in chemical biology and pharmacological research, rat sarcoma (RAS) gene mutations play a critical driving factor in several fatal cancers. Despite the great progress of RAS subtype-specific inhibitors, rapid acquired drug resistance could limit their further clinical applications. Proteolysis targeting chimera (PROTAC) has emerged as a powerful tool to handle "undruggable" targets and exhibited significant therapeutic benefit for the combat of drug resistance. Owing to unique molecular mechanism and binding kinetics, PROTAC is expected to become a feasible strategy to break the bottleneck of classical RAS inhibitors. This review aims to discuss the current advances of RAS inhibitors and especially focus on PROTAC strategy targeting RAS mutations and their downstream effectors for relevant cancer treatment.

Unveiling the State Transition Mechanisms of Ras Proteins through Enhanced Sampling and QM/MM Simulations

In cells, wild-type RasGTP complexes exist in two distinct states: active State 2 and inactive State 1. These complexes regulate their functions by transitioning between the two states. However, the mechanisms underlying this state transition have not been clearly elucidated. To address this, we conducted a detailed simulation study to characterize the energetics of the stable states involved in the state transitions of the HRasGTP complex, specifically from State 2 to State 1. This was achieved by employing multiscale quantum mechanics/molecular mechanics and enhanced sampling molecular dynamics methods. Based on the simulation results, we constructed the two-dimensional free energy landscapes that provide crucial information about the conformational changes of the HRasGTP complex from State 2 to State 1. Furthermore, we also explored the conformational changes from the intermediate state to the product state during guanosine triphosphate hydrolysis. This study on the conformational changes involved in the HRas state transitions serves as a valuable reference for understanding the corresponding events of both KRas and NRas as well.

Conformational States of the GDP- and GTP-Bound HRAS Affected by A59E and K117R: An Exploration from Gaussian Accelerated Molecular Dynamics

The HRAS protein is considered a critical target for drug development in cancers. It is vital for effective drug development to understand the effects of mutations on the binding of GTP and GDP to HRAS. We conducted Gaussian accelerated molecular dynamics (GaMD) simulations and free energy landscape (FEL) calculations to investigate the impacts of two mutations (A59E and K117R) on GTP and GDP binding and the conformational states of the switch domain. Our findings demonstrate that these mutations not only modify the flexibility of the switch domains, but also affect the correlated motions of these domains. Furthermore, the mutations significantly disrupt the dynamic behavior of the switch domains, leading to a conformational change in HRAS. Additionally, these mutations significantly impact the switch domain's interactions, including their hydrogen bonding with ligands and electrostatic interactions with magnesium ions. Since the switch domains are crucial for the binding of HRAS to effectors, any alterations in their interactions or conformational states will undoubtedly disrupt the activity of HRAS. This research provides valuable information for the design of drugs targeting HRAS.

Molecular mechanism of regulation of RhoA GTPase by phosphorylation of RhoGDI

Rho-specific guanine nucleotide dissociation inhibitors (RhoGDIs) play a crucial role in the regulation of Rho family GTPases. They act as negative regulators that prevent the activation of Rho GTPases by forming complexes with the inactive GDP-bound state of GTPase. Release of Rho GTPase from the RhoGDI-bound complex is necessary for Rho GTPase activation. Biochemical studies provide evidence of a "phosphorylation code," where phosphorylation of some specific residues of RhoGDI selectively releases its GTPase partner (RhoA, Rac1, Cdc42, etc.). This work attempts to understand the molecular mechanism behind this specific phosphorylation-induced reduction in binding affinity. Using several microseconds long atomistic molecular dynamics simulations of the wild-type and phosphorylated states of the RhoA-RhoGDI complex, we propose a molecular-interaction-based mechanistic model for the dissociation of the complex. Phosphorylation induces major structural changes, particularly in the positively charged polybasic region (PBR) of RhoA and the negatively charged N-terminal region of RhoGDI that contribute most to the binding affinity. Molecular mechanics Poisson-Boltzmann surface area binding energy calculations show a significant weakening of interaction on phosphorylation at the RhoA-specific site of RhoGDI. In contrast, phosphorylation at a Rac1-specific site does not affect the overall binding affinity significantly, which confirms the presence of a phosphorylation code. RhoA-specific phosphorylation leads to a reduction in the number of contacts between the PBR of RhoA and the N-terminal region of RhoGDI, which manifests a reduction of the binding affinity. Using hydrogen bond occupancy analysis and energetic perturbation network, we propose a mechanistic model for the allosteric response, i.e., long-range signal propagation from the site of phosphorylation to the PBR and buried geranylgeranyl group in the form of rearrangement and rewiring of hydrogen bonds and salt bridges. Our results highlight the crucial role of specific electrostatic interactions in manifestation of the phosphorylation code.

Nucleic Acid-Based Approaches to Tackle KRAS Mutant Cancers

Activating mutations in KRAS are highly relevant to various cancers, driving persistent efforts toward the development of drugs that can effectively inhibit KRAS activity. Previously, KRAS was considered 'undruggable'; however, the recent advances in our understanding of RNA and nucleic acid chemistry and delivery formulations have sparked a paradigm shift in the approach to KRAS inhibition. We are currently witnessing a large wave of next-generation drugs for KRAS mutant cancers-nucleic acid-based therapeutics. In this review, we discuss the current progress in targeting KRAS mutant tumors and outline significant developments in nucleic acid-based strategies. We delve into their mechanisms of action, address existing challenges, and offer insights into the current clinical trial status of these approaches. We aim to provide a thorough understanding of the potential of nucleic acid-based strategies in the field of KRAS mutant cancer therapeutics.

The current landscape of using direct inhibitors to target KRASG12C-mutated NSCLC

Mutation in KRAS protooncogene represents one of the most common genetic alterations in NSCLC and has posed a great therapeutic challenge over the past ~ 40 years since its discovery. However, the pioneer work from Shokat's lab in 2013 has led to a recent wave of direct KRASG12C inhibitors that utilize the switch II pocket identified. Notably, two of the inhibitors have recently received US FDA approval for their use in the treatment of KRASG12C mutant NSCLC. Despite this success, there remains the challenge of combating the resistance that cell lines, xenografts, and patients have exhibited while treated with KRASG12C inhibitors. This review discusses the varying mechanisms of resistance that limit long-lasting effective treatment of those direct inhibitors and highlights several novel therapeutic approaches including a new class of KRASG12C (ON) inhibitors, combinational therapies across the same and different pathways, and combination with immunotherapy/chemotherapy as possible solutions to the pressing question of adaptive resistance.

Discovery of Hit Compounds Targeting the P4 Allosteric Site of K-RAS, Identified through Ensemble-Based Virtual Screening

Mutants of Ras are oncogenic drivers of a large number of human tumors. Despite being recognized as an attractive target for the treatment of cancer, the high affinity for its substrate tagged the protein as undruggable for a few years. The identification of cryptic pockets on the protein surface gave the opportunity to identify molecules capable of acting as allosteric modulators. Several molecules were disclosed in recent years, with sotorasib and adagrasib already approved for clinical use. The present study makes use of computational methods to characterize eight prospective allosteric pockets (P1-P8) in K-Ras, four of which (P1-P4) were previously characterized in the literature. The present study also describes the results of a virtual screening study focused on the discovery of hit compounds, binders of the P4 site that can be considered as peptidomimetics of a fragment of the SOS αI helix, a guanine exchange factor of Ras. After a detailed description of the computational procedure followed, we disclose five hit compounds, prospective binders of the P4 allosteric site that exhibit an inhibitory capability higher than 30% in a cell proliferation assay at 50 μM.

Wavelet coherence phase analysis decodes the universal switching mechanism of Ras GTPase superfamily

The Ras superfamily of GTPases regulate critical cellular processes by shuttling between GTP-bound ON and GDP-bound OFF states. This switching mechanism is attributed to the conformational changes in two loops, SWI and SWII, upon GTP binding and hydrolysis. Since these conformational changes vary across the Ras superfamily, there is no generic parameter to define their functional states. A unique wavelet coherence (WC) analysis-based approach developed here shows that the structural changes in switch regions could be mapped onto the wavelet coherence phase couplings (WPCs). Thus, WPCs could serve as unique parameters to define their functional states. Disentanglement of WPCs in oncogenic GTPases shows how breakdown of structural allostery leads to their aberrant function. These observations stand out even for simulated ensemble of switch region conformers. Overall, for the first time, we show that WPCs could unravel the latent structural deviations in Ras proteins to decode their universal switching mechanism.

Impacts of Mutations in the P-Loop on Conformational Alterations of KRAS Investigated with Gaussian Accelerated Molecular Dynamics Simulations

G12 mutations heavily affect conformational transformation and activity of KRAS. In this study, Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the GDP-bound wild-type (WT), G12A, G12D, and G12R KRAS to probe mutation-mediated impacts on conformational alterations of KRAS. The results indicate that three G12 mutations obviously affect the structural flexibility and internal dynamics of the switch domains. The analyses of the free energy landscapes (FELs) suggest that three G12 mutations induce more conformational states of KRAS and lead to more disordered switch domains. The principal component analysis shows that three G12 mutations change concerted motions and dynamics behavior of the switch domains. The switch domains mostly overlap with the binding region of KRAS to its effectors. Thus, the high disorder states and concerted motion changes of the switch domains induced by G12 mutations affect the activity of KRAS. The analysis of interaction network of GDP with KRAS signifies that the instability in the interactions of GDP and magnesium ion with the switch domain SW1 drives the high disordered state of the switch domains. This work is expected to provide theoretical aids for understanding the function of KRAS.

A complete allosteric map of a GTPase switch in its native cellular network

Allosteric regulation is central to protein function in cellular networks. A fundamental open question is whether cellular regulation of allosteric proteins occurs only at a few defined positions or at many sites distributed throughout the structure. Here, we probe the regulation of GTPases-protein switches that control signaling through regulated conformational cycling-at residue-level resolution by deep mutagenesis in the native biological network. For the GTPase Gsp1/Ran, we find that 28% of the 4,315 assayed mutations show pronounced gain-of-function responses. Twenty of the sixty positions enriched for gain-of-function mutations are outside the canonical GTPase active site switch regions. Kinetic analysis shows that these distal sites are allosterically coupled to the active site. We conclude that the GTPase switch mechanism is broadly sensitive to cellular allosteric regulation. Our systematic discovery of new regulatory sites provides a functional map to interrogate and target GTPases controlling many essential biological processes.

Binding modes of GDP, GTP and GNP to NRAS deciphered by using Gaussian accelerated molecular dynamics simulations

Probing binding modes of GDP, GTP and GNP to NRAS are of significance for understanding the regulation mechanism on the activity of RAS proteins. Four separate Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the apo, GDP-, GTP- and GNP-bound NRAS. Dynamics analyses suggest that binding of three ligands highly affects conformational states of the switch domains from NRAS, which disturbs binding of NRAS to its effectors. The analyses of free energy landscapes (FELs) indicate that binding of GDP, GTP and GNP induces more energetic states of NRAS compared to the apo NRAS but the presence of GNP makes the switch domains more ordered than binding of GDP and GNP. The information of interaction networks of ligands with NRAS reveals that the π-π interaction of residue F28 and the salt bridge interactions of K16 and D119 with ligands stabilize binding of GDP, GTP and GNP to NRAS. Meanwhile magnesium ion plays a bridge role in interactions of ligands with NRAS, which is favourable for associations of GDP, GTP and GNP with NRAS. This work is expected to provide useful information for deeply understanding the function and activity of NRAS.

KRAS Inhibitor that Simultaneously Inhibits Nucleotide Exchange Activity and Effector Engagement

We describe a small molecule ligand ACA-14 (2-hydroxy-5-{[(2-phenylcyclopropyl) carbonyl] amino} benzoic acid) as an initial lead for the development of direct inhibitors of KRAS, a notoriously difficult anticancer drug target. We show that the compound binds to KRAS near the switch regions with affinities in the low micromolar range and exerts different effects on KRAS interactions with binding partners. Specifically, ACA-14 impedes the interaction of KRAS with its effector Raf and reduces both intrinsic and SOS-mediated nucleotide exchange rates. Likely as a result of these effects, ACA-14 inhibits signal transduction through the MAPK pathway in cells expressing mutant KRAS and inhibits the growth of pancreatic and colon cancer cells harboring mutant KRAS. We thus propose compound ACA-14 as a useful initial lead for the development of broad-acting inhibitors that target multiple KRAS mutants and simultaneously deplete the fraction of GTP-loaded KRAS while abrogating the effector-binding ability of the already GTP-loaded fraction.

Synthesis, Biological Evaluation, and Binding Mode of a New Class of Oncogenic K-Ras4b Inhibitors

Ras proteins are implicated in some of the most common life-threatening cancers. Despite intense research during the past three decades, progress towards small-molecule inhibitors of mutant Ras proteins still has been limited. Only recently has significant progress been made, in particular with ligands for binding sites located in the switch II and between the switch I and switch II region of K-Ras4B. However, the structural diversity of inhibitors identified for those sites to date is narrow. Herein, we show that hydrazones and oxime ethers of specific bis(het)aryl ketones represent structurally variable chemotypes for new GDP/GTP-exchange inhibitors with significant cellular activity.

Mechanisms of isoform-specific residue influence on GTP-bound HRas, KRas, and NRas

HRas, KRas, and NRas are GTPases with a common set of effectors that control many cell-signaling pathways, including proliferation through Raf kinase. Their G-domains are nearly identical in sequence, with a few isoform-specific residues that have an effect on dynamics and biochemical properties. Here, we use accelerated molecular dynamics (aMD) simulations consistent with solution x-ray scattering experiments to elucidate mechanisms through which isoform-specific residues associated with each Ras isoform affects functionally important regions connected to the active site. HRas-specific residues cluster in loop 8 to stabilize the nucleotide-binding pocket, while NRas-specific residues on helix 3 directly affect the conformations of switch I and switch II. KRas, the most globally flexible of the isoforms, shows greatest fluctuations in the switch regions enhanced by a KRas-specific residue in loop 7 and a highly dynamic loop 8 region. The analysis of isoform-specific residue effects on Ras proteins is supported by NMR experiments and is consistent with previously published biochemical data.

Deciphering Conformational Changes of the GDP-Bound NRAS Induced by Mutations G13D, Q61R, and C118S through Gaussian Accelerated Molecular Dynamic Simulations

The conformational changes in switch domains significantly affect the activity of NRAS. Gaussian-accelerated molecular dynamics (GaMD) simulations of three separate replicas were performed to decipher the effects of G13D, Q16R, and C118S on the conformational transformation of the GDP-bound NRAS. The analyses of root-mean-square fluctuations and dynamics cross-correlation maps indicated that the structural flexibility and motion modes of the switch domains involved in the binding of NRAS to effectors are highly altered by the G13D, Q61R, and C118Smutations. The free energy landscapes (FELs) suggested that mutations induce more energetic states in NRAS than the GDP-bound WT NRAS and lead to high disorder in the switch domains. The FELs also indicated that the different numbers of sodium ions entering the GDP binding regions compensate for the changes in electrostatic environments caused by mutations, especially for G13D. The GDP–residue interactions revealed that the disorder in the switch domains was attributable to the unstable hydrogen bonds between GDP and two residues, V29 and D30. This work is expected to provide information on the energetic basis and dynamics of conformational changes in switch domains that can aid in deeply understanding the target roles of NRAS in anticancer treatment.

Structural insights into Ras regulation by SIN1

Both the mTORC2 and Ras-ERK pathways respond to growth factor stimulation and play critical roles in cell growth and proliferation, disarray of these pathways leads to many diseases, especially cancer. These two signaling pathways crosstalk at many levels; recently it's become clear that the SIN1 component of mTORC2 could interact with Ras family small GTPases, but how these two proteins interact at the molecular level and the functional outcomes of this interaction remain to be addressed. In this work we determined the high-resolution structure of Ras-SIN1 complexes and revealed the detailed interaction mechanism. We also showed that Ras-SIN1 association inhibits insulin-induced ERK activation. Insights from this work could improve our understanding of the disease-causing mechanism of errant mTORC2 or Ras proteins.

Over the years it has been established that SIN1, a key component of mTORC2, could interact with Ras family small GTPases through its Ras-binding domain (RBD). The physical association of Ras and SIN1/mTORC2 could potentially affect both mTORC2 and Ras-ERK pathways. To decipher the precise molecular mechanism of this interaction, we determined the high-resolution structures of HRas/KRas-SIN1 RBD complexes, showing the detailed interaction interface. Mutation of critical interface residues abolished Ras-SIN1 interaction and in SIN1 knockout cells we demonstrated that Ras-SIN1 association promotes SGK1 activity but inhibits insulin-induced ERK activation. With structural comparison and competition fluorescence resonance energy transfer (FRET) assays we showed that HRas-SIN1 RBD association is much weaker than HRas-Raf1 RBD but is slightly stronger than HRas-PI3K RBD interaction, providing a possible explanation for the different outcome of insulin or EGF stimulation. We also found that SIN1 isoform lacking the PH domain binds stronger to Ras than other longer isoforms and the PH domain appears to have an inhibitory effect on Ras-SIN1 binding. In addition, we uncovered a Ras dimerization interface that could be critical for Ras oligomerization. Our results advance our understanding of Ras-SIN1 association and crosstalk between growth factor-stimulated pathways.

Delineating The RAS Conformational Landscape

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

Classification of KRAS-Activating Mutations and the Implications for Therapeutic Intervention

Members of the family of RAS proto-oncogenes, discovered just over 40 years ago, were among the first cancer-initiating genes to be discovered. Of the three RAS family members, KRAS is the most frequently mutated in human cancers. Despite intensive biological and biochemical study of RAS proteins over the past four decades, we are only now starting to devise therapeutic strategies to target their oncogenic properties. Here, we highlight the distinct biochemical properties of common and rare KRAS alleles, enabling their classification into functional subtypes. We also discuss the implications of this functional classification for potential therapeutic avenues targeting mutant subtypes.

SIGNIFICANCE

Efforts in the recent past to inhibit KRAS oncogenicity have focused on kinases that function in downstream signal transduction cascades, although preclinical successes have not translated to patients with KRAS-mutant cancer. Recently, clinically effective covalent inhibitors of KRASG12C have been developed, establishing two principles that form a foundation for future efforts. First, KRAS is druggable. Second, each mutant form of KRAS is likely to have properties that make it uniquely druggable.

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.

Conformations and binding pockets of HRas and its guanine nucleotide exchange factors complexes in the guanosine triphosphate exchange process

The human Son of Sevenless (SOS) activates the signal-transduction protein Ras by forming the complex SOS·Ras and accelerating the guanosine triphosphate (GTP) exchange in Ras. Inhibition of SOS·Ras could regulate the function of Ras in cells and has emerged as an effective strategy for battling Ras related cancers. A key factor to the success of this approach is to understand the conformational change of Ras during the GTP exchange process. In this study, we perform an extensive molecular dynamics simulation to characterize the specific conformations of Ras without and with guanine nucleotide exchange factors (GEFs) of SOS, especially for the substates of State 1 of HRasGTP∙Mg²⁺ . The potent binding pockets on the surfaces of the RasGDP∙Mg²⁺ , the S1.1 and S1.2 substates in State 1 of RasGTP∙Mg²⁺ and the ternary complexes with SOS are predicted, including the binding sites of other domains of SOS. These findings help to obtain a more thorough understanding of Ras functions in the GTP cycling process and provide a structural foundation for future drug design.

From Distinct to Differential Conformational Dynamics to Map Allosteric Communication Pathways in Proteins

Initiation of biological processes involving protein-ligand binding, transient protein-protein interactions, or amino acid modifications alters the conformational dynamics of proteins. Accompanying these biological processes are ensuing coupled atomic level conformational changes within the proteins. These conformational changes collectively connect multiple amino acid residues at distal allosteric, binding, and/or active sites. Local changes due to, for example, binding of a regulatory ligand at an allosteric site initiate the allosteric regulation. The allosteric signal propagates throughout the protein structure, causing changes at distal sites, activating, deactivating, or modifying the function of the protein. Hence, dynamical responses within protein structures to stimuli contain critical information on protein function. In this Perspective, we examine the description of allosteric regulation from protein dynamical responses and associated alternative and emerging computational approaches to map allosteric communication pathways between distal sites in proteins at the atomic level.

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.

Q61 mutant-mediated dynamics changes of the GTP-KRAS complex probed by Gaussian accelerated molecular dynamics and free energy landscapes

Understanding the molecular mechanism of the GTP-KRAS binding is significant for improving the target roles of KRAS in cancer treatment. In this work, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations were applied to decode the effect of Q61A, Q61H and Q61L on the activity of KRAS. Dynamics analyses based on MR-GaMD trajectory reveal that motion modes and dynamics behavior of the switch domain in KRAS are heavily affected by the three Q61 mutants. Information of free energy landscapes (FELs) shows that Q61A, Q61H and Q61L induce structural disorder of the switch domain and disturb the activity of KRAS. Analysis of the interaction network uncovers that the decrease in the stability of hydrogen bonding interactions (HBIs) of GTP with residues V29 and D30 induced by Q61A, Q61H and Q61L is responsible for the structural disorder of the switch-I and that in the occupancy of the hydrogen bond between GTP and residue G60 leads to the structural disorder of the switch-II. Thus, the high disorder of the switch domain caused by three current Q61 mutants produces a significant effect on binding of KRAS to its effectors. This work is expected to provide useful information for further understanding function and target roles of KRAS in anti-cancer drug development.

Structural communication between the GTPase Sec4p and its activator Sec2p: Determinants of GEF activity and early deformations to nucleotide release

Ras GTPases are molecular switches that cycle between OFF and ON states depending on the bound nucleotide (i.e. GDP-bound and GTP-bound, respectively).

The Rab GTPase, Sec4p, plays regulatory roles in multiple steps of intracellular vesicle trafficking. Nucleotide release is catalyzed by the Guanine Nucleotide Exchange Factor (GEF) Sec2p.

Here, the integration of structural information with molecular dynamics (MD) simulations addressed a number of questions concerning the intrinsic and stimulated dynamics of Sec2p and Sec4p as well as the chain of structural deformations leading to GEF-assisted activation of the Rab GTPase.

Sec2p holds an intrinsic ability to adopt the conformation found in the crystallographic complexes with Sec4p, thus suggesting that the latter selects and shifts the conformational equilibrium towards a pre-existing bound-like conformation of Sec2p.

The anchoring of Sec4p to a suitable conformation of Sec2p favors the Sec2p-assisted pulling on itself of the α1/switch 1 (SWI) loop and of SWI, which loose any contact with GDP. Those deformations of Sec4p would occur earlier. Formation of the final Sec2p-Sec4p hydrophobic interface, accomplishes later. Disruption of the nucleotide cage would cause firstly loss of interactions with the guanine ring and secondly loss of interactions with the phosphates.

The ease in sampling the energy landscape and adopting a bound-like conformation likely favors the catalyzing ability of GEFs for Ras GTPases.

Targeted Long-Read Sequencing Decodes the Transcriptional Atlas of the Founding RAS Gene Family Members

The complicity of human RAS proteins in cancer is a well-documented fact, both due to the mutational hyperactivation of these GTPases and the overexpression of the genes encoding these proteins. Thus, it can be easily assumed that the study of RAS genes at the transcriptional and post-transcriptional level is of the utmost importance. Although previous research has shed some light on the basic mechanisms by which GTPases are involved in tumorigenesis, limited information is known regarding the transcriptional profile of the genes encoding these proteins. The present study highlights for the first time the wide spectrum of the mRNAs generated by the three most significant RAS genes (KRAS, NRAS and HRAS), providing an in-depth analysis of the splicing events and exon/intron boundaries. The implementation of a versatile, targeted nanopore-sequencing approach led to the identification of 39 novel RAS mRNA transcript variants and to the elucidation of their expression profiles in a broad panel of human cell lines. Although the present work unveiled multiple hidden aspects of the RAS gene family, further study is required to unravel the biological function of all the novel alternative transcript variants, as well as the putative protein isoforms.

Past and Future Strategies to Inhibit Membrane Localization of the KRAS Oncogene

KRAS is one of the most studied oncogenes. It is well known that KRAS undergoes post-translational modifications at its C-terminal end. These modifications are essential for its membrane location and activity. Despite significant efforts made in the past three decades to target the mechanisms involved in its membrane localization, no therapies have been approved and taken into the clinic. However, many studies have recently reintroduced interest in the development of KRAS inhibitors, either by directly targeting KRAS or indirectly through the inhibition of critical steps involved in post-translational KRAS modifications. In this review, we summarize the approaches that have been applied over the years to inhibit the membrane localization of KRAS in cancer and propose a new anti-KRAS strategy that could be used in clinic.

Allosteric regulation of binding specificity of HVEM for CD160 and BTLA ligands upon G89F mutation

Molecular interactions mediated by engagement of the Herpes virus entry mediator (HVEM) with members of TNF and Ig superfamily generate distinct signals in T cell activation pathways that modulate inflammatory and inhibitory responses. HVEM interacts with CD160 and B and T lymphocyte attenuator (BTLA), both members of the immunoglobulin (Ig) superfamily, which share a common binding site that is unique from that of LIGHT, a TNF ligand. BTLA or CD160 engagement with HVEM deliver inhibitory or stimulatory signals to the host immune response in a context dependent fashion, whereas HVEM engagement with LIGHT results in pro-inflammatory responses. We identified a mutation in human HVEM, G89F, which directly interferes with the human LIGHT interaction, but interestingly, also differentially modulates the binding of human BTLA and CD160 via an apparent allosteric mechanism involving recognition surfaces remote from the site of the mutation. Specifically, the G89F mutation enhances binding of CD160, while decreasing that of BTLA to HVEM in cell-based assays. Molecular dynamics simulations for wild-type and G89F mutant HVEM, bound to different sets of ligands, were performed to define the molecular basis of this unexpected allosteric effect. These results were leveraged to design additional human HVEM mutants with altered binding specificities.

Conformational Fluctuations in GTP-Bound K-Ras: A Metadynamics Perspective with Harmonic Linear Discriminant Analysis

Biomacromolecules often undergo significant conformational rearrangements during function. In proteins, these motions typically consist in nontrivial, concerted rearrangement of multiple flexible regions. Mechanistic, thermodynamics, and kinetic predictions can be obtained via molecular dynamics simulations, provided that the simulation time is at least comparable to the relevant time scale of the process of interest. Because of the substantial computational cost, however, plain MD simulations often have difficulty in obtaining sufficient statistics for converged estimates, requiring the use of more-advanced techniques. Central in many enhanced sampling methods is the definition of a small set of relevant degrees of freedom (collective variables) that are able to describe the transitions between different metastable states of the system. The harmonic linear discriminant analysis (HLDA) has been shown to be useful for constructing low-dimensional collective variables in various complex systems. Here, we apply HLDA to study the free-energy landscape of a monomeric protein around its native state. More precisely, we study the K-Ras protein bound to GTP, focusing on two flexible loops and on the region associated with oncogenic mutations. We perform microsecond-long biased simulations on the wild type and on G12C, G12D, G12 V mutants, describe the resulting free-energy landscapes, and compare our predictions with previous experimental and computational studies. The fast interconversion between open and closed macroscopic states and their similar thermodynamic stabilities are observed. The mutation-induced effects include the alternations of the relative stabilities of different conformational states and the introduction of many microscopic metastable states. Together, our results demonstrate the applicability of the HLDA-based protocol for the conformational sampling of multiple flexible regions in folded proteins.

RAS Dimers: The Novice Couple at the RAS-ERK Pathway Ball

Signals conveyed through the RAS-ERK pathway constitute a pivotal regulatory element in cancer-related cellular processes. Recently, RAS dimerization has been proposed as a key step in the relay of RAS signals, critically contributing to RAF activation. RAS clustering at plasma membrane microdomains and endomembranes facilitates RAS dimerization in response to stimulation, promoting RAF dimerization and subsequent activation. Remarkably, inhibiting RAS dimerization forestalls tumorigenesis in cellular and animal models. Thus, the pharmacological disruption of RAS dimers has emerged as an additional target for cancer researchers in the quest for a means to curtail aberrant RAS activity.

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

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

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.

Conformational Features of Ras: Key Hydrogen-Bonding Interactions of Gln61 in the Intermediate State during GTP Hydrolysis

The Ras protein is one of the most important drug targets for battling cancers. To effectively design novel drugs of Ras, we characterize here its conformational ensembles for the hydrolysis intermediate state RasGDP·Pi and the product state RasGDP by extensive replica-exchange molecular dynamics simulations. Several substates for RasGDP·Pi have been identified, while structural analyses have revealed an unrecognized hydrogen-bonding network that stabilizes the hydrolysis intermediate state. More interestingly, Gln61, which is involved in numerous oncogenic mutations, was found to be engaged in this hydrogen-bonding network, adopting a specific conformation that always points to Pi in contrast to that in the RasGTP state. The simulations also reveal that RasGDP has more than one substate, suggesting a conformational selection mechanism for the interaction between Ras and the guanine nucleotide exchange factors (GEFs). These findings offer new opportunities for the drug design of Ras by stabilizing the hydrolysis intermediate or disrupting its interaction with the GEFs.

K-Ras prenylation as a potential anticancer target

KRAS is one of the most commonly mutated oncogene and a negative predictive factor for a number of targeted therapies. Therefore, the development of targeting strategies against mutant KRAS is urgently needed. One potential strategy involves disruption of K-Ras membrane localization, which is necessary for its proper function. In this review, we summarize the current data about the importance of membrane-anchorage of K-Ras and provide a critical evaluation of this targeting paradigm focusing mainly on prenylation inhibition. Additionally, we performed a RAS mutation-specific analysis of prenylation-related drug sensitivity data from a publicly available database (https://depmap.org/repurposing/) of three classes of prenylation inhibitors: statins, N-bisphosphonates, and farnesyl-transferase inhibitors. We observed significant differences in sensitivity to N-bisphosphonates and farnesyl-transferase inhibitors depending on KRAS mutational status and tissue of origin. These observations emphasize the importance of factors affecting efficacy of prenylation inhibition, like distinct features of different KRAS mutations, tissue-specific mutational patterns, K-Ras turnover, and changes in regulation of prenylation process. Finally, we enlist the factors that might be responsible for the large discrepancy between the outcomes in preclinical and clinical studies including methodological pitfalls, the incomplete understanding of K-Ras protein turnover, and the variation of KRAS dependency in KRAS mutant tumors.

Conformational transformation of switch domains in GDP/K-Ras induced by G13 mutants: An investigation through Gaussian accelerated molecular dynamics simulations and principal component analysis

Mutations in K-Ras are involved in a large number of all human cancers, thus, K-Ras is regarded as a promising target for anticancer drug design. Understanding the target roles of K-Ras is important for providing insights on the molecular mechanism underlying the conformational transformation of the switch domains in K-Ras due to mutations. In this study, multiple replica Gaussian accelerated molecular (MR-GaMD) simulations and principal component analysis (PCA) were applied to probe the effect of G13A, G13D and G13I mutations on conformational transformations of the switch domains in GDP-associated K-Ras. The results suggest that G13A, G13D and G13I enhance the structural flexibility of the switch domains, change the correlated motion modes of the switch domains and strengthen the total motion strength of K-Ras compared with the wild-type (WT) K-Ras. Free energy landscape analyses not only show that the switch domains of the GDP-bound inactive K-Ras mainly exist as a closed state but also indicate that mutations evidently alter the free energy profile of K-Ras and affect the conformational transformation of the switch domains between the closed and open states. Analyses of hydrophobic interaction contacts and hydrogen bonding interactions show that the mutations scarcely change the interaction network of GDP with K-Ras and only disturb the interaction of GDP with the switch (SW1). In summary, two newly introduced mutations, G13A and G13I, play similar adjustment roles in the conformational transformations of two switch domains to G13D and are possibly utilized to tune the activity of K-Ras and the binding of guanine nucleotide exchange factors.

Structural bioinformatics enhances mechanistic interpretation of genomic variation, demonstrated through the analyses of 935 distinct RAS family mutations

MOTIVATION

Protein-coding genetic alterations are frequently observed in Clinical Genetics, but the high yield of variants of uncertain significance remains a limitation in decision making. RAS-family GTPases are cancer drivers, but only 54 variants, across all family members, fall within well-known hotspots. However, extensive sequencing has identified 881 non-hotspot variants for which significance remains to be investigated.

RESULTS

Here, we evaluate 935 missense variants from seven RAS genes, observed in cancer, RASopathies and the healthy adult population. We characterized hotspot variants, previously studied experimentally, using 63 sequence- and 3D structure-based scores, chosen by their breadth of biophysical properties. Applying scores that display best correlation with experimental measures, we report new valuable mechanistic inferences for both hot-spot and non-hotspot variants. Moreover, we demonstrate that 3D scores have little-to-no correlation with those based on DNA sequence, which are commonly used in Clinical Genetics. Thus, combined, these new knowledge bear significant relevance.

AVAILABILITY AND IMPLEMENTATION

All genomic and 3D scores, and markdown for generating figures, are provided in our supplemental data.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Monoubiquitination of KRAS at Lysine104 and Lysine147 Modulates Its Dynamics and Interaction with Partner Proteins

KRAS, a 21 kDa guanine nucleotide-binding protein that functions as a molecular switch, plays a key role in regulating cellular growth. Dysregulation of this key signaling node leads to uncontrolled cell growth, a hallmark of cancer cells. KRAS undergoes post-translational modification by monoubiquitination at various locations, including at lysine104 (K104) and lysine147 (K147). Previous studies have suggested that K104 stabilizes helix-2/helix-3 interactions and K147 is involved in nucleotide binding. However, the impact of monoubiquitination at these residues on the overall structure, dynamics, or function of KRAS is not fully understood. In this study, we examined KRAS monoubiquitination at these sites using data from extensive (12 μs aggregate time) molecular dynamics simulations complemented by nuclear magnetic resonance spectroscopy data. We found that ubiquitin forms dynamic nonspecific interactions with various regions of KRAS and that ubiquitination at both sites modulates conformational fluctuations. In both cases, ubiquitin samples a broad range of conformational space and does not form long-lasting noncovalent contacts with KRAS but it adopts several preferred orientations relative to KRAS. To examine the functional impact of these preferred orientations, we performed a systematic comparison of the dominant configurations of the ubiquitin/KRAS simulated complex with experimental structures of KRAS bound to regulatory and effector proteins as well as a model membrane. Results from these analyses suggest that conformational selection and population shift may minimize the deleterious effects of KRAS ubiquitination at K104 and K147 on binding to some but not all interaction partners. Our findings thus provide new insights into the steric effects of ubiquitin and suggest a potential avenue for therapeutic targeting.

Mutation-Induced Impacts on the Switch Transformations of the GDP- and GTP-Bound K-Ras: Insights from Multiple Replica Gaussian Accelerated Molecular Dynamics and Free Energy Analysis

Mutations yield significant effect on the structural flexibility of two switch domains, SW1 and SW2, in K-Ras, which is considered as an important target of anticancer drug design. To unveil a molecular mechanism with regard to mutation-mediated tuning on the activity of K-Ras, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations followed by analysis of free energy landscapes (FELs) are performed on the GDP- and GTP-bound wild-type (WT), G12V, and D33E K-Ras. The results suggest that G12V and D33E not only evidently change the flexibility of SW1 and SW2 but also greatly affect correlated motions of SW1 and SW2 separately relative to the P-loop and SW1, which exerts a certain tuning on the activity of K-Ras. The information stemming from the analyses of FELs reveals that the conformations of SW1 and SW2 are in high disorders in the GDP- and GTP-associated WT and mutated K-Ras, possibly producing significant effect on binding of guanine nucleotide exchange factors or effectors to K-Ras. The interaction networks of GDP and GTP with K-Ras are identified and the results uncover that the instability in hydrogen-bonding interactions of SW1 with GDP and GTP is mostly responsible for conformational disorder of SW1 and SW2 as well as tunes the activity of oncogenic K-Ras.

Approaches to inhibiting oncogenic K-Ras

Activating somatic K-Ras mutations are associated with >15% all human tumors and up to 90% of specific tumor types such as pancreatic cancer. Successfully inhibiting abnormal K-Ras signaling would therefore be a game changer in cancer therapy. However, K-Ras has long been considered an undruggable target for various reasons. This view is now changing by the discovery of allosteric inhibitors that directly target K-Ras and inhibit its functions, and by the identification of new mechanisms to dislodge it from the plasma membrane and thereby abrogate its cellular activities. In this review, we will discuss recent progresses and challenges to inhibiting aberrant K-Ras functions by these two approaches. We will also provide a broad overview of other approaches such as inhibition of K-Ras effectors, and offer a brief perspective on the way forward.

Stabilization of the Hinge Region of Human E-selectin Enhances Binding Affinity to Ligands Under Force

Introduction

E-selectin is a member of the selectin family of cell adhesion molecules expressed on the plasma membrane of inflamed endothelium and facilitates initial leukocyte tethering and subsequent cell rolling during the early stages of the inflammatory response via binding to glycoproteins expressing sialyl Lewis^X and sialyl Lewis^A (sLe^X/A). Existing crystal structures of the extracellular lectin/EGF-like domain of E-selectin complexed with sLe^X have revealed that E-selectin can exist in two conformation states, a low affinity (bent) conformation, and a high affinity (extended) conformation. The differentiating characteristic of the two conformations is the interdomain angle between the lectin and the EGF-like domain.

Methods

Using molecular dynamics (MD) simulations we observed that in the absence of tensile force E-selectin undergoes spontaneous switching between the two conformational states at equilibrium. A single amino acid substitution at residue 2 (serine to tyrosine) on the lectin domain favors the extended conformation.

Results

Steered molecular dynamics (SMD) simulations of E-selectin and PSGL-1 in conjunction with experimental cell adhesion assays show a longer binding lifetime of E-selectin (S2Y) to PSGL-1 compared to wildtype protein.

Conclusions

The findings in this study advance our understanding into how the structural makeup of E-selectin allosterically influences its adhesive dynamics.

Structure of an inactive conformation of GTP-bound RhoA GTPase

By using ³¹P NMR, we present evidence that the Rho family GTPase RhoA, similar to Ras GTPases, exists in an equilibrium of conformations when bound to GTP. High-resolution crystal structures of RhoA bound to the GTP analog GMPPNP and to GDP show that they display a similar overall inactive conformation. In contrast to the previously reported crystal structures of GTP analog-bound forms of two RhoA dominantly active mutants (G14V and Q63L), GMPPNP-bound RhoA assumes an open conformation in the Switch I loop with a previously unseen interaction between the γ-phosphate and Pro36, instead of the canonical Thr37. Molecular dynamics simulations found that the oncogenic RhoA^G14V mutant displays a reduced flexibility in the Switch regions, consistent with a crystal structure of GDP-bound RhoA^G14V. Thus, GDP- and GTP-bound RhoA can present similar inactive conformations, and the molecular dynamics in the Switch regions are likely to have a role in RhoA activation.

Unveiling conformational dynamics changes of H-Ras induced by mutations based on accelerated molecular dynamics

Uncovering molecular basis with regard to the conformational change of two switches I and II in the GppNHp (GNP)-bound H-Ras is highly significant for the understanding of Ras signaling. For this purpose, accelerated molecular dynamics (aMD) simulations and principal component (PC) analysis are integrated to probe the effect of mutations G12V, T35S and Q61K on conformational transformation between two switches of the GNP-bound H-Ras. The RMSF and cross-correlation analyses suggest that three mutations exert a vital effect on the flexibility and internal dynamics of two switches in the GNP-bound H-Ras. The results stemming from PC analysis indicate that two switches in the GNP-bound WT H-Ras tend to form a closed state in most conformations, while those in the GNP-bound mutated H-Ras display transformation between different states. This conclusion is further supported by free energy landscapes constructed by using the distances of residues 12 away from 35 and 35 away from 61 as reaction coordinates and different experimental studies. Interaction scanning is performed on aMD trajectories and the information shows that conformational transformations of two switches I and II induced by mutations extremely affect the GNP-residue interactions. Meanwhile, the scanning results also signify that residues G15, A18, F28, K117, A146 and K147 form stable contacts with GNP, while residues D30, E31, Y32, D33, P34 and E62 in two switches I and II produce unstable contacts with GNP. This study not only reveals dynamic behavior changes of two switches in H-Ras induced by mutations, but also unveils general principles and mechanisms with regard to functional conformational changes of H-Ras.

Natural Products Attenuating Biosynthesis, Processing, and Activity of Ras Oncoproteins: State of the Art and Future Perspectives

RAS genes encode signaling proteins, which, in mammalian cells, act as molecular switches regulating critical cellular processes as proliferation, growth, differentiation, survival, motility, and metabolism in response to specific stimuli. Deregulation of Ras functions has a high impact on human health: gain-of-function point mutations in RAS genes are found in some developmental disorders and thirty percent of all human cancers, including the deadliest. For this reason, the pathogenic Ras variants represent important clinical targets against which to develop novel, effective, and possibly selective pharmacological inhibitors. Natural products represent a virtually unlimited resource of structurally different compounds from which one could draw on for this purpose, given the improvements in isolation and screening of active molecules from complex sources. After a summary of Ras proteins molecular and regulatory features and Ras-dependent pathways relevant for drug development, we point out the most promising inhibitory approaches, the known druggable sites of wild-type and oncogenic Ras mutants, and describe the known natural compounds capable of attenuating Ras signaling. Finally, we highlight critical issues and perspectives for the future selection of potential Ras inhibitors from natural sources.

Mechanisms of Ras Membrane Organization and Signaling: Ras Rocks Again

Ras is the most frequently mutated oncogene and recent drug development efforts have spurred significant new research interest. Here we review progress toward understanding how Ras functions in nanoscale, proteo-lipid signaling complexes on the plasma membrane, called nanoclusters. We discuss how G-domain reorientation is plausibly linked to Ras-nanoclustering and -dimerization. We then look at how these mechanistic features could cooperate in the engagement and activation of RAF by Ras. Moreover, we show how this structural information can be integrated with microscopy data that provide nanoscale resolution in cell biological experiments. Synthesizing the available data, we propose to distinguish between two types of Ras nanoclusters, an active, immobile RAF-dependent type and an inactive/neutral membrane anchor-dependent. We conclude that it is possible that Ras reorientation enables dynamic Ras dimerization while the whole Ras/RAF complex transits into an active state. These transient di/oligomer interfaces of Ras may be amenable to pharmacological intervention. We close by highlighting a number of open questions including whether all effectors form active nanoclusters and whether there is an isoform specific composition of Ras nanocluster.

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

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

The Bio3D packages for structural bioinformatics

Bio3D is a family of R packages for the analysis of biomolecular sequence, structure, and dynamics. Major functionality includes biomolecular database searching and retrieval, sequence and structure conservation analysis, ensemble normal mode analysis, protein structure and correlation network analysis, principal component, and related multivariate analysis methods. Here, we review recent package developments, including a new underlying segregation into separate packages for distinct analysis, and introduce a new method for structure analysis named ensemble difference distance matrix analysis (eDDM). The eDDM approach calculates and compares atomic distance matrices across large sets of homologous atomic structures to help identify the residue wise determinants underlying specific functional processes. An eDDM workflow is detailed along with an example application to a large protein family. As a new member of the Bio3D family, the Bio3D-eddm package supports both experimental and theoretical simulation-generated structures, is integrated with other methods for dissecting sequence-structure-function relationships, and can be used in a highly automated and reproducible manner. Bio3D is distributed as an integrated set of platform independent open source R packages available from: http://thegrantlab.org/bio3d/.

Drugging all RAS isoforms with one pocket

Activating mutations in the three human RAS genes, KRAS, NRAS and HRAS, are among the most common oncogenic drivers in human cancers. Covalent KRAS^G12C inhibitors, which bind to the switch II pocket in the 'off state' of KRAS, represent the first direct KRAS drugs that entered human clinical trials. However, the remaining 85% of non-KRAS^G12C-driven cancers remain undrugged as do NRAS and HRAS and no drugs targeting the 'on state' have been discovered so far. The switch I/II pocket is a second pocket for which the nanomolar inhibitor BI-2852 has been discovered. Here, we elucidate inhibitor binding modes in KRAS, NRAS and HRAS on and off and discuss future strategies to drug all RAS isoforms with this one pocket.