The mechanism of activation of monomeric B-Raf V600E

An overview of RAF kinases and their inhibitors (2019-2023)

Protein kinases (PKs) including RAF, perform a principal role in regulating countless cellular events such as cell growth, differentiation, and angiogenesis. Overexpression and mutation of RAF kinases are significant contributors to the development and spread of cancer. Therefore, RAF kinase inhibitors show promising outcomes as anti-cancer small molecules by suppressing the expression of RAF protein, blocking RAS/RAF interaction, or inhibiting RAF enzymes. Currently, there are insufficient reports about approving drugs with minimal degree of toxicity. Therefore, it is an urgent need to develop new RAF kinase inhibitors correlated with increased anticancer activity and lower cytotoxicity. This review outlines reported RAF kinase inhibitors for cancer treatment in patents and literature from 2019 to 2023. It highlights the available inhibitors by shedding light on their chemical structures, biochemical profiles, and current status. Additionally, we highlighted the hinge region-binding moiety of the reported compounds by showing the hydrogen bond patterns of representative inhibitors with the hinge region for each class. In recent years, RAF kinase inhibitors have gained considerable attention in cancer research and drug development due to their potential to be studied under clinical trials and their demonstration of various degrees of efficacy and safety profiles across different cancer types. However, addressing challenges related to drug resistance and safety represents a major avenue for the optimization and enhancement of RAF kinase inhibitors. Strategies to overcome such obstacles were discussed such as developing novel pan-RAF inhibitors, RAF dimer inhibitors, and combination treatments.

Representing mutations for predicting cancer drug response

MOTIVATION

Predicting cancer drug response requires a comprehensive assessment of many mutations present across a tumor genome. While current drug response models generally use a binary mutated/unmutated indicator for each gene, not all mutations in a gene are equivalent.

RESULTS

Here, we construct and evaluate a series of predictive models based on leading methods for quantitative mutation scoring. Such methods include VEST4 and CADD, which score the impact of a mutation on gene function, and CHASMplus, which scores the likelihood a mutation drives cancer. The resulting predictive models capture cellular responses to dabrafenib, which targets BRAF-V600 mutations, whereas models based on binary mutation status do not. Performance improvements generalize to other drugs, extending genetic indications for PIK3CA, ERBB2, EGFR, PARP1, and ABL1 inhibitors. Introducing quantitative mutation features in drug response models increases performance and mechanistic understanding.

AVAILABILITY AND IMPLEMENTATION

Code and example datasets are available at https://github.com/pgwall/qms.

Targeting the RAS upstream and downstream signaling pathway for cancer treatment

Cancer often involves the overactivation of RAS/RAF/MEK/ERK (MAPK) and PI3K-Akt-mTOR pathways due to mutations in genes like RAS, RAF, PTEN, and PIK3CA. Various strategies are employed to address the overactivation of these pathways, among which targeted therapy emerges as a promising approach. Directly targeting specific proteins, leads to encouraging results in cancer treatment. For instance, RTK inhibitors such as imatinib and afatinib selectively target these receptors, hindering ligand binding and reducing signaling initiation. These inhibitors have shown potent efficacy against Non-Small Cell Lung Cancer. Other inhibitors, like lonafarnib targeting Farnesyltransferase and GGTI 2418 targeting geranylgeranyl Transferase, disrupt post-translational modifications of proteins. Additionally, inhibition of proteins like SOS, SH2 domain, and Ras demonstrate promising anti-tumor activity both in vivo and in vitro. Targeting downstream components with RAF inhibitors such as vemurafenib, dabrafenib, and sorafenib, along with MEK inhibitors like trametinib and binimetinib, has shown promising outcomes in treating cancers with BRAF-V600E mutations, including myeloma, colorectal, and thyroid cancers. Furthermore, inhibitors of PI3K (e.g., apitolisib, copanlisib), AKT (e.g., ipatasertib, perifosine), and mTOR (e.g., sirolimus, temsirolimus) exhibit promising efficacy against various cancers such as Invasive Breast Cancer, Lymphoma, Neoplasms, and Hematological malignancies. This review offers an overview of small molecule inhibitors targeting specific proteins within the RAS upstream and downstream signaling pathways in cancer.

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.

PreMode predicts mode-of-action of missense variants by deep graph representation learning of protein sequence and structural context

Accurate prediction of the functional impact of missense variants is important for disease gene discovery, clinical genetic diagnostics, therapeutic strategies, and protein engineering. Previous efforts have focused on predicting a binary pathogenicity classification, but the functional impact of missense variants is multi-dimensional. Pathogenic missense variants in the same gene may act through different modes of action (i.e., gain/loss-of-function) by affecting different aspects of protein function. They may result in distinct clinical conditions that require different treatments. We developed a new method, PreMode, to perform gene-specific mode-of-action predictions. PreMode models effects of coding sequence variants using SE(3)-equivariant graph neural networks on protein sequences and structures. Using the largest-to-date set of missense variants with known modes of action, we showed that PreMode reached state-of-the-art performance in multiple types of mode-of-action predictions by efficient transfer-learning. Additionally, PreMode's prediction of G/LoF variants in a kinase is consistent with inactive-active conformation transition energy changes. Finally, we show that PreMode enables efficient study design of deep mutational scans and optimization in protein engineering.

Mechanism of Dimer Selectivity and Binding Cooperativity of BRAF Inhibitors

Aberrant signaling of BRAFV600E is a major cancer driver. Current FDA-approved RAF inhibitors selectively inhibit the monomeric BRAFV600E and suffer from tumor resistance. Recently, dimer-selective and equipotent RAF inhibitors have been developed; however, the mechanism of dimer selectivity is poorly understood. Here, we report extensive molecular dynamics (MD) simulations of the monomeric and dimeric BRAFV600E in the apo form or in complex with one or two dimer-selective (PHI1) or equipotent (LY3009120) inhibitor(s). The simulations uncovered the unprecedented details of the remarkable allostery in BRAFV600E dimerization and inhibitor binding. Specifically, dimerization retrains and shifts the αC helix inward and increases the flexibility of the DFG motif; dimer compatibility is due to the promotion of the αC-in conformation, which is stabilized by a hydrogen bond formation between the inhibitor and the αC Glu501. A more stable hydrogen bond further restrains and shifts the αC helix inward, which incurs a larger entropic penalty that disfavors monomer binding. This mechanism led us to propose an empirical way based on the co-crystal structure to assess the dimer selectivity of a BRAFV600E inhibitor. Simulations also revealed that the positive cooperativity of PHI1 is due to its ability to preorganize the αC and DFG conformation in the opposite protomer, priming it for binding the second inhibitor. The atomically detailed view of the interplay between BRAF dimerization and inhibitor allostery as well as cooperativity has implications for understanding kinase signaling and contributes to the design of protomer selective RAF inhibitors.

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.

Cancer and Autism: How PTEN Mutations Degrade Function at the Membrane and Isoform Expression in the Human Brain

Mutations causing loss of PTEN lipid phosphatase activity can promote cancer, benign tumors (PHTS), and neurodevelopmental disorders (NDDs). Exactly how they preferentially trigger distinct phenotypic outcomes has been puzzling. Here, we demonstrate that PTEN mutations differentially allosterically bias P loop dynamics and its connection to the catalytic site, affecting catalytic activity. NDD-related mutations are likely to sample conformations of the functional wild-type state, while sampled conformations for the strong, cancer-related driver mutation hotspots favor catalysis-primed conformations, suggesting that NDD mutations are likely to be weaker, and our large-scale simulations show why. Prenatal PTEN isoform expression data suggest exons 5 and 7, which harbor NDD mutations, as cancer-risk carriers. Since cancer requires more than a single mutation, our conformational and genomic analysis helps discover how same protein mutations can foster different clinical manifestations, articulates a role for co-occurring background latent driver mutations, and uncovers relationships of splicing isoform expression to life expectancy.

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.

Allosteric coupling asymmetry mediates paradoxical activation of BRAF by type II inhibitors

The type II class of RAF inhibitors currently in clinical trials paradoxically activate BRAF at subsaturating concentrations. Activation is mediated by induction of BRAF dimers, but why activation rather than inhibition occurs remains unclear. Using biophysical methods tracking BRAF dimerization and conformation we built an allosteric model of inhibitor-induced dimerization that resolves the allosteric contributions of inhibitor binding to the two active sites of the dimer, revealing key differences between type I and type II RAF inhibitors. For type II inhibitors the allosteric coupling between inhibitor binding and BRAF dimerization is distributed asymmetrically across the two dimer binding sites, with binding to the first site dominating the allostery. This asymmetry results in efficient and selective induction of dimers with one inhibited and one catalytically active subunit. Our allosteric models quantitatively account for paradoxical activation data measured for 11 RAF inhibitors. Unlike type II inhibitors, type I inhibitors lack allosteric asymmetry and do not activate BRAF homodimers. Finally, NMR data reveal that BRAF homodimers are dynamically asymmetric with only one of the subunits locked in the active αC-in state. This provides a structural mechanism for how binding of only a single αC-in inhibitor molecule can induce potent BRAF dimerization and activation.

Challenges and Opportunities in the Crusade of BRAF Inhibitors: From 2002 to 2022

Serine/threonine-protein kinase B-Raf (BRAF; RAF = rapidly accelerated fibrosarcoma) plays an important role in the mitogen-activated protein kinase (MAPK) signaling cascade. Somatic mutations in the BRAF gene were first discovered in 2002 by Davies et al., which was a major breakthrough in cancer research. Subsequently, three different classes of BRAF mutants have been discovered. This class includes class I monomeric mutants (BRAFV600), class II BRAF homodimer mutants (non-V600), and class III BRAF heterodimers (non-V600). Cancers caused by these include melanoma, thyroid cancer, ovarian cancer, colorectal cancer, nonsmall cell lung cancer, and others. In this study, we have highlighted the major binding pockets in BRAF protein, their active and inactive conformations with inhibitors, and BRAF dimerization and its importance in paradoxical activation and BRAF mutation. We have discussed the first-, second-, and third-generation drugs approved by the Food and Drug Administration and drugs under clinical trials with all four different binding approaches with DFG-IN/OUT and αC-IN/OUT for BRAF protein. We have investigated particular aspects and difficulties with all three generations of inhibitors. Finally, this study has also covered recent developments in synthetic BRAF inhibitors (from their discovery in 2002 to 2022), their unique properties, and importance in inhibiting BRAF mutants.

Molecular analyses of the C-terminal CRAF variants associated with cardiomyopathy reveal their opposing impacts on the active conformation of the kinase domain

The germline mutations in the C-terminus of CRAF kinase, particularly L603, and S612T/L613V, are associated with congenital heart disorders, for example, dilated cardiomyopathy (DCM) and hypertrophic cardiomyopathy (HCM). The experimental data suggest that genetic alternation at position 603 impairs, while those at positions 612/613 enhance the CRAF kinase activity. However, the underlying mechanistic details by which these mutations increase or decrease kinase activity remain elusive. Therefore, we applied molecular dynamic simulation to investigate the impacts of these point mutations on the conformation of the CRAF kinase domain. The results revealed that the substitution of Leucine 603 for proline transits the kinase domain to a state that exhibits the molecular hallmarks of an inactive kinase, for example, a closed activation loop, 'αC-helix out' conformation and a distorted regulatory hydrophobic spine. However, two HCM-associated variants (S612T and L613V) show features of an active conformation, such as an open activation loop conformation, 'αC-helix in', the assembly of the hydrophobic spine, and more surface-exposed catalytic residues of phosphoryl transfer reaction. Overall, our study provides a mechanistic basis for the contradictory effects of the CRAF variants associated with HCM and DCM.

How PTEN mutations degrade function at the membrane and life expectancy of carriers of mutations in the human brain

PTEN dysfunction, caused by loss of lipid phosphatase activity or deletion, promotes pathologies, cancer, benign tumors, and neurodevelopmental disorders (NDDs). Despite efforts, exactly how the mutations trigger distinct phenotypic outcomes, cancer or NDD, has been puzzling. It has also been unclear how to distinguish between mutations harbored by isoforms, are they cancer or NDDs-related. Here we address both. We demonstrate that PTEN mutations differentially allosterically bias P-loop dynamics and its connection to the catalytic site, affecting catalytic activity. NDD-related mutations are likely to sample conformations present in the wild-type, while sampled conformations sheltering cancer-related hotspots favor catalysis-prone conformations, suggesting that NDD mutations are weaker. Analysis of isoform expression data indicates that if the transcript has NDD-related mutations, alone or in combination with cancer hotspots, there is high prenatal expression. If no mutations within the measured days, low expression levels. Cancer mutations promote stronger signaling and cell proliferation; NDDs' are weaker, influencing brain cell differentiation. Further, exon 5 is impacted by NDD or non-NDD mutations, while exon 7 is exclusively impacted by NDD mutations. Our comprehensive conformational and genomic analysis helps discover how same allele mutations can foster different clinical manifestations and uncovers correlations of splicing isoform expression to life expectancy.

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.

Insights into pralsetinib resistance to the non-gatekeeper RET kinase G810C mutation through molecular dynamics simulations

OBJECTIVE

RET (rearranged during transfection) kinase, as a transmembrane receptor tyrosine kinase, is a therapeutic target for several human cancer such as non-small cell lung cancer (NSCLC) and thyroid cancer. Pralsetinib is a recently approved drug for the treatment of RET-driven NSCLC and thyroid cancers. A single point mutation G810C at the C-lobe of the RET kinase causes pralsetinib resistance to this non-gatekeeper variant. However, the detailed mechanism remains poorly understood.

METHODS

Here, multiple microsecond molecular dynamics (MD) simulations, molecular mechanics/generalized born surface area (MM/GBSA) binding free energy calculations, and community network analysis were performed to reveal the mechanism of pralsetinib resistance to the RET G810C mutant.

RESULTS

The simulations showed that the G810C mutation had a minor effect on the overall conformational dynamics of the RET kinase domain. Energetic analysis suggested that the G810C mutation reduced the binding affinity of pralsetinib to the mutant. Per-residue energy contribution and structural analyses revealed that the hydrogen bonding interactions between pralsetinib and the hinge residues Glu805 and Ala807 were disrupted in the G810C mutant, which were responsible for the decreased binding affinity of pralsetinib to the mutant.

CONCLUSIONS

The obtained results may provide understanding of the mechanism of pralsetinib resistance to the non-gatekeeper RET G810C mutant.

Mechanistic insights into the role of calcium in the allosteric regulation of the calmodulin-regulated death-associated protein kinase

Calcium (Ca²⁺) signaling plays an important role in the regulation of many cellular functions. Ca²⁺-binding protein calmodulin (CaM) serves as a primary effector of calcium function. Ca²⁺/CaM binds to the death-associated protein kinase 1 (DAPK1) to regulate intracellular signaling pathways. However, the mechanism underlying the influence of Ca²⁺ on the conformational dynamics of the DAPK1-CaM interactions is still unclear. Here, we performed large-scale molecular dynamics (MD) simulations of the DAPK1-CaM complex in the Ca²⁺-bound and-unbound states to reveal the importance of Ca²⁺. MD simulations revealed that removal of Ca²⁺ increased the anti-correlated inter-domain motions between DAPK1 and CaM, which weakened the DAPK1-CaM interactions. Binding free energy calculations validated the decreased DAPK1-CaM interactions in the Ca²⁺-unbound state. Structural analysis further revealed that Ca²⁺ removal caused the significant conformational changes at the DAPK1-CaM interface, especially the helices α1, α2, α4, α6, and α7 from the CaM and the basic loop and the phosphate-binding loop from the DAPK1. These results may be useful to understand the biological role of Ca²⁺ in physiological processes.

Establishment and characterization of urothelial carcinoma cell lines with and without BRAF mutation (V595E) in dogs

Each 5 urothelial carcinoma (UC) cell lines with and without the v-Raf murine sarcoma virus oncogene homolog B (BRAF) gene mutation (V595E) were established and examined V595E-related tumorigenic characteristics in dogs. No typical morphological features were observed in cloned cells with and without V595E. The cell proliferation of both cloned cells showed logarithmic growth curve and those doubling time were 24.9 ± 4.1 h in V595E ( +) and 29.3 ± 11.3 h in V595E ( -). On the growth curve of xenotransplanted tumor in severe combined immunodeficiency mice, 3 out of 5 V595E ( +) and 2 out of 5 V595E ( -) cloned cells revealed gradually and remarkably increasing curve, indicating clearly tumorigenicity. The xenotransplanted tumors with V595E ( +) showed typical features of UC, such as solid proliferation of pleomorphic tumor cells, formation of papillary structure, and glandular structure. Additionally, various vascular formation was observed, probably indicating an advanced growth phase of UC. In mitogen-activated protein kinase (MAPK) signaling pathway, cytoplasmic phosphorylated-BRAF (pBRAF) and cytoplasmic and nuclear phosphorylated-ERK1/2 (pERK1/2) were detected in all 4 tumors with V595E ( +), whereas only cytoplasmic and nuclear pERK1/2 was detected in tumors with V595E ( -). Since V595E can directly activate MAPK signaling pathway, coincidence of V595E with pBRAF (phosphor Thr598/Ser601) indicates acquired resistance to BRAF inhibitors. These established UC cell lines, especially V595E ( +) cell lines, are useful tool for understanding pathophysiological states and controlling therapeutic manners of UC in dogs.

Mechanistic Insights into the Protection Effect of Argonaute-RNA Complex on the HCV Genome

While host miRNA usually plays an antiviral role, the relentless tides of viral evolution have carved out a mechanism to recruit host miRNA as a viral protector. By complementing miR-122 at the 5' end of the genome, the hepatitis C virus (HCV) gene can form a complex with Argonaute 2 (Ago2) protein to protect the 5' end of HCV RNA from exonucleolytic attacks. Experiments showed that the disruption of the stem-loop 1(SL1) structure and the 9th nucleotide (T9) of HCV site 1 RNA could enhance the affinity of the Ago2 protein to the HCV site 1 RNA (target RNA). However, the underlying mechanism of how the conformation and dynamics of the Ago2: miRNA: target RNA complex is affected by the SL1 and T9 remains unclear. To address this, we performed large-scale molecular dynamics simulations on the AGO2-miRNA complex binding with the WT target, T9-abasic target and SL1-disruption target, respectively. The results revealed that the T9 and SL1 structures could induce the departing motion of the PAZ, PIWI and N domains, propping up the mouth of the central groove which accommodates the target RNA, causing the instability of the target RNA and disrupting the Ago2 binding. The coordinated motion among the PAZ, PIWI and N domains were also weakened by the T9 and SL1 structures. Moreover, we proposed a new model wherein the Ago2 protein could adopt a more constraint conformation with the proximity and more correlated motions of the PAZ, N and PIWI domains to protect the target RNA from dissociation. These findings reveal the mechanism of the Ago2-miRNA complex's protective effect on the HCV genome at the atomic level, which will offer guidance for the design of drugs to confront the protection effect and engineering of Ago2 as a gene-regulation tool.

Coincidence of v-raf murine sarcoma viral oncogene homolog B mutation (V595E) with phosphorylated v-raf murine sarcoma viral oncogene homolog B in urothelial carcinoma in dogs

Expression of phosphorylated v-raf murine sarcoma viral oncogene homolog B (pBRAF) and phosphorylated extracellular signal-regulated kinase1/2 (pERK1/2) were investigated in urothelial carcinoma (UC) in dogs with or without the BRAF gene mutation (V595E). Among the 10 cases of UC with V595E (-), cytoplasmic immunoreactivity against pBRAF of neoplastic cells was reported in 8, with 7 displaying moderate reactivity and 1 displaying intense reactivity. Nuclear immunoreactivity against pBRAF was detected in 5 cases; however, these reactivities were non-specific, due to pBRAF being limited in the cytoplasm. In addition, positive cytoplasmic immunoreactivity against pERK1/2 of neoplastic cells was detected in 7 cases and nuclear immunoreactivity against ERK1/2 was detected in 6 cases. Among the 13 cases of UC with V595E (+), cytoplasmic immunoreactivity against pBRAF of neoplastic cells was detected in all 13 cases and nuclear immunoreactivity against pBRAF was detected in 10 cases; however, the nuclear immunoreactivity was non-specific. Cytoplasmic immunoreactivity against pERK1/2 of neoplastic cells was detected in all 13 cases and nuclear immunoreactivity against pERK1/2 was also detected in all cases. As nuclear pERK1/2 indicates a progressive signaling process in the mitogen-activated protein kinase pathway, V595E (+) UC might be in its growing stage. Probable phosphorylated sites of pBRAF at Thr598/Ser601, detected in this study, are major and essential sites of the upstream rat sarcoma viral oncogene homolog (RAS) signaling pathway. In human cancers, the BRAF mutation never coincides with oncogenic RAS. To our knowledge, this is the first report on the simultaneous occurrence of the BRAF mutation (V595E) and pBRAF expression (at Thr598/Ser601) in dogs with UC with V595E (+).

Allostery: Allosteric Cancer Drivers and Innovative Allosteric Drugs

Here, we discuss the principles of allosteric activating mutations, propagation downstream of the signals that they prompt, and allosteric drugs, with examples from the Ras signaling network. We focus on Abl kinase where mutations shift the landscape toward the active, imatinib binding-incompetent conformation, likely resulting in the high affinity ATP outcompeting drug binding. Recent pharmacological innovation extends to allosteric inhibitor (GNF-5)-linked PROTAC, targeting Bcr-Abl1 myristoylation site, and broadly, allosteric heterobifunctional degraders that destroy targets, rather than inhibiting them. Designed chemical linkers in bifunctional degraders can connect the allosteric ligand that binds the target protein and the E3 ubiquitin ligase warhead anchor. The physical properties and favored conformational state of the engineered linker can precisely coordinate the distance and orientation between the target and the recruited E3. Allosteric PROTACs, noncompetitive molecular glues, and bitopic ligands, with covalent links of allosteric ligands and orthosteric warheads, increase the effective local concentration of productively oriented and placed ligands. Through covalent chemical or peptide linkers, allosteric drugs can collaborate with competitive drugs, degrader anchors, or other molecules of choice, driving innovative drug discovery.

Structure of the RAF1-HSP90-CDC37 complex reveals the basis of RAF1 regulation

RAF kinases are RAS-activated enzymes that initiate signaling through the MAPK cascade to control cellular proliferation, differentiation, and survival. Here, we describe the structure of the full-length RAF1 protein in complex with HSP90 and CDC37 obtained by cryoelectron microscopy. The reconstruction reveals a RAF1 kinase with an unfolded N-lobe separated from its C-lobe. The hydrophobic core of the N-lobe is trapped in the HSP90 dimer, while CDC37 wraps around the chaperone and interacts with the N- and C-lobes of the kinase. The structure indicates how CDC37 can discriminate between the different members of the RAF family. Our structural analysis also reveals that the folded RAF1 assembles with 14-3-3 dimers, suggesting that after folding RAF1 follows a similar activation as B-RAF. Finally, disruption of the interaction between CDC37 and the DFG segment of RAF1 unveils potential vulnerabilities in attempting the pharmacological degradation of RAF1 for therapeutic purposes.

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

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

Mechanistic Insights into the Mechanism of Inhibitor Selectivity toward the Dark Kinase STK17B against Its High Homology STK17A

As a member of the death-associated protein kinase (DAPK) family, STK17B plays an important role in the regulation of cellular apoptosis and has been considered as a promising drug target for hepatocellular carcinoma. However, the highly conserved ATP-binding site of protein kinases represents a challenge to design selective inhibitors for a specific DAPK isoform. In this study, molecular docking, multiple large-scale molecular dynamics (MD) simulations, and binding free energy calculations were performed to decipher the molecular mechanism of the binding selectivity of PKIS43 toward STK17B against its high homology STK17A. MD simulations revealed that STK17A underwent a significant conformational arrangement of the activation loop compared to STK17B. The binding free energy predictions suggested that the driving force to control the binding selectivity of PKIS43 was derived from the difference in the protein–ligand electrostatic interactions. Furthermore, the per-residue free energy decomposition unveiled that the energy contribution from Arg41 at the phosphate-binding loop of STK17B was the determinant factor responsible for the binding specificity of PKIS43. This study may provide useful information for the rational design of novel and potent selective inhibitors toward STK17B.

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.

The structural basis of BCR-ABL recruitment of GRB2 in chronic myelogenous leukemia

BCR-ABL drives chronic myeloid leukemia (CML). BCR binding to GRB2 transduces signaling via the Ras/MAPK pathway. Despite considerable data confirming the binding, molecular-level understanding of exactly how the two proteins interact, and especially, what are the determinants of the specificity of the SH2^GRB2 domain-pBCR recognition are still open questions. Yet, this is vastly important for understanding binding selectivity, and for predicting the phosphorylated receptors, or peptides, that are likely to bind. Here, we uncover these determinants and ascertain to what extent they relate to the affinity of the interaction. Toward this end, we modeled the complexes of the phosphorylated BCR (pBCR) and SH2^GRB2 and other pY/Y-peptide-SH2 complexes and compared their specificity and affinity. We observed that pBCR's ₁₇₆FpYVNV₁₈₀ motif is favorable and specific to SH2^GRB2, similar to pEGFR, but not other complexes. SH2^GRB2 contains two binding pockets: pY-binding recognition pocket triggers binding, and the specificity pocket whose interaction is governed by N179 in pBCR and W121 in SH2^GRB2. Our proposed motif with optimal affinity to SH2^GRB2 is E/D-pY-E/V-N-I/L. Collectively, we provide the structural basis of BCR-ABL recruitment of GRB2, outline its specificity hallmarks, and delineate a blueprint for prediction of BCR-binding scaffolds and for therapeutic peptide design.

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.

On the development of B-Raf inhibitors acting through innovative mechanisms

B-Raf is a protein kinase participating to the regulation of many biological processes in cells. Several studies have demonstrated that this protein is frequently upregulated in human cancers, especially when it bears activating mutations. In the last years, few ATP-competitive inhibitors of B-Raf have been marketed for the treatment of melanoma and are currently under clinical evaluation on a variety of other types of cancer. Although the introduction of drugs targeting B-Raf has provided significant advances in cancer treatment, responses to ATP-competitive inhibitors remain limited, mainly due to selectivity issues, side effects, narrow therapeutic windows, and the insurgence of drug resistance.

Impressive research efforts have been made so far towards the identification of novel ATP-competitive modulators with improved efficacy against cancers driven by mutant Raf monomers and dimers, some of them showing good promises. However, several limitations could still be envisioned for these compounds, according to literature data. Besides, increased attentions have arisen around approaches based on the design of allosteric modulators, polypharmacology, proteolysis targeting chimeras (PROTACs) and drug repurposing for the targeting of B-Raf proteins. The design of compounds acting through such innovative mechanisms is rather challenging. However, valuable therapeutic opportunities can be envisioned on these drugs, as they act through innovative mechanisms in which limitations typically observed for approved ATP-competitive B-Raf inhibitors are less prone to emerge. In this article, current approaches adopted for the design of non-ATP competitive inhibitors targeting B-Raf are described, discussing also on the possibilities, ligands acting through such innovative mechanisms could provide for the obtainment of more effective therapies.

Novel potential oncogenic and druggable mutations of FGFRs recur in the kinase domain across cancer types

Fibroblast growth factor receptors (FGFRs) are recurrently altered by single nucleotide variants (SNVs) in many human cancers. The prevalence of SNVs in FGFRs depends on the cancer type. In some tumors, such as the urothelial carcinoma, mutations of FGFRs occur at very high frequency (up to 60%). Many characterized mutations occur in the extracellular or transmembrane domains, while fewer known mutations are found in the kinase domain. In this study, we performed a bioinformatics analysis to identify novel putative cancer driver or therapeutically actionable mutations of the kinase domain of FGFRs. To pinpoint those mutations that may be clinically relevant, we exploited the recurrence of alterations on analogous amino acid residues within the kinase domain (PK_Tyr_Ser-Thr) of different kinases as a predictor of functional impact. By exploiting MutationAligner and LowMACA bioinformatics resources, we highlighted novel uncharacterized mutations of FGFRs which recur in other protein kinases. By revealing unanticipated correspondence with known variants, we were able to infer their functional effects, as alterations clustering on similar residues in analogous proteins have a high probability to elicit similar effects. As FGFRs represent an important class of oncogenes and drug targets, our study opens the way for further studies to validate their driver and/or actionable nature and, in the long term, for a more efficacious application of precision oncology.

Kincore: a web resource for structural classification of protein kinases and their inhibitors

The active form of kinases is shared across different family members, as are several commonly observed inactive forms. We previously performed a clustering of the conformation of the activation loop of all protein kinase structures in the Protein Data Bank (PDB) into eight classes based on the dihedral angles that place the Phe side chain of the DFG motif at the N-terminus of the activation loop. Our clusters are strongly associated with the placement of the activation loop, the C-helix, and other structural elements of kinases. We present Kincore, a web resource providing access to our conformational assignments for kinase structures in the PDB. While other available databases provide conformational states or drug type but not both, KinCore includes the conformational state and the inhibitor type (Type 1, 1.5, 2, 3, allosteric) for each kinase chain. The user can query and browse the database using these attributes or determine the conformational labels of a kinase structure using the web server or a standalone program. The database and labeled structure files can be downloaded from the server. Kincore will help in understanding the conformational dynamics of these proteins and guide development of inhibitors targeting specific states. Kincore is available at http://dunbrack.fccc.edu/kincore.

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 mechanism of Raf activation through dimerization

Raf, a threonine/serine kinase in the Raf/MEK/ERK pathway, regulates cell proliferation. Raf's full activation requires dimerization. Aberrant activation through dimerization is an important therapeutic target. Despite its clinical importance, fundamental questions, such as how the side-to-side dimerization promotes the OFF-to-ON transition of Raf's kinase domain and how the fully activated ON-state kinase domain is stabilized in the dimer for Raf signaling, remain unanswered. Herein, we decipher an atomic-level mechanism of Raf activation through dimerization, clarifying this enigma. The mechanism reveals that the replacement of intramolecular π–π stacking by intermolecular π–π stacking at the dimer interface releases the structural constraint of the αC-helix, promoting the OFF-to-ON transition. During the transition, the inhibitory hydrophobic interactions were disrupted, making the phosphorylation sites in A-loop approach the HRD motif for cis-autophosphorylation. Once fully activated, the ON-state kinase domain can be stabilized by a newly identified functional N-terminal basic (NtB) motif in the dimer for Raf signaling. This work provides atomic level insight into critical steps in Raf activation and outlines a new venue for drug discovery against Raf dimerization.

We decipher an atomic-level mechanism of Raf activation through dimerization, revealing that the disruption of intramolecular π–π stacking at the dimer interface promotes the OFF-to-ON transition.

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.