All-atomic molecular dynamic studies of human CDK8: insight into the A-loop, point mutations and binding with its partner CycC

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.

Unraveling Bacterial Single-Stranded Sequence Specificities: Insights from Molecular Dynamics and MMPBSA Analysis of Oligonucleotide Probes

We utilized molecular dynamics (MD) simulations and Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) free energy calculations to investigate the specificity of two oligonucleotide probes, namely probe B and probe D, in detecting single-stranded DNA (ssDNA) within three bacteria families: Enterobacteriaceae, Pasteurellaceae, and Vibrionaceae. Due to the limited understanding of molecular mechanisms in the previous research, we have extended the discussion to focus specifically on investigating the binding process of bacteria-probe DNA duplexes, with an emphasis on analyzing the binding free energy. The role of electrostatic contributions in the specificity between the oligonucleotide probes and the bacterial ssDNAs was investigated and found to be crucial. Our calculations yielded results that were highly consistent with the experimental data. Through our study, we have successfully exhibited the benefits of utilizing in-silico approaches as a powerful virtual-screening tool, particularly in research areas that demand a thorough comprehension of molecular interactions.

Introducing mirror-image discrimination capability to the TSR-based method for capturing stereo geometry and understanding hierarchical structure relationships of protein receptor family

We have developed a Triangular Spatial Relationship (TSR)-based computational method for protein structure comparison and motif discovery that is both sequence and structure alignment-free. A protein 3D structure is modeled by all possible triangles that are constructed with every three Cα atoms of amino acids as vertices. Every triangle is represented using an integer (a key). The keys are calculated by a rule-based formula which is a function of a representative length, a representative angle, and the vertex labels associated with amino acids. A 3D structure is thereby represented by a vector of integers (TSR keys). Global or local structure comparisons are achieved by computing all keys or a set of keys, respectively. Many enzymatic reactions and notable marketed drugs are highly stereospecific. Thus, in this paper, we propose a modified key calculation formula by including a mechanism for discriminating mirror-image keys to capture stereo geometry. We assign a positive or a negative sign to the integers representing mirror-image keys. Applying the new key calculation function provides the ability to further discriminate mirror-image keys that were previously considered identical. As the result, applying the mirror-image discrimination capability (i) significantly increases the number of distinct keys; (ii) decreases the number of common keys; (iii) decreases structural similarity; (iv) increases the opportunity to identify specific keys for each type of the receptors. The specific keys identified in this study for the cases of without (not applying) and with (applying) mirror-image discrimination can be considered as the structure signatures that exclusively belong to a certain type of receptors. Applying mirror-image discrimination introduces stereospecificity to keys for allowing more precise modeling of ligand - target interactions. The development of mirror-image TSR keys of Cα atom, in conjunction with the integration of Cα TSR keys with all-atom TSR keys for amino acids and drugs, will lead to a new and promising computational method for aiding drug design and discovery.

Cdk8 attenuates lipogenesis by inhibiting SREBP-dependent transcription in Drosophila

Fine-tuning of lipogenic gene expression is important for the maintenance of long-term homeostasis of intracellular lipids. The SREBP family of transcription factors are master regulators that control the transcription of lipogenic and cholesterogenic genes, but the mechanisms modulating SREBP-dependent transcription are still not fully understood. We previously reported that CDK8, a subunit of the transcription co-factor Mediator complex, phosphorylates SREBP at a conserved threonine residue. Here, using Drosophila as a model system, we observed that the phosphodeficient SREBP proteins (SREBP-Thr390Ala) were more stable and more potent in stimulating the expression of lipogenic genes and promoting lipogenesis in vivo than wild-type SREBP. In addition, starvation blocked the effects of wild-type SREBP-induced lipogenic gene transcription, whereas phosphodeficient SREBP was resistant to this effect. Furthermore, our biochemical analyses identified six highly conserved amino acid residues in the N-terminus disordered region of SREBP that are required for its interactions with both Cdk8 and the MED15 subunit of the small Mediator complex. These results support that the concerted actions of Cdk8 and MED15 are essential for the tight regulation of SREBP-dependent transcription. This article has an associated First Person interview with the first author of the paper.

A Computational Study on the Interaction of NSP10 and NSP14: Unraveling the RNA Synthesis Proofreading Mechanism in SARS-CoV-2, SARS-CoV, and MERS-CoV

The interaction of exoribonuclease (ExoN) nonstructural protein (NSP14) with NSP10 co-factors is crucial for high-fidelity proofreading activity of coronavirus replication and transcription. Proofreading function is critical for maintaining the large genomes to ensure replication proficiency; therefore, while maintaining the viral replication fitness, quick resistance has been reported to the nucleotide analogue (NA) drugs. Therefore, targeting the NSP14 and NSP10 interacting interface with small molecules or peptides could be a better strategy to obstruct replication processes of coronaviruses (CoVs). A comparative study on the binding mechanism of NSP10 with the NSP14 ExoN domain of SARS-CoV-2, SARS-CoV, MERS-CoV, and four SARS-CoV-2 NSP14mutant complexes has been carried out. Protein-protein interaction (PPI) dynamics, per-residue binding free energy (BFE) analyses, and the identification of interface hotspot residues have been studied using molecular dynamics simulations and various computational tools. The BFE of the SARS-CoV NSP14-NSP10 complex was higher when compared to novel SARS-CoV-2 and MERS. However, SARS-CoV-2 NSP14mutant systems display a higher BFE as compared to the wild type (WT) but lower than SARS-CoV and MERS. Despite the high BFE, the SARS-CoV NSP14-NSP10 complex appears to be structurally more flexible in many regions especially the catalytic site, which is not seen in SARS-CoV-2 and its mutant or MERS complexes. The significantly high residue energy contribution of key interface residues and hotspots reveals that the high binding energy between NSP14 and NSP10 may enhance the functional activity of the proofreading complex, as the NSP10-NSP14 interaction is essential in maintaining the stability of the ExoN domain for the replicative fitness of CoVs. The factors discussed for SARS-CoV-2 complexes may be responsible for NSP14 ExoN having a high replication proficiency, significantly leading to the evolution of new variants of SARS-CoV-2. The NSP14 residues V66, T69, D126, and I201and eight residues of NSP10 (L16, F19, V21, V42, M44, H80, K93, and F96) are identified as common hotspots. Overall, the interface area, hotspot locations, bonded/nonbonded contacts, and energies between NSP14 and NSP10 may pave a way in designing potential inhibitors to disrupt NSP14-NSP10 interactions of CoVs especially SARS-CoV-2.

Development of a TSR-Based Method for Protein 3-D Structural Comparison With Its Applications to Protein Classification and Motif Discovery

Development of protein 3-D structural comparison methods is important in understanding protein functions. At the same time, developing such a method is very challenging. In the last 40 years, ever since the development of the first automated structural method, ~200 papers were published using different representations of structures. The existing methods can be divided into five categories: sequence-, distance-, secondary structure-, geometry-based, and network-based structural comparisons. Each has its uniqueness, but also limitations. We have developed a novel method where the 3-D structure of a protein is modeled using the concept of Triangular Spatial Relationship (TSR), where triangles are constructed with the C_α atoms of a protein as vertices. Every triangle is represented using an integer, which we denote as “key,” A key is computed using the length, angle, and vertex labels based on a rule-based formula, which ensures assignment of the same key to identical TSRs across proteins. A structure is thereby represented by a vector of integers. Our method is able to accurately quantify similarity of structure or substructure by matching numbers of identical keys between two proteins. The uniqueness of our method includes: (i) a unique way to represent structures to avoid performing structural superimposition; (ii) use of triangles to represent substructures as it is the simplest primitive to capture shape; (iii) complex structure comparison is achieved by matching integers corresponding to multiple TSRs. Every substructure of one protein is compared to every other substructure in a different protein. The method is used in the studies of proteases and kinases because they play essential roles in cell signaling, and a majority of these constitute drug targets. The new motifs or substructures we identified specifically for proteases and kinases provide a deeper insight into their structural relations. Furthermore, the method provides a unique way to study protein conformational changes. In addition, the results from CATH and SCOP data sets clearly demonstrate that our method can distinguish alpha helices from beta pleated sheets and vice versa. Our method has the potential to be developed into a powerful tool for efficient structure-BLAST search and comparison, just as BLAST is for sequence search and alignment.

Post-translational Modifications of Serine/Threonine and Histidine Kinases and Their Roles in Signal Transductions in Synechocystis Sp. PCC 6803

Cyanobacterium Synechocystis sp. PCC 6803, a popular model organism for researches in photosynthesis and biofuel production, contains plant-like photosynthetic machineries which significantly contribute to global carbon fixation. There are 12 eukaryotic-type Ser/Thr kinases (SpkA-L) and 49 His kinases (Hik1-49) of two-component systems in the genome of Synechocystis sp. PCC 6803. They are the key regulators in sensing and transmitting stimuli including light- and glucose-mediate signal transduction. Proteomic studies were able to identify all the kinases. The majority of kinases no matter whether they have a predicted transmembrane domain were identified in the membrane fractions. Six Ser/Thr kinases (SpkA-D, F and G) and ten His kinases (Hik4, 12, 14, 21, 26-27, 29, 36, 43, and 46) were identified to have one or more of the three types of post-translational modifications: phosphorylation, acetylation, and thiol oxidation. Interestingly, SpkG has the phosphorylatable threonine residue that was aligned with the phosphorylated threonine residue in the activation loop of human CDK7, demonstrating conserved phosphorylation between cyanobacterial and human kinases. Transcriptomics and proteomics revealed differential expression of the kinases in heterotrophic and photoheterotrophic compared with photoautotrophic conditions, indicating their roles in regulating the growth modes of cyanobacteria. In summary, this review focuses on the discussions on post-transcriptional modifications, transcriptomic, and proteomic studies of Ser/Thr and His kinases. This together with our published review in 2019 present a complete story of an overview of sequences, domain architectures, and biochemical and physiological functions of cyanobacterial kinases with adequate details in the context of high throughput systems. We also emphasize the importance of discovering upstream molecules and substrates to understand the exact functions of the kinases in vivo. As an attempt, a model is proposed in which Hik31, His33, Sll1334, and IcfG are hypothesized to be critical for switching between autotrophic and heterotrophic modes based on the results from the phenotypes of the gene knockout strains combined with their post-translational modifications, and gene expression profiles.

All-Atomic Molecular Dynamic Studies of Human and Drosophila CDK8: Insights into Their Kinase Domains, the LXXLL Motifs, and Drug Binding Site

Cyclin-dependent kinase 8 (CDK8) and its regulatory partner Cyclin C (CycC) play conserved roles in modulating RNA polymerase II (Pol II)-dependent gene expression. To understand the structure and function relations of CDK8, we analyzed the structures of human and Drosophila CDK8 proteins using molecular dynamics simulations, combined with functional analyses in Drosophila. Specifically, we evaluated the structural differences between hCDK8 and dCDK8 to predict the effects of the LXXLL motif mutation (AQKAA), the P154L mutations, and drug binding on local structures of the CDK8 proteins. First, we have observed that both the LXXLL motif and the kinase activity of CDK8 are required for the normal larval-to-pupal transition in Drosophila. Second, our molecular dynamic analyses have revealed that hCDK8 has higher hydrogen bond occupation of His149-Asp151 and Asp151-Asn156 than dCDK8. Third, the substructure of Asp282, Phe283, Arg285, Thr287 and Cys291 can distinguish human and Drosophila CDK8 structures. In addition, there are two hydrogen bonds in the LXXLL motif: a lower occupation between L312 and L315, and a relatively higher occupation between L312 and L316. Human CDK8 has higher hydrogen bond occupation between L312 and L316 than dCDK8. Moreover, L312, L315 and L316 in the LXXLL motif of CDK8 have the specific pattern of hydrogen bonds and geometries, which could be crucial for the binding to nuclear receptors. Furthermore, the P154L mutation dramatically decreases the hydrogen bond between L312 and L315 in hCDK8, but not in dCDK8. The mutations of P154L and AQKAA modestly alter the local structures around residues 154. Finally, we identified the inhibitor-induced conformational changes of hCDK8, and our results suggest a structural difference in the drug-binding site between hCDK8 and dCDK8. Taken together, these results provide the structural insights into the roles of the LXXLL motif and the kinase activity of CDK8 in vivo.

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.

Evaluation of drug candidature: In silico ADMET, binding interactions with CDK7 and normal cell line studies of potentially anti-breast cancer enamidines

We have recently explored novel class of potentially anti-breast cancer active enamidines in which four molecules 4a-c and 4h showed higher anticancer activity compared to standard drug doxorubicin. As a part of extension of this work, we have further evaluated in silico cheminformatic studies on bioactivity prediction of synthesized series of enamidines using mole information. The normal cell line study of four lead compounds 4a-c and 4h against African green monkey kidney vero strain further revealed that the compounds complemented good selectivity in inhibition of cancer cells. The in silico bioactivity and molecular docking studies also revealed that the compounds have significant interactions with the drug targets. The results reveal that enamidine moieties are vital for anti-breast cancer activity as they possess excellent drug-like characteristics, being potentially good inhibitors of cyclin dependent kinases7 (CDK7).

Sequences, Domain Architectures, and Biological Functions of the Serine/Threonine and Histidine Kinases in Synechocystis sp. PCC 6803

The cyanobacterium Synechocystis sp. PCC 6803 (hereafter Synechocystis) is a photoautotrophic prokaryote with plant-like photosynthetic machineries which significantly contribute to global carbon fixation and atmospheric oxygen production. Because of the relatively short cell doubling time, small size of the genome, and the ease for genetic manipulation, Synechocystis is a popular model organism for studies including photosynthesis and biofuel production. The cyanobacterium contains 12 eukaryotic type Ser/Thr kinases (SpkA-L) and 49 histidine kinases (Hik1-47 and Sll1334 and Sll5060 are named as Hik48 and Hik49, respectively, in this review) of the two-component system. All SpkA-L kinases have a eukaryotic kinase DFG signature in their A-loops. Based on the types of the kinase domains, the Spks can be separated into three groups: one group contains SpkA and SpkG which are related to human kinases, while SpkH-L are in another group that is distinct from human kinases. The third group contains SpkB-F which are between the first two groups. Four histidine kinases (Hiks17, 36, 45, and 48) lack a clear histidine kinase domain, and the conserved phosphorylatable histidine residue could not be identified for six histidine kinases (Hiks11, 18, 29, 37, 39, and 43) even though they have clear histidine kinase domains. Each of the remaining 39 has a histidine kinase domain with the conserved histidine residue. Eight hybrid histidine kinases contain one or two receiver domains, and they all, except Hik25 (Slr0222), have the conserved phosphorylatable aspartate. The disruptants of all kinases except hik13 and hik15 have been generated, and the majority of them have modest or no obvious phenotypes, indicating other kinases could functionally compensate the loss of a particular kinase. This review presents a comprehensive discussion including a spectrum of sequence, domain architecture, in vivo function, and proteomics investigations of Ser/Thr and histidine kinases. Understanding the sequences, domain architectures, and biology of the kinases will help to integrate "omic" data to clarify their exact biochemical functions.

A molecular dynamics investigation of CDK8/CycC and ligand binding: conformational flexibility and implication in drug discovery

Abnormal activity of cyclin-dependent kinase 8 (CDK8) along with its partner protein cyclin C (CycC) is a common feature of many diseases including colorectal cancer. Using molecular dynamics (MD) simulations, this study determined the dynamics of the CDK8-CycC system and we obtained detailed breakdowns of binding energy contributions for four type-I and five type-II CDK8 inhibitors. We revealed system motions and conformational changes that will affect ligand binding, confirmed the essentialness of CycC for inclusion in future computational studies, and provide guidance in development of CDK8 binders. We employed unbiased all-atom MD simulations for 500 ns on twelve CDK8-CycC systems, including apoproteins and protein-ligand complexes, then performed principal component analysis (PCA) and measured the RMSF of key regions to identify protein dynamics. Binding pocket volume analysis identified conformational changes that accompany ligand binding. Next, H-bond analysis, residue-wise interaction calculations, and MM/PBSA were performed to characterize protein-ligand interactions and find the binding energy. We discovered that CycC is vital for maintaining a proper conformation of CDK8 to facilitate ligand binding and that the system exhibits motion that should be carefully considered in future computational work. Surprisingly, we found that motion of the activation loop did not affect ligand binding. Type-I and type-II ligand binding is driven by van der Waals interactions, but electrostatic energy and entropic penalties affect type-II binding as well. Binding of both ligand types affects protein flexibility. Based on this we provide suggestions for development of tighter-binding CDK8 inhibitors and offer insight that can aid future computational studies.

Recent Developments and Applications of the MMPBSA Method

The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.

In silico study of porphyrin-anthraquinone hybrids as CDK2 inhibitor

Cyclin-Dependent Kinases (CDKs) are known to play crucial roles in controlling cell cycle progression of eukaryotic cell and inhibition of their activity has long been considered as potential strategy in anti-cancer drug research. In the present work, a series of porphyrin-anthraquinone hybrids bearing meso-substituents, i.e. either pyridine or pyrazole rings were designed and computationally evaluated for their Cyclin Dependent Kinase-2 (CDK2) inhibitory activity using molecular docking, molecular dynamics simulation, and binding free energy calculation. The molecular docking simulation revealed that all six porphyrin hybrids were able to bind to ATP-binding site of CDK2 and interacted with key residues constituted the active cavity of CDK2, while molecular dynamics simulation indicated that all porphyrins bound to CDK2 were stable for 6ns. The binding free energies predicted by MM-PBSA method showed that most compounds exhibited higher affinity than that of native ligand (4-anilinoquinazoline, DTQ) and the affinity of mono-H₂PyP-AQ was about three times better than that of DTQ, indicating its potential to be advanced as a new CDK2 inhibitor.

Experimental and molecular dynamics studies showed that CBP KIX mutation affects the stability of CBP:c-Myb complex

The coactivators CBP (CREBBP) and its paralog p300 (EP300), two conserved multi-domain proteins in eukaryotic organisms, regulate gene expression in part by binding DNA-binding transcription factors. It was previously reported that the CBP/p300 KIX domain mutant (Y650A, A654Q, and Y658A) altered both c-Myb-dependent gene activation and repression, and that mice with these three point mutations had reduced numbers of platelets, B cells, T cells, and red blood cells. Here, our transient transfection assays demonstrated that mouse embryonic fibroblast cells containing the same mutations in the KIX domain and without a wild-type allele of either CBP or p300, showed decreased c-Myb-mediated transcription. Dr. Wright's group solved a 3-D structure of the mouse CBP:c-Myb complex using NMR. To take advantage of the experimental structure and function data and improved theoretical calculation methods, we performed MD simulations of CBP KIX, CBP KIX with the mutations, and c-Myb, as well as binding energy analysis for both the wild-type and mutant complexes. The binding between CBP and c-Myb is mainly mediated by a shallow hydrophobic groove in the center where the side-chain of Leu302 of c-Myb plays an essential role and two salt bridges at the two ends. We found that the KIX mutations slightly decreased stability of the CBP:c-Myb complex as demonstrated by higher binding energy calculated using either MM/PBSA or MM/GBSA methods. More specifically, the KIX mutations affected the two salt bridges between CBP and c-Myb (CBP-R646 and c-Myb-E306; CBP-E665 and c-Myb-R294). Our studies also revealed differing dynamics of the hydrogen bonds between CBP-R646 and c-Myb-E306 and between CBP-E665 and c-Myb-R294 caused by the CBP KIX mutations. In the wild-type CBP:c-Myb complex, both of the hydrogen bonds stayed relatively stable. In contrast, in the mutant CBP:c-Myb complex, hydrogen bonds between R646 and E306 showed an increasing trend followed by a decreasing trend, and hydrogen bonds of the E665:R294 pair exhibited a fast decreasing trend over time during MD simulations. In addition, our data showed that the KIX mutations attenuate CBP's hydrophobic interaction with Leu302 of c-Myb. Furthermore, our 500-ns MD simulations showed that CBP KIX with the mutations has a slightly lower potential energy than wild-type CBP. The CBP KIX structures with or without its interacting protein c-Myb are different for both wild-type and mutant CBP KIX, and this is likewise the case for c-Myb with or without CBP, suggesting that the presence of an interacting protein influences the structure of a protein. Taken together, these analyses will improve our understanding of the exact functions of CBP and its interaction with c-Myb.

CDK8-Cyclin C Mediates Nutritional Regulation of Developmental Transitions through the Ecdysone Receptor in Drosophila

The steroid hormone ecdysone and its receptor (EcR) play critical roles in orchestrating developmental transitions in arthropods. However, the mechanism by which EcR integrates nutritional and developmental cues to correctly activate transcription remains poorly understood. Here, we show that EcR-dependent transcription, and thus, developmental timing in Drosophila, is regulated by CDK8 and its regulatory partner Cyclin C (CycC), and the level of CDK8 is affected by nutrient availability. We observed that cdk8 and cycC mutants resemble EcR mutants and EcR-target genes are systematically down-regulated in both mutants. Indeed, the ability of the EcR-Ultraspiracle (USP) heterodimer to bind to polytene chromosomes and the promoters of EcR target genes is also diminished. Mass spectrometry analysis of proteins that co-immunoprecipitate with EcR and USP identified multiple Mediator subunits, including CDK8 and CycC. Consistently, CDK8-CycC interacts with EcR-USP in vivo; in particular, CDK8 and Med14 can directly interact with the AF1 domain of EcR. These results suggest that CDK8-CycC may serve as transcriptional cofactors for EcR-dependent transcription. During the larval–pupal transition, the levels of CDK8 protein positively correlate with EcR and USP levels, but inversely correlate with the activity of sterol regulatory element binding protein (SREBP), the master regulator of intracellular lipid homeostasis. Likewise, starvation of early third instar larvae precociously increases the levels of CDK8, EcR and USP, yet down-regulates SREBP activity. Conversely, refeeding the starved larvae strongly reduces CDK8 levels but increases SREBP activity. Importantly, these changes correlate with the timing for the larval–pupal transition. Taken together, these results suggest that CDK8-CycC links nutrient intake to developmental transitions (EcR activity) and fat metabolism (SREBP activity) during the larval–pupal transition.

During the larval-pupal transition in Drosophila, CDK8-CycC helps to link nutrient intake to development by activating ecdysone receptor-dependent transcription and to fat metabolism by inhibiting SREBP-activated gene expression.

Arthropods are estimated to account for over 80% of animal species on earth. Characterized by their rigid exoskeletons, juvenile arthropods must periodically shed their thick outer cuticles by molting in order to grow. The steroid hormone ecdysone plays an essential role in regulating the timing of developmental transitions, but exactly how ecdysone and its receptor EcR activates transcription correctly after integrating nutritional and developmental cues remains unknown. Our developmental genetic analyses of two Drosophila mutants, cdk8 and cycC, show that they are lethal during the prepupal stage, with aberrant accumulation of fat and a severely delayed larval–pupal transition. As we have reported previously, CDK8-CycC inhibits fat accumulation by directly inactivating SREBP, a master transcription factor that controls the expression of lipogenic genes, which explains the abnormal fat accumulation in the cdk8 and cycC mutants. We find that CDK8 and CycC are required for EcR to bind to its target genes, serving as transcriptional cofactors for EcR-dependent gene expression. The expression of EcR target genes is compromised in cdk8 and cycC mutants and underpins the retarded pupariation phenotype. Starvation of feeding larvae precociously up-regulates CDK8 and EcR, prematurely down-regulates SREBP activity, and leads to early pupariation, whereas re-feeding starved larvae has opposite effects. Taken together, these results suggest that CDK8 and CycC play important roles in coordinating nutrition intake with fat metabolism by directly inhibiting SREBP-dependent gene expression and regulating developmental timing by activating EcR-dependent transcription in Drosophila.

Mediator kinase module and human tumorigenesis

Mediator is a conserved multi-subunit signal processor through which regulatory informatiosn conveyed by gene-specific transcription factors is transduced to RNA Polymerase II (Pol II). In humans, MED13, MED12, CDK8 and Cyclin C (CycC) comprise a four-subunit "kinase" module that exists in variable association with a 26-subunit Mediator core. Genetic and biochemical studies have established the Mediator kinase module as a major ingress of developmental and oncogenic signaling through Mediator, and much of its function in signal-dependent gene regulation derives from its resident CDK8 kinase activity. For example, CDK8-targeted substrate phosphorylation impacts transcription factor half-life, Pol II activity and chromatin chemistry and functional status. Recent structural and biochemical studies have revealed a precise network of physical and functional subunit interactions required for proper kinase module activity. Accordingly, pathologic change in this activity through altered expression or mutation of constituent kinase module subunits can have profound consequences for altered signaling and tumor formation. Herein, we review the structural organization, biological function and oncogenic potential of the Mediator kinase module. We focus principally on tumor-associated alterations in kinase module subunits for which mechanistic relationships as opposed to strictly correlative associations are established. These considerations point to an emerging picture of the Mediator kinase module as an oncogenic unit, one in which pathogenic activation/deactivation through component change drives tumor formation through perturbation of signal-dependent gene regulation. It follows that therapeutic strategies to combat CDK8-driven tumors will involve targeted modulation of CDK8 activity or pharmacologic manipulation of dysregulated CDK8-dependent signaling pathways.