Mechanism of activation and the rewired network: New drug design concepts

Unveiling hidden connections in omics data via pyPARAGON: an integrative hybrid approach for disease network construction

Network inference or reconstruction algorithms play an integral role in successfully analyzing and identifying causal relationships between omics hits for detecting dysregulated and altered signaling components in various contexts, encompassing disease states and drug perturbations. However, accurate representation of signaling networks and identification of context-specific interactions within sparse omics datasets in complex interactomes pose significant challenges in integrative approaches. To address these challenges, we present pyPARAGON (PAgeRAnk-flux on Graphlet-guided network for multi-Omic data integratioN), a novel tool that combines network propagation with graphlets. pyPARAGON enhances accuracy and minimizes the inclusion of nonspecific interactions in signaling networks by utilizing network rather than relying on pairwise connections among proteins. Through comprehensive evaluations on benchmark signaling pathways, we demonstrate that pyPARAGON outperforms state-of-the-art approaches in node propagation and edge inference. Furthermore, pyPARAGON exhibits promising performance in discovering cancer driver networks. Notably, we demonstrate its utility in network-based stratification of patient tumors by integrating phosphoproteomic data from 105 breast cancer tumors with the interactome and demonstrating tumor-specific signaling pathways. Overall, pyPARAGON is a novel tool for analyzing and integrating multi-omic data in the context of signaling networks. pyPARAGON is available at https://github.com/netlab-ku/pyPARAGON.

SHP2-EGFR States in Dephosphorylation Can Inform Selective SHP2 Inhibitors, Dampening RasGAP Action

SHP2 is a positive regulator of the EGFR-dependent Ras/MAPK pathway. It dephosphorylates a regulatory phosphorylation site in EGFR that serves as the binding site to RasGAP (RASA1 or p120RasGAP). RASA1 is activated by binding to the EGFR phosphate group. Active RASA1 deactivates Ras by hydrolyzing Ras-bound GTP to GDP. Thus, SHP2 dephosphorylation of EGFR effectively prevents RASA1-mediated deactivation of Ras, thereby stimulating proliferation. Despite knowledge of this vital regulation in cell life, mechanistic in-depth structural understanding of the involvement of SHP2, EGFR, and RASA1 in the Ras/MAPK pathway has largely remained elusive. Here we elucidate the interactions, the factors influencing EGFR's recruitment of RASA1, and SHP2's recognition of the substrate site in EGFR. We reveal that RASA1 specifically interacts with the DEpY992LIP motif in EGFR featuring a proline residue at the +3 position C-terminal to pY primarily through its nSH2 domain. This interaction is strengthened by the robust attraction of two acidic residues, E991 and D990, of EGFR to two basic residues in the BC-loop near the pY-binding pocket of RASA1's nSH2. In the stable precatalytic state of SHP2 with EGFR (DADEpY992LIPQ), the E-loop of SHP2's active site favors the interaction with the (-2)-position D990 and (-4)-position D988 N-terminal to pY992 in EGFR, while the pY-loop constrains the (+4)-position Q996 C-terminal to pY992. These specific interactions not only provide a structural basis for identifying negative regulatory sites in other RTKs but can inform selective, high-affinity active-site SHP2 inhibitors tailored for SHP2 mutants.

CDK2 and CDK4: Cell Cycle Functions Evolve Distinct, Catalysis-Competent Conformations, Offering Drug Targets

Cyclin-dependent kinases (CDKs), particularly CDK4 and CDK2, are crucial for cell cycle progression from the Gap 1 (G1) to the Synthesis (S) phase by phosphorylating targets such as the Retinoblastoma Protein (Rb). CDK4, paired with cyclin-D, operates in the long G1 phase, while CDK2 with cyclin-E, manages the brief G1-to-S transition, enabling DNA replication. Aberrant CDK signaling leads to uncontrolled cell proliferation, which is a hallmark of cancer. Exactly how they accomplish their catalytic phosphorylation actions with distinct efficiencies poses the fundamental, albeit overlooked question. Here we combined available experimental data and modeling of the active complexes to establish their conformational functional landscapes to explain how the two cyclin/CDK complexes differentially populate their catalytically competent states for cell cycle progression. Our premise is that CDK catalytic efficiencies could be more important for cell cycle progression than the cyclin-CDK biochemical binding specificity and that efficiency is likely the prime determinant of cell cycle progression. We observe that CDK4 is more dynamic than CDK2 in the ATP binding site, the regulatory spine, and the interaction with its cyclin partner. The N-terminus of cyclin-D acts as an allosteric regulator of the activation loop and the ATP-binding site in CDK4. Integrated with a suite of experimental data, we suggest that the CDK4 complex is less capable of remaining in the active catalytically competent conformation, and may have a lower catalytic efficiency than CDK2, befitting their cell cycle time scales, and point to critical residues and motifs that drive their differences. Our mechanistic landscape may apply broadly to kinases, and we propose two drug design strategies: (i) allosteric Inhibition by conformational stabilization for targeting allosteric CDK4 regulation by cyclin-D, and (ii) dynamic entropy-optimized targeting which leverages the dynamic, entropic aspects of CDK4 to optimize drug binding efficacy.

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.

Neurodevelopmental disorders, like cancer, are connected to impaired chromatin remodelers, PI3K/mTOR, and PAK1-regulated MAPK

Neurodevelopmental disorders (NDDs) and cancer share proteins, pathways, and mutations. Their clinical symptoms are different. However, individuals with NDDs have higher probabilities of eventually developing cancer. Here, we review the literature and ask how the shared features can lead to different medical conditions and why having an NDD first can increase the chances of malignancy. To explore these vital questions, we focus on dysregulated PI3K/mTOR, a major brain cell growth pathway in differentiation, and MAPK, a critical pathway in proliferation, a hallmark of cancer. Differentiation is governed by chromatin organization, making aberrant chromatin remodelers highly likely agents in NDDs. Dysregulated chromatin organization and accessibility influence the lineage of specific cell brain types at specific embryonic development stages. PAK1, with pivotal roles in brain development and in cancer, also regulates MAPK. We review, clarify, and connect dysregulated pathways with dysregulated proliferation and differentiation in cancer and NDDs and highlight PAK1 role in brain development and MAPK regulation. Exactly how PAK1 activation controls brain development, and why specific chromatin remodeler components, e.g., BAF170 encoded by SMARCC2 in autism, await clarification.

Pan-cancer clinical impact of latent drivers from double mutations

Here, we discover potential 'latent driver' mutations in cancer genomes. Latent drivers have low frequencies and minor observable translational potential. As such, to date they have escaped identification. Their discovery is important, since when paired in cis, latent driver mutations can drive cancer. Our comprehensive statistical analysis of the pan-cancer mutation profiles of ~60,000 tumor sequences from the TCGA and AACR-GENIE cohorts identifies significantly co-occurring potential latent drivers. We observe 155 same gene double mutations of which 140 individual components are cataloged as latent drivers. Evaluation of cell lines and patient-derived xenograft response data to drug treatment indicate that in certain genes double mutations may have a prominent role in increasing oncogenic activity, hence obtaining a better drug response, as in PIK3CA. Taken together, our comprehensive analyses indicate that same-gene double mutations are exceedingly rare phenomena but are a signature for some cancer types, e.g., breast, and lung cancers. The relative rarity of doublets can be explained by the likelihood of strong signals resulting in oncogene-induced senescence, and by doublets consisting of non-identical single residue components populating the background mutational load, thus not identified.

A New View of Activating Mutations in Cancer

A vast effort has been invested in the identification of driver mutations of cancer. However, recent studies and observations call into question whether the activating mutations or the signal strength are the major determinant of tumor development. The data argue that signal strength determines cell fate, not the mutation that initiated it. In addition to activating mutations, factors that can impact signaling strength include (i) homeostatic mechanisms that can block or enhance the signal, (ii) the types and locations of additional mutations, and (iii) the expression levels of specific isoforms of genes and regulators of proteins in the pathway. Because signal levels are largely decided by chromatin structure, they vary across cell types, states, and time windows. A strong activating mutation can be restricted by low expression, whereas a weaker mutation can be strengthened by high expression. Strong signals can be associated with cell proliferation, but too strong a signal may result in oncogene-induced senescence. Beyond cancer, moderate signal strength in embryonic neural cells may be associated with neurodevelopmental disorders, and moderate signals in aging may be associated with neurodegenerative diseases, like Alzheimer's disease. The challenge for improving patient outcomes therefore lies in determining signaling thresholds and predicting signal strength.

Allosteric regulation of autoinhibition and activation of c-Abl

c-Abl, a non-receptor tyrosine kinase, regulates cell growth and survival in healthy cells and causes chronic myeloid leukemia (CML) when fused by Bcr. Its activity is blocked in the assembled inactive state, where the SH3 and SH2 domains dock into the kinase domain, reducing its conformational flexibility, resulting in the autoinhibited state. It is active in an extended 'open' conformation. Allostery governs the transitions between the autoinhibited and active states. Even though experiments revealed the structural hallmarks of the two states, a detailed grasp of the determinants of c-Abl autoinhibition and activation at the atomic level, which may help innovative drug discovery, is still lacking. Here, using extensive molecular dynamics simulations, we decipher exactly how these determinants regulate it. Our simulations confirm and extend experimental data that the myristoyl group serves as the switch for c-Abl inhibition/activation. Its dissociation from the kinase domain promotes the SH2-SH3 release, initiating c-Abl activation. We show that the precise SH2/N-lobe interaction is required for full activation of c-Abl. It stabilizes a catalysis-favored conformation, priming it for catalytic action. Bcr-Abl allosteric drugs elegantly mimic the endogenous myristoyl-mediated autoinhibition state of c-Abl 1b. Allosteric activating mutations shift the ensemble to the active state, blocking ATP-competitive drugs. Allosteric drugs alter the active-site conformation, shifting the ensemble to re-favor ATP-competitive drugs. Our work provides a complete mechanism of c-Abl activation and insights into critical parameters controlling at the atomic level c-Abl inactivation, leading us to propose possible strategies to counter reemergence of drug resistance.

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 the P-Loop Conformation in the Determination of Inhibitor Selectivity Toward the Hepatocellular Carcinoma-Associated Dark Kinase STK17B

As a member of the death-associated protein kinase family of serine/threonine kinases, the STK17B has been associated with diverse diseases such as hepatocellular carcinoma. However, the conformational dynamics of the phosphate-binding loop (P-loop) in the determination of inhibitor selectivity profile to the STK17B are less understood. Here, a multi-microsecond length molecular dynamics (MD) simulation of STK17B in the three different states (ligand-free, ADP-bound, and ligand-bound states) was carried out to uncover the conformational plasticity of the P-loop. Together with the analyses of principal component analysis, cross-correlation and generalized correlation motions, secondary structural analysis, and community network analysis, the conformational dynamics of the P-loop in the different states were revealed, in which the P-loop flipped into the ADP-binding site upon the inhibitor binding and interacted with the inhibitor and the C-lobe, strengthened the communication between the N- and C-lobes. These resulting interactions contributed to inhibitor selectivity profile to the STK17B. Our results may advance our understanding of kinase inhibitor selectivity and offer possible implications for the design of highly selective inhibitors for other protein kinases.

The mechanism of activation of MEK1 by B-Raf and KSR1

MEK1 interactions with B-Raf and KSR1 are key steps in Ras/Raf/MEK/ERK signaling. Despite this, vital mechanistic details of how these execute signal transduction are still enigmatic. Among these is why, despite B-Raf and KSR1 kinase domains similarity, the B-Raf/MEK1 and KSR1/MEK1 complexes have distinct contributions to MEK1 activation, and broadly, what is KSR1's role. Our molecular dynamics simulations clarify these still unresolved ambiguities. Our results reveal that the proline-rich (P-rich) loop of MEK1 plays a decisive role in MEK1 activation loop (A-loop) phosphorylation. In the inactive B-Raf/MEK1 heterodimer, the collapsed A-loop of B-Raf interacts with the P-rich loop and A-loop of MEK1, minimizing MEK1 A-loop fluctuation and preventing it from phosphorylation. In the active B-Raf/MEK1 heterodimer, the P-rich loop moves in concert with the A-loop of B-Raf as it extends. This reduces the number of residues interacting with MEK1 A-loop, allowing increased A-loop fluctuation, and bringing Ser222 closer to ATP for phosphorylation. B-Raf αG-helix Arg662 promotes MEK1 activation by orienting Ser218 towards ATP. In KSR1/MEK1, the KSR1 αG-helix has Ala826 in place of B-Raf Arg662. This difference results in much fewer interactions between KSR1 αG-helix and MEK1 A-loop, thus a more flexible A-loop. We postulate that if KSR1 were to adopt an active configuration with an extended A-loop as seen in other protein kinases, then the MEK1 P-rich loop would extend in a similar manner, as seen in the active B-Raf/MEK1 heterodimer. This would result in highly flexible MEK1 A-loop, and KSR1 functioning as an active, B-Raf-like, kinase.

Neurodevelopmental disorders, immunity, and cancer are connected

Immunity could be viewed as the common factor in neurodevelopmental disorders and cancer. The immune and nervous systems coevolve as the embryo develops. Immunity can release cytokines that activate MAPK signaling in neural cells. In specific embryonic brain cell types, dysregulated signaling that results from germline or embryonic mutations can promote changes in chromatin organization and gene accessibility, and thus expression levels of essential genes in neurodevelopment. In cancer, dysregulated signaling can emerge from sporadic somatic mutations during human life. Neurodevelopmental disorders and cancer share similarities. In neurodevelopmental disorders, immunity, and cancer, there appears an almost invariable involvement of small GTPases (e.g., Ras, RhoA, and Rac) and their pathways. TLRs, IL-1, GIT1, and FGFR signaling pathways, all can be dysregulated in neurodevelopmental disorders and cancer. Although there are signaling similarities, decisive differentiating factors are timing windows, and cell type specific perturbation levels, pointing to chromatin reorganization. Finally, we discuss drug discovery.

Allosteric pluripotency: challenges and opportunities

Allosteric pluripotency arises when the functional response of an allosteric receptor to an allosteric stimulus depends on additional allosteric modulators. Here, we discuss allosteric pluripotency as observed in the prototypical Protein Kinase A (PKA) as well as in other signaling systems, from typical multidomain signaling proteins to bacterial enzymes. We identify key drivers of pluripotent allostery and illustrate how hypothesizing allosteric pluripotency may solve apparent discrepancies currently present in the literature regarding the dual nature of known allosteric modulators. We also outline the implications of allosteric pluripotency for cellular signaling and allosteric drug design, and analyze the challenges and opportunities opened by the pluripotent nature of allostery.

Regulatory spine RS3 residue of protein kinases: a lipophilic bystander or a decisive element in the small-molecule kinase inhibitor binding?

In recent years, protein kinases have been one of the most pursued drug targets. These determined efforts have resulted in ever increasing numbers of small-molecule kinase inhibitors reaching to the market, offering novel treatment options for patients with distinct diseases. One essential component related to the activation and normal functionality of a protein kinase is the regulatory spine (R-spine). The R-spine is formed of four conserved residues named as RS1–RS4. One of these residues, RS3, located in the C-terminal part of αC-helix, is usually accessible for the inhibitors from the ATP-binding cavity as its side chain is lining the hydrophobic back pocket in many protein kinases. Although the role of RS3 has been well acknowledged in protein kinase function, this residue has not been actively considered in inhibitor design, even though many small-molecule kinase inhibitors display interactions to this residue. In this minireview, we will cover the current knowledge of RS3, its relationship with the gatekeeper, and the role of RS3 in kinase inhibitor interactions. Finally, we comment on the future perspectives how this residue could be utilized in the kinase inhibitor design.

Targeting Histone Deacetylases: Opportunities for Cancer Treatment and Chemoprevention

The dysregulation of gene expression is a critical event involved in all steps of tumorigenesis. Aberrant histone and non-histone acetylation modifications of gene expression due to the abnormal activation of histone deacetylases (HDAC) have been reported in hematologic and solid types of cancer. In this sense, the cancer-associated epigenetic alterations are promising targets for anticancer therapy and chemoprevention. HDAC inhibitors (HDACi) induce histone hyperacetylation within target proteins, altering cell cycle and proliferation, cell differentiation, and the regulation of cell death programs. Over the last three decades, an increasing number of synthetic and naturally derived compounds, such as dietary-derived products, have been demonstrated to act as HDACi and have provided biological and molecular insights with regard to the role of HDAC in cancer. The first part of this review is focused on the biological roles of the Zinc-dependent HDAC family in malignant diseases. Accordingly, the small-molecules and natural products such as HDACi are described in terms of cancer therapy and chemoprevention. Furthermore, structural considerations are included to improve the HDACi selectivity and combinatory potential with other specific targeting agents in bifunctional inhibitors and proteolysis targeting chimeras. Additionally, clinical trials that combine HDACi with current therapies are discussed, which may open new avenues in terms of the feasibility of HDACi’s future clinical applications in precision cancer therapies.