Epiblastin A Induces Reprogramming of Epiblast Stem Cells Into Embryonic Stem Cells by Inhibition of Casein Kinase 1

Past, Present, and Future Perspectives on Computer-Aided Drug Design Methodologies

The application of computational approaches in drug discovery has been consolidated in the last decades. These families of techniques are usually grouped under the common name of "computer-aided drug design" (CADD), and they now constitute one of the pillars in the pharmaceutical discovery pipelines in many academic and industrial environments. Their implementation has been demonstrated to tremendously improve the speed of the early discovery steps, allowing for the proficient and rational choice of proper compounds for a desired therapeutic need among the extreme vastness of the drug-like chemical space. Moreover, the application of CADD approaches allows the rationalization of biochemical and interactive processes of pharmaceutical interest at the molecular level. Because of this, computational tools are now extensively used also in the field of rational 3D design and optimization of chemical entities starting from the structural information of the targets, which can be experimentally resolved or can also be obtained with other computer-based techniques. In this work, we revised the state-of-the-art computer-aided drug design methods, focusing on their application in different scenarios of pharmaceutical and biological interest, not only highlighting their great potential and their benefits, but also discussing their actual limitations and eventual weaknesses. This work can be considered a brief overview of computational methods for drug discovery.

Thermal Titration Molecular Dynamics (TTMD): Not Your Usual Post-Docking Refinement

Molecular docking is one of the most widely used computational approaches in the field of rational drug design, thanks to its favorable balance between the rapidity of execution and the accuracy of provided results. Although very efficient in exploring the conformational degrees of freedom available to the ligand, docking programs can sometimes suffer from inaccurate scoring and ranking of generated poses. To address this issue, several post-docking filters and refinement protocols have been proposed throughout the years, including pharmacophore models and molecular dynamics simulations. In this work, we present the first application of Thermal Titration Molecular Dynamics (TTMD), a recently developed method for the qualitative estimation of protein-ligand unbinding kinetics, to the refinement of docking results. TTMD evaluates the conservation of the native binding mode throughout a series of molecular dynamics simulations performed at progressively increasing temperatures with a scoring function based on protein-ligand interaction fingerprints. The protocol was successfully applied to retrieve the native-like binding pose among a set of decoy poses of drug-like ligands generated on four different pharmaceutically relevant biological targets, including casein kinase 1δ, casein kinase 2, pyruvate dehydrogenase kinase 2, and SARS-CoV-2 main protease.

Qualitative Estimation of Protein-Ligand Complex Stability through Thermal Titration Molecular Dynamics Simulations

The prediction of ligand efficacy has long been linked to thermodynamic properties such as the equilibrium dissociation constant, which considers both the association and the dissociation rates of a defined protein-ligand complex. In the last 15 years, there has been a paradigm shift, with an increased interest in the determination of kinetic properties such as the drug-target residence time since they better correlate with ligand efficacy compared to other parameters. In this article, we present thermal titration molecular dynamics (TTMD), an alternative computational method that combines a series of molecular dynamics simulations performed at progressively increasing temperatures with a scoring function based on protein-ligand interaction fingerprints for the qualitative estimation of protein-ligand-binding stability. The protocol has been applied to four different pharmaceutically relevant test cases, including protein kinase CK1δ, protein kinase CK2, pyruvate dehydrogenase kinase 2, and SARS-CoV-2 main protease, on a variety of ligands with different sizes, structures, and experimentally determined affinity values. In all four cases, TTMD was successfully able to distinguish between high-affinity compounds (low nanomolar range) and low-affinity ones (micromolar), proving to be a useful screening tool for the prioritization of compounds in a drug discovery campaign.

Circadian clock modulating small molecules repurposing as inhibitors of SARS-CoV-2 Mpro for pharmacological interventions in COVID-19 pandemic

The COVID-19 pandemic caused by SARS-CoV-2 is a global health emergency warranting the development of targeted treatment. The main protease M^pro is considered as a key drug target in coronavirus infections because of its vital role in the proteolytic processing of two essential polyproteins required for the replication and transcription of viral RNA. Targeting and inhibiting the M^pro activity represents a valid approach to prevent the SARS-CoV-2 replication and spread. Based on the structure-assisted drug designing, here we report a circadian clock-modulating small molecule “SRT2183” as a potent inhibitor of M^pro to block the replication of SARS-CoV-2. The findings are expected to pave the way for the development of therapeutics for COVID-19.

Cooperative genetic networks drive embryonic stem cell transition from naïve to formative pluripotency

In the mammalian embryo, epiblast cells must exit the naïve state and acquire formative pluripotency. This cell state transition is recapitulated by mouse embryonic stem cells (ESCs), which undergo pluripotency progression in defined conditions in vitro. However, our understanding of the molecular cascades and gene networks involved in the exit from naïve pluripotency remains fragmentary. Here, we employed a combination of genetic screens in haploid ESCs, CRISPR/Cas9 gene disruption, large‐scale transcriptomics and computational systems biology to delineate the regulatory circuits governing naïve state exit. Transcriptome profiles for 73 ESC lines deficient for regulators of the exit from naïve pluripotency predominantly manifest delays on the trajectory from naïve to formative epiblast. We find that gene networks operative in ESCs are also active during transition from pre‐ to post‐implantation epiblast in utero. We identified 496 naïve state‐associated genes tightly connected to the in vivo epiblast state transition and largely conserved in primate embryos. Integrated analysis of mutant transcriptomes revealed funnelling of multiple gene activities into discrete regulatory modules. Finally, we delineate how intersections with signalling pathways direct this pivotal mammalian cell state transition.

The phosphorylation switch that regulates ticking of the circadian clock

In our 24/7 well-lit world, it's easy to skip or delay sleep to work, study, and play. However, our circadian rhythms are not easily fooled; the consequences of jet lag and shift work are many and severe, including metabolic, mood, and malignant disorders. The internal clock that keeps track of time has at its heart the reversible phosphorylation of the PERIOD proteins, regulated by isoforms of casein kinase 1 (CK1). In-depth biochemical, genetic, and structural studies of these kinases, their mutants, and their splice variants have combined over the past several years to provide a robust understanding of how the core clock is regulated by a phosphoswitch whereby phosphorylation of a stabilizing site on PER blocks phosphorylation of a distant phosphodegron. The recent structure of a circadian mutant form of CK1 implicates an internal activation loop switch that regulates this phosphoswitch and points to new approaches to regulation of the clock.

Phosphorylation of NANOG by casein kinase I regulates embryonic stem cell self-renewal

The self-renewal efficiency of mouse embryonic stem cells (ESCs) is determined by the concentration of the transcription factor NANOG. While NANOG binds thousands of sites in chromatin, the regulatory systems that control DNA binding are poorly characterised. Here, we show that NANOG is phosphorylated by casein kinase I, and identify target residues. Phosphomimetic substitutions at phosphorylation sites within the homeodomain (S130 and S131) have site-specific functional effects. Phosphomimetic substitution of S130 abolishes DNA binding by NANOG and eliminates LIF-independent self-renewal. In contrast, phosphomimetic substitution of S131 enhances LIF-independent self-renewal, without influencing DNA binding. Modelling the DNA-homeodomain complex explains the disparate effects of these phosphomimetic substitutions. These results indicate how phosphorylation may influence NANOG homeodomain interactions that underpin ESC self-renewal.

Simultaneous Control of Endogenous and User-Defined Genetic Pathways Using Unique ecDHFR Pharmacological Chaperones

Destabilizing domains (DDs), such as a mutated form of Escherichia coli dihydrofolate reductase (ecDHFR), confer instability and promote protein degradation. However, when combined with small-molecule stabilizers (e.g., the antibiotic trimethoprim), DDs allow positive regulation of fusion protein abundance. Using a combinatorial screening approach, we identified and validated 17 unique 2,4-diaminopyrimidine/triazine-based ecDHFR DD stabilizers, at least 15 of which were ineffective antibiotics against E. coli and S. aureus. Identified stabilizers functioned in vivo to control an ecDHFR DD-firefly luciferase in the mouse eye and/or the liver. Next, stabilizers were leveraged to perform synergistic dual functions in vitro (HeLa cell death sensitization) and in vivo (repression of ocular inflammation) by stabilizing a user-defined ecDHFR DD while also controlling endogenous signaling pathways. Thus, these newly identified pharmacological chaperones allow for simultaneous control of compound-specific endogenous and user-defined genetic pathways, the combination of which may provide synergistic effects in complex biological scenarios.

Kinase and GPCR polypharmacological approach for the identification of efficient anticancer medicines

Discovery of an anticancer medicine using a single target protein has often been unsuccessful due to the complexity of pathogenic mechanisms as well as the presence of redundant signaling pathways. In this work, we attempted to find promising anticancer drug candidates by simultaneously targeting casein kinase 1 delta (CK1δ) and muscarinic acetylcholine receptor M3 (M3R). Through the structure-based virtual screening and de novo design with the modified potential function for protein-ligand binding, a series of benzo[4,5]imidazo[1,2-a][1,3,5]triazine-2-amine (BITA) derivatives were identified as CK1δ inhibitors and also as M3R antagonists. The biochemical potencies of these bifunctional molecules reached the nanomolar and low-micromolar levels with respect to CK1δ and M3R, respectively. A common interaction feature in the calculated CK1δ-inhibitor and M3R-antagonist complexes is that the BITA moiety is well-stabilized in the orthosteric site of M3R and the hinge region of CK1δ through the establishment of the three hydrogen bonds and the hydrophobic contacts in the vicinity. The computational and experimental results found in this work exemplify the efficiency of kinase and GPCR polypharmacology in developing anticancer medicines.

3,4-Dibromo-7-Azaindole Modulates Arabidopsis Circadian Clock by Inhibiting Casein Kinase 1 Activity

The circadian clock is a timekeeping system for regulation of numerous biological daily rhythms. One characteristic of the circadian clock is that period length remains relatively constant in spite of environmental fluctuations, such as temperature change. Here, using the curated collection of in-house small molecule chemical library (ITbM chemical library), we show that small molecule 3,4-dibromo-7-azaindole (B-AZ) lengthened the circadian period of Arabidopsis thaliana (Arabidopsis). B-AZ has not previously been reported to have any biological and biochemical activities. Target identification can elucidate the mode of action of small molecules, but we were unable to make a molecular probe of B-AZ for target identification. Instead, we performed other analysis, gene expression profiling that potentially reveals mode of action of molecules. Short-term treatment of B-AZ decreased the expression of four dawn- and morning-phased clock-associated genes, CIRCADIAN CLOCK-ASSOCIATED 1 (CCA1), LATE ELONGATED HYPOCOTYL (LHY), PSEUDO-RESPONSE REGULATOR 9 (PRR9) and PRR7. Consistently, amounts of PRR5 and TIMING OF CAB EXPRESSION 1 (TOC1) proteins, transcriptional repressors of CCA1, LHY, PRR9 and PRR7 were increased upon B-AZ treatment. B-AZ inhibited Casein Kinase 1 family (CK1) that phosphorylates PRR5 and TOC1 for targeted degradation. A docking study and molecular dynamics simulation suggested that B-AZ interacts with the ATP-binding pocket of human CK1 delta, whose amino acid sequences are highly similar to those of Arabidopsis CK1. B-AZ-induced period-lengthening effect was attenuated in prr5 toc1 mutants. Collectively, this study provides a novel and simple structure CK1 inhibitor that modulates circadian clock via accumulation of PRR5 and TOC1.

Discovering New Casein Kinase 1d Inhibitors with an Innovative Molecular Dynamics Enabled Virtual Screening Workflow

The value of including protein flexibility in structure-based drug design (SBDD) is widely documented, and currently, molecular dynamics (MD) simulations represent a powerful tool to investigate protein dynamics. Yet, the inclusion of MD-derived information in pre-existing SBDD workflows is still far from trivial. We recently published an integrated MD-FLAP (Fingerprints for Ligands and Proteins) approach combining MD, clustering and Linear Discriminant Analysis (LDA) for enhancing accuracy, efficacy, and for protein conformational selection in virtual screening (VS) campaigns. Here we prospectively applied the MD-FLAP workflow to discover novel chemotypes inhibiting the Casein Kinase 1 delta (CSNK1D) enzyme. We first obtained a VS model able to separate active from inactive compounds, with a global AUC of 0.9 and a partial ROC enrichment at 0.5% of 0.18, and use it to mine the internal Pfizer screening database. Seven active molecules sharing a phenyl-indazole scaffold, not yet reported among CSNK1D inhibitors, were found. The most potent inhibitor showed an IC₅₀ of 134 nM.

Reducing Limitation in Probe Design: The Development of a Diazirine-Compatible Suzuki-Miyaura Cross Coupling Reaction

Access to high quality photoaffinity probe molecules is often constrained by synthetic limitations related to diazirine installation. A survey of recently published photoaffinity probe syntheses identified the Suzuki-Miyaura (S-M) coupling reaction, ubiquitous in drug discovery, as being underutilized to incorporate diazirines. To test whether advances in modern cross-coupling catalysis might enable efficient S-M couplings tolerant of the diazirine moiety, a fragment-based screening approach was employed. A model S-M coupling reaction was screened under various conditions in the presence of an aromatic diazirine fragment. This screen identified reaction conditions that gave good yields of S-M coupling product while minimally perturbing the diazirine reporter fragment. These conditions were found to be highly scalable and exhibited broad scope when applied to a chemistry informer library of 24 pharmaceutically relevant aryl boron pinacol esters. Furthermore, these conditions were used to synthesize a known diazirine-containing probe molecule with improved synthetic efficiency.

The Cell Painting Assay as a Screening Tool for the Discovery of Bioactivities in New Chemical Matter

Multiparametric phenotypic screening based on cellular morphology interrogates many biological pathways simultaneously and is therefore a valuable screening tool for the discovery of new biological activities. The cell painting assay stains various cellular features using six different dyes in one well. By automated image analysis, hundreds of parameters are calculated from the images which deliver a phenotypic profile of the cell. It has been shown that compounds with similar modes of action deliver similar phenotypic profiles. Using a reference set of compounds with known modes of action, it is possible to assign probable modes of action to new compounds and to discover compounds with potentially new modes of action.Here we describe the cell painting assay as a screening tool using a hit identification workflow which has been implemented using open-source software.

Casein kinase 1α: biological mechanisms and theranostic potential

Casein kinase 1α (CK1α) is a multifunctional protein belonging to the CK1 protein family that is conserved in eukaryotes from yeast to humans. It regulates signaling pathways related to membrane trafficking, cell cycle progression, chromosome segregation, apoptosis, autophagy, cell metabolism, and differentiation in development, circadian rhythm, and the immune response as well as neurodegeneration and cancer. Given its involvement in diverse cellular, physiological, and pathological processes, CK1α is a promising therapeutic target. In this review, we summarize what is known of the biological functions of CK1α, and provide an overview of existing challenges and potential opportunities for advancing theranostics.

From Otic Induction to Hair Cell Production: Pax2EGFP Cell Line Illuminates Key Stages of Development in Mouse Inner Ear Organoid Model

Producing hair cells of the inner ear is the major goal of ongoing research that combines advances in developmental and stem cell biology. The recent advent of an inner ear organoid protocol-resulting in three-dimensional stem cell-derived tissues resembling vestibular sensory epithelia-has sparked interest in applications such as regeneration, drug discovery, and disease modeling. In this study, we adapted this protocol for a novel mouse embryonic stem cell line with a fluorescent reporter for Pax2 expression. We used Pax2^EGFP/+ organoid formation to model otic induction, the pivotal developmental event when preplacodal tissue adopts otic fate. We found upregulation of Pax2 and activation of ERK downstream of fibroblast growth factor signaling in organoid formation as in embryonic inner ear development. Pax2 expression was evident from the EGFP reporter beginning at the vesicle formation stage and persisting through generation of the sensory epithelium. The native ventralizing signal sonic hedgehog was largely absent from the cell aggregates as otic vesicles began to form, confirming the dorsal vestibular organoid fate. Nonetheless, cochlear- or vestibular-like neurons appeared to delaminate from the derived otic vesicles and formed synaptic contacts with hair cells in the organoids. Cell lines with transcriptional reporters such as Pax2^EGFP/+ facilitate direct evaluation of morphological changes during organoid production, a major asset when establishing and validating the culture protocol.

Chemical-Induced Naive Pluripotency

The capacity to harness the properties of pluripotent stem cells has a wide-reaching impact on regenerative medicine, drug discovery, and also basic science. Two recent publications by Ursu et al. and Illich et al. demonstrate that inhibition of a CK1 isoform efficiently induces naive pluripotency in epiblast stem cells.

Distinct Signaling Requirements for the Establishment of ESC Pluripotency in Late-Stage EpiSCs

It has previously been reported that mouse epiblast stem cell (EpiSC) lines comprise heterogeneous cell populations that are functionally equivalent to cells of either early- or late-stage postimplantation development. So far, the establishment of the embryonic stem cell (ESC) pluripotency gene regulatory network through the widely known chemical inhibition of MEK and GSK3beta has been impractical in late-stage EpiSCs. Here, we show that chemical inhibition of casein kinase 1alpha (CK1alpha) induces the conversion of recalcitrant late-stage EpiSCs into ESC pluripotency. CK1alpha inhibition directly results in the simultaneous activation of the WNT signaling pathway, together with inhibition of the TGFbeta/SMAD2 signaling pathway, mediating the rewiring of the gene regulatory network in favor of an ESC-like state. Our findings uncover a molecular mechanism that links CK1alpha to ESC pluripotency through the direct modulation of WNT and TGFbeta signaling.

•

Inhibition of CK1alpha induces ESC conversion in EpiSCs recalcitrant to 2i/LIF

•

The ESC conversion acts via WNT activation and TGFbeta/SMAD2 inhibition

•

MEK inhibition stabilizes the conversion and restores germline competence

•

CK1 inhibition promotes activation and maintenance of the pluripotency network