Targeting RAS phosphorylation in cancer therapy: Mechanisms and modulators

Molecular Mechanism of Phosphorylation-Mediated Impacts on the Conformation Dynamics of GTP-Bound KRAS Probed by GaMD Trajectory-Based Deep Learning

The phosphorylation of different sites produces a significant effect on the conformational dynamics of KRAS. Gaussian accelerated molecular dynamics (GaMD) simulations were combined with deep learning (DL) to explore the molecular mechanism of the phosphorylation-mediated effect on conformational dynamics of the GTP-bound KRAS. The DL finds that the switch domains are involved in obvious differences in conformation contacts and suggests that the switch domains play a key role in the function of KRAS. The analyses of free energy landscapes (FELs) reveal that the phosphorylation of pY32, pY64, and pY137 leads to more disordered states of the switch domains than the wild-type (WT) KRAS and induces conformational transformations between the closed and open states. The results from principal component analysis (PCA) indicate that principal motions PC1 and PC2 are responsible for the closed and open states of the phosphorylated KRAS. Interaction networks were analyzed and the results verify that the phosphorylation alters interactions of GTP and magnesium ion Mg2+ with the switch domains. It is concluded that the phosphorylation pY32, pY64, and pY137 tune the activity of KRAS through changing conformational dynamics and interactions of the switch domains. We anticipated that this work could provide theoretical aids for deeply understanding the function of KRAS.

The effect of phosphorylation on the conformational dynamics and allostery of the association of death-associated protein kinase with calmodulin

Protein phosphorylation plays an important role in the signal transduction and is capable of regulation of cell activity. The death-associated protein kinase 1 (DAPK1), as a Ser/Thr kinase, interacts with calmodulin (CaM) to regulate apoptotic and autophagic signaling. Autophosphorylation of DAPK1 at Ser308 located at the autoregulatory domain (ARD) blocks CaM binding and inhibits kinase catalytic activity. However, the mechanism underlying the influence of Ser308 phosphorylation (pS308) on the DAPK1 activity remains unclear. Here, we performed multiple, microsecond length molecular dynamics (MD) simulations, the molecular mechanics generalized Born/surface area (MM-GBSA) binding free energy calculations, principal component analysis, and dynamic cross-correlation analysis to unravel the conformational dynamics and allostery of the DAPK1 - CaM interaction triggered by the pS308 at the ARD. MD simulations showed that pS308 affected the conformational stability of the DAPK1 - CaM complex. Further energetic and structural exploration revealed that pS308 weakened the association of the phosphorylated DAPK1 to CaM, which lowered the susceptibility of DAPK1 to be activated by CaM. This result can provide mechanistic insights into the molecular underpinning through which the DAPK1 kinase activity is modulated by the auto-phosphorylation.Communicated by Ramaswamy H. Sarma.

STK19 is a DNA/RNA-binding protein critical for DNA damage repair and cell proliferation

STK19 was originally identified as a manganese-dependent serine/threonine-specific protein kinase, but its function has been highly debated. Here, the crystal structure of STK19 revealed that it does not contain a kinase domain, but three intimately packed winged helix (WH) domains. The third WH domain mediated homodimerization and double-stranded DNA binding, both being important for its nuclear localization. STK19 participated in the nucleotide excision repair (NER) and mismatch repair (MMR) pathways by recruiting damage repair factors such as RPA2 and PCNA. STK19 also bound double-stranded RNA through the DNA-binding interface and regulated the expression levels of many mRNAs. Furthermore, STK19 knockdown cells exhibited very slow cell proliferation, which cannot be rescued by dimerization or DNA-binding mutants. Therefore, this work concludes that STK19 is highly unlikely to be a kinase but a DNA/RNA-binding protein critical for DNA damage repair (DDR) and cell proliferation. To prevent further confusions, we renamed this protein as TWH19 (Tandem Winged Helix protein formerly known as STK19).

Role of c-Src in Carcinogenesis and Drug Resistance

The aberrant transformation of normal cells into cancer cells, known as carcinogenesis, is a complex process involving numerous genetic and molecular alterations in response to innate and environmental stimuli. The Src family kinases (SFK) are key components of signaling pathways implicated in carcinogenesis, with c-Src and its oncogenic counterpart v-Src often playing a significant role. The discovery of c-Src represents a compelling narrative highlighting groundbreaking discoveries and valuable insights into the molecular mechanisms underlying carcinogenesis. Upon oncogenic activation, c-Src activates multiple downstream signaling pathways, including the PI3K-AKT pathway, the Ras-MAPK pathway, the JAK-STAT3 pathway, and the FAK/Paxillin pathway, which are important for cell proliferation, survival, migration, invasion, metastasis, and drug resistance. In this review, we delve into the discovery of c-Src and v-Src, the structure of c-Src, and the molecular mechanisms that activate c-Src. We also focus on the various signaling pathways that c-Src employs to promote oncogenesis and resistance to chemotherapy drugs as well as molecularly targeted agents.

Elevated Levels of Mislocalised, Constitutive Ras Signalling Can Drive Quiescence by Uncoupling Cell-Cycle Regulation from Metabolic Homeostasis

The small GTPase Ras plays an important role in connecting external and internal signalling cues to cell fate in eukaryotic cells. As such, the loss of RAS regulation, localisation, or expression level can drive changes in cell behaviour and fate. Post-translational modifications and expression levels are crucial to ensure Ras localisation, regulation, function, and cell fate, exemplified by RAS mutations and gene duplications that are common in many cancers. Here, we reveal that excessive production of yeast Ras2, in which the phosphorylation-regulated serine at position 225 is replaced with alanine or glutamate, leads to its mislocalisation and constitutive activation. Rather than inducing cell death, as has been widely reported to be a consequence of constitutive Ras2 signalling in yeast, the overexpression of RAS2S225A or RAS2S225E alleles leads to slow growth, a loss of respiration, reduced stress response, and a state of quiescence. These effects are mediated via cAMP/PKA signalling and transcriptional changes, suggesting that quiescence is promoted by an uncoupling of cell-cycle regulation from metabolic homeostasis. The quiescent cell fate induced by the overexpression of RAS2S225A or RAS2S225E could be rescued by the deletion of CUP9, a suppressor of the dipeptide transporter Ptr2, or the addition of peptone, implying that a loss of metabolic control, or a failure to pass a metabolic checkpoint, is central to this altered cell fate. Our data suggest that the combination of an increased RAS2 copy number and a dominant active mutation that leads to its mislocalisation can result in growth arrest and add weight to the possibility that approaches to retarget RAS signalling could be employed to develop new therapies.

Overcoming phenotypic switching: targeting protein-protein interactions in cancer

Alternative protein-protein interactions (PPIs) arising from mutations or post-translational modifications (PTMs), termed phenotypic switching (PS), are critical for the transmission of alternative pathogenic signals and are particularly significant in cancer. In recent years, PPIs have emerged as promising targets for rational drug design, primarily because their high specificity facilitates targeting of disease-related signaling pathways. However, obstacles exist at the molecular level that arise from the properties of the interaction interfaces and the propensity of small molecule drugs to interact with more than one cleft surface. The difficulty in identifying small molecules that act as activators or inhibitors to counteract the biological effects of mutations raises issues that have not been encountered before. For example, small molecules can bind tightly but may not act as drugs or bind to multiple sites (interaction promiscuity). Another reason is the absence of significant clefts on protein surfaces; if a pocket is present, it may be too small, or its geometry may prevent binding. PS, which arises from oncogenic (alternative) signaling, causes drug resistance and forms the basis for the systemic robustness of tumors. In this review, the properties of PPI interfaces relevant to the design and development of targeting drugs are examined. In addition, the interactions between three tyrosine kinase inhibitors (TKIs) employed as drugs are discussed. Finally, potential novel targets of one of these drugs were identified in silico.

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.

On the linkage of thermodynamics and pathogenicity

This review outlines the effect of disease-causing mutations on proteins' thermodynamics. Two major thermodynamics quantities, which are essential for structural integrity, the folding and binding free energy changes caused by missense mutations, are considered. It is emphasized that disease effects in case of complex diseases may originate from several mutations over several genes, while monogenic diseases are caused by mutation is a single gene. Nevertheless, in both cases it is shown that pathogenic mutations cause larger perturbations of the above-mentioned thermodynamics quantities as compared with the benign mutations. Recent works demonstrating the effect of pathogenic mutations on the above-mentioned thermodynamics quantities, as well as on structural dynamics and allosteric pathways, are reviewed.

Setting sail: Maneuvering SHP2 activity and its effects in cancer

Since the discovery of tyrosine phosphorylation being a critical modulator of cancer signaling, proteins regulating phosphotyrosine levels in cells have fast become targets of therapeutic intervention. The nonreceptor protein tyrosine phosphatase (PTP) coded by the PTPN11 gene "SHP2" integrates phosphotyrosine signaling from growth factor receptors into the RAS/RAF/ERK pathway and is centrally positioned in processes regulating cell development and oncogenic transformation. Dysregulation of SHP2 expression or activity is linked to tumorigenesis and developmental defects. Even as a compelling anti-cancer target, SHP2 was considered "undruggable" for a long time owing to its conserved catalytic PTP domain that evaded drug development. Recently, SHP2 has risen from the "undruggable curse" with the discovery of small molecules that manipulate its intrinsic allostery for effective inhibition. SHP2's unique domain arrangement and conformation(s) allow for a truly novel paradigm of inhibitor development relying on skillful targeting of noncatalytic sites on proteins. In this review we summarize the biological functions, signaling properties, structural attributes, allostery and inhibitors of SHP2.

Mutant Kras and mTOR crosstalk drives hepatocellular carcinoma development via PEG3/STAT3/BEX2 signaling

Background & Aims: Abnormal activation of mTOR through loss of tuberous sclerosis complex (Tsc) frequently occurs in hepatocellular carcinoma (HCC). Mutant Kras could induce aggressive HCCs. Here, we aim to identify the predictive or prognostic biomarkers for HCC patients with Kras mutant and mTOR hyperactivation, and to provide potential therapeutic approaches for this subtype of HCCs. Methods: We generated transgenic mice in which hepatocytic mTOR was hyperactivated through Tsc1 insufficiency with or without oncogenic Kras^G12D. Bioinformatics and gain- or loss-of-function studies were used to illustrate the mechanisms underlying oncogenic pathway alterations. Transcriptional profiling was used to identify biomarker for the subtype of HCC. The therapeutic efficacy of targeting mTOR was tested in a liver orthotropic homogeneous murine model. Results: Oncogenic Kras^G12D facilitated mTOR activation via the Mek/Erk/ROS axis, leading to HCC tumorigenesis and metastasis. Inhibition of Mek/Erk enhanced the anticancer effect of mTOR inhibitor via reduction of mTOR activity. Paternally expressed 3 (PEG3) was responsible for Kras/Erk- and mTOR-driven HCC. Elevated PEG3 protein interacted with STAT3 and promoted its transcriptional activity, resulting in the upregulation of proliferation- and metastasis-related proteins. Targeting mTOR significantly inhibited these actions in vitro and in vivo. Moreover, in clinical samples, PEG3 was identified as a new poor prognostic marker for HCC patients with Kras/Erk and mTOR hyperactivation. Conclusion: These findings reveal the underlying mechanism of hepatocytic Kras/Erk-driven mTOR activation and its downstream targets (PEG3 and STAT3) in HCC, identify PEG3 as a new prognostic biomarker for HCC with Kras/Erk and mTOR hyperactivation, and provide a potential therapeutic strategy for this subset of HCC patients.

Pan-KRAS inhibitors suppress proliferation through feedback regulation in pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is currently one of the most lethal cancers worldwide. Several basic studies have confirmed that Kirsten rat sarcoma virus (KRAS) is a key driver gene for the occurrence of PDAC, and KRAS mutations have also been found in most patients in clinical studies. In this study, two pan-KRAS inhibitors, BI-2852 and BAY-293, were chosen as chemical probes to investigate their antitumor potency in PDAC. Their inhibitory effects on KRAS activation were validated in vitro and their antiproliferative potency in PDAC cell lines were profiled, with half-maximal inhibitory concentration (IC50) values of approximately 1 μM, demonstrating the therapeutic potential of pan-KRAS inhibitors in the treatment of PDAC. However, feedback regulation in the KRAS pathway weakened inhibitor activity, which was observed by a 50 times difference in BAY-293 from in vitro activity. Furthermore, pan-KRAS inhibitors effectively inhibited cell proliferation in 3D organoids cultured from PDAC patient samples; however, there were some variations between individuals. These results provide a sufficient theoretical foundation for KRAS as a clinical therapeutic target and for the application of pan-KRAS inhibitors in the treatment of PDAC, with important scientific significance in translational medicine.

Computer-Aided Drug Design Boosts RAS Inhibitor Discovery

The Rat Sarcoma (RAS) family (NRAS, HRAS, and KRAS) is endowed with GTPase activity to regulate various signaling pathways in ubiquitous animal cells. As proto-oncogenes, RAS mutations can maintain activation, leading to the growth and proliferation of abnormal cells and the development of a variety of human cancers. For the fight against tumors, the discovery of RAS-targeted drugs is of high significance. On the one hand, the structural properties of the RAS protein make it difficult to find inhibitors specifically targeted to it. On the other hand, targeting other molecules in the RAS signaling pathway often leads to severe tissue toxicities due to the lack of disease specificity. However, computer-aided drug design (CADD) can help solve the above problems. As an interdisciplinary approach that combines computational biology with medicinal chemistry, CADD has brought a variety of advances and numerous benefits to drug design, such as the rapid identification of new targets and discovery of new drugs. Based on an overview of RAS features and the history of inhibitor discovery, this review provides insight into the application of mainstream CADD methods to RAS drug design.

Markov State Models and Molecular Dynamics Simulations Reveal the Conformational Transition of the Intrinsically Disordered Hypervariable Region of K-Ras4B to the Ordered Conformation

K-Ras4B, the most frequently mutated Ras isoform in human tumors, plays a vital part in cell growth, differentiation, and survival. Its tail, the C-terminal hypervariable region (HVR), is involved in anchoring K-Ras4B at the cellular plasma membrane and in isoform-specific protein-protein interactions and signaling. In the inactive guanosine diphosphate-bound state, the intrinsically disordered HVR interacts with the catalytic domain at the effector-binding region, rendering K-Ras4B in its autoinhibited state. Activation releases the HVR from the catalytic domain, with its ensemble favoring an ordered α-helical structure. The large-scale conformational transition of the HVR from the intrinsically disordered to the ordered conformation remains poorly understood. Here, we deploy a computational scheme that integrates a transition path-generation algorithm, extensive molecular dynamics simulation, and Markov state model analysis to investigate the conformational landscape of the HVR transition pathway. Our findings reveal a stepwise pathway for the HVR transition and uncover several key conformational substates along the transition pathway. Importantly, key interactions between the HVR and the catalytic domain are unraveled, highlighting the pathogenesis of K-Ras4B mild mutations in several congenital developmental anomaly syndromes. Together, these findings provide a deeper understanding of the HVR transition mechanism and the regulation of K-Ras4B activity at an atomic level.

Mechanistic insights into the clinical Y96D mutation with acquired resistance to AMG510 in the KRASG12C

Special oncogenic mutations in the RAS proteins lead to the aberrant activation of RAS and its downstream signaling pathways. AMG510, the first approval drug for KRAS, covalently binds to the mutated cysteine 12 of KRASG12C protein and has shown promising antitumor activity in clinical trials. Recent studies have reported that the clinically acquired Y96D mutation could severely affect the effectiveness of AMG510. However, the underlying mechanism of the drug-resistance remains unclear. To address this, we performed multiple microsecond molecular dynamics simulations on the KRASG12C−AMG510 and KRASG12C/Y96D−AMG510 complexes at the atomic level. The direct interaction between the residue 96 and AMG510 was impaired owing to the Y96D mutation. Moreover, the mutation yielded higher flexibility and more coupled motion of the switch II and α3-helix, which led to the departing motion of the switch II and α3-helix. The resulting departing motion impaired the interaction between the switch II and α3-helix and subsequently induced the opening and loosening of the AMG510 binding pocket, which further disrupted the interaction between the key residues in the pocket and AMG510 and induced an increased solvent exposure of AMG510. These findings reveal the resistance mechanism of AMG510 to KRASG12C/Y96D, which will help to offer guidance for the development of KRAS targeted drugs to overcome acquired resistance.

Understanding gilteritinib resistance to FLT3-F691L mutation through an integrated computational strategy

FMS-like tyrosine kinase 3 (FLT3) serves as an important drug target for acute myeloid leukemia (AML), and gene mutations of FLT3 have been closely associated with AML patients with an incidence rate of ~ 30%. However, the mechanism of the clinically relevant F691L gatekeeper mutation conferred resistance to the drug gilteritinib remained poorly understood. In this study, multiple microsecond molecular dynamics (MD) simulations, end-point free energy calculations, and dynamic correlated and network analyses were performed to investigate the molecular basis of gilteritinib resistance to the FLT3-F691L mutation. The simulations revealed that the resistant mutation largely induced the conformational changes of the activation loop (A-loop), the phosphate-binding loop, and the helix αC of the FLT3 protein. The binding abilities of the gilteritinib to the wild-type and the F691L mutant were different through the binding free energy prediction. The simulation results further indicated that the driving force to determine the binding affinity of gilteritinib was derived from the differences in the energy terms of electrostatic and van der Waals interactions. Moreover, the per-residue free energy decomposition suggested that the four residues (Phe803, Gly831, Leu832, and Ala833) located at the A-loop of FLT3 had a significant impact on the binding affinity of gilteritinib to the F691L mutant. This study may provide useful information for the design of novel FLT3 inhibitors specially targeting the F691L gatekeeper mutant.

Natural Products in Preventing Tumor Drug Resistance and Related Signaling Pathways

Drug resistance is still an obstacle in cancer therapy, leading to the failure of tumor treatment. The emergence of tumor drug resistance has always been a main concern of oncologists. Therefore, overcoming tumor drug resistance and looking for new strategies for tumor treatment is a major focus in the field of tumor research. Natural products serve as effective substances against drug resistance because of their diverse chemical structures and pharmacological effects. We reviewed the signaling pathways involved in the development of tumor drug resistance, including Epidermal growth factor receptor (EGFR), Renin-angiotensin system (Ras), Phosphatidylinositol-3-kinase/protein kinase B (PI3K/Akt), Wnt, Notch, Transforming growth factor-beta (TGF-β), and their specific signaling pathway inhibitors derived from natural products. This can provide new ideas for the prevention of drug resistance in cancer therapy.

Autopromotion of K-Ras4B Feedback Activation Through an SOS-Mediated Long-Range Allosteric Effect

The Ras-specific guanine nucleotide exchange factors Son of Sevenless (SOS) regulates Ras activation by converting inactive GDP-bound to active GTP-bound states. The catalytic activity of Ras is further allosterically regulated by GTP−Ras bound to a distal site through a positive feedback loop. To address the mechanism underlying the long-range allosteric activation of the catalytic K-Ras4B by an additional allosteric GTP–Ras through SOS, we employed molecular dynamics simulation of the K-Ras4B^G13D•SOS^cat complex with and without an allosteric GTP-bound K-Ras4B^G13D. We found that the binding of an allosteric GTP−K-Ras4B^G13D enhanced the affinity between the catalytic K-Ras4B^G13D and SOS^cat, forming a more stable conformational state. The peeling away of the switch I from the nucleotide binding site facilitated the dissociation of GDP, thereby contributing to the increased nucleotide exchange rate. The community networks further showed stronger edge connection upon allosteric GTP−K-Ras4B^G13D binding, which represented an increased interaction between catalytic K-Ras4B^G13D and SOS^cat. Moreover, GTP−K-Ras4B^G13D binding transmitted allosteric signaling pathways though the Cdc25 domain of SOS that enhanced the allosteric regulatory from the K-Ras4B^G13D allosteric site to the catalytic site. This study may provide an in-depth mechanism for abnormal activation and allosteric regulation of K-Ras4B^G13D.

Mechanistic Insights Into Co-Administration of Allosteric and Orthosteric Drugs to Overcome Drug-Resistance in T315I BCR-ABL1

Chronic myeloid leukemia (CML) is a myeloproliferative neoplasm, driven by the BCR-ABL1 fusion oncoprotein. The discovery of orthosteric BCR-ABL1 tyrosine kinase inhibitors (TKIs) targeting its active ATP-binding pocket, such as first-generation Imatinib and second-generation Nilotinib (NIL), has profoundly revolutionized the therapeutic landscape of CML. However, currently targeted therapeutics still face considerable challenges with the inevitable emergence of drug-resistant mutations within BCR-ABL1. One of the most common resistant mutations in BCR-ABL1 is the T315I gatekeeper mutation, which confers resistance to most current TKIs in use. To resolve such conundrum, co-administration of orthosteric TKIs and allosteric drugs offers a novel paradigm to tackle drug resistance. Remarkably, previous studies have confirmed that the dual targeting BCR-ABL1 utilizing orthosteric TKI NIL and allosteric inhibitor ABL001 resulted in eradication of the CML xenograft tumors, exhibiting promising therapeutic potential. Previous studies have demonstrated the cooperated mechanism of two drugs. However, the conformational landscapes of synergistic effects remain unclear, hampering future efforts in optimizations and improvements. Hence, extensive large-scale molecular dynamics (MD) simulations of wide type (WT), WT-NIL, T315I, T315I-NIL, T315I-ABL001 and T315I-ABL001-NIL systems were carried out in an attempt to address such question. Simulation data revealed that the dynamic landscape of NIL-bound BCR-ABL1 was significantly reshaped upon ABL001 binding, as it shifted from an active conformation towards an inactive conformation. The community network of allosteric signaling was analyzed to elucidate the atomistic overview of allosteric regulation within BCR-ABL1. Moreover, binding free energy analysis unveiled that the affinity of NIL to BCR-ABL1 increased by the induction of ABL001, which led to its favorable binding and the release of drug resistance. The findings uncovered the in-depth structural mechanisms underpinning dual-targeting towards T315I BCR-ABL1 to overcome its drug resistance and will offer guidance for the rational design of next generations of BCR-ABL1 modulators and future combinatory therapeutic regimens.

Harnessing Reversed Allosteric Communication: A Novel Strategy for Allosteric Drug Discovery

Allostery is a fundamental and extensive mechanism of intramolecular signal transmission. Allosteric drugs possess several unique pharmacological advantages over traditional orthosteric drugs, including greater selectivity, better physicochemical properties, and lower off-target toxicity. However, owing to the complexity of allosteric regulation, experimental approaches for the development of allosteric modulators are traditionally serendipitous. Recently, the reversed allosteric communication theory has been proposed, providing a feasible tool for the unbiased detection of allosteric sites. Herein, we review the latest research on the reversed allosteric communication effect using the examples of sirtuin 6, epidermal growth factor receptor, 3-phosphoinositide-dependent protein kinase 1, and Related to A and C kinases (RAC) serine/threonine protein kinase B and recapitulate the methodologies of reversed allosteric communication strategy. The novel reversed allosteric communication strategy greatly expands the horizon of allosteric site identification and allosteric mechanism exploration and is expected to accelerate an end-to-end framework for drug discovery.