The mechanism of PI3Kα activation at the atomic level

Free energy landscape of the PI3Kα C-terminal activation

The gene PIK3CA, encoding the catalytic subunit p110α of PI3Kα, is the second most frequently mutated gene in cancer, with the highest frequency oncogenic mutants occurring in the C-terminus of the kinase domain. The C-terminus has a dual function in regulating the kinase, playing a putative auto-inhibitory role for kinase activity and being absolutely essential for binding to the cell membrane. However, the molecular mechanisms by which these C-terminal oncogenic mutations cause PI3Kα overactivation remain unclear. To understand how a spectrum of C-terminal mutations of PI3Kα alter kinase activity compared to the WT, we perform unbiased and biased Molecular Dynamics simulations of several C-terminal mutants and report the free energy landscapes for the C-terminal "closed-to-open" transition in the WT, H1047R, G1049R, M1043L and N1068KLKR mutants. Results are consistent with HDX-MS experimental data and provide a molecular explanation why H1047R and G1049R reorient the C-terminus with a different mechanism compared to the WT and M1043L and N1068KLKR mutants. Moreover, we show that in the H1047R mutant, the cavity, where the allosteric ligands STX-478 and RLY-2608 bind, is more accessible contrary to the WT. This study provides insights into the molecular mechanisms underlying activation of oncogenic PI3Kα by C-terminal mutations and represents a valuable resource for continued efforts in the development of mutant selective inhibitors as therapeutics.

Quantum biochemistry description of PI3Kα enzyme bound to selective inhibitors

The PI3K class I is composed of four PI3K isoforms that serve as regulatory enzymes governing cellular metabolism, proliferation, and survival. The hyperactivation of PI3Kα is observed in various types of cancer and is linked to poor prognosis. Unfortunately, the development inhibitors selectively targeting one of the isoforms remains challenging, with only few agents in clinical use. The main difficulty arises from the high conservation among residues at the ATP-binding pocket across isoforms, which also serves as target pocket for inhibitors. In this work, molecular dynamics and quantum calculations were performed to investigate the molecular features guiding the binding of selective inhibitors, alpelisib and GDC-0326, into the ATP-binding pocket of PI3Kα. While molecular dynamics allowed crystallographic coordinates to relax, the interaction eergy between each amino acid residues and inhibitors was obtained by combining the Molecular Fractionation with Conjugated Caps scheme with Density Functional Theory calculations. In addition, the atomic charge of ligands in the bound and unbound (free) was calculated. Results indicated that the most relevant residues for the binding of alpelisib are Ile932, Glu859, Val851, Val850, Tyr836, Met922, Ile800, and Ile848, while the most important residues for the binding of GDC-0326 are Ile848, Ile800, Ile932, Gln859, Glu849, and Met922. In addition, residues Trp780, Ile800, Tyr836, Ile848, Gln859 Val850, Val851, Ile932 and Met922 are common hotspots for both inhibitors. Overall, the results from this work contribute to improving the understanding of the molecular mechanisms controlling selectivity and highlight important interactions to be considered during the rational design of new agents.Communicated by Ramaswamy H. Sarma.

Ras, RhoA, and vascular pharmacology in neurodevelopment and aging

Small GTPases Ras, Rac, and RhoA are crucial regulators of cellular functions. They also act in dysregulated cell proliferation and transformation. Multiple publications have focused on illuminating their roles and mechanisms, including in immune system pathologies. Their functions in neurology-related diseases, neurodegeneration and neurodevelopment, are also emerging, as well as their potential as pharmacological targets in both pathologies. Observations increasingly suggest that these pathologies may relate to activation (or suppression) of signaling by members of the Ras superfamily, especially Ras, Rho, and Rac isoforms, and components of their signaling pathways. Germline (or embryonic) mutations that they harbor are responsible for neurodevelopmental disorders, such as RASopathies, autism spectrum disorder, and dilated cardiomyopathy. In aging, they promote neurodegenerative diseases, with Rho GTPase featuring in their pharmacology, as in the case of Alzheimer's disease (AD). Significantly, drugs with observed anti-AD activity, particularly those involved in cardiovascular systems, are associated with the RhoA signaling, as well as cerebral vasculature in brain development and aging. This leads us to suggest that anti-AD drugs could inform neurodevelopmental disorders, including pediatric low-grade gliomas pharmacology. Neurodevelopmental disorders associated with RhoA, like autism, are also connected with vascular systems, thus could be targets of vascular system-connected drugs.

Capturing Autoinhibited PDK1 Reveals the Linker's Regulatory Role, Informing Innovative Inhibitor Design

PDK1 is crucial for PI3K/AKT/mTOR and Ras/MAPK cancer signaling. It phosphorylates AKT in a PIP3-dependent but S6K, SGK, and RSK kinases in a PIP3-independent manner. Unlike its substrates, its autoinhibited monomeric state has been unclear, likely due to its low population time, and phosphorylation in the absence of PIP3 has been puzzling too. Here, guided by experimental data, we constructed models and performed all-atom molecular dynamics simulations. In the autoinhibited PDK1 conformation that resembles autoinhibited AKT, binding of the linker between the kinase and PH domains to the PIF-binding pocket promotes the formation of the Glu130-Lys111 salt bridge and weakens the association of the kinase domain with the PH domain, shifting the population from the autoinhibited state to states accessible to the membrane and its kinase substrates. The interaction of the substrates' hydrophobic motif and the PDK1 PIF-binding pocket facilitates the release of the autoinhibition even in the absence of PIP3. Phosphorylation of the serine-rich motif within the linker further attenuates the association of the PH domain with the kinase domain. These suggest that while the monomeric autoinhibited state is relatively stable, it can readily shift to its active, catalysis-prone state to phosphorylate its diverse substrates. Our findings reveal the PDK1 activation mechanism and discover the regulatory role of PDK1's linker, which lead to two innovative linker-based inhibitor strategies: (i) locking the autoinhibited PDK1 through optimization of the interactions of AKT inhibitors with the PH domain of PDK1 and (ii) analogs (small molecules or peptidomimetics) that mimic the linker interactions with the PIF-binding pocket.

Hsp90 Promotes Gastric Cancer Cell Metastasis and Stemness by Regulating the Regional Distribution of Glycolysis-Related Metabolic Enzymes in the Cytoplasm

Heat-shock protein 90 (Hsp90) plays a crucial role in tumorigenesis and tumor progression; however, its mechanism of action in gastric cancer (GC) remains unclear. Here, the role of Hsp90 in GC metabolism is the focus of this research. High expression of Hsp90 in GC tissues can interact with glycolysis, collectively affecting prognosis in clinical samples. Both in vitro and in vivo experiments demonstrate that Hsp90 is able to regulate the migration and stemness properties of GC cells. Metabolic phenotype analyses indicate that Hsp90 influences glycolytic metabolism. Mechanistically, Hsp90 interacts with glycolysis-related enzymes, forming multienzyme complexes to enhance glycolysis efficiency and yield. Additionally, Hsp90 binds to cytoskeleton-related proteins, regulating the regional distribution of glycolytic enzymes at the cell margin and lamellar pseudopods. This effect could lead to a local increase in efficient energy supply from glycolysis, further promoting epithelial-mesenchymal transition (EMT) and metastasis. In summary, Hsp90, through its interaction with metabolic enzymes related to glycolysis, forms multi-enzyme complexes and regulates regional distribution of glycolysis by dynamic cytoskeletal adjustments, thereby promoting the migration and stemness of GC cells. These conclusions also support the potential for a combined targeted approach involving Hsp90, glycolysis, and the cytoskeleton in clinical therapy.

The value of protein allostery in rational anticancer drug design: an update

INTRODUCTION

Allosteric drugs are advantageous. However, they still face hurdles, including identification of allosteric sites that will effectively alter the active site. Current strategies largely focus on identifying pockets away from the active sites into which the allosteric ligand will dock and do not account for exactly how the active site is altered. Favorable allosteric inhibitors dock into sites that are nearby the active sites and follow nature, mimicking diverse allosteric regulation strategies.

AREAS COVERED

The following article underscores the immense significance of allostery in drug design, describes current allosteric strategies, and especially offers a direction going forward. The article concludes with the authors' expert perspectives on the subject.

EXPERT OPINION

To select a productive venue in allosteric inhibitor development, we should learn from nature. Currently, useful strategies follow this route. Consider, for example, the mechanisms exploited in relieving autoinhibition and in harnessing allosteric degraders. Mimicking compensatory, or rescue mutations may also fall into such a thesis, as can molecular glues that capture features of scaffolding proteins. Capturing nature and creatively tailoring its mimicry can continue to innovate allosteric drug discovery.

A deep learning-based theoretical protocol to identify potentially isoform-selective PI3Kα inhibitors

Phosphoinositide 3-kinase alpha (PI3Kα) is one of the most frequently dysregulated kinases known for their pivotal role in many oncogenic diseases. While the side effects linked to existing drugs against PI3Kα-induced cancers provide an avenue for further research, the significant structural conservation among PI3Ks makes it extremely difficult to develop new isoform-selective PI3Kα inhibitors. Embracing this challenge, we herein designed a hybrid protocol by integrating machine learning (ML) with in silico drug-designing strategies. A deep learning classification model was developed and trained on the physicochemical descriptors data of known PI3Kα inhibitors and used as a screening filter for a database of small molecules. This approach led us to the prediction of 662 compounds showcasing appropriate features to be considered as PI3Kα inhibitors. Subsequently, a multiphase molecular docking was applied to further characterize the predicted hits in terms of their binding affinities and binding modes in the targeted cavity of the PI3Kα. As a result, a total of 12 compounds were identified whereas the best poses highlighted the efficiency of these ligands in maintaining interactions with the crucial residues of the protein to be targeted for the inhibition of associated activity. Notably, potential activity of compound 12 in counteracting PI3Kα function was found in a previous in vitro study. Following the drug-likeness and pharmacokinetic characterizations, six compounds (compounds 1, 2, 3, 6, 7, and 11) with suitable ADME-T profiles and promising bioavailability were selected. The mechanistic studies in dynamic mode further endorsed the potential of identified hits in blocking the ATP-binding site of the receptor with higher binding affinities than the native inhibitor, alpelisib (BYL-719), particularly the compounds 1, 2, and 11. These outcomes support the reliability of the developed classification model and the devised computational strategy for identifying new isoform-selective drug candidates for PI3Kα inhibition.

Cryo-EM structures reveal two allosteric inhibition modes of PI3KαH1047R involving a re-shaping of the activation loop

PI3Kα is a lipid kinase that phosphorylates PIP2 and generates PIP3. The hyperactive PI3Kα mutation, H1047R, accounts for about 14% of breast cancer, making it a highly attractive target for drug discovery. Here, we report the cryo-EM structures of PI3KαH1047R bound to two different allosteric inhibitors QR-7909 and QR-8557 at a global resolution of 2.7 Å and 3.0 Å, respectively. The structures reveal two distinct binding pockets on the opposite sides of the activation loop. Structural and MD simulation analyses show that the allosteric binding of QR-7909 and QR-8557 inhibit PI3KαH1047R hyper-activity by reducing the fluctuation and mobility of the activation loop. Our work provides a strong rational basis for a further optimization and development of highly selective drug candidates to treat PI3KαH1047R-driven cancers.

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.

A p85 isoform switch enhances PI3K activation on endosomes by a MAP4- and PI3P-dependent mechanism

Phosphatidylinositol 3-kinase α (PI3Kα) is a heterodimer of p110α catalytic and p85 adaptor subunits that is activated by agonist-stimulated receptor tyrosine kinases. Although p85α recruits p110α to activated receptors on membranes, p85α loss, which occurs commonly in cancer, paradoxically promotes agonist-stimulated PI3K/Akt signaling. p110α localizes to microtubules via microtubule-associated protein 4 (MAP4), facilitating its interaction with activated receptor kinases on endosomes to initiate PI3K/Akt signaling. Here, we demonstrate that in response to agonist stimulation and p85α knockdown, the residual p110α, coupled predominantly to p85β, exhibits enhanced recruitment with receptor tyrosine kinases to endosomes. Moreover, the p110α C2 domain binds PI3-phosphate, and this interaction is also required to recruit p110α to endosomes and for PI3K/Akt signaling. Stable knockdown of p85α, which mimics the reduced p85α levels observed in cancer, enhances cell growth and tumorsphere formation, and these effects are abrogated by MAP4 or p85β knockdown, underscoring their role in the tumor-promoting activity of p85α loss.

Structural insights into the activation mechanism of phosphoinositide 3-kinase alpha

Phosphoinositide 3-kinases (PI3Ks) are lipid kinases known to regulate important cellular functions by phosphorylating the inositol ring of inositol-phospholipids (PtdIns) at 3' position. The PI3Kα is a heterodimer and the activation of the catalytic subunit (p110α) is regulated by its regulatory subunit (p85α). The current work deals with studying the activation mechanism of the PI3Kα using multi micro-second molecular dynamic simulations. Structural changes involved in activation mechanism is studied by gradually releasing the inhibitory effects of different domains of regulatory subunit namely, n-terminal SH2 (nSH2) and inter SH2 (iSH2). The observation shows that even in the presence of n-terminal and inter SH2 domain (niSH2) of regulatory subunit, the catalytic domain has some intrinsic activation activity and the presence of c-terminal SH2 (cSH2) domain may be required for complete inhibition. The release of nSH2 domain leads to loss of interactions between iSH2 domain (regulatory subunit) and C2 and kinase domain (catalytic subunit). The study shows that early events in the activation mechanism involve the movement of the ABD domain of the catalytic subunit along with the linker region between ABD and RBD region which may lead to movement of ABD closer to the CLobe of the kinase domain. This movement is essentially as it triggers the rearrangement of CLobe especially the catalytic loop and activation loop which bring catalytic important residues closer to ATP and PIP2(phosphatidylinositol-4,5-bisphosphate). Water mediated interaction analysis reveal that water may be playing an important role in the transfer of phosphate from ATP to PIP2. The study shows that initial signal for release of inhibitory effect of the regulatory subunit might be propagated through the linker region between ABD and RBD through allosteric effect to different regions of the protein. These understanding of early events during the activation mechanism may help in the design of better therapeutic targeting PI3K.

Direct K-Ras Inhibitors to Treat Cancers: Progress, New Insights, and Approaches to Treat Resistance

Here we discuss approaches to K-Ras inhibition and drug resistance scenarios. A breakthrough offered a covalent drug against K-RasG12C. Subsequent innovations harnessed same-allele drug combinations, as well as cotargeting K-RasG12C with a companion drug to upstream regulators or downstream kinases. However, primary, adaptive, and acquired resistance inevitably emerge. The preexisting mutation load can explain how even exceedingly rare mutations with unobservable effects can promote drug resistance, seeding growth of insensitive cell clones, and proliferation. Statistics confirm the expectation that most resistance-related mutations are in cis, pointing to the high probability of cooperative, same-allele effects. In addition to targeted Ras inhibitors and drug combinations, bifunctional molecules and innovative tri-complex inhibitors to target Ras mutants are also under development. Since the identities and potential contributions of preexisting and evolving mutations are unknown, selecting a pharmacologic combination is taxing. Collectively, our broad review outlines considerations and provides new insights into pharmacology and resistance.

The allosteric mechanism of mTOR activation can inform bitopic inhibitor optimization

mTOR serine/threonine kinase is a cornerstone in the PI3K/AKT/mTOR pathway. Yet, the detailed mechanism of activation of its catalytic core is still unresolved, likely due to mTOR complexes' complexity. Its dysregulation was implicated in cancer and neurodevelopmental disorders. Using extensive molecular dynamics (MD) simulations and compiled published experimental data, we determine exactly how mTOR's inherent motifs can control the conformational changes in the kinase domain, thus kinase activity. We also chronicle the critical regulation by the unstructured negative regulator domain (NRD). When positioned inside the catalytic cleft (NRD IN state), mTOR tends to adopt a deep and closed catalytic cleft. This is primarily due to the direct interaction with the FKBP-rapamycin binding (FRB) domain which restricts it, preventing substrate access. Conversely, when outside the catalytic cleft (NRD OUT state), mTOR favors an open conformation, exposing the substrate-binding site on the FRB domain. We further show how an oncogenic mutation (L2427R) promotes shifting the mTOR ensemble toward the catalysis-favored state. Collectively, we extend mTOR's "active-site restriction" mechanism and clarify mutation action. In particular, our mechanism suggests that RMC-5552 (RMC-6272) bitopic inhibitors may benefit from adjustment of the (PEG8) linker length when targeting certain mTOR variants. In the cryo-EM mTOR/RMC-5552 structure, the distance between the allosteric and orthosteric inhibitors is ∼22.7 Å. With a closed catalytic cleft, this linker bridges the sites. However, in our activation mechanism, in the open cleft it expands to ∼24.7 Å, offering what we believe to be the first direct example of how discovering an activation mechanism can potentially increase the affinity of inhibitors targeting mutants.

Oliveira AS, Mulholland AJ, Serapian SA, Colombo G. Decrypting Allostery in Membrane-Bound K-Ras4B Using Complementary In Silico Approaches Based on Unbiased Molecular Dynamics Simulations

Protein functions are dynamically regulated by allostery, which enables conformational communication even between faraway residues, and expresses itself in many forms, akin to different "languages": allosteric control pathways predominating in an unperturbed protein are often unintuitively reshaped whenever biochemical perturbations arise (e.g., mutations). To accurately model allostery, unbiased molecular dynamics (MD) simulations require integration with a reliable method able to, e.g., detect incipient allosteric changes or likely perturbation pathways; this is because allostery can operate at longer time scales than those accessible by plain MD. Such methods are typically applied singularly, but we here argue their joint application─as a "multilingual" approach─could work significantly better. We successfully prove this through unbiased MD simulations (∼100 μs) of the widely studied, allosterically active oncotarget K-Ras4B, solvated and embedded in a phospholipid membrane, from which we decrypt allostery using four showcase "languages": Distance Fluctuation analysis and the Shortest Path Map capture allosteric hotspots at equilibrium; Anisotropic Thermal Diffusion and Dynamical Non-Equilibrium MD simulations assess perturbations upon, respectively, either superheating or hydrolyzing the GTP that oncogenically activates K-Ras4B. Chosen "languages" work synergistically, providing an articulate, mutually coherent, experimentally consistent picture of K-Ras4B allostery, whereby distinct traits emerge at equilibrium and upon GTP cleavage. At equilibrium, combined evidence confirms prominent allosteric communication from the membrane-embedded hypervariable region, through a hub comprising helix α5 and sheet β5, and up to the active site, encompassing allosteric "switches" I and II (marginally), and two proposed pockets. Upon GTP cleavage, allosteric perturbations mostly accumulate on the switches and documented interfaces.

SHP2 clinical phenotype, cancer, or RASopathies, can be predicted by mutant conformational propensities

SHP2 phosphatase promotes full activation of the RTK-dependent Ras/MAPK pathway. Its mutations can drive cancer and RASopathies, a group of neurodevelopmental disorders (NDDs). Here we ask how same residue mutations in SHP2 can lead to both cancer and NDD phenotypes, and whether we can predict what the outcome will be. We collected and analyzed mutation data from the literature and cancer databases and performed molecular dynamics simulations of SHP2 mutants. We show that both cancer and Noonan syndrome (NS, a RASopathy) mutations favor catalysis-prone conformations. As to cancer versus RASopathies, we demonstrate that cancer mutations are more likely to accelerate SHP2 activation than the NS mutations at the same genomic loci, in line with NMR data for K-Ras4B more aggressive mutations. The compiled experimental data and dynamic features of SHP2 mutants lead us to propose that different from strong oncogenic mutations, SHP2 activation by NS mutations is less likely to induce a transition of the ensemble from the SHP2 inactive state to the active state. Strong signaling promotes cell proliferation, a hallmark of cancer. Weak, or moderate signals are associated with differentiation. In embryonic neural cells, dysregulated differentiation is connected to NDDs. Our innovative work offers structural guidelines for identifying and correlating mutations with clinical outcomes, and an explanation for why bearers of RASopathy mutations may have a higher probability of cancer. Finally, we propose a drug strategy against SHP2 variants-promoting cancer and RASopathies.

p55γ degrades RIP3 via MG53 to suppress ischaemia-induced myocardial necroptosis and mediates cardioprotection of preconditioning

AIMS

Regulated necrosis (necroptosis) and apoptosis are important biological features of myocardial infarction, ischaemia-reperfusion (I/R) injury, and heart failure. However, the molecular mechanisms underlying myocardial necroptosis remain elusive. Ischaemic preconditioning (IPC) is the most powerful intrinsic cardioprotection against myocardial I/R injury. In this study, we aimed to determine whether IPC suppresses I/R-induced necroptosis and the underlying molecular mechanisms.

METHODS AND RESULTS

We generated p55γ transgenic and knockout mice and used ligation of left anterior descending coronary artery to produce an in vivo I/R model. The effects of p55γ and its downstream molecules were subsequently identified using mass spectroscopy and co-immunoprecipitation and pulldown assays. We found that p55γ expression was down-regulated in failing human myocardium caused by coronary heart disease as well as in I/R mouse hearts. Cardiac-specific p55γ overexpression ameliorated the I/R-induced necroptosis. In striking contrast, p55γ deficiency (p55γ-/-) and cardiac-specific deletion of p55γ (p55γc-KO) worsened I/R-induced injury. IPC up-regulated p55γ expression in vitro and in vivo. Using reporter and chromatin immunoprecipitation assays, we found that Hif1α transcriptionally regulated p55γ expression and mediated the cardioprotection of IPC. IPC-mediated suppression of necroptosis was attenuated in p55γ-/- and p55γc-KO hearts. Mechanistically, p55γ overexpression decreased the protein levels of RIP3 rather than the mRNA levels, while p55γ deficiency increased the protein abundance of RIP3. IPC attenuated the I/R-induced up-regulation of RIP3, which was abolished in p55γ-deficient mice. Up-regulation of RIP3 attenuated the p55γ- or IPC-induced inhibition of necroptosis in vivo. Importantly, p55γ directly bound and degraded RIP3 in a ubiquitin-dependent manner. We identified MG53 as the E3 ligase that mediated the p55γ-induced degradation of RIP3. In addition, we also found that p55γ activated the RISK pathway during IPC.

CONCLUSIONS

Our findings reveal that activation of the MG53-RIP3 signal pathway by p55γ protects the heart against I/R-induced necroptosis and underlies IPC-induced cardioprotection.

Protein conformational ensembles in function: roles and mechanisms

The sequence-structure-function paradigm has dominated twentieth century molecular biology. The paradigm tacitly stipulated that for each sequence there exists a single, well-organized protein structure. Yet, to sustain cell life, function requires (i) that there be more than a single structure, (ii) that there be switching between the structures, and (iii) that the structures be incompletely organized. These fundamental tenets called for an updated sequence-conformational ensemble-function paradigm. The powerful energy landscape idea, which is the foundation of modernized molecular biology, imported the conformational ensemble framework from physics and chemistry. This framework embraces the recognition that proteins are dynamic and are always interconverting between conformational states with varying energies. The more stable the conformation the more populated it is. The changes in the populations of the states are required for cell life. As an example, in vivo, under physiological conditions, wild type kinases commonly populate their more stable "closed", inactive, conformations. However, there are minor populations of the "open", ligand-free states. Upon their stabilization, e.g., by high affinity interactions or mutations, their ensembles shift to occupy the active states. Here we discuss the role of conformational propensities in function. We provide multiple examples of diverse systems, including protein kinases, lipid kinases, and Ras GTPases, discuss diverse conformational mechanisms, and provide a broad outlook on protein ensembles in the cell. We propose that the number of molecules in the active state (inactive for repressors), determine protein function, and that the dynamic, relative conformational propensities, rather than the rigid structures, are the hallmark of cell life.

Praeruptorin C alleviates cognitive impairment in type 2 diabetic mice through restoring PI3K/AKT/GSK3β pathway

Diabetic encephalopathy is a common consequence of diabetes mellitus that causes cognitive dysfunction and neuropsychiatric disorders. Praeruptorin C (Pra-C) from the traditional Chinese medicinal herb Peucedanum praeruptorum Dunn. is a potential antioxidant and neuroprotective agent. This study was conducted to investigate the molecular mechanisms underlying the effect of Pra-C on diabetic cognitive impairment. A novel object recognition test and the Morris water maze test were performed to assess the behavioral performance of mice. Electrophysiological recordings were made to monitor synaptic plasticity in the hippocampus. A protein-protein interaction network of putative Pra-C targets was constructed, and molecular docking simulations were performed to predict the potential mechanisms of the action of Pra-C. Protein expression levels were detected by western blotting. Pra-C administration significantly lowered body weight and fasting blood glucose levels and alleviated learning and memory deficits in type 2 diabetic mice. Network pharmacology and molecular docking results suggested that Pra-C affects the PI3K/AKT/GSK3β signaling pathway. Western blot analysis confirmed significant increases in phosphorylated PI3K, AKT, and GSK3β levels in vivo and in vitro upon Pra-C administration. Pra-C alleviated cognitive impairment in type 2 diabetic mice by activating PI3K/AKT/GSK3β pathway.

Strategy toward Kinase-Selective Drug Discovery

Kinase drug selectivity is the ground challenge in cancer research. Due to the structurally similar kinase drug pockets, off-target inhibitor toxicity has been a major cause for clinical trial failures. The pockets are similar but not identical. Here, we describe a transformation invariant protocol to identify distinct geometric features in the drug pocket that can distinguish one kinase from all others. We integrate available experimental structures with the artificial intelligence-based structural kinome, performing a kinome-wide structural bioinformatic analysis to establish the structural principles of kinase drug selectivity. We generate the structural landscape from the experimental kinase-ligand complexes and propose a binary network that encapsulates the information. The results show that all kinases contain binary units that are shared by less than seven other kinases in the kinome. 331 kinases contain unique binary units that may distinguish them from all others. The structural features encoded by these binary units in the network represent the inhibitor-accessible geometric space that may capture the kinome-wide selectivity. Our proposed binary network with the unsupervised clustering can serve as a general structural bioinformatic protocol for extracting the distinguishing structural features for any protein from their families. We apply the binary network to epidermal growth factor receptor tyrosine kinase inhibitor selectivity by targeting the gate area and the AKT1 serine/threonine kinase selectivity by binding to the αC-helix region and the allosteric pocket. Finally, we develop the cross-platform software, KDS (Kinase Drug Selectivity), for customized visualization and analysis of the binary networks in the human kinome (https://github.com/CBIIT/KDS).

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.

Addressing the Reciprocal Crosstalk between the AR and the PI3K/AKT/mTOR Signaling Pathways for Prostate Cancer Treatment

The reduction in androgen synthesis and the blockade of the androgen receptor (AR) function by chemical castration and AR signaling inhibitors represent the main treatment lines for the initial stages of prostate cancer. Unfortunately, resistance mechanisms ultimately develop due to alterations in the AR pathway, such as gene amplification or mutations, and also the emergence of alternative pathways that render the tumor less or, more rarely, completely independent of androgen activation. An essential oncogenic axis activated in prostate cancer is the phosphatidylinositol-3-kinase (PI3K)/AKT/mammalian target of rapamycin (mTOR) pathway, as evidenced by the frequent alterations of the negative regulator phosphatase and tensin homolog (PTEN) and by the activating mutations in PI3K subunits. Additionally, crosstalk and reciprocal feedback loops between androgen signaling and the PI3K/AKT/mTOR signaling cascade that activate pro-survival signals and play an essential role in disease recurrence and progression have been evidenced. Inhibitors addressing different players of the PI3K/AKT/mTOR pathway have been evaluated in the clinic. Only a limited benefit has been reported in prostate cancer up to now due to the associated side effects, so novel combination approaches and biomarkers predictive of patient response are urgently needed. Here, we reviewed recent data on the crosstalk between AR signaling and the PI3K/AKT/mTOR pathway, the selective inhibitors identified, and the most advanced clinical studies, with a focus on combination treatments. A deeper understanding of the complex molecular mechanisms involved in disease progression and treatment resistance is essential to further guide therapeutic approaches with improved outcomes.

The allosteric effect of the upper half of SENP1 contributes to its substrate selectivity for SUMO1 over SUMO2

SUMOylation regulates various cellular process and SENP1 (SUMO-specific protease 1) serves as a SUMO (small ubiquitin-related modifier) specific protease that participates in the SUMO cycle. Given its extensive influences on metabolic activities, SENP1 has gained more and more attentions in clinical treatments. However, there remains a question on why does the SENP1 prefer to process SUMO1 rather than SUMO2. Here, we performed molecular dynamics simulations of SENP1-SUMO1, SENP1-SUMO2, and apo SENP1 systems and observed distinct conformational dynamics in the upper half of the clamp and the three loops in the catalytic center of the SENP1. Principal component analysis revealed that the most prominent canonical variable represented the spatial distribution of the upper half of the clamp, while the openness of the cleft was closely related to the catalytic ability of SENP1. Further analysis of the SENP1-SUMO interactions revealed that the extensive and strong interactions between the SENP1 and SUMO1 were both in the interface of the upper half region and the catalytic center. Dynamic cross-correlation matrix analysis demonstrated that the inter-residue correlations in the SUMO1 system was much stronger, especially in the two essential regions belonging to the upper and lower half of cleft. Based on these observations, we proposed an allosteric propagation model and further testified it using the community analysis. These results revealed the propagation pathway of allosteric communication that contributed to the substrate discrimination of SENP1 upon SUMO1 and SUMO2.Communicated by Ramaswamy H. Sarma.

Allosteric Boost by TAB1 on the TAK1 Kinase Favorably Sculpts the Thermodynamic Landscape of Activation

The intricate mechanisms of allosteric regulation in kinases are of general interest to the scientific community for potential therapeutic implications. However, the diversity among kinases and their regulatory routes requires a case-by-case study to widen the repertoire of known mechanisms. The present study achieves this by understanding TAK1 kinase activation by TAB1 as a model phenomenon for the first time. Despite the known capacity of TAK1 to switch between its inactive ("DFG-out") and active-like ("DFG-in") conformations, the questionable role of TAB1 in offering an energetic favor to this has been addressed here using sequential combination of enhanced sampling methods like targeted molecular dynamics (TMD) and Gaussian accelerated molecular dynamics (GaMD). It reveals how a minimal domain of TAB1 sufficiently acts like a "catalytic gear" by favorably sculpting TAK1's thermodynamic landscape (potential of mean force in 2D) that accelerates "in"-"out" conformational switching of the conserved DFG motif. Standard molecular dynamics simulations (∼5 μs) reveal that TAB1 fascinatingly exploits the "lever-like" αF helix of TAK1 kinase domain to remotely propel the DFG motif via subtle helical "unfolding-folding" modifications within the kinase activation loop. The presence of two charged residues on terminal poles of αF helix imparts it, with this unique "lever-like" utility, and this turns out to be one important signature of co-evolution between TAK1 and TAB1. The entire mechanism of TAB1's impact transduction, which is found to be analogous to the moves in the popular "Chinese checker" game, gives a clear proof of the "dynamics-driven allostery" concept in kinases. The findings further benchmark TAK1's known autophosphorylation capacity. A novel insight into kinase allostery is thus provided, which potentiates investigation of similar capacities in other kinases.

Targeting RAS mutants in malignancies: successes, failures, and reasons for hope

RAS genes are the most frequently mutated oncogenes and play critical roles in the development and progression of malignancies. The mutation, isoform (KRAS, HRAS, and NRAS), position, and type of substitution vary depending on the tissue types. Despite decades of developing RAS-targeted therapies, only small subsets of these inhibitors are clinically effective, such as the allele-specific inhibitors against KRASG12C . Targeting the remaining RAS mutants would require further experimental elucidation of RAS signal transduction, RAS-altered metabolism, and the associated immune microenvironment. This study reviews the mechanisms and efficacy of novel targeted therapies for different RAS mutants, including KRAS allele-specific inhibitors, combination therapies, immunotherapies, and metabolism-associated therapies.

Higher-order interactions of Bcr-Abl can broaden chronic myeloid leukemia (CML) drug repertoire

Bcr-Abl, a nonreceptor tyrosine kinase, is associated with leukemias, especially chronic myeloid leukemia (CML). Deletion of Abl's N-terminal region, to which myristoyl is linked, renders the Bcr-Abl fusion oncoprotein constitutively active. The substitution of Abl's N-terminal region by Bcr enables Bcr-Abl oligomerization. Oligomerization is critical: it promotes clustering on the membrane, which is essential for potent MAPK signaling and cell proliferation. Here we decipher the Bcr-Abl specific, step-by-step oligomerization process, identify a specific packing surface, determine exactly how the process is structured and identify its key elements. Bcr's coiled coil (CC) domain at the N-terminal controls Bcr-Abl oligomerization. Crystallography validated oligomerization via Bcr-Abl dimerization between two Bcr CC domains, with tetramerization via tight packing between two binary assemblies. However, the structural principles guiding Bcr CC domain oligomerization are unknown, hindering mechanistic understanding and drugs exploiting it. Using molecular dynamics (MD) simulations, we determine that the binary complex of the Bcr CC domain serves as a basic unit in the quaternary complex providing a specific surface for dimer-dimer packing and higher-order oligomerization. We discover that the small α1-helix is the key. In the binary assembly, the helix forms interchain aromatic dimeric packing, and in the quaternary assembly, it contributes to the specific dimer-dimer packing. Our mechanism is supported by the experimental literature. It offers the key elements controlling this process which can expand the drug discovery strategy, including by Bcr CC-derived peptides, and candidate residues for small covalent drugs, toward quenching oligomerization, supplementing competitive and allosteric tyrosine kinase inhibitors.

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.

The orchestrated signaling by PI3Kα and PTEN at the membrane interface

The oncogene PI3Kα and the tumor suppressor PTEN represent two antagonistic enzymatic activities that regulate the interconversion of the phosphoinositide lipids PI(4,5)P2 and PI(3,4,5)P3 in membranes. As such, they are defining components of phosphoinositide-based cellular signaling and membrane trafficking pathways that regulate cell survival, growth, and proliferation, and are often deregulated in cancer. In this review, we highlight aspects of PI3Kα and PTEN interplay at the intersection of signaling and membrane trafficking. We also discuss the mechanisms of PI3Kα- and PTEN- membrane interaction and catalytic activation, which are fundamental for our understanding of the structural and allosteric implications on signaling at the membrane interface and may aid current efforts in pharmacological targeting of these proteins.

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.

The Role of Conformational Dynamics and Allostery in the Control of Distinct Efficacies of Agonists to the Glucocorticoid Receptor

Glucocorticoid receptor (GR) regulates various cellular functions. Given its broad influence on metabolic activities, it has been the target of drug discovery for decades. However, how drugs induce conformational changes in GR has remained elusive. Herein, we used five GR agonists (dex, AZ938, pred, cor, and dibC) with different efficacies to investigate which aspect of the ligand induced the differences in efficacy. We performed molecular dynamics simulations on the five systems (dex-, AZ938-, pred-, cor-, and dibC-bound systems) and observed a distinct discrepancy in the conformation of the cofactor TIF2. Moreover, we discovered ligand-induced differences regarding the level of conformational changes posed by the binding of cofactor TIF2 and identified a pair of essential residues D590 and T39. We further found a positive correlation between the efficacies of ligands and the interaction of the two binding pockets' domains, where D590 and T739 were involved, implying their significance in the participation of allosteric communication. Using community network analysis, two essential communities containing D590 and T739 were identified with their connectivity correlating to the efficacy of ligands. The potential communication pathways between these two residues were revealed. These results revealed the underlying mechanism of allosteric communication between the ligand-binding and cofactor-binding pockets and identified a pair of important residues in the allosteric communication pathway, which can serve as a guide for future drug discovery.

Aspirin sensitivity of PIK3CA-mutated Colorectal Cancer: potential mechanisms revisited

PIK3CA mutations are amongst the most prevalent somatic mutations in cancer and are associated with resistance to first-line treatment along with low survival rates in a variety of malignancies. There is evidence that patients carrying PIK3CA mutations may benefit from treatment with acetylsalicylic acid, commonly known as aspirin, particularly in the setting of colorectal cancer. In this regard, it has been clarified that Class IA Phosphatidylinositol 3-kinases (PI3K), whose catalytic subunit p110α is encoded by the PIK3CA gene, are involved in signal transduction that regulates cell cycle, cell growth, and metabolism and, if disturbed, induces carcinogenic effects. Although PI3K is associated with pro-inflammatory cyclooxygenase-2 (COX-2) expression and signaling, and COX-2 is among the best-studied targets of aspirin, the mechanisms behind this clinically relevant phenomenon are still unclear. Indeed, there is further evidence that the protective, anti-carcinogenic effect of aspirin in this setting may be mediated in a COX-independent manner. However, until now the understanding of aspirin's prostaglandin-independent mode of action is poor. This review will provide an overview of the current literature on this topic and aims to analyze possible mechanisms and targets behind the aspirin sensitivity of PIK3CA-mutated cancers.

Insights into the Allosteric Effect of SENP1 Q597A Mutation on the Hydrolytic Reaction of SUMO1 via an Integrated Computational Study

Small ubiquitin-related modifier (SUMO)-specific protease 1 (SENP1) is a cysteine protease that catalyzes the cleavage of the C-terminus of SUMO1 for the processing of SUMO precursors and deSUMOylation of target proteins. SENP1 is considered to be a promising target for the treatment of hepatocellular carcinoma (HCC) and prostate cancer. SENP1 Gln597 is located at the unstructured loop connecting the helices α4 to α5. The Q597A mutation of SENP1 allosterically disrupts the hydrolytic reaction of SUMO1 through an unknown mechanism. Here, extensive multiple replicates of microsecond molecular dynamics (MD) simulations, coupled with principal component analysis, dynamic cross-correlation analysis, community network analysis, and binding free energy calculations, were performed to elucidate the detailed mechanism. Our MD simulations showed that the Q597A mutation induced marked dynamic conformational changes in SENP1, especially in the unstructured loop connecting the helices α4 to α5 which the mutation site occupies. Moreover, the Q597A mutation caused conformational changes to catalytic Cys603 and His533 at the active site, which might impair the catalytic activity of SENP1 in processing SUMO1. Moreover, binding free energy calculations revealed that the Q597A mutation had a minor effect on the binding affinity of SUMO1 to SENP1. Together, these results may broaden our understanding of the allosteric modulation of the SENP1-SUMO1 complex.

Kinases on Double Duty: A Review of UniProtKB Annotated Bifunctionality within the Kinome

Phosphorylation facilitates the regulation of all fundamental biological processes, which has triggered extensive research of protein kinases and their roles in human health and disease. In addition to their phosphotransferase activity, certain kinases have evolved to adopt additional catalytic functions, while others have completely lost all catalytic activity. We searched the Universal Protein Resource Knowledgebase (UniProtKB) database for bifunctional protein kinases and focused on kinases that are critical for bacterial and human cellular homeostasis. These kinases engage in diverse functional roles, ranging from environmental sensing and metabolic regulation to immune-host defense and cell cycle control. Herein, we describe their dual catalytic activities and how they contribute to disease pathogenesis.

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.

Open Structural Data in Precision Medicine

Three-dimensional protein structural data at the molecular level are pivotal for successful precision medicine. Such data are crucial not only for discovering drugs that act to block the active site of the target mutant protein but also for clarifying to the patient and the clinician how the mutations harbored by the patient work. The relative paucity of structural data reflects their cost, challenges in their interpretation, and lack of clinical guidelines for their utilization. Rapid technological advancements in experimental high-resolution structural determination increasingly generate structures. Computationally, modeling algorithms, including molecular dynamics simulations, are becoming more powerful, as are compute-intensive hardware, particularly graphics processing units, overlapping with the inception of the exascale era. Accessible, freely available, and detailed structural and dynamical data can be merged with big data to powerfully transform personalized pharmacology. Here we review protein and emerging genome high-resolution data, along with means, applications, and examples underscoring their usefulness in precision medicine. Expected final online publication date for the Annual Review of Biomedical Data Science, Volume 5 is August 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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

Precision oncology benefits from effective early phase drug discovery decisions. Recently, drugging inactive protein conformations has shown impressive successes, raising the cardinal questions of which targets can profit and what are the principles of the active/inactive protein pharmacology. Cancer driver mutations have been established to mimic the protein activation mechanism. We suggest that the decision whether to target an inactive (or active) conformation should largely rest on the protein mechanism of activation. We next discuss the recent identification of double (multiple) same-allele driver mutations and their impact on cell proliferation and suggest that like single driver mutations, double drivers also mimic the mechanism of activation. We further suggest that the structural perturbations of double (multiple) in cis mutations may reveal new surfaces/pockets for drug design. Finally, we underscore the preeminent role of the cellular network which is deregulated in cancer. Our structure-based review and outlook updates the traditional Mechanism of Action, informs decisions, and calls attention to the intrinsic activation mechanism of the target protein and the rewired tumor-specific network, ushering innovative considerations in precision medicine.

How can same-gene mutations promote both cancer and developmental disorders?

The question of how same-gene mutations can drive both cancer and neurodevelopmental disorders has been puzzling. It has also been puzzling why those with neurodevelopmental disorders have a high risk of cancer. Ras, MEK, PI3K, PTEN, and SHP2 are among the oncogenic proteins that can harbor mutations that encode diseases other than cancer. Understanding why some of their mutations can promote cancer, whereas others promote neurodevelopmental diseases, and why even the same mutations may promote both phenotypes, has important clinical ramifications. Here, we review the literature and address these tantalizing questions. We propose that cell type–specific expression of the mutant protein, and of other proteins in the respective pathway, timing of activation (during embryonic development or sporadic emergence), and the absolute number of molecules that the mutations activate, alone or in combination, are pivotal in determining the pathological phenotypes—cancer and (or) developmental disorders.

Therapeutic implications of activating noncanonical PIK3CA mutations in head and neck squamous cell carcinoma

Alpelisib selectively inhibits the p110α catalytic subunit of PI3Kα and is approved for treatment of breast cancers harboring canonical PIK3CA mutations. In head and neck squamous cell carcinoma (HNSCC), 63% of PIK3CA mutations occur at canonical hotspots. The oncogenic role of the remaining 37% of PIK3CA noncanonical mutations is incompletely understood. We report a patient with HNSCC with a noncanonical PIK3CA mutation (Q75E) who exhibited a durable (12 months) response to alpelisib in a phase II clinical trial. Characterization of all 32 noncanonical PIK3CA mutations found in HNSCC using several functional and phenotypic assays revealed that the majority (69%) were activating, including Q75E. The oncogenic impact of these mutations was validated in 4 cellular models, demonstrating that their activity was lineage independent. Further, alpelisib exhibited antitumor effects in a xenograft derived from a patient with HNSCC containing an activating noncanonical PIK3CA mutation. Structural analyses revealed plausible mechanisms for the functional phenotypes of the majority of the noncanonical PIK3CA mutations. Collectively, these findings highlight the importance of characterizing the function of noncanonical PIK3CA mutations and suggest that patients with HNSCC whose tumors harbor activating noncanonical PIK3CA mutations may benefit from treatment with PI3Kα inhibitors.

Activation pathway of a G protein-coupled receptor uncovers conformational intermediates as targets for allosteric drug design

G protein-coupled receptors (GPCRs) are the most common proteins targeted by approved drugs. A complete mechanistic elucidation of large-scale conformational transitions underlying the activation mechanisms of GPCRs is of critical importance for therapeutic drug development. Here, we apply a combined computational and experimental framework integrating extensive molecular dynamics simulations, Markov state models, site-directed mutagenesis, and conformational biosensors to investigate the conformational landscape of the angiotensin II (AngII) type 1 receptor (AT1 receptor) - a prototypical class A GPCR-activation. Our findings suggest a synergistic transition mechanism for AT1 receptor activation. A key intermediate state is identified in the activation pathway, which possesses a cryptic binding site within the intracellular region of the receptor. Mutation of this cryptic site prevents activation of the downstream G protein signaling and β-arrestin-mediated pathways by the endogenous AngII octapeptide agonist, suggesting an allosteric regulatory mechanism. Together, these findings provide a deeper understanding of AT1 receptor activation at an atomic level and suggest avenues for the design of allosteric AT1 receptor modulators with a broad range of applications in GPCR biology, biophysics, and medicinal chemistry.

The Importance of Being PI3K in the RAS Signaling Network

Ras proteins are essential mediators of a multitude of cellular processes, and its deregulation is frequently associated with cancer appearance, progression, and metastasis. Ras-driven cancers are usually aggressive and difficult to treat. Although the recent Food and Drug Administration (FDA) approval of the first Ras G12C inhibitor is an important milestone, only a small percentage of patients will benefit from it. A better understanding of the context in which Ras operates in different tumor types and the outcomes mediated by each effector pathway may help to identify additional strategies and targets to treat Ras-driven tumors. Evidence emerging in recent years suggests that both oncogenic Ras signaling in tumor cells and non-oncogenic Ras signaling in stromal cells play an essential role in cancer. PI3K is one of the main Ras effectors, regulating important cellular processes such as cell viability or resistance to therapy or angiogenesis upon oncogenic Ras activation. In this review, we will summarize recent advances in the understanding of Ras-dependent activation of PI3K both in physiological conditions and cancer, with a focus on how this signaling pathway contributes to the formation of a tumor stroma that promotes tumor cell proliferation, migration, and spread.

The mechanism of activation of monomeric B-Raf V600E

Oncogenic mutations in the serine/threonine kinase B-Raf, particularly the V600E mutation, are frequent in cancer, making it a major drug target. Although much is known about B-Raf's active and inactive states, questions remain about the mechanism by which the protein changes between these two states. Here, we utilize molecular dynamics to investigate both wild-type and V600E B-Raf to gain mechanistic insights into the impact of the Val to Glu mutation. The results show that the wild-type and mutant follow similar activation pathways involving an extension of the activation loop and an inward motion of the αC-helix. The V600E mutation, however, destabilizes the inactive state by disrupting hydrophobic interactions present in the wild-type structure while the active state is stabilized through the formation of a salt bridge between Glu600 and Lys507. Additionally, when the activation loop is extended, the αC-helix is able to move between an inward and outward orientation as long as the DFG motif adopts a specific orientation. In that orientation Phe595 rotates away from the αC-helix, allowing the formation of a salt bridge between Lys483 and Glu501. These mechanistic insights have implications for the development of new Raf inhibitors.

The structural basis of Akt PH domain interaction with calmodulin

Akt plays a key role in the Ras/PI3K/Akt/mTOR signaling pathway. In breast cancer, Akt translocation to the plasma membrane is enabled by the interaction of its pleckstrin homology domain (PHD) with calmodulin (CaM). At the membrane, the conformational change promoted by PIP3 releases CaM and facilitates Thr308 and Ser473 phosphorylation and activation. Here, using modeling and molecular dynamics simulations, we aim to figure out how CaM interacts with Akt's PHD at the atomic level. Our simulations show that CaM-PHD interaction is thermodynamically stable and involves a β-strand rather than an α-helix, in agreement with NMR data, and that electrostatic and hydrophobic interactions are critical. The PHD interacts with CaM lobes; however, multiple modes are possible. IP4, the polar head of PIP3, weakens the CaM-PHD interaction, implicating the release mechanism at the plasma membrane. Recently, we unraveled the mechanism of PI3Kα activation at the atomistic level and the structural basis for Ras role in the activation. Here, our atomistic structural data clarify the mechanism of how CaM interacts, delivers, and releases Akt-the next node in the Ras/PI3K pathway-at the plasma membrane.

Mechanism of allosteric activation of SIRT6 revealed by the action of rationally designed activators

The recent discovery of activator compounds binding to an allosteric site on the NAD⁺-dependent protein lysine deacetylase, sirtuin 6 (SIRT6) has attracted interest and presents a pharmaceutical target for aging-related and cancer diseases. However, the mechanism underlying allosteric activation of SIRT6 by the activator MDL-801 remains largely elusive because no major conformational changes are observed upon activator binding. By combining molecular dynamics simulations with biochemical and kinetic analyses of wild-type SIRT6 and its variant M136A, we show that conformational rotation of 2-methyl-4-fluoro-5-bromo substituent on the right phenyl ring (R-ring) of MDL-801, which uncovers previously unseen hydrophobic interactions, contributes to increased activating deacetylation activity of SIRT6. This hypothesis is further supported by the two newly synthesized MDL-801 derivatives through the removal of the 5-Br atom on the R-ring (MDL-801-D1) or the restraint of the rotation of the R-ring (MDL-801-D2). We further propose that the 5-Br atom serves as an allosteric driver that controls the ligand allosteric efficacy. Our study highlights the effect of allosteric enzyme catalytic activity by activator binding and provides a rational approach for enhancing deacetylation activity.

Deciphering the selective binding mechanisms of anaplastic lymphoma kinase-derived neuroblastoma tumor neoepitopes to human leukocyte antigen

Neuroblastoma (NB), as a metastatic form of solid tumor, has a high fatality rate found in early childhood. The two anaplastic lymphoma kinase (ALK) neoepitopes nonamer and decamer used in cancer immunotherapy against NB cancer can selectively bind to the human leukocyte antigen (HLA-B*15:01) groove with high affinities, whereas the native self-peptide is unable to interact with the HLA-B*15:01. Here, we performed molecular dynamics (MD) simulations and subsequent molecular mechanics-generalized born surface area (MM-GBSA) binding free energy calculations to explore the selective binding mechanisms of nonamer and decamer to the HLA-B*15:01 against the self-peptide. MD simulations revealed the significant conformational dynamics of the self-peptide in the HLA-B*15:01 groove against the nonamer and decamer. Binding free energy calculations showed that the binding affinities of HLA-B*15:01-neoepitope complexes were followed in the order decamer > nonamer > self-peptide. Detailed analysis of HLA-B*15:01-neoepitope structural complexes showed that compared to the nonamer, the self-peptide tended to move outward to the solvent, whereas the decamer moved deep to the HLA-B*15:01 groove. These different dynamic observations of the ALK neoepitopes can explain the distinct binding affinities of self-peptide, nonamer, and decamer to the HLA-B*15:01. The results may be useful for the design of more selective ALK neoepitopes.

The mechanism of full activation of tumor suppressor PTEN at the phosphoinositide-enriched membrane

Tumor suppressor PTEN, the second most highly mutated protein in cancer, dephosphorylates signaling lipid PIP₃ produced by PI3Ks. Excess PIP₃ promotes cell proliferation. The mechanism at the membrane of this pivotal phosphatase is unknown hindering drug discovery. Exploiting explicit solvent simulations, we tracked full-length PTEN trafficking from the cytosol to the membrane. We observed its interaction with membranes composed of zwitterionic phosphatidylcholine, anionic phosphatidylserine, and phosphoinositides, including signaling lipids PIP₂ and PIP₃. We tracked its moving away from the zwitterionic and getting absorbed onto anionic membrane that harbors PIP₃. We followed it localizing on microdomains enriched in signaling lipids, as PI3K does, and observed PIP₃ allosterically unfolding the N-terminal PIP₂ binding domain, positioning it favorably for the polybasic motif interaction with PIP₂. Finally, we determined PTEN catalytic action at the membrane, all in line with experimental observations, deciphering the mechanisms of how PTEN anchors to the membrane and restrains cancer.

•

PTEN localizes on membrane microdomains enriched in phosphoinositides, as PI3K does

•

Full PTEN activation requires both signaling lipids, PIP₂ and PIP₃

•

Strong salt bridge interactions sustain stable PTEN membrane localization

•

Substrate-induced P loop conformational change implicates PTEN catalytic activity

Burn Ointment Promotes Cutaneous Wound Healing by Modulating the PI3K/AKT/mTOR Signaling Pathway

Burn ointment (BO) is a clinically useful medicine for the treatment of burns and scalds. However, there is no enough scientific evidence to report the effect of BO on wound healing and its analgesic and anti-inflammatory efficacy. The aim of this work was to evaluate the anti-inflammatory and analgesic efficacy of BO and to reveal the potential wound healing properties and related mechanisms of BO. In this work, the content of active ingredients of BO was determined by high-performance liquid chromatography (HPLC). Two animal models of inflammation were used to study its anti-inflammatory activity, and a hot plate method was used to evaluate its analgesic effect. In addition, mouse incision and rat burn models were used to investigate the effect of BO on the anti-inflammatory and wound healing mechanisms. The results showed that BO was safe for topical application, and BO could significantly inhibit auricular swelling in mice and paw swelling in rats and significantly prolong the latency period of paw licking in the hot plate experiment in mice. It can also accelerate wound healing and repair scars by promoting the formation of new epithelial tissues in rat burn models. In addition, BO significantly downregulated the serum level of TNF-α and significantly increased the serum levels of VEGF and TGF-β1. Also, BO promoted the expression of collagen I and increased the ratio in p-PI3K/PI3K, p-AKT/AKT, and p-mTOR/mTOR pathways. Our results demonstrate the safety and efficacy of BO and suggest that activation of the PI3K/AKT/mTOR signaling pathway may play an important role in the promotion of wound healing by BO.

Inhibition of Nonfunctional Ras

Intuitively, functional states should be targeted; not nonfunctional ones. So why could drugging the inactive K-Ras4BG12Cwork-but drugging the inactive kinase will likely not? The reason is the distinct oncogenic mechanisms. Kinase driver mutations work by stabilizing the active state and/or destabilizing the inactive state. Either way, oncogenic kinases are mostly in the active state. Ras driver mutations work by quelling its deactivation mechanisms, GTP hydrolysis, and nucleotide exchange. Covalent inhibitors that bind to the inactive GDP-bound K-Ras4BG12C conformation can thus work. By contrast, in kinases, allosteric inhibitors work by altering the active-site conformation to favor orthosteric drugs. From the translational standpoint this distinction is vital: it expedites effective pharmaceutical development and extends the drug classification based on the mechanism of action. Collectively, here we postulate that drug action relates to blocking the mechanism of activation, not to whether the protein is in the active or inactive state.

PI3K Driver Mutations: A Biophysical Membrane-Centric Perspective

Ras activates its effectors at the membrane. Active PI3Kα and its associated kinases/phosphatases assemble at membrane regions enriched in signaling lipids. In contrast, the Raf kinase domain extends into the cytoplasm and its assembly is away from the crowded membrane surface. Our structural membrane-centric outlook underscores the spatiotemporal principles of membrane and signaling lipids, which helps clarify PI3Kα activation. Here we focus on mechanisms of activation driven by PI3Kα driver mutations, spotlighting the PI3Kα double (multiple) activating mutations. Single mutations can be potent, but double mutations are stronger: their combination is specific, a single strong driver cannot fully activate PI3K, and two weak drivers may or may not do so. In contrast, two strong drivers may successfully activate PI3K, where one, for example, H1047R, modulates membrane interactions facilitating substrate binding at the active site (km) and the other, for example, E542K and E545K, reduces the transition state barrier (ka), releasing autoinhibition by nSH2. Although mostly unidentified, weak drivers are expected to be common, so we ask here how common double mutations are likely to be and why PI3Kα with double mutations responds effectively to inhibitors. We provide a structural view of hotspot and weak driver mutations in PI3Kα activation, explain their mechanisms, compare these with mechanisms of Raf activation, and point to targeting cell-specific, chromatin-accessible, and parallel (or redundant) pathways to thwart the expected emergence of drug resistance. Collectively, our biophysical outlook delineates activation and highlights the challenges of drug resistance.