PIP4K2C inhibition reverses autophagic flux impairment induced by SARS-CoV-2.

In search for broad-spectrum antivirals, we discovered a small molecule inhibitor, RMC-113, that potently suppresses the replication of multiple RNA viruses including SARS-CoV-2 in human lung organoids. We demonstrated selective dual inhibition of the lipid kinases PIP4K2C and PIKfyve by RMC-113 and target engagement by its clickable analog. Advanced lipidomics revealed alteration of SARS-CoV-2-induced phosphoinositide signature by RMC-113 and linked its antiviral effect with functional PIP4K2C and PIKfyve inhibition. We discovered PIP4K2C's roles in SARS-CoV-2 entry, RNA replication, and assembly/egress, validating it as a druggable antiviral target. Integrating proteomics, single-cell transcriptomics, and functional assays revealed that PIP4K2C binds SARS-CoV-2 nonstructural protein 6 and regulates virus-induced impairment of autophagic flux. Reversing this autophagic flux impairment is a mechanism of antiviral action of RMC-113. These findings reveal virus-induced autophagy regulation via PIP4K2C, an understudied kinase, and propose dual inhibition of PIP4K2C and PIKfyve as a candidate strategy to combat emerging viruses.

In Silico Evaluation of the Thr58-Associated Conserved Water with KRAS Switch-II Pocket Binders.

The KRAS switch-II pocket (SII-P) has proven to be one of the most successful tools for targeting KRAS with small molecules to date. This has been demonstrated with several KRAS(G12C)-targeting covalent inhibitors, already resulting in two FDA-approved drugs. Several earlier-stage compounds have also been reported to engage KRAS SII-P with other position 12 mutants, including G12D, G12S, and G12R. A highly conserved water molecule exists in the KRAS SII-P, linking Thr58 of switch-II and Gly10 of ß1 sheet. This conserved water is also present in the cocrystal structures of most of the disclosed small-molecule inhibitors but is only displaced by a handful of SII-P binders. Here, we evaluated the conserved water molecule energetics by the WaterMap for the SII-P binders with publicly disclosed structures and studied the water behavior in the presence of selected inhibitors by microsecond timescale molecular dynamics (MD) simulations using two water models (total simulation time of 120 µs). Our data revealed the high-energy nature of this hydration site when coexisting with an SII-P binder and that there is a preference for a single isolated hydration site in this location within the most advanced compounds. Furthermore, water displacement was only achieved with a few disclosed compounds and was suboptimal, as for instance a cyanomethyl group as a water displacer appears to introduce repulsion with the native conformation of Thr58. These results suggested that this conserved water should be considered more central when designing new inhibitors, especially in the design of noncovalent inhibitors targeting the SII-P.

Discovery and Development of First-in-Class ACKR3/CXCR7 Superagonists for Platelet Degranulation Modulation.

The atypical chemokine receptor 3 (ACKR3), formerly known as CXC-chemokine receptor 7 (CXCR7), has been postulated to regulate platelet function and thrombus formation. Herein, we report the discovery and development of first-in-class ACKR3 agonists, which demonstrated superagonistic properties with Emax values of up to 160% compared to the endogenous reference ligand CXCL12 in a ß-arrestin recruitment assay. Initial in silico screening using an ACKR3 homology model identified two hits, C10 (EC50 19.1 µM) and C11 (EC50 = 11.4 µM). Based on these hits, extensive structure-activity relationship studies were conducted by synthesis and testing of derivatives. It resulted in the identification of the novel thiadiazolopyrimidinone-based compounds 26 (LN5972, EC50 = 3.4 µM) and 27 (LN6023, EC50 = 3.5 µM). These compounds are selective for ACKR3 versus CXCR4 and show metabolic stability. In a platelet degranulation assay, these agonists effectively reduced P-selectin expression by up to 97%, suggesting potential candidates for the treatment of platelet-mediated thrombosis.

Discrepancy in interactions and conformational dynamics of pregnane X receptor (PXR) bound to an agonist and a novel competitive antagonist.

Pregnane X receptor (PXR) is a nuclear receptor with an essential role in regulating drug metabolism genes. While the mechanism of action for ligand-mediated PXR agonism is well-examined, its ligand-mediated inhibition or antagonism is poorly understood. Here we employ microsecond timescale all-atom molecular dynamics (MD) simulations to investigate how our newly identified dual kinase and PXR inhibitor, compound 100, acts as a competitive PXR antagonist and not as a full agonist. We study the PXR ligand binding domain conformational changes associated with compound 100 and compare the results to the full agonist SR12813, in presence and absence of the coactivator. Furthermore, we complement our research by experimentally disclosing the effect of eight key-residue mutations on PXR activation. Finally, simulations of P2X4 inhibitor (BAY-1797) in complex with PXR, which shares an identical structural moiety with compound 100, provide further insights to ligand-induced PXR behaviour. Our MD data suggests ligand-specific influence on conformations of different PXR-LBD regions, including α6 region, αAF-2, α1-α2', ß1'-α3 and ß1-ß1' loop. Our results provide important insights on conformational behaviour of PXR and offers guidance how to alleviate PXR agonism or to promote PXR antagonism.

Target Hopping from Protein Kinases to PXR: Identification of Small-Molecule Protein Kinase Inhibitors as Selective Modulators of Pregnane X Receptor from TüKIC Library.

Small-molecule protein kinase inhibitors are used for the treatment of cancer, but off-target effects hinder their clinical use. Especially off-target activation of the pregnane X receptor (PXR) has to be considered, as it not only governs drug metabolism and elimination, but also can promote tumor growth and cancer drug resistance. Consequently, PXR antagonism has been proposed for improving cancer drug therapy. Here we aimed to identify small-molecule kinase inhibitors of the Tübingen Kinase Inhibitor Collection (TüKIC) compound library that would act also as PXR antagonists. By a combination of in silico screen and confirmatory cellular reporter gene assays, we identified four novel PXR antagonists and a structurally related agonist with a common phenylaminobenzosuberone scaffold. Further characterization using biochemical ligand binding and cellular protein interaction assays classified the novel compounds as mixed competitive/noncompetitive, passive antagonists, which bind PXR directly and disrupt its interaction with coregulatory proteins. Expression analysis of prototypical PXR target genes ABCB1 and CYP3A4 in LS174T colorectal cancer cells and HepaRG hepatocytes revealed novel antagonists as selective receptor modulators, which showed gene- and tissue-specific effects. These results demonstrate the possibility of dual PXR and protein kinase inhibitors, which might represent added value in cancer therapy.

Regulatory spine RS3 residue of protein kinases: a lipophilic bystander or a decisive element in the small-molecule kinase inhibitor binding?

In recent years, protein kinases have been one of the most pursued drug targets. These determined efforts have resulted in ever increasing numbers of small-molecule kinase inhibitors reaching to the market, offering novel treatment options for patients with distinct diseases. One essential component related to the activation and normal functionality of a protein kinase is the regulatory spine (R-spine). The R-spine is formed of four conserved residues named as RS1-RS4. One of these residues, RS3, located in the C-terminal part of αC-helix, is usually accessible for the inhibitors from the ATP-binding cavity as its side chain is lining the hydrophobic back pocket in many protein kinases. Although the role of RS3 has been well acknowledged in protein kinase function, this residue has not been actively considered in inhibitor design, even though many small-molecule kinase inhibitors display interactions to this residue. In this minireview, we will cover the current knowledge of RS3, its relationship with the gatekeeper, and the role of RS3 in kinase inhibitor interactions. Finally, we comment on the future perspectives how this residue could be utilized in the kinase inhibitor design.

Decisive role of water and protein dynamics in residence time of p38α MAP kinase inhibitors.

Target residence time plays a crucial role in the pharmacological activity of small molecule inhibitors. Little is known, however, about the underlying causes of inhibitor residence time at the molecular level, which complicates drug optimization processes. Here, we employ all-atom molecular dynamics simulations (~400 µs in total) to gain insight into the binding modes of two structurally similar p38α MAPK inhibitors (type I and type I½) with short and long residence times that otherwise show nearly identical inhibitory activities in the low nanomolar IC50 range. Our results highlight the importance of protein conformational stability and solvent exposure, buried surface area of the ligand and binding site resolvation energy for residence time. These findings are further confirmed by simulations with a structurally diverse short residence time inhibitor SB203580. In summary, our data provide guidance in compound design when aiming for inhibitors with improved target residence time.

Addressing a Trapped High-Energy Water: Design and Synthesis of Highly Potent Pyrimidoindole-Based Glycogen Synthase Kinase-3ß Inhibitors.

In small molecule binding, water is not a passive bystander but rather takes an active role in the binding site, which may be decisive for the potency of the inhibitor. Here, by addressing a high-energy water, we improved the IC50 value of our co-crystallized glycogen synthase kinase-3ß (GSK-3ß) inhibitor by nearly two orders of magnitude. Surprisingly, our results demonstrate that this high-energy water was not displaced by our potent inhibitor (S)-3-(3-((7-ethynyl-9H-pyrimido[4,5-b]indol-4-yl)(methyl)amino)piperidin-1-yl)propanenitrile ((S)-15, IC50 value of 6 nM). Instead, only a subtle shift in the location of this water molecule resulted in a dramatic decrease in the energy of this high-energy hydration site, as shown by the WaterMap analysis combined with microsecond timescale molecular dynamics simulations. (S)-15 demonstrated both a favorable kinome selectivity profile and target engagement in a cellular environment and reduced GSK-3 autophosphorylation in neuronal SH-SY5Y cells. Overall, our findings highlight that even a slight adjustment in the location of a high-energy water can be decisive for ligand binding.

Identification and characterization of novel splice variants of human farnesoid X receptor.

Farnesoid X receptor (FXR, NR1H4) is a ligand-activated nuclear receptor, which regulates bile acid, lipid and glucose metabolism. Due to these functions, FXR has been investigated as a potential drug target for the treatment of liver diseases, such as primary biliary cholangitis and non-alcoholic steatohepatitis. Based on the previously described four splice variants, it has been suggested that alternative promoter usage and splicing may have an impact on total FXR activity as a result of encoding functionally diverse variants. Here we aimed for a systematic analysis of human hepatic FXR splice variants. In addition to the previously described FXRα1-4, we identified four novel splice variants (FXRα5-8) in human hepatocytes, which resulted from previously undetected exon skipping events. These newly identified isoforms displayed diminished DNA binding and impaired transactivation activities. Isoform FXRα5, which suppressed the transactivation activity of the functional isoform FXRα2, was further characterized as deficient in heterodimerization, coactivator recruitment and ligand binding. These findings were further supported by molecular dynamics simulations, which offered an explanation for the behavior of this isoform on the molecular level. FXRα5 exhibited low uniform expression levels in nearly all human tissues. Our systematic analysis of FXR splice variants in human hepatocytes resulted in the identification of four novel FXR isoforms, which all proved to be functionally deficient, but one novel variant, FXRα5, also displayed dominant negative activity. The possible associations with and roles of these novel isoforms in human liver diseases require further investigation.

Design and synthesis of novel fluorescently labeled analogs of vemurafenib targeting MKK4.

The mitogen-activated protein kinase kinase 4 (MKK4) plays a key role in liver regeneration and is under investigation as a target for stimulating hepatocytes to increased proliferation. Therefore, new small molecules inhibiting MKK4 may represent a promising approach for treating acute and chronic liver diseases. Fluorescently labeled compounds are useful tools for high-throughput screenings of large compound libraries. Here we utilized the azaindole-based scaffold of FDA-approved BRAF inhibitor vemurafenib 1, which displays off-target activity on MKK4, as a starting point in our fluorescent compound design. Chemical variation of the scaffold and optimization led to a selection of fluorescent 5-TAMRA derivatives which possess high binding affinities on MKK4. Compound 45 represents a suitable tool compound for Fluorescence polarization assays to identify new small-molecule inhibitors of MKK4.

The autoinhibited state of MKK4: Phosphorylation, putative dimerization and R134W mutant studied by molecular dynamics simulations.

Protein kinases are crucial components of the cell-signalling machinery that orchestrate and convey messages to their downstream targets. Most often, kinases are activated upon a phosphorylation to their activation loop, which will shift the kinase into the active conformation. The Dual specificity mitogen-activated protein kinase kinase 4 (MKK4) exists in a unique conformation in its inactive unphosphorylated state, where its activation segment appears in a stable α-helical conformation. However, the precise role of this unique conformational state of MKK4 is unknown. Here, by all-atom molecular dynamics simulations (MD simulations), we show that this inactive state is unstable as monomer even when unphosphorylated and that the phosphorylation of the activation segment further destabilizes the autoinhibited α-helix. The specific phosphorylation pattern of the activation segment has also a unique influence on MKK4 dynamics. Furthermore, we observed that this specific inactive state is stable as a dimer, which becomes destabilized upon phosphorylation. Finally, we noticed that the most frequent MKK4 mutation observed in cancer, R134W, which role has not been disclosed to date, contributes to the dimer stability. Based on these data we postulate that MKK4 occurs as a dimer in its inactive autoinhibited state, providing an additional layer for its activity regulation.

Discovery and Evaluation of Enantiopure 9H-pyrimido[4,5-b]indoles as Nanomolar GSK-3ß Inhibitors with Improved Metabolic Stability.

Glycogen synthase kinase-3ß (GSK-3ß) is a potential target in the field of Alzheimer's disease drug discovery. We recently reported a new class of 9H-pyrimido[4,5-b]indole-based GSK-3ß inhibitors, of which 3-(3-((7-chloro-9H-pyrimido[4,5-b]indol-4-yl)(methyl)amino)piperidin-1-yl)propanenitrile (1) demonstrated promising inhibitory potency. However, this compound underwent rapid degradation by human liver microsomes. Starting from 1, we prepared a series of amide-based derivatives and studied their structure-activity relationships against GSK-3ß supported by 1 µs molecular dynamics simulations. The biological potency of this series was substantially enhanced by identifying the eutomer configuration at the stereocenter. Moreover, the introduction of an amide bond proved to be an effective strategy to eliminate the metabolic hotspot. The most potent compounds, (R)-3-(3-((7-chloro-9H-pyrimido[4,5-b]indol-4-yl)(methyl)amino)piperidin-1-yl)-3-oxopropanenitrile ((R)-2) and (R)-1-(3-((7-bromo-9Hpyrimido[4,5-b]indol-4-yl)(methyl)amino)piperidin-1-yl)propan-1-one ((R)-28), exhibited IC50 values of 480 nM and 360 nM, respectively, and displayed improved metabolic stability. Their favorable biological profile is complemented by minimal cytotoxicity and neuroprotective properties.

KRAS(G12C)-AMG 510 interaction dynamics revealed by all-atom molecular dynamics simulations.

The first KRAS(G12C) targeting inhibitor in clinical development, AMG 510, has shown promising antitumor activity in clinical trials. On the molecular level, however, the interaction dynamics of this covalently bound drug-protein complex has been undetermined. Here, we disclose the interaction dynamics of the KRAS(G12C)-AMG 510 complex by long timescale all-atom molecular dynamics (MD) simulations (total of 75 µs). Moreover, we investigated the influence of the recently reported post-translational modification (PTM) of KRAS' N-terminus, removal of initiator methionine (iMet1) with acetylation of Thr2, to this complex. Our results demonstrate that AMG 510 does not entrap KRAS into a single conformation, as one would expect based on the crystal structure, but rather into an ensemble of conformations. AMG 510 binding is extremely stable regardless of highly dynamic interface of KRAS' switches. Overall, KRAS(G12C)-AMG 510 complex partially mimic the native dynamics of GDP bound KRAS; however, AMG 510 stabilizes the α3-helix region. N-terminally modified KRAS displays similar interaction dynamics with AMG 510 as when Met1 is present, but this PTM appears to stabilize ß2-ß3-loop. These results provide novel conformational insights on the molecular level to KRAS(G12C)-AMG 510 interactions and dynamics, providing new perspectives to RAS related drug discovery.

Bioisosteric Replacement of Arylamide-Linked Spine Residues with N-Acylhydrazones and Selenophenes as a Design Strategy to Novel Dibenzosuberone Derivatives as Type I 1/2 p38α MAP Kinase Inhibitors.

The recent disclosure of type I 1/2 inhibitors for p38α MAPK demonstrated how the stabilization of the R-spine can be used as a strategy to greatly increase the target residence time (TRT) of inhibitors. Herein, for the first time, we describe N-acylhydrazone and selenophene residues as spine motifs, yielding metabolically stable inhibitors with high potency on enzymatic, NanoBRET, and whole blood assays, improved metabolic stability, and prolonged TRT.

The current understanding of KRAS protein structure and dynamics.

One of the most common drivers in human cancer is the mutant KRAS protein. Not so long ago KRAS was considered as an undruggable oncoprotein. After a long struggle, however, we finally see some light at the end of the tunnel as promising KRAS targeted therapies are in or approaching clinical trials. In recent years, together with the promising progress in RAS drug discovery, our understanding of KRAS has increased tremendously. This progress has been accompanied with a resurgence of publicly available KRAS structures, which were limited to nine structures less than ten years ago. Furthermore, the ever-increasing computational capacity has made biologically relevant timescales accessible, enabling molecular dynamics (MD) simulations to study the dynamics of KRAS protein in more detail at the atomistic level. In this minireview, my aim is to provide the reader an overview of the publicly available KRAS structural data, insights to conformational dynamics revealed by experiments and what we have learned from MD simulations. Also, I will discuss limitations of the current data and provide suggestions for future research related to KRAS, which would fill out the existing gaps in our knowledge and provide guidance in deciphering this enigmatic oncoprotein.

Pyridinylimidazoles as GSK3ß Inhibitors: The Impact of Tautomerism on Compound Activity via Water Networks.

Glycogen synthase kinase-3ß (GSK3ß) is involved in many pathological conditions and represents an attractive drug target. We previously reported dual GSK3ß/p38α mitogen-activated protein kinase inhibitors and identified N-(4-(4-(4-fluorophenyl)-2-methyl-1H-imidazol-5-yl)pyridin-2-yl)cyclopropanecarboxamide (1) as a potent dual inhibitor of both target kinases. In this study, we aimed to design selective GSK3ß inhibitors based on our pyridinylimidazole scaffold. Our efforts resulted in several novel and potent GSK3ß inhibitors with IC50 values in the low nanomolar range. 5-(2-(Cyclopropanecarboxamido)pyridin-4-yl)-4-cyclopropyl-1H-imidazole-2-carboxamide (6g) displayed very good kinase selectivity as well as metabolical stability and inhibited GSK3ß activity in neuronal SH-SY5Y cells. Interestingly, we observed the importance of the 2-methylimidazole's tautomeric state for the compound activity. Finally, we reveal how this crucial tautomerism effect is surmounted by imidazole-2-carboxamides, which are able to stabilize the binding via enhanced water network interactions, regardless of their tautomeric state.

Design, Synthesis and Biological Evaluation of 7-Chloro-9H-pyrimido[4,5-b]indole-based Glycogen Synthase Kinase-3ß Inhibitors.

Glycogen synthase kinase-3ß (GSK-3ß) represents a relevant drug target for the treatment of neurodegenerative pathologies including Alzheimer's disease. We herein report on the optimization of a novel class of GSK-3ß inhibitors based on the tofacitinib-derived screen hit 3-((3R,4R)-3-((7-chloro-9H-pyrimido[4,5-b]indol-4-yl)(methyl)amino)-4-methylpiperidin-1-yl)-3-oxopropanenitrile (1). We synthesized a series of 19 novel 7-chloro-9H-pyrimido[4,5-b]indole-based derivatives and studied their structure-activity relationships with focus on the cyanoacetyl piperidine moiety. We unveiled the crucial role of the nitrile group and its importance for the activity of this compound series. A successful rigidization approach afforded 3-(3aRS,7aSR)-(1-(7-chloro-9H-pyrimido[4,5-b]indol-4-yl)octahydro-6H-pyrrolo[2,3-c]pyridin-6-yl)-propanenitrile (24), which displayed an IC50 value of 130 nM on GSK-3ß and was further characterized by its metabolic stability. Finally, we disclosed the putative binding modes of the most potent inhibitors within the ATP binding site of GSK-3ß by 1 µs molecular dynamics simulations.

Pyridinylimidazoles as dual glycogen synthase kinase 3ß/p38α mitogen-activated protein kinase inhibitors.

Compounds simultaneously inhibiting two targets that are involved in the progression of the same complex disease may exhibit additive or even synergistic therapeutic effects. Here we unveil 2,4,5-trisubstituted imidazoles as dual inhibitors of p38α mitogen-activated protein kinase and glycogen synthase kinase 3ß (GSK3ß). Both enzymes are potential therapeutic targets for neurodegenerative disorders, like Alzheimer's disease. A set of 39 compounds was synthesized and evaluated in kinase activity assays for their ability to inhibit both target kinases. Among the synthesized compounds, potent dual-target-directed inhibitors showing IC50 values down to the low double-digit nanomolar range, were identified. One of the best balanced dual inhibitors presented in here is N-(4-(2-ethyl-4-(4-fluorophenyl)-1H-imidazol-5-yl)pyridin-2-yl)cyclopropanecarboxamide (20c) (p38α, IC50 = 16 nM; GSK3ß, IC50 = 35 nM) featuring an excellent metabolic stability and an appreciable isoform selectivity over the closely related GSK3α. Our findings were rationalized by computational docking studies based on previously published X-ray structures.

Assessment of mutation probabilities of KRAS G12 missense mutants and their long-timescale dynamics by atomistic molecular simulations and Markov state modeling.

A mutated KRAS protein is frequently observed in human cancers. Traditionally, the oncogenic properties of KRAS missense mutants at position 12 (G12X) have been considered as equal. Here, by assessing the probabilities of occurrence of all KRAS G12X mutations and KRAS dynamics we show that this assumption does not hold true. Instead, our findings revealed an outstanding mutational bias. We conducted a thorough mutational analysis of KRAS G12X mutations and assessed to what extent the observed mutation frequencies follow a random distribution. Unique tissue-specific frequencies are displayed with specific mutations, especially with G12R, which cannot be explained by random probabilities. To clarify the underlying causes for the nonrandom probabilities, we conducted extensive atomistic molecular dynamics simulations (170 µs) to study the differences of G12X mutations on a molecular level. The simulations revealed an allosteric hydrophobic signaling network in KRAS, and that protein dynamics is altered among the G12X mutants and as such differs from the wild-type and is mutation-specific. The shift in long-timescale conformational dynamics was confirmed with Markov state modeling. A G12X mutation was found to modify KRAS dynamics in an allosteric way, which is especially manifested in the switch regions that are responsible for the effector protein binding. The findings provide a basis to understand better the oncogenic properties of KRAS G12X mutants and the consequences of the observed nonrandom frequencies of specific G12X mutations.

Binding Affinity via Docking: Fact and Fiction.

In 1982, Kuntz et al. published an article with the title "A Geometric Approach to Macromolecule-Ligand Interactions", where they described a method "to explore geometrically feasible alignment of ligands and receptors of known structure". Since then, small molecule docking has been employed as a fast way to estimate the binding pose of a given compound within a specific target protein and also to predict binding affinity. Remarkably, the first docking method suggested by Kuntz and colleagues aimed to predict binding poses but very little was specified about binding affinity. This raises the question as to whether docking is the right tool to estimate binding affinity. The short answer is no, and this has been concluded in several comprehensive analyses. However, in this opinion paper we discuss several critical aspects that need to be reconsidered before a reliable binding affinity prediction through docking is realistic. These are not the only issues that need to be considered, but they are perhaps the most critical ones. We also consider that in spite of the huge efforts to enhance scoring functions, the accuracy of binding affinity predictions is perhaps only as good as it was 10⁻20 years ago. There are several underlying reasons for this poor performance and these are analyzed. In particular, we focus on the role of the solvent (water), the poor description of H-bonding and the lack of the systems' true dynamics. We hope to provide readers with potential insights and tools to overcome the challenging issues related to binding affinity prediction via docking.