In Silico Study of the Acquired Resistance Caused by the Secondary Mutations of KRAS G12C Protein Using Long Time Molecular Dynamics Simulation and Markov State Model Analysis

Pathways and mechanism of MRTX1133 binding to KRAS G12D elucidated by molecular dynamics simulations and Markov state models

KRAS G12D is the most common oncogenic mutation identified in several types of cancer. Therefore, design of inhibitors targeting KRAS G12D represents a promising strategy for anticancer therapy. MRTX1133 is a highly potent inhibitor (approximate experiment Kd ≈ 0.0002 nM) of KRAS G12D and is currently in Phase 1/2 study, however, pathways of the compound binding to KRAS G12D has remained unknown, and the mechanism underlying the complicated dynamic process are challenging to capture experimentally, which hinder the structure-based anti-cancer drug design. Here, using MRTX1133 as a probe, unbiased molecular dynamics (MD) was used to simulate the process of MRTX1133 spontaneously binding to KRAS G12D. In six of 42 independent MD simulation (a total of 99 μs), MRTX1133 was observed to successfully associate with KRAS G12D. The kinetically metastable states refer to the potential pathways of MRTX1133 binding to KRAS G12D were revealed by Markov state models (MSM) analysis. Additionally, 8 key residues that are essential for MRTX1133 recognition and tight binding at the preferred low energy states were identified by MM/GBSA analysis. In sum, this study provides a new perspective on understanding the pathways and mechanism of MRTX1133 binding to KRAS G12D.

Circular RNAs in the KRAS pathway: Emerging players in cancer progression

Circular RNAs (circRNAs) have been recognized as key components in the intricate regulatory network of the KRAS pathway across various cancers. The KRAS pathway, a central signalling cascade crucial in tumorigenesis, has gained substantial emphasis as a possible therapeutic target. CircRNAs, a subgroup of non-coding RNAs known for their closed circular arrangement, play diverse roles in gene regulation, contributing to the intricate landscape of cancer biology. This review consolidates existing knowledge on circRNAs within the framework of the KRAS pathway, emphasizing their multifaceted functions in cancer progression. Notable circRNAs, such as Circ_GLG1 and circITGA7, have been identified as pivotal regulators in colorectal cancer (CRC), influencing KRAS expression and the Ras signaling pathway. Aside from their significance in gene regulation, circRNAs contribute to immune evasion, apoptosis, and drug tolerance within KRAS-driven cancers, adding complexity to the intricate interplay. While our comprehension of circRNAs in the KRAS pathway is evolving, challenges such as the diverse landscape of KRAS mutant tumors and the necessity for synergistic combination therapies persist. Integrating cutting-edge technologies, including deep learning-based prediction methods, holds the potential for unveiling disease-associated circRNAs and identifying novel therapeutic targets. Sustained research efforts are crucial to comprehensively unravel the molecular mechanisms governing the intricate interplay between circRNAs and the KRAS pathway, offering insights that could potentially revolutionize cancer diagnostics and treatment strategies.

Assessing the mechanism of fast-cycling cancer-associated mutations of Rac1 small Rho GTPase

Rho-GTPases proteins function as molecular switches alternating from an active to an inactive state upon Guanosine triphosphate (GTP) binding and hydrolysis to Guanosine diphosphate (GDP). Among them, Rac subfamily regulates cell dynamics, being overexpressed in distinct cancer types. Notably, these proteins are object of frequent cancer-associated mutations at Pro29 (P29S, P29L, and P29Q). To assess the impact of these mutations on Rac1 structure and function, we performed extensive all-atom molecular dynamics simulations on wild-type (wt) and oncogenic isoforms of this protein in GDP- and GTP-bound states. Our results unprecedentedly elucidate that P29Q/S-induced structural and dynamical perturbations of Rac1 core domain weaken the binding of the catalytic site Mg2+ ion, and reduce the GDP residence time within protein, enhancing the GDP/GTP exchange rate and Rac1 activity. This broadens our knowledge of the role of cancer-associated mutations on small GTPases mechanism supplying valuable information for future drug discovery efforts targeting specific Rac1 isoforms.

A novel intranasal peptide vaccine inhibits non-small cell lung cancer with KRAS mutation

KRAS mutations occur commonly in the lung and can lead to the development of non-small cell lung cancer (NSCLC). While the mutated KRAS protein is a neoantigen, it usually does not generate an effective anti-tumor immune response on mucosal/epithelial surfaces. Despite this, mutated KRAS remains a potential target for immunotherapy since immune targeting of this protein in animal models has been effective at eliminating tumor cells. We attempted to develop a KRAS vaccine using mutated and wild-type KRAS peptides in combination with a nanoemulsion (NE) adjuvant. The efficacy of this approach was tested in an inducible mutant KRAS-mouse lung tumor model. Animals were immunized intranasally using NE with KRAS peptides. These animals had decreased CD4+FoxP3+ T cells in both lymph nodes and spleen. Immunized animals also showed higher IFN-γ and IL-17a levels to mutated KRAS that were produced by CD8+ T cells and enhancement in KRAS-specific Th1 and Th17 responses that persisted for 3 months after the last vaccination. Importantly, the immunized animals had significantly decreased tumor incidence compared to control animals. In conclusion, a mucosal approach to KRAS vaccination demonstrated the ability to induce local KRAS-specific immune responses in the lung and resulted in reduced tumor incidence.

In silico design of a lipid-like compound targeting KRAS4B-G12D through non-covalent bonds

One of the most common drivers in human cancer is the peripheral membrane protein KRAS4B, able to promote oncogenic signalling. To signal, oncogenic KRAS4B not only requires a sufficient nucleotide exchange, but also needs to recruit effectors by exposing its effector-binding sites while anchoring to the phospholipid bilayer where KRAS4B-mediated signalling events occur. The enzyme phosphodiesterase-δ plays an important role in sequestering KRAS4B from the cytoplasm and targeting it to cellular membranes of different cell species. In this work, we present an in silico design of a lipid-like compound that has the remarkable feature of being able to target both an oncogenic KRAS4B-G12D mutant and the phosphodiesterase-δ enzyme. This double action is accomplished by adding a lipid tail (analogous to the farnesyl group of the KRAS4B protein) to an previously known active compound (2H-1,2,4-benzothiadiazine, 3,4-dihydro-,1,1-dioxide). The proposed lipid-like molecule was found to lock KRAS4B-G12D in its GDP-bound state by adjusting the effector-binding domain to be blocked by the interface of the lipid bilayer. Meanwhile, it can tune GTP-bound KRAS4B-G12D to shift from the active orientation state to the inactive state. The proposed compound is also observed to stably accommodate itself in the prenyl-binding pocket of phosphodiesterase-δ, which impairs KRAS4B enrichment at the lipid bilayer, potentially reducing the proliferation of KRAS4B inside the cytoplasm and its anchoring at the bilayer. In conclusion, we report a potential inhibitor of KRAS4B-G12D with a lipid tail attached to a specific warhead, a compound which has not yet been considered for drugs targeting RAS mutants. Our work provides new ways to target KRAS4B-G12D and can also foster drug discovery efforts for the targeting of oncogenes of the RAS family and beyond.

A Nexus between Genetic and Non-Genetic Mechanisms Guides KRAS Inhibitor Resistance in Lung Cancer

Several studies in the last few years have determined that, in contrast to the prevailing dogma that drug resistance is simply due to Darwinian evolution-the selection of mutant clones in response to drug treatment-non-genetic changes can also lead to drug resistance whereby tolerant, reversible phenotypes are eventually relinquished by resistant, irreversible phenotypes. Here, using KRAS as a paradigm, we illustrate how this nexus between genetic and non-genetic mechanisms enables cancer cells to evade the harmful effects of drug treatment. We discuss how the conformational dynamics of the KRAS molecule, that includes intrinsically disordered regions, is influenced by the binding of the targeted therapies contributing to conformational noise and how this noise impacts the interaction of KRAS with partner proteins to rewire the protein interaction network. Thus, in response to drug treatment, reversible drug-tolerant phenotypes emerge via non-genetic mechanisms that eventually enable the emergence of irreversible resistant clones via genetic mutations. Furthermore, we also discuss the recent data demonstrating how combination therapy can help alleviate KRAS drug resistance in lung cancer, and how new treatment strategies based on evolutionary principles may help minimize or even preclude the emergence of drug resistance.

Classification of GTP-dependent K-Ras4B active and inactive conformational states

Classifying reliably active and inactive molecular conformations of wildtype (WT) and mutated oncogenic proteins is a key, ongoing challenge in molecular cancer studies. Here, we probe the GTP-bound K-Ras4B conformational dynamics using long-time atomistic molecular dynamics (MD) simulations. We extract and analyze the detailed underlying free energy landscape of WT K-Ras4B. We use two key reaction coordinates, labeled d₁ and d₂ (i.e., distances coordinating the P_β atom of the GTP ligand with two key residues, T35 and G60), shown to correlate closely with activities of WT and mutated K-Ras4B. However, our new K-Ras4B conformational kinetics study reveals a more complex network of equilibrium Markovian states. We show that a new reaction coordinate is required to account for the orientation of acidic K-Ras4B sidechains such as D38 with respect to the interface with binding effector RAF1 and rationalize the activation/inactivation propensities and the corresponding molecular binding mechanisms. We use this understanding to unveil how a relatively conservative mutation (i.e., D33E, in the switch I region) can lead to significantly different activation propensities compared with WT K-Ras4B. Our study sheds new light on the ability of residues near the K-Ras4B-RAF1 interface to modulate the network of salt bridges at the binding interface with the RAF1 downstream effector and, thus, to influence the underlying GTP-dependent activation/inactivation mechanism. Altogether, our hybrid MD-docking modeling approach enables the development of new in silico methods for quantitative assessment of activation propensity changes (e.g., due to mutations or local binding environment). It also unveils the underlying molecular mechanisms and facilitates the rational design of new cancer drugs.