Structure-kinetic relationship study of CDK8/CycC specific compounds

Pseudoirreversible inhibition elicits persistent efficacy of a sphingosine 1-phosphate receptor 1 antagonist

Sphingosine 1-phosphate receptor 1 (S1PR1), a G protein-coupled receptor, is required for lymphocyte trafficking, and is a promising therapeutic target in inflammatory diseases. Here, we synthesize a competitive S1PR1 antagonist, KSI-6666, that effectively suppresses pathogenic inflammation. Metadynamics simulations suggest that the interaction of KSI-6666 with a methionine residue Met124 in the ligand-binding pocket of S1PR1 may inhibit the dissociation of KSI-6666 from S1PR1. Consistently, in vitro functional and mutational analyses reveal that KSI-6666 causes pseudoirreversible inhibition of S1PR1, dependent on the Met124 of the protein and substituents on the distal benzene ring of KSI-6666. Moreover, in vivo study suggests that this pseudoirreversible inhibition is responsible for the persistent activity of KSI-6666.

Structural basis of the human transcriptional Mediator complex modulated by its dissociable Kinase module

The eukaryotic Mediator, comprising a large Core (cMED) and a dissociable CDK8 kinase module (CKM), regulates RNA Polymerase II (Pol II)-dependent transcription. cMED recruits Pol II and promotes pre-initiation complex (PIC) formation in a manner inhibited by the CKM, which is also implicated in post-initiation control of gene expression. Herein we report cryo-electron microscopy structures of the human complete Mediator and its CKM, which explains the basis for CKM inhibition of cMED-activated transcription. The CKM binds to cMED through an intrinsically disordered region (IDR) in MED13 and HEAT repeats in MED12. The CKM inhibits transcription by allocating its MED13 IDR to occlude binding of Pol II and MED26 to cMED and further obstructing cMED-PIC assembly through steric hindrance with TFIIH and the +1 nucleosome. Notably, MED12 binds to the cMED Hook, positioning CDK8 downstream of the transcription start site, which sheds new light on its stimulatory function in post-initiation events.

Highlighting the Major Role of Cyclin C in Cyclin-Dependent Kinase 8 Activity through Molecular Dynamics Simulations

Dysregulation of cyclin-dependent kinase 8 (CDK8) activity has been associated with many diseases, including colorectal and breast cancer. As usual in the CDK family, the activity of CDK8 is controlled by a regulatory protein called cyclin C (CycC). But, while human CDK family members are generally activated in two steps, that is, the binding of the cyclin to CDK and the phosphorylation of a residue in the CDK activation loop, CDK8 does not require the phosphorylation step to be active. Another peculiarity of CDK8 is its ability to be associated with CycC while adopting an inactive form. These specificities raise the question of the role of CycC in the complex CDK8-CycC, which appears to be more complex than the other members of the CDK family. Through molecular dynamics (MD) simulations and binding free energy calculations, we investigated the effect of CycC on the structure and dynamics of CDK8. In a second step, we particularly focused our investigation on the structural and molecular basis of the protein-protein interaction between the two partners by finely analyzing the energetic contribution of residues and simulating the transition between the active and the inactive form. We found that CycC has a stabilizing effect on CDK8, and we identified specific interaction hotspots within its interaction surface compared to other human CDK/Cyc pairs. Targeting these specific interaction hotspots could be a promising approach in terms of specificity to effectively disrupt the interaction between CDK8. The simulation of the conformational transition from the inactive to the active form of CDK8 suggests that the residue Glu99 of CycC is involved in the orientation of three conserved arginines of CDK8. Thus, this residue may assume the role of the missing phosphorylation step in the activation mechanism of CDK8. In a more general view, these results point to the importance of keeping the CycC in computational studies when studying the human CDK8 protein in both the active and the inactive form.

Hit discovery of potential CDK8 inhibitors and analysis of amino acid mutations for cancer therapy through computer-aided drug discovery

Cyclin-dependent kinase 8 (CDK8) has emerged as a promising target for inhibiting cancer cell function, intensifying efforts towards the development of CDK8 inhibitors as potential cancer therapeutics. Mutations in CDK8, a protein kinase, are also implicated as a primary factor associated with tumor formation. In this study, we identified potential inhibitors through virtual screening for CDK8 and single amino acid mutations in CDK8, namely D173A (Aspartate 173 mutate to Alanine), D189N (Aspartate 189 mutate to Asparagine), T196A (Threonine 196 mutate to Alanine) and T196D (Threonine 196 mutate to Aspartate). Four databases (CHEMBEL, ZINC, MCULE, and MolPort) containing 65,209,131 molecules have been searched to identify new inhibitors for CDK8 and its single mutations. In the first step, structure-based pharmacophore modeling in the Pharmit server was used to select the compounds to know the inhibitors. Then molecules with better predicted drug-like molecule properties were selected. The final filter used to select more effective inhibitors among the previously selected molecules was molecular docking. Finally, 13 hits for CDK8, 11 hits for D173A, 11 hits for D189N, 15 hits for T196A, and 12 hits for T196D were considered potential inhibitors. A majority of the virtual screening hits exhibited satisfactorily predict pharmacokinetic characteristics and toxicity properties.

Discovery of Potent Antimalarial Type II Kinase Inhibitors with Selectivity over Human Kinases

While progress has been made in the effort to eradicate malaria, the disease remains a significant threat to global health. Acquired resistance to frontline treatments is emerging in Africa, urging a need for the development of novel antimalarial agents. Repurposing human kinase inhibitors provides a potential expedited route given the availability of a diverse array of kinase-targeting drugs that are approved or in clinical trials. Phenotypic screening of a library of type II human kinase inhibitors identified compound 1 as a lead antimalarial, which was initially developed to target human ephrin type A receptor 2 (EphA2). Here, we report a structure-activity relationship study and lead optimization of compound 1, which led to compound 33, with improved antimalarial activity and selectivity.

Unveiling the noncanonical activation mechanism of CDKs: insights from recent structural studies

The Cyclin-dependent kinases (CDKs) play crucial roles in a range of essential cellular processes. While the classical two-step activation mechanism is generally applicable to cell cycle-related CDKs, both CDK7 and CDK8, involved in transcriptional regulation, adopt distinct mechanisms for kinase activation. In both cases, binding to their respective cyclin partners results in only partial activity, while their full activation requires the presence of an additional subunit. Recent structural studies of these two noncanonical kinases have provided unprecedented insights into their activation mechanisms, enabling us to understand how the third subunit coordinates the T-loop stabilization and enhances kinase activity. In this review, we summarize the structure and function of CDK7 and CDK8 within their respective functional complexes, while also describing their noncanonical activation mechanisms. These insights open new avenues for targeted drug discovery and potential therapeutic interventions in various diseases related to CDK7 and CDK8.

Estimation of Drug-Target Residence Time by Targeted Molecular Dynamics Simulations

Drug-target residence time has emerged as a key selection factor in drug discovery since the binding duration of a drug molecule to its protein target can significantly impact its in vivo efficacy. The challenge in studying the residence time, in early drug discovery stages, lies in how to cost-effectively determine the residence time for the systematic assessment of compounds. Currently, there is still a lack of computational protocols to quickly estimate such a measure, particularly for large and flexible protein targets and drugs. Here, we report an efficient computational protocol, based on targeted molecular dynamics, to rank drug candidates by their residence time and to obtain insights into ligand-target dissociation mechanisms. The method was assessed on a dataset of 10 arylpyrazole inhibitors of CDK8, a large, flexible, and clinically important target, for which the experimental residence time of the inhibitors ranges from minutes to hours. The compounds were correctly ranked according to their estimated residence time scores compared to their experimental values. The analysis of protein-ligand interactions along the dissociation trajectories highlighted the favorable contribution of hydrophobic contacts to residence time and revealed key residues that strongly affect compound residence time.

Structure-Kinetic Relationship Studies for the Development of Long Residence Time LpxC Inhibitors

UDP-3-O-(R-3-hydroxymyristoyl)-N-acetylglucosamine deacetylase (LpxC) is a promising drug target in Gram-negative bacteria. Previously, we described a correlation between the residence time of inhibitors on Pseudomonas aeruginosa LpxC (paLpxC) and the post-antibiotic effect (PAE) caused by the inhibitors on the growth of P. aeruginosa. Given that drugs with prolonged activity following compound removal may have advantages in dosing regimens, we have explored the structure-kinetic relationship for paLpxC inhibition by analogues of the pyridone methylsulfone PF5081090 (1) originally developed by Pfizer. Several analogues have longer residence times on paLpxC than 1 (41 min) including PT913, which has a residence time of 124 min. PT913 also has a PAE of 4 h, extending the original correlation observed between residence time and PAE. Collectively, the studies provide a platform for the rational modulation of paLpxC inhibitor residence time and the potential development of antibacterial agents that cause prolonged suppression of bacterial growth.

Quantitative Structure-Activity Relationship (QSAR) Study Predicts Small-Molecule Binding to RNA Structure

The diversity of RNA structural elements and their documented role in human diseases make RNA an attractive therapeutic target. However, progress in drug discovery and development has been hindered by challenges in the determination of high-resolution RNA structures and a limited understanding of the parameters that drive RNA recognition by small molecules, including a lack of validated quantitative structure-activity relationships (QSARs). Herein, we develop QSAR models that quantitatively predict both thermodynamic- and kinetic-based binding parameters of small molecules and the HIV-1 transactivation response (TAR) RNA model system. Small molecules bearing diverse scaffolds were screened against TAR using surface plasmon resonance. Multiple linear regression (MLR) combined with feature selection afforded robust models that allowed direct interpretation of the properties critical for both binding strength and kinetic rate constants. These models were validated with new molecules, and their accurate performance was confirmed via comparison to ensemble tree methods, supporting the general applicability of this platform.

Lessons Learned from Past Cyclin-Dependent Kinase Drug Discovery Efforts

Inhibition of cyclin-dependent kinases (CDKs) has become an effective therapeutic strategy for treating various diseases, especially cancer. Over almost three decades, although great efforts have been made to discover CDK inhibitors, many of which have entered clinical trials, only four CDK inhibitors have been approved. In the process of CDK inhibitor development, many difficulties and misunderstandings have hampered their discovery and clinical applications, which mainly include inadequate understanding of the biological functions of CDKs, less attention paid to pan- and multi-CDK inhibitors, nonideal isoform selectivity of developed selective CDK inhibitors, overlooking the metabolic stability of early discovered CDK inhibitors, no effective resistance solutions, and a lack of available combination therapy and effective biomarkers for CDK therapies. After reviewing the mechanisms of CDKs and the research progress of CDK inhibitors, this perspective summarizes and discusses these difficulties or lessons, hoping to facilitate the successful discovery of more useful CDK inhibitors.

A Selective and Orally Bioavailable Quinoline-6-Carbonitrile-Based Inhibitor of CDK8/19 Mediator Kinase with Tumor-Enriched Pharmacokinetics

Senexins are potent and selective quinazoline inhibitors of CDK8/19 Mediator kinases. To improve their potency and metabolic stability, quinoline-based derivatives were designed through a structure-guided strategy based on the simulated drug-target docking model of Senexin A and Senexin B. A library of quinoline-Senexin derivatives was synthesized to explore the structure-activity relationship (SAR). An optimized compound 20a (Senexin C) exhibits potent CDK8/19 inhibitory activity with high selectivity. Senexin C is more metabolically stable and provides a more sustained inhibition of CDK8/19-dependent cellular gene expression when compared with the prototype inhibitor Senexin B. In vivo pharmacokinetic (PK) and pharmacodynamic (PD) evaluation using a novel tumor-based PD assay showed good oral bioavailability of Senexin C with a strong tumor-enrichment PK profile and tumor-PD marker responses. Senexin C inhibits MV4-11 leukemia growth in a systemic in vivo model with good tolerability.

Computational studies to identify the common type-I and type-II inhibitors against the CDK8 enzyme

In this study, multicomplex-based pharmacophore modeling was conducted on the structural proteome of the two states of CDK8 protein, that is, DMG-in and out. Three pharmacophores having six, five, and four features were selected as the representative models to conduct the virtual screening process using the prepared drug-like natural product database. The screened candidates were subjected to molecular docking studies on DMG-in (5XS2) and out (4F6U) conformation of the CDK8 protein. Subsequently, the common four docked candidates of 5XS2 and 4F6U were selected to perform the molecular dynamics simulation studies. Apart from one of the complexes of DMG-in (5XS2-UNPD163102), all other complexes displayed stable dynamic behavior. The interaction and stability studies of the docked complexes were compared with the references selected from the two conformations (DMG-in and out) of the protein. The current work leads to the identification of three common DMG-in and out hits with diverse scaffolds which can be employed as the initial leads for the design of the novel CDK8 inhibitors.

Unraveling binding mechanism and kinetics of macrocyclic Gαq protein inhibitors

G proteins represent intracellular switches that transduce signals relayed from G protein-coupled receptors. The structurally related macrocyclic depsipeptides FR900359 (FR) and YM-254890 (YM) are potent, selective inhibitors of the Gα_q protein family. We recently discovered that radiolabeled FR and YM display strongly divergent residence times, which translates into significantly longer antiasthmatic effects of FR. The present study is aimed at investigating the molecular basis for this observed disparity. Based on docking studies, we mutated amino acid residues of the Gα_q protein predicted to interact with FR or YM, and recombinantly expressed the mutated Gα_q proteins in cells in which the native Gα_q proteins had been knocked out by CRISPR-Cas9. Both radioligands showed similar association kinetics, and their binding followed a conformational selection mechanism, which was rationalized by molecular dynamics simulation studies. Several mutations of amino acid residues near the putative binding site of the "lipophilic anchors" of FR, especially those predicted to interact with the isopropyl group present in FR but not in YM, led to dramatically accelerated dissociation kinetics. Our data indicate that the long residence time of FR depends on lipophilic interactions within its binding site. The observed structure-kinetic relationships point to a complex binding mechanism of FR, which likely involves snap-lock- or dowel-like conformational changes of either ligand or protein, or both. These experimental data will be useful for the design of compounds with a desired residence time, a parameter that has now been recognized to be of utmost importance in drug development.

SAR study of small molecule inhibitors of the programmed cell death-1/programmed cell death-ligand 1 interaction

The development of small molecule inhibitors of programmed cell death-1/programmed cell death-ligand 1 (PD-1/PD-L1) has drawn research interest for the treatment of cancer. Recently, we reported the discovery of a novel dimeric core small molecule PD-1/PD-L1 inhibitor. In an effort to discover more potent inhibitors, we further explored the dimeric core scaffold. Our investigations of the structure-activity-relationship revealed that introduction of lipophilic substituents onto one of the di-alkoxylated phenyl rings improved binding affinities to PD-L1, and inhibitory activities of PD-1/PD-L1 in cellular assays. Furthermore, conversion of the ether linker part to an olefin linker not only improved binding affinity but also led to slow dissociation binding kinetics. We also explored more potent, as well as downsized, scaffolds. Compounds bearing a linear chain in place of one of the di-alkoxylated phenyl rings exhibited good binding affinity with improved ligand efficiency (LE). Representative compounds demonstrated potent inhibitory activities of PD-1/PD-L1 in the submicromolar range in cellular assays as well as cellular function in the mixed lymphocyte reaction (MLR) assay with efficacy comparable to anti-PD-1 antibody. Our results provide applicable information for the design of more potent inhibitors targeting PD-1/PD-L1 pathway.

Novel irreversible covalent BTK inhibitors discovered using DNA-encoded chemistry

Libraries of DNA-Encoded small molecules created using combinatorial chemistry and synthetic oligonucleotides are being applied to drug discovery projects across the pharmaceutical industry. The majority of reported projects describe the discovery of reversible, i.e. non-covalent, target modulators. We synthesized multiple DNA-encoded chemical libraries terminated in electrophiles and then used them to discover covalent irreversible inhibitors and report the successful discovery of acrylamide- and epoxide-terminated Bruton's Tyrosine Kinase (BTK) inhibitors. We also demonstrate their selectivity, potency and covalent cysteine engagement using a range of techniques including X-ray crystallography, thermal transition shift assay, reporter displacement assay and intact protein complex mass spectrometry. The epoxide BTK inhibitors described here are the first ever reported to utilize this electrophile for this target.

Potent and orally bioavailable CDK8 inhibitors: Design, synthesis, structure-activity relationship analysis and biological evaluation

CDK8 regulates transcription either by phosphorylation of transcription factors or, as part of a four-subunit kinase module, through a reversible association of the kinase module with the Mediator complex, a highly conserved transcriptional coactivator. Deregulation of CDK8 has been found in various types of human cancer, while the role of CDK8 in supressing anti-cancer response of natural killer cells is being understood. Currently, CDK8-targeting cancer drugs are highly sought-after. Herein we detail the discovery of a series of novel pyridine-derived CDK8 inhibitors. Medicinal chemistry optimisation gave rise to 38 (AU1-100), a potent CDK8 inhibitor with oral bioavailability. The compound inhibited the proliferation of MV4-11 acute myeloid leukaemia cells with the kinase activity of cellular CDK8 dampened. No systemic toxicology was observed in the mice treated with 38. These results warrant further pre-clinical studies of 38 as an anti-cancer agent.

Structure and noncanonical Cdk8 activation mechanism within an Argonaute-containing Mediator kinase module

The Cdk8 kinase module (CKM) in Mediator, comprising Med13, Med12, CycC, and Cdk8, regulates RNA polymerase II transcription through kinase-dependent and -independent functions. Numerous pathogenic mutations causative for neurodevelopmental disorders and cancer congregate in CKM subunits. However, the structure of the intact CKM and the mechanism by which Cdk8 is non-canonically activated and functionally affected by oncogenic CKM alterations are poorly understood. Here, we report a cryo-electron microscopy structure of Saccharomyces cerevisiae CKM that redefines prior CKM structural models and explains the mechanism of Med12-dependent Cdk8 activation. Med12 interacts extensively with CycC and activates Cdk8 by stabilizing its activation (T-)loop through conserved Med12 residues recurrently mutated in human tumors. Unexpectedly, Med13 has a characteristic Argonaute-like bi-lobal architecture. These findings not only provide a structural basis for understanding CKM function and pathological dysfunction, but also further impute a previously unknown regulatory mechanism of Mediator in transcriptional modulation through its Med13 Argonaute-like features.

Mediator's Kinase Module: A Modular Regulator of Cell Fate

Selective gene expression is crucial in maintaining the self-renewing and multipotent properties of stem cells. Mediator is a large, evolutionarily conserved, multi-subunit protein complex that modulates gene expression by relaying signals from cell type-specific transcription factors to RNA polymerase II. In humans, this complex consists of 30 subunits arranged in four modules. One critical module of the Mediator complex is the kinase module consisting of four subunits: MED12, MED13, CDK8, and CCNC. The kinase module exists in variable association with the 26-subunit Mediator core and affects transcription through phosphorylation of transcription factors and by controlling Mediator structure and function. Many studies have shown the kinase module to be a key player in the maintenance of stem cells that is distinct from a general role in transcription. Genetic studies have revealed that dysregulation of this kinase subunit contributes to the development of many human diseases. In this review, we discuss the importance of the Mediator kinase module by examining how this module functions with the more recently identified transcriptional super-enhancers, how changes in the kinase module and its activity can lead to the development of human disease, and the role of this unique module in directing and maintaining cell state. As we look to use stem cells to understand human development and treat human disease through both cell-based therapies and tissue engineering, we need to remain aware of the on-going research and address critical gaps in knowledge related to the molecular mechanisms that control cell fate.

Angel or Devil ? - CDK8 as the new drug target

Cyclin-dependent kinase 8 (CDK8) plays an momentous role in transcription regulation by forming kinase module or transcription factor phosphorylation. A large number of evidences have identified CDK8 as an important factor in cancer occurrence and development. In addition, CDK8 also participates in the regulation of cancer cell stress response to radiotherapy and chemotherapy, assists tumor cell invasion, metastasis, and drug resistance. Therefore, CDK8 is regarded as a promising target for cancer therapy. Most studies in recent years supported the role of CDK8 as a carcinogen, however, under certain conditions, CDK8 exists as a tumor suppressor. The functional diversity of CDK8 and its exceptional role in different types of cancer have aroused great interest from scientists but even more controversy during the discovery of CDK8 inhibitors. In addition, CDK8 appears to be an effective target for inflammation diseases and immune system disorders. Therefore, we summarized the research results of CDK8, involving physiological/pathogenic mechanisms and the development status of compounds targeting CDK8, provide a reference for the feasibility evaluation of CDK8 as a therapeutic target, and guidance for researchers who are involved in this field for the first time.

Pyrazolyl-Ureas as Interesting Scaffold in Medicinal Chemistry

The pyrazole nucleus has long been known as a privileged scaffold in the synthesis of biologically active compounds. Within the numerous pyrazole derivatives developed as potential drugs, this review is focused on molecules characterized by a urea function directly linked to the pyrazole nucleus in a different position. In the last 20 years, the interest of numerous researchers has been especially attracted by pyrazolyl-ureas showing a wide spectrum of biological activities, ranging from the antipathogenic activities (bacteria, plasmodium, toxoplasma, and others) to the anticarcinogenic activities. In particular, in the anticancer field, pyrazolyl-ureas have been shown to interact at the intracellular level on many pathways, in particular on different kinases such as Src, p38-MAPK, TrKa, and others. In addition, some of them evidenced an antiangiogenic potential that deserves to be explored. This review therefore summarizes all these biological data (from 2000 to date), including patented compounds.

Design, synthesis, and biological evaluation of new pyrazoloquinazoline derivatives as dual COX-2/5-LOX inhibitors

A new series of pyrazoloquinazoline derivatives equipped with different chalcones was designed, synthesized, and identified through ¹ H nuclear magnetic resonance (NMR), ¹³ C NMR, and infrared spectroscopic techniques. Our design strategy of the quinazolinone-privileged scaffold as a new scaffold was based on merging pharmacophores previously reported to exhibit cyclooxygenase-2 (COX-2)/5-lipoxygenase (5-LOX) inhibitory activity. All the newly synthesized derivatives were biologically evaluated for COX and 5-LOX inhibitory activity and COX-2 selectivity, using celecoxib and zileuton as reference drugs, as they exhibited promising anti-inflammatory activity. Compound 3j was found to be the most promising derivative, with IC₅₀ values of 667 and 47 nM against COX-1 and COX-2, respectively, which are superior to that of celecoxib (IC₅₀ value against COX-2 = 95 nM), showing an SI of 14.2 that was much better than celecoxib. Compounds 3f and 3h exhibited COX-1 inhibition, with IC₅₀ values of 1,485 and 684 nM, respectively. The synthesized compounds showed a significant inhibitory activity against 5-LOX, with IC₅₀ values ranging from 0.6 to 4.3 µM, where compounds 3f and 3h were found to be the most potent derivatives, with IC₅₀ values of 0.6 and 1.0 µM, respectively, in comparison with that of zileuton (IC₅₀  = 0.8 µM). These promising derivatives, 3f, 3h, and 3j, were further investigated in vivo for anti-inflammatory, gastric ulcerogenic effects, and prostaglandin production (PGE2) in rat serum. The molecular docking studies concerning the binding sites of COX-2 and 5-LOX revealed similar orientation, compared with reported inhibitors, which encouraged us to design new leads targeting COX-2 and 5-LOX as dual inhibitors, as a new avenue in anti-inflammatory therapy.

Synthesis and biological evaluation of small molecule modulators of CDK8/Cyclin C complex with phenylaminoquinoline scaffold

Background

CDK8/CycC complex has kinase activity towards the carboxyterminal domain of RNA polymerase II, and contributes to the regulation of transcription via association with the mediator complex. Different human malignancies, mainly colorectal and gastric cancers, were produced as a result of overexpression of CDK8/CycC in the mediator complex. Therefore, CDK8/CycC complex represents as a cancer oncogene and it has become a potential target for developing CDK8/CycC modulators.

Methods

A series of nine 4-phenylaminoquinoline scaffold-based compounds 5a-i was synthesized, and biologically evaluated as potential CDK8/CycC complex inhibitors.

Results

The scaffold substituent effects on the intrinsic inhibitory activity toward CDK8/CycC complex are addressed trying to present a novel outlook of CDK8/CycC Complex inhibitors with 4-phenylaminoquinoline scaffold in cancer therapy. The secondary benzenesulfonamide analogues proved to be the most potent compounds in suppressing CDK8/CycC enzyme, whereas, their primary benzenesulfonamide analogues showed inferior activity. Moreover, the benzene reversed sulfonamide analogues were totally inactive.

Discussion

The titled scaffold showed promising inhibitory activity data and there is a crucial role of un/substituted sulfonamido group for CDK8/CycC complex inhibitory activity. Compound 5d showed submicromolar potency against CDK8/CycC (IC₅₀ = 0.639 µM) and it can be used for further investigations and to design another larger library of phenylaminoquinoline scaffold-based analogues in order to establish detailed SARs.

Transient States and Barriers from Molecular Simulations and the Milestoning Theory: Kinetics in Ligand-Protein Recognition and Compound Design

This study presents a novel computational approach to study molecular recognition and binding kinetics for drug-like compounds dissociating from a flexible protein system. The intermediates and their free energy profile during ligand association and dissociation processes control ligand-protein binding kinetics and bring a more complete picture of ligand-protein binding. The method applied the milestoning theory to extract kinetics and thermodynamics information from running short classical molecular dynamics (MD) simulations for frames from a given dissociation path. High-dimensional ligand-protein motions (3N-6 degrees of freedom) during ligand dissociation were reduced by use of principal component modes for assigning more than 100 milestones, and classical MD runs were allowed to travel multiple milestones to efficiently obtain ensemble distribution of initial structures for MD simulations and estimate the transition time and rate during ligand traveling between milestones. We used five pyrazolourea ligands and cyclin-dependent kinase 8 with cyclin C (CDK8/CycC) as our model system as well as metadynamics and a pathway search method to sample dissociation pathways. With our strategy, we constructed the free energy profile for highly mobile biomolecular systems. The computed binding free energy and residence time correctly ranked the pyrazolourea ligand series, in agreement with experimental data. Guided by a barrier of a ligand passing an αC helix and activation loop, we introduced one hydroxyl group to parent compounds to design our ligands with increased residence time and validated our prediction by experiments. This work provides a novel and robust approach to investigate dissociation kinetics of large and flexible systems for understanding unbinding mechanisms and designing new small-molecule drugs with desired binding kinetics.

Selective and Potent CDK8/19 Inhibitors Enhance NK-Cell Activity and Promote Tumor Surveillance

Natural killer (NK) cells play a pivotal role in controlling cancer. Multiple extracellular receptors and internal signaling nodes tightly regulate NK activation. Cyclin-dependent kinases of the mediator complex (CDK8 and CDK19) were described as a signaling intermediates in NK cells. Here, we report for the first time the development and use of CDK8/19 inhibitors to suppress phosphorylation of STAT1^S727 in NK cells and to augment the production of the cytolytic molecules perforin and granzyme B (GZMB). Functionally, this resulted in enhanced NK-cell-mediated lysis of primary leukemia cells. Treatment with the CDK8/19 inhibitor BI-1347 increased the response rate and survival of mice bearing melanoma and breast cancer xenografts. In addition, CDK8/19 inhibition augmented the antitumoral activity of anti-PD-1 antibody and SMAC mimetic therapy, both agents that promote T-cell-mediated antitumor immunity. Treatment with the SMAC mimetic compound BI-8382 resulted in an increased number of NK cells infiltrating EMT6 tumors. Combination of the CDK8/19 inhibitor BI-1347, which augments the amount of degranulation enzymes, with the SMAC mimetic BI-8382 resulted in increased survival of mice carrying the EMT6 breast cancer model. The observed survival benefit was dependent on an intermittent treatment schedule of BI-1347, suggesting the importance of circumventing a hyporesponsive state of NK cells. These results suggest that CDK8/19 inhibitors can be combined with modulators of the adaptive immune system to inhibit the growth of solid tumors, independent of their activity on cancer cells, but rather through promoting NK-cell function.

A precisely positioned MED12 activation helix stimulates CDK8 kinase activity

The Mediator kinase module regulates eukaryotic transcription by phosphorylating transcription-related targets and by modulating the association of Mediator and RNA polymerase II. The activity of its catalytic core, cyclin-dependent kinase 8 (CDK8), is controlled by Cyclin C and regulatory subunit MED12, with its deregulation contributing to numerous malignancies. Here, we combine in vitro biochemistry, cross-linking coupled to mass spectrometry, and in vivo studies to describe the binding location of the N-terminal segment of MED12 on the CDK8/Cyclin C complex and to gain mechanistic insights into the activation of CDK8 by MED12. Our data demonstrate that the N-terminal portion of MED12 wraps around CDK8, whereby it positions an "activation helix" close to the T-loop of CDK8 for its activation. Intriguingly, mutations in the activation helix that are frequently found in cancers do not diminish the affinity of MED12 for CDK8, yet likely alter the exact positioning of the activation helix. Furthermore, we find the transcriptome-wide gene-expression changes in human cells that result from a mutation in the MED12 activation helix to correlate with deregulated genes in breast and colon cancer. Finally, functional assays in the presence of kinase inhibitors reveal that binding of MED12 remodels the active site of CDK8 and thereby precludes the inhibition of ternary CDK8 complexes by type II kinase inhibitors. Taken together, our results not only allow us to propose a revised model of how CDK8 activity is regulated by MED12, but also offer a path forward in developing small molecules that target CDK8 in its MED12-bound form.

[Ligand-receptor binding kinetics in drug design]

Traditionally, the thermodynamic values of affinity are considered as the main criterion for the development of new drugs. Usually, these values for drugs are measured <i>in vitro</i> at steady concentrations of the receptor and ligand, which are differed from <i>in vivo</i> environment. Recent studies have shown that the kinetics of the process of drug binding to its receptor make significant contribution in the drug effectiveness. This has increased attention in characterizing and predicting the rate constants of association and dissociation of the receptor ligand at the stage of preclinical studies of drug candidates. A drug with a long residence time can determine ligand-receptor selectivity (kinetic selectivity), maintain pharmacological activity of the drug at its low concentration in vivo. The paper discusses the theoretical basis of protein-ligand binding, molecular determinants that control the kinetics of the drug-receptor binding. Understanding the molecular features underlying the kinetics of receptor-ligand binding will contribute to the rational design of drugs with desired properties.

Accelerating Rare Dissociative Processes in Biomolecules Using Selectively Scaled MD Simulations

Molecular dynamics (MD) simulations can be a powerful tool for modeling complex dissociative processes such as ligand unbinding. However, many biologically relevant dissociative processes occur on timescales that far exceed the timescales of typical MD simulations. Here, we implement and apply an enhanced sampling method in which specific energy terms in the potential energy function are selectively "scaled" to accelerate dissociative events during simulations. Using ligand unbinding as an example of a complex dissociative process, we selectively scaled up ligand-water interactions in an attempt to increase the rate of ligand unbinding. After applying our selectively scaled MD (ssMD) approach to several cyclin-dependent kinase-inhibitor complexes, we discovered that we could accelerate ligand unbinding, thereby allowing, in some cases, unbinding events to occur within as little as 2 ns. Moreover, we found that we could make realistic estimates of the initial unbinding times (τ_unbind^sim) as well as the accompanying change in free energy (ΔG^sim) of the inhibitors from our ssMD simulation data. To accomplish this, we employed a previously described Kramers'-based rate extrapolation method and a newly described free energy extrapolation method. Because our ssMD approach is general, it should find utility as an easy-to-deploy, enhanced sampling method for modeling other dissociative processes.

Recent Progress towards Chemically-Specific Coarse-Grained Simulation Models with Consistent Dynamical Properties

Coarse-grained (CG) models can provide computationally efficient and conceptually simple characterizations of soft matter systems. While generic models probe the underlying physics governing an entire family of free-energy landscapes, bottom-up CG models are systematically constructed from a higher-resolution model to retain a high level of chemical specificity. The removal of degrees of freedom from the system modifies the relationship between the relative time scales of distinct dynamical processes through both a loss of friction and a “smoothing” of the free-energy landscape. While these effects typically result in faster dynamics, decreasing the computational expense of the model, they also obscure the connection to the true dynamics of the system. The lack of consistent dynamics is a serious limitation for CG models, which not only prevents quantitatively accurate predictions of dynamical observables but can also lead to qualitatively incorrect descriptions of the characteristic dynamical processes. With many methods available for optimizing the structural and thermodynamic properties of chemically-specific CG models, recent years have seen a stark increase in investigations addressing the accurate description of dynamical properties generated from CG simulations. In this review, we present an overview of these efforts, ranging from bottom-up parameterizations of generalized Langevin equations to refinements of the CG force field based on a Markov state modeling framework. We aim to make connections between seemingly disparate approaches, while laying out some of the major challenges as well as potential directions for future efforts.

De Novo Missense Substitutions in the Gene Encoding CDK8, a Regulator of the Mediator Complex, Cause a Syndromic Developmental Disorder

The Mediator is an evolutionarily conserved, multi-subunit complex that regulates multiple steps of transcription. Mediator activity is regulated by the reversible association of a four-subunit module comprising CDK8 or CDK19 kinases, together with cyclin C, MED12 or MED12L, and MED13 or MED13L. Mutations in MED12, MED13, and MED13L were previously identified in syndromic developmental disorders with overlapping phenotypes. Here, we report CDK8 mutations (located at 13q12.13) that cause a phenotypically related disorder. Using whole-exome or whole-genome sequencing, and by international collaboration, we identified eight different heterozygous missense CDK8 substitutions, including 10 shown to have arisen de novo, in 12 unrelated subjects; a recurrent mutation, c.185C>T (p.Ser62Leu), was present in five individuals. All predicted substitutions localize to the ATP-binding pocket of the kinase domain. Affected individuals have overlapping phenotypes characterized by hypotonia, mild to moderate intellectual disability, behavioral disorders, and variable facial dysmorphism. Congenital heart disease occurred in six subjects; additional features present in multiple individuals included agenesis of the corpus callosum, ano-rectal malformations, seizures, and hearing or visual impairments. To evaluate the functional impact of the mutations, we measured phosphorylation at STAT1-Ser727, a known CDK8 substrate, in a CDK8 and CDK19 CRISPR double-knockout cell line transfected with wild-type (WT) or mutant CDK8 constructs. These experiments demonstrated a reduction in STAT1 phosphorylation by all mutants, in most cases to a similar extent as in a kinase-dead control. We conclude that missense mutations in CDK8 cause a developmental disorder that has phenotypic similarity to syndromes associated with mutations in other subunits of the Mediator kinase module, indicating probable overlap in pathogenic mechanisms.

Structural insights into the functional diversity of the CDK-cyclin family

Since their characterization as conserved modules that regulate progression through the eukaryotic cell cycle, cyclin-dependent protein kinases (CDKs) in higher eukaryotic cells are now also emerging as significant regulators of transcription, metabolism and cell differentiation. The cyclins, though originally characterized as CDK partners, also have CDK-independent roles that include the regulation of DNA damage repair and transcriptional programmes that direct cell differentiation, apoptosis and metabolic flux. This review compares the structures of the members of the CDK and cyclin families determined by X-ray crystallography, and considers what mechanistic insights they provide to guide functional studies and distinguish CDK- and cyclin-specific activities. Aberrant CDK activity is a hallmark of a number of diseases, and structural studies can provide important insights to identify novel routes to therapy.

Cyclin C: The Story of a Non-Cycling Cyclin

The class I cyclin family is a well-studied group of structurally conserved proteins that interact with their associated cyclin-dependent kinases (Cdks) to regulate different stages of cell cycle progression depending on their oscillating expression levels. However, the role of class II cyclins, which primarily act as transcription factors and whose expression remains constant throughout the cell cycle, is less well understood. As a classic example of a transcriptional cyclin, cyclin C forms a regulatory sub-complex with its partner kinase Cdk8 and two accessory subunits Med12 and Med13 called the Cdk8-dependent kinase module (CKM). The CKM reversibly associates with the multi-subunit transcriptional coactivator complex, the Mediator, to modulate RNA polymerase II-dependent transcription. Apart from its transcriptional regulatory function, recent research has revealed a novel signaling role for cyclin C at the mitochondria. Upon oxidative stress, cyclin C leaves the nucleus and directly activates the guanosine 5’-triphosphatase (GTPase) Drp1, or Dnm1 in yeast, to induce mitochondrial fragmentation. Importantly, cyclin C-induced mitochondrial fission was found to increase sensitivity of both mammalian and yeast cells to apoptosis. Here, we review and discuss the biology of cyclin C, focusing mainly on its transcriptional and non-transcriptional roles in tumor promotion or suppression.

Recent development of CDK inhibitors: An overview of CDK/inhibitor co-crystal structures

The cyclin-dependent protein kinases (CDKs) are protein-serine/threonine kinases that display crucial effects in regulation of cell cycle and transcription. While the excessive expression of CDKs is intimate related to the development of diseases including cancers, which provides opportunities for disease treatment. A large number of small molecules are explored targeting CDKs. CDK/inhibitor co-crystal structures play an important role during the exploration of inhibitors. So far nine kinds of CDK/inhibitor co-crystals have been determined, they account for the highest proportion among the Protein Data Bank (PDB) deposited crystal structures. Herein, we review main co-crystals of CDKs in complex with corresponding inhibitors reported in recent years, focusing our attention on the binding models and the pharmacological activities of inhibitors.

CDK8 as a therapeutic target for cancers and recent developments in discovery of CDK8 inhibitors

Cyclin-dependent kinases 8 (CDK8) regulates transcriptional process via associating with the mediator complex or phosphorylating transcription factors (TF). Overexpression of CDK8 has been observed in various cancers. It mediates aberrant activation of Wnt/β-catenin signaling pathway, which is initially recognized and best studied in colorectal cancer (CRC). CDK8 acts as an oncogene and represents a potential target for developing novel CDK8 inhibitors in cancer therapeutics. However, other study has revealed its contrary role. The function of CDK8 is context dependent. Even so, a variety of potent and selective CDK8 inhibitors have been discovered after crystal structures were resolved in two states (active or inactive). In this review, we summarize co-crystal structures, biological mechanisms, dysregulation in cancers and recent progress in the field of CDK8 inhibitors, trying to offer an outlook of CDK8 inhibitors in cancer therapy in future.

Discovery of CDK8/CycC Ligands with a New Virtual Screening Tool

Selective inhibition of cyclin-dependent kinase 8 and cyclin C (CDK8/CycC) has been suggested as a promising strategy for decreasing mitogenic signals in cancer cells with reduced toxicity toward normal cells. We developed a novel virtual screening protocol for drug development and applied it to the discovery of new CDK8/CycC type II ligands, which is likely to achieve long residence time and specificity. We first analyzed the binding thermodynamics of 11 published pyrazolourea ligands using molecular dynamics simulations and a free-energy calculation method, VM2, and extracted the key binding information to assist virtual screening. The urea moiety was found to be the critical structural contributor of the reference ligands. Starting with the urea moiety, we conducted substructure-based searches with our newly developed superposition and single-point energy evaluation method, followed by free-energy calculations, and singled out three purchasable compounds for bioassay testing. The ranking from the experimental results is completely consistent with the predicted rankings. A potent drug-like compound was found to have a K_d value of 42.5 nm, which is similar to those of the most potent reference ligands; this provided a good starting point for further improvement. This study shows that our novel virtual screening protocol is an accurate and efficient tool for drug development.

Binding Kinetics Survey of the Drugged Kinome

Target residence time is emerging as an important optimization parameter in drug discovery, yet target and off-target engagement dynamics have not been clearly linked to the clinical performance of drugs. Here we developed high-throughput binding kinetics assays to characterize the interactions of 270 protein kinase inhibitors with 40 clinically relevant targets. Analysis of the results revealed that on-rates are better correlated with affinity than off-rates and that the fraction of slowly dissociating drug-target complexes increases from early/preclinical to late stage and FDA-approved compounds, suggesting distinct contributions by each parameter to clinical success. Combining binding parameters with PK/ADME properties, we illustrate in silico and in cells how kinetic selectivity could be exploited as an optimization strategy. Furthermore, using bio- and chemoinformatics we uncovered structural features influencing rate constants. Our results underscore the value of binding kinetics information in rational drug design and provide a resource for future studies on this subject.

KinITC-One Method Supports both Thermodynamic and Kinetic SARs as Exemplified on FimH Antagonists

Affinity data, such as dissociation constants (K_D ) or inhibitory concentrations (IC₅₀ ), are widely used in drug discovery. However, these parameters describe an equilibrium state, which is often not established in vivo due to pharmacokinetic effects and they are therefore not necessarily sufficient for evaluating drug efficacy. More accurate indicators for pharmacological activity are the kinetics of binding processes, as they shed light on the rate of formation of protein-ligand complexes and their half-life. Nonetheless, although highly desirable for medicinal chemistry programs, studies on structure-kinetic relationships (SKR) are still rare. With the recently introduced analytical tool kinITC this situation may change, since not only thermodynamic but also kinetic information of the binding process can be deduced from isothermal titration calorimetry (ITC) experiments. Using kinITC, ITC data of 29 mannosides binding to the bacterial adhesin FimH were re-analyzed to make their binding kinetics accessible. To validate these kinetic data, surface plasmon resonance (SPR) experiments were conducted. The kinetic analysis by kinITC revealed that the nanomolar affinities of the FimH antagonists arise from both (i) an optimized interaction between protein and ligand in the bound state (reduced off-rate constant k_off ) and (ii) a stabilization of the transition state or a destabilization of the unbound state (increased on-rate constant k_on ). Based on congeneric ligand modifications and structural input from co-crystal structures, a strong relationship between the formed hydrogen-bond network and k_off could be concluded, whereas electrostatic interactions and conformational restrictions upon binding were found to have mainly an impact on k_on .

Structure-kinetic relationships that control the residence time of drug-target complexes: insights from molecular structure and dynamics

Time-dependent target occupancy is a function of both the thermodynamics and kinetics of drug-target interactions. However, while the optimization of thermodynamic affinity through approaches such as structure-based drug design is now relatively straight forward, less is understood about the molecular interactions that control the kinetics of drug complex formation and breakdown since this depends on both the ground and transition state energies on the binding reaction coordinate. In this opinion we highlight several recent examples that shed light on current approaches that are elucidating the factors that control the life-time of the drug-target complex.

Dynamics of human protein kinase Aurora A linked to drug selectivity

Protein kinases are major drug targets, but the development of highly-selective inhibitors has been challenging due to the similarity of their active sites. The observation of distinct structural states of the fully-conserved Asp-Phe-Gly (DFG) loop has put the concept of conformational selection for the DFG-state at the center of kinase drug discovery. Recently, it was shown that Gleevec selectivity for the Tyr-kinase Abl was instead rooted in conformational changes after drug binding. Here, we investigate whether protein dynamics after binding is a more general paradigm for drug selectivity by characterizing the binding of several approved drugs to the Ser/Thr-kinase Aurora A. Using a combination of biophysical techniques, we propose a universal drug-binding mechanism, that rationalizes selectivity, affinity and long on-target residence time for kinase inhibitors. These new concepts, where protein dynamics in the drug-bound state plays the crucial role, can be applied to inhibitor design of targets outside the kinome.

Protein kinases are a family of enzymes found in all living organisms. These enzymes help to control many biological processes, including cell division. When particular protein kinases do not work correctly, cells may start to divide uncontrollably, which can lead to cancer. One example is the kinase Aurora A, which is over-active in many common human cancers. As a result, researchers are currently trying to design drugs that reduce the activity of Aurora A in the hope that these could form new anticancer treatments.

In general, drugs are designed to be as specific in their action as possible to reduce the risk of harmful side effects to the patient. Designing a drug that affects a single protein kinase, however, is difficult because there are hundreds of different kinases in the body, all with similar structures. Because drugs often work by binding to specific structural features, a drug that targets one protein kinase can often alter the activity of a large number of others too.

Gleevec is a successful anti-leukemia drug that specifically works on one target kinase, producing minimal side effects. It was recently discovered that the drug works through a phenomenon called ‘induced fit’. This means that after the drug binds it causes a change in the enzyme’s overall shape that alters the activity of the enzyme. The shape change is complex, and so even small structural differences can change the effect of a particular drug.

Do other drugs that target other protein kinases also produce induced fit effects? To find out, Pitsawong, Buosi, Otten, Agafonov et al. studied how three anti-cancer drugs interact with Aurora A: two drugs specifically designed to switch off Aurora A, and Gleevec (which does not target Aurora A).

The two drugs that specifically target Aurora A were thought to work by targeting one structural feature of the enzyme. However, the biochemical and biophysical experiments performed by Pitsawong et al. revealed that these drugs instead work through an induced fit effect. By contrast, Gleevec did not trigger an induced fit on Aurora A and so bound less tightly to it.

In light of these results, Pitsawong et al. suggest that future efforts to design drugs that target protein kinases should focus on exploiting the induced fit process. This will require more research into the structure of particular kinases.

Halogen-Aromatic π Interactions Modulate Inhibitor Residence Times

Prolonged drug residence times may result in longer-lasting drug efficacy, improved pharmacodynamic properties, and "kinetic selectivity" over off-targets with high drug dissociation rates. However, few strategies have been elaborated to rationally modulate drug residence time and thereby to integrate this key property into the drug development process. Herein, we show that the interaction between a halogen moiety on an inhibitor and an aromatic residue in the target protein can significantly increase inhibitor residence time. By using the interaction of the serine/threonine kinase haspin with 5-iodotubercidin (5-iTU) derivatives as a model for an archetypal active-state (type I) kinase-inhibitor binding mode, we demonstrate that inhibitor residence times markedly increase with the size and polarizability of the halogen atom. The halogen-aromatic π interactions in the haspin-inhibitor complexes were characterized by means of kinetic, thermodynamic, and structural measurements along with binding-energy calculations.

Recent advances in CDK inhibitors for cancer therapy

Inhibition of CDKs is an attractive approach to cancer therapy due to their vital role in cell growth and transcription. Pan-CDK inhibitors have shown some clinical benefit, and trials are ongoing. Selective CDK4 and CDK6 inhibitors have been licensed for the treatment of hormone responsive, RB-positive breast cancer in combination with antihormonal agents. Selective inhibitors of CDKs 5, 7, 8, 9 and 12 have been identified across a range of chemotypes.