How ATP-Competitive Inhibitors Allosterically Modulate Tyrosine Kinases That Contain a Src-like Regulatory Architecture

Conformation selection by ATP-competitive inhibitors and allosteric communication in ERK2

Activation of the extracellular signal-regulated kinase-2 (ERK2) by phosphorylation has been shown to involve changes in protein dynamics, as determined by hydrogen-deuterium exchange mass spectrometry (HDX-MS) and NMR relaxation dispersion measurements. These can be described by a global exchange between two conformational states of the active kinase, named 'L' and 'R,' where R is associated with a catalytically productive ATP-binding mode. An ATP-competitive ERK1/2 inhibitor, Vertex-11e, has properties of conformation selection for the R-state, revealing movements of the activation loop that are allosterically coupled to the kinase active site. However, the features of inhibitors important for R-state selection are unknown. Here, we survey a panel of ATP-competitive ERK inhibitors using HDX-MS and NMR and identify 14 new molecules with properties of R-state selection. They reveal effects propagated to distal regions in the P+1 and helix αF segments surrounding the activation loop, as well as helix αL16. Crystal structures of inhibitor complexes with ERK2 reveal systematic shifts in the Gly loop and helix αC, mediated by a Tyr-Tyr ring stacking interaction and the conserved Lys-Glu salt bridge. The findings suggest a model for the R-state involving small movements in the N-lobe that promote compactness within the kinase active site and alter mobility surrounding the activation loop. Such properties of conformation selection might be exploited to modulate the protein docking interface used by ERK substrates and effectors.

Discovery of lupus nephritis targeted inhibitors based on De novo molecular design: comprehensive application of vinardo scoring, ADMET analysis, and molecular dynamics simulation

Lupus Nephritis (LN) is an autoimmune disease affecting the kidneys, and conventional drug studies have limitations due to its imprecise and complex pathogenesis. Therefore, the aim of this study was to design a novel Lupus Nephritis-targeted drug with good clinical due potential, high potency and selectivity by computer-assisted approach.NIK belongs to the serine/threonine protein kinase, which is gaining attention as a drug target for Lupus Nephritis. we used bioinformatics, homology modelling and sequence comparison analysis, small molecule ab initio design, ADMET analysis, molecular docking, molecular dynamics simulation, and MM/PBSA analysis to design and explore the selectivity and efficiency of a novel Lupus Nephritis-targeting drug, ClImYnib, and a classical NIK inhibitor, NIK SMI1. We used bioinformatics techniques to determine the correlation between lupus nephritis and the NF-κB signaling pathway. De novo drugs design was used to create a NIK-targeted inhibitor, ClImYnib, with lower toxicity, after which we used molecular dynamics to simulate NIK SMI1 against ClImYnib, and the simulation results showed that ClImYnib had better selectivity and efficiency. Our research delves into the molecular mechanism of protein ligands, and we have designed and validated an excellent NIK inhibitor using multiple computational simulation methods. More importantly, it provides an idea of target designing small molecules.Communicated by Ramaswamy H. Sarma.

Profiling of drug resistance in Src kinase at scale uncovers a regulatory network coupling autoinhibition and catalytic domain dynamics

Kinase inhibitors are effective cancer therapies, but resistance often limits clinical efficacy. Despite the cataloging of numerous resistance mutations, our understanding of kinase inhibitor resistance is still incomplete. Here, we comprehensively profiled the resistance of ∼3,500 Src tyrosine kinase mutants to four different ATP-competitive inhibitors. We found that ATP-competitive inhibitor resistance mutations are distributed throughout Src's catalytic domain. In addition to inhibitor contact residues, residues that participate in regulating Src's phosphotransferase activity were prone to the development of resistance. Unexpectedly, we found that a resistance-prone cluster of residues located on the top face of the N-terminal lobe of Src's catalytic domain contributes to autoinhibition by reducing catalytic domain dynamics, and mutations in this cluster led to resistance by lowering inhibitor affinity and promoting kinase hyperactivation. Together, our studies demonstrate how drug resistance profiling can be used to define potential resistance pathways and uncover new mechanisms of kinase regulation.

Conformation Selection by ATP-competitive Inhibitors and Allosteric Communication in ERK2

Activation of the extracellular signal regulated kinase-2 (ERK2) by phosphorylation has been shown to involve changes in protein dynamics, as determined by hydrogen-deuterium exchange mass spectrometry (HDX-MS) and NMR relaxation dispersion measurements. These can be described by a global exchange between two conformational states of the active kinase, named "L" and "R", where R is associated with a catalytically productive ATP-binding mode. An ATP-competitive ERK1/2 inhibitor, Vertex-11e, has properties of conformation selection for the R-state, revealing movements of the activation loop that are allosterically coupled to the kinase active site. However, the features of inhibitors important for R-state selection are unknown. Here we survey a panel of ATP-competitive ERK inhibitors using HDX-MS and NMR and identify 14 new molecules with properties of R-state selection. They reveal effects propagated to distal regions in the P+1 and helix αF segments surrounding the activation loop, as well as helix αL16. Crystal structures of inhibitor complexes with ERK2 reveal systematic shifts in the Gly loop and helix αC, mediated by a Tyr-Tyr ring stacking interaction and the conserved Lys-Glu salt bridge. The findings suggest a model for the R-state involving small movements in the N-lobe that promote compactness within the kinase active site and alter mobility surrounding the activation loop. Such properties of conformation selection might be exploited to modulate the protein docking interface used by ERK substrates and effectors.

The Importance of the Pyrazole Scaffold in the Design of Protein Kinases Inhibitors as Targeted Anticancer Therapies

The altered activation or overexpression of protein kinases (PKs) is a major subject of research in oncology and their inhibition using small molecules, protein kinases inhibitors (PKI) is the best available option for the cure of cancer. The pyrazole ring is extensively employed in the field of medicinal chemistry and drug development strategies, playing a vital role as a fundamental framework in the structure of various PKIs. This scaffold holds major importance and is considered a privileged structure based on its synthetic accessibility, drug-like properties, and its versatile bioisosteric replacement function. It has proven to play a key role in many PKI, such as the inhibitors of Akt, Aurora kinases, MAPK, B-raf, JAK, Bcr-Abl, c-Met, PDGFR, FGFRT, and RET. Of the 74 small molecule PKI approved by the US FDA, 8 contain a pyrazole ring: Avapritinib, Asciminib, Crizotinib, Encorafenib, Erdafitinib, Pralsetinib, Pirtobrutinib, and Ruxolitinib. The focus of this review is on the importance of the unfused pyrazole ring within the clinically tested PKI and on the additional required elements of their chemical structures. Related important pyrazole fused scaffolds like indazole, pyrrolo[1,2-b]pyrazole, pyrazolo[4,3-b]pyridine, pyrazolo[1,5-a]pyrimidine, or pyrazolo[3,4-d]pyrimidine are beyond the subject of this work.

Molecular determinants of Hsp90 dependence of Src kinase revealed by deep mutational scanning

Hsp90 is a molecular chaperone involved in the refolding and activation of numerous protein substrates referred to as clients. While the molecular determinants of Hsp90 client specificity are poorly understood and limited to a handful of client proteins, strong clients are thought to be destabilized and conformationally extended. Here, we measured the phosphotransferase activity of 3929 variants of the tyrosine kinase Src in both the presence and absence of an Hsp90 inhibitor. We identified 84 previously unknown functionally dependent client variants. Unexpectedly, many destabilized or extended variants were not functionally dependent on Hsp90. Instead, functionally dependent client variants were clustered in the αF pocket and β1-β2 strand regions of Src, which have yet to be described in driving Hsp90 dependence. Hsp90 dependence was also strongly correlated with kinase activity. We found that a combination of activation, global extension, and general conformational flexibility, primarily induced by variants at the αF pocket and β1-β2 strands, was necessary to render Src functionally dependent on Hsp90. Moreover, the degree of activation and flexibility required to transform Src into a functionally dependent client varied with variant location, suggesting that a combination of regulatory domain disengagement and catalytic domain flexibility are required for chaperone dependence. Thus, by studying the chaperone dependence of a massive number of variants, we highlight factors driving Hsp90 client specificity and propose a model of chaperone-kinase interactions.

Pharmacological and In Silico Analysis of Oat Avenanthramides as EGFR Inhibitors: Effects on EGF-Induced Lung Cancer Cell Growth and Migration

Avena sativa L. is a wholegrain cereal and an important edible crop. Oats possesses high nutritional and health promoting values and contains high levels of bioactive compounds, including a group of phenolic amides, named avenanthramides (Avns), exerting antioxidant, anti-inflammatory, and anticancer activities. Epidermal growth factor receptor (EGFR) represents one of the most known oncogenes and it is frequently up-regulated or mutated in human cancers. The oncogenic effects of EGFR include enhanced cell growth, angiogenesis, and metastasis, and down-regulation or inhibition of EGFR signaling has therapeutic benefit. Front-line EGFR tyrosine kinase inhibitor therapy is the standard therapy for patients with EGFR-mutated lung cancer. However, the clinical effects of EGFR inhibition may be lost after a few months of treatment due to the onset of resistance. Here, we showed the anticancer activity of Avns, focusing on EGFR activation and signaling pathway. Lung cancer cellular models have been used to evaluate the activity of Avns on tumor growth, migration, EMT, and anoikis induced by EGF. In addition, docking and molecular dynamics simulations showed that the Avns bind with high affinity to a region in the vicinity of αC-helix and the DGF motif of EGFR, jeopardizing the target biological function. Altogether, our results reveal a new pharmacological activity of Avns as EGFR tyrosine kinase inhibitors.

Conformation-tunable ATP-competitive kinase inhibitors

Small molecule kinase inhibitors have shown immense clinical utility for diverse indications. While >60 kinase inhibitors have been approved (and many more in clinical trials), it remains unclear whether the clinical efficacy of a kinase inhibitor is solely dependent on enzymatic inhibition, or whether non-catalytic functions play a role in the efficacy of some kinase inhibitors. Here, we designed and synthesized a series of pyrazolopyrimidine kinase inhibitors that modulate the global kinase conformation of c-Src kinase. Expanding upon our findings from the pyrazolopyrimidine inhibitor series, we designed, synthesized, and evaluated three pair of conformation-selective kinase inhibitors, each with a unique hinge-binding scaffold. We profiled each pair of kinase inhibitors across 468 kinases and identified 38 kinases that could be studied using these pair of conformation-selective inhibitors. We also explore the binding of conformation-selective kinase inhibitors to mutant kinases of EGFR, FLT3, and KIT. Together, these studies yield important insight into the design of conformation-tunable kinase inhibitors and provide a toolset of compounds to study the role of protein conformation on kinase signaling.

Pseudokinases repurpose flexibility signatures associated with the protein kinase fold for noncatalytic roles

The bilobal protein kinase-like fold in pseudokinases lack one or more catalytic residues, conserved in canonical protein kinases, and are considered enzymatically deficient. Tertiary structures of pseudokinases reveal that their loops topologically equivalent to activation segments of kinases adopt contracted configurations, which is typically extended in active conformation of kinases. Herein, anisotropic network model based normal mode analysis (NMA) was conducted on 51 active conformation structures of protein kinases and 26 crystal structures of pseudokinases. Our observations indicate that although backbone fluctuation profiles are similar for individual kinase-pseudokinase families, low intensity mean square fluctuations in pseudo-activation segment and other sub-structures impart rigidity to pseudokinases. Analyses of collective motions from functional modes reveal that pseudokinases, compared to active kinases, undergo distinct conformational transitions using the same structural fold. All-atom NMA of protein kinase-pseudokinase pairs from each family, sharing high amino acid sequence identities, yielded distinct community clusters, partitioned by residues exhibiting highly correlated fluctuations. It appears that atomic fluctuations from equivalent activation segments guide community membership and network topologies for respective kinase and pseudokinase. Our findings indicate that such adaptations in backbone and side-chain fluctuations render pseudokinases competent for catalysis-independent roles.

TYROSINE KINASES: COMPLEX MOLECULAR SYSTEMS CHALLENGING COMPUTATIONAL METHODOLOGIES

Classical molecular dynamics (MD) simulations based on atomic models play an increasingly important role in a wide range of applications in physics, biology, and chemistry. Nonetheless, generating genuine knowledge about biological systems using MD simulations remains challenging. Protein tyrosine kinases are important cellular signaling enzymes that regulate cell growth, proliferation, metabolism, differentiation, and migration. Due to the large conformational changes and long timescales involved in their function, these kinases present particularly challenging problems to modern computational and theoretical frameworks aimed at elucidating the dynamics of complex biomolecular systems. Markov state models have achieved limited success in tackling the broader conformational ensemble and biased methods are often employed to examine specific long timescale events. Recent advances in machine learning continue to push the limitations of current methodologies and provide notable improvements when integrated with the existing frameworks. A broad perspective is drawn from a critical review of recent studies.

Making the invisible visible: Toward structural characterization of allosteric states, interaction networks, and allosteric regulatory mechanisms in protein kinases

Despite the established view of protein kinases as dynamic and versatile allosteric regulatory machines, our knowledge of allosteric functional states, allosteric interaction networks, and the intrinsic folding energy landscapes is surprisingly limited. We discuss the latest developments in structural characterization of allosteric molecular events underlying protein kinase dynamics and functions using structural, biophysical, and computational biology approaches. The recent studies highlighted progress in making the invisible aspects of protein kinase 'life' visible, including the determination of hidden allosteric states and mapping of allosteric energy landscapes, discovery of new mechanisms underlying ligand-induced modulation of allosteric activity, evolutionary adaptation of kinase allostery, and characterization of allosteric interaction networks as the intrinsic driver of kinase adaptability and signal transmission in the regulatory assemblies.

Reining in BTK: Interdomain Interactions and Their Importance in the Regulatory Control of BTK

Since Dr. Ogden Bruton's 1952 paper describing the first human primary immunodeficiency disease, the peripheral membrane binding signaling protein, aptly named Bruton's tyrosine kinase (BTK), has been the target of intense study. Dr. Bruton's description of agammaglobulinemia set the stage for ultimately understanding key signaling steps emanating from the B cell receptor. BTK is a multidomain tyrosine kinase and in the decades since Dr. Bruton's discovery it has become clear that genetic defects in the regulatory domains or the catalytic domain can lead to immunodeficiency. This finding underscores the intricate regulatory mechanisms within the BTK protein that maintain appropriate levels of signaling both in the resting B cell and during an immune challenge. In recent decades, BTK has become a target for clinical intervention in treating B cell malignancies. The survival reliance of B cell malignancies on B cell receptor signaling has allowed small molecules that target BTK to become essential tools in treating patients with hematological malignancies. The first-in-class Ibrutinib and more selective second-generation inhibitors all target the active site of the multidomain BTK protein. Therapeutic interventions targeting BTK have been successful but are plagued by resistance mutations that render drug treatment ineffective for some patients. This review will examine the molecular mechanisms that drive drug resistance, the long-range conformational effects of active site inhibitors on the BTK regulatory apparatus, and emerging opportunities to allosterically target the BTK kinase to improve therapeutic interventions using combination therapies.

Differential impact of BTK active site inhibitors on the conformational state of full-length BTK

Bruton’s tyrosine kinase (BTK) is targeted in the treatment of B-cell disorders including leukemias and lymphomas. Currently approved BTK inhibitors, including Ibrutinib, a first-in-class covalent inhibitor of BTK, bind directly to the kinase active site. While effective at blocking the catalytic activity of BTK, consequences of drug binding on the global conformation of full-length BTK are unknown. Here, we uncover a range of conformational effects in full-length BTK induced by a panel of active site inhibitors, including large-scale shifts in the conformational equilibria of the regulatory domains. Additionally, we find that a remote Ibrutinib resistance mutation, T316A in the BTK SH2 domain, drives spurious BTK activity by destabilizing the compact autoinhibitory conformation of full-length BTK, shifting the conformational ensemble away from the autoinhibited form. Future development of BTK inhibitors will need to consider long-range allosteric consequences of inhibitor binding, including the emerging application of these BTK inhibitors in treating COVID-19.

Treatments for blood cancers, such as leukemia and lymphoma, rely heavily on chemotherapy, using drugs that target a vulnerable aspect of the cancer cells. B-cells, a type of white blood cell that produces antibodies, require a protein called Bruton’s tyrosine kinase, or BTK for short, to survive. The drug ibrutinib (Imbruvica) is used to treat B-cell cancers by blocking BTK. The BTK protein consists of several regions. One of them, known as the kinase domain, is responsible for its activity as an enzyme (which allows it to modify other proteins by adding a ‘tag’ known as a phosphate group). The other regions of BTK, known as regulatory modules, control this activity. In BTK’s inactive form, the regulatory modules attach to the kinase domain, blocking the regulatory modules from interacting with other proteins. When BTK is activated, it changes its conformation so the regulatory regions detach and become available for interactions with other proteins, at the same time exposing the active kinase domain.

Ibrutinib and other BTK drugs in development bind to the kinase domain to block its activity. However, it is not known how this binding affects the regulatory modules. Previous efforts to study how drugs bind to BTK have used a version of the protein that only had the kinase domain, instead of the full-length protein.

Now, Joseph et al. have studied full-length BTK and how it binds to five different drugs. The results reveal that ibrutinib and another drug called dasatinib both indirectly disrupt the normal position of the regulatory domains pushing BTK toward a conformation that resembles the activated state. By contrast, the three other compounds studied do not affect the inactive structure. Joseph et al. also examined a mutation in BTK that confers resistance against ibrutinib. This mutation increases the activity of BTK by disrupting the inactive structure, leading to B cells surviving better.

Understanding how drug resistance mechanisms can work will lead to better drug treatment strategies for cancer. BTK is also a target in other diseases such as allergies or asthma and even COVID-19. If interactions between partner proteins and the regulatory domain are important in these diseases, then they may be better treated with drugs that maintain the regulatory modules in their inactive state. This research will help to design drugs that are better able to control BTK activity.

Parallel Chemoselective Profiling for Mapping Protein Structure

Solution-based structural techniques complement high-resolution structural data by providing insight into the oft-missed links between protein structure and dynamics. Here, we present Parallel Chemoselective Profiling, a solution-based structural method for characterizing protein structure and dynamics. Our method utilizes deep mutational scanning saturation mutagenesis data to install amino acid residues with specific chemistries at defined positions on the solvent-exposed surface of a protein. Differences in the extent of labeling of installed mutant residues are quantified using targeted mass spectrometry, reporting on each residue's local environment and structural dynamics. Using our method, we studied how conformation-selective, ATP-competitive inhibitors affect the local and global structure and dynamics of full-length Src kinase. Our results highlight how parallel chemoselective profiling can be used to study a dynamic multi-domain protein, and suggest that our method will be a useful addition to the relatively small toolkit of existing protein footprinting techniques.