MISATO: machine learning dataset of protein-ligand complexes for structure-based drug discovery

Large language models have greatly enhanced our ability to understand biology and chemistry, yet robust methods for structure-based drug discovery, quantum chemistry and structural biology are still sparse. Precise biomolecule-ligand interaction datasets are urgently needed for large language models. To address this, we present MISATO, a dataset that combines quantum mechanical properties of small molecules and associated molecular dynamics simulations of ~20,000 experimental protein-ligand complexes with extensive validation of experimental data. Starting from the existing experimental structures, semi-empirical quantum mechanics was used to systematically refine these structures. A large collection of molecular dynamics traces of protein-ligand complexes in explicit water is included, accumulating over 170 μs. We give examples of machine learning (ML) baseline models proving an improvement of accuracy by employing our data. An easy entry point for ML experts is provided to enable the next generation of drug discovery artificial intelligence models.

SARS-CoV-2 Mpro oligomerization as a potential target for therapy

The main protease (Mpro) of SARS-CoV-2 is critical in the virus's replication cycle, facilitating the maturation of polyproteins into functional units. Due to its conservation across taxa, Mpro is a promising target for broad-spectrum antiviral drugs. Targeting Mpro with small molecule inhibitors, such as nirmatrelvir combined with ritonavir (Paxlovid™), which the FDA has approved for post-exposure treatment and prophylaxis, can effectively interrupt the replication process of the virus. A key aspect of Mpro's function is its ability to form a functional dimer. However, the mechanics of dimerization and its influence on proteolytic activity remain less understood. In this study, we utilized biochemical, structural, and molecular modelling approaches to explore Mpro dimerization. We evaluated critical residues, specifically Arg4 and Arg298, that are essential for dimerization. Our results show that changes in the oligomerization state of Mpro directly affect its enzymatic activity and dimerization propensity. We discovered a synergistic relationship influencing dimer formation, involving both intra- and intermolecular interactions. These findings highlight the potential for developing allosteric inhibitors targeting Mpro, offering promising new directions for therapeutic strategies.

Modulation of peroxisomal import by the PEX13 SH3 domain and a proximal FxxxF binding motif

Import of proteins into peroxisomes depends on PEX5, PEX13 and PEX14. By combining biochemical methods and structural biology, we show that the C-terminal SH3 domain of PEX13 mediates intramolecular interactions with a proximal FxxxF motif. The SH3 domain also binds WxxxF peptide motifs in the import receptor PEX5, demonstrating evolutionary conservation of such interactions from yeast to human. Strikingly, intramolecular interaction of the PEX13 FxxxF motif regulates binding of PEX5 WxxxF/Y motifs to the PEX13 SH3 domain. Crystal structures reveal how FxxxF and WxxxF/Y motifs are recognized by a non-canonical surface on the SH3 domain. The PEX13 FxxxF motif also mediates binding to PEX14. Surprisingly, the potential PxxP binding surface of the SH3 domain does not recognize PEX14 PxxP motifs, distinct from its yeast ortholog. Our data show that the dynamic network of PEX13 interactions with PEX5 and PEX14, mediated by diaromatic peptide motifs, modulates peroxisomal matrix import.

Despite the odds: formation of the SARS-CoV-2 methylation complex

Coronaviruses modify their single-stranded RNA genome with a methylated cap during replication to mimic the eukaryotic mRNAs. The capping process is initiated by several nonstructural proteins (nsp) encoded in the viral genome. The methylation is performed by two methyltransferases, nsp14 and nsp16, while nsp10 acts as a co-factor to both. Additionally, nsp14 carries an exonuclease domain which operates in the proofreading system during RNA replication of the viral genome. Both nsp14 and nsp16 were reported to independently bind nsp10, but the available structural information suggests that the concomitant interaction between these three proteins would be impossible due to steric clashes. Here, we show that nsp14, nsp10, and nsp16 can form a heterotrimer complex upon significant allosteric change. This interaction is expected to encourage the formation of mature capped viral mRNA, modulating nsp14's exonuclease activity, and protecting the viral RNA. Our findings show that nsp14 is amenable to allosteric regulation and may serve as a novel target for therapeutic approaches.

Epstein-Barr virus-driven B cell lymphoma mediated by a direct LMP1-TRAF6 complex

Epstein-Barr virus (EBV) latent membrane protein 1 (LMP1) drives viral B cell transformation and oncogenesis. LMP1's transforming activity depends on its C-terminal activation region 2 (CTAR2), which induces NF-κB and JNK by engaging TNF receptor-associated factor 6 (TRAF6). The mechanism of TRAF6 recruitment to LMP1 and its role in LMP1 signalling remains elusive. Here we demonstrate that TRAF6 interacts directly with a viral TRAF6 binding motif within CTAR2. Functional and NMR studies supported by molecular modeling provide insight into the architecture of the LMP1-TRAF6 complex, which differs from that of CD40-TRAF6. The direct recruitment of TRAF6 to LMP1 is essential for NF-κB activation by CTAR2 and the survival of LMP1-driven lymphoma. Disruption of the LMP1-TRAF6 complex by inhibitory peptides interferes with the survival of EBV-transformed B cells. In this work, we identify LMP1-TRAF6 as a critical virus-host interface and validate this interaction as a potential therapeutic target in EBV-associated cancer.

Discovery of Inhibitory Fragments That Selectively Target Spire2-FMN2 Interaction

Here, we report the fragment-based drug discovery of potent and selective fragments that disrupt the Spire2-FMN2 but not the Spire1-FMN2 interaction. Hit fragments were identified in a differential scanning fluorimetry-based screen of an in-house library of 755 compounds and subsequently validated in multiple orthogonal biophysical assays, including fluorescence polarization, microscale thermophoresis, and 1H-15N HSQC nuclear magnetic resonance. Extensive structure-activity relationships combined with molecular docking followed by chemical optimization led to the discovery of compound 13, which exhibits micromolar potency and high ligand efficiency (LE = 0.38). Therefore, this fragment represents a validated starting point for the future development of selective chemical probes targeting the Spire2-FMN2 interaction.

Binding mechanism and biological effects of flavone DYRK1A inhibitors for the design of new antidiabetics

The selective inhibition of kinases from the diabetic kinome is known to promote the regeneration of beta cells and provide an opportunity for the curative treatment of diabetes. The effect can be achieved by carefully tailoring the selectivity of inhibitor toward a particular kinase, especially DYRK1A, previously associated with Down syndrome and Alzheimer's disease. Recently DYRK1A inhibition has been shown to promote both insulin secretion and beta cells proliferation. Here, we show that commonly available flavones are effective inhibitors of DYRK1A. The observed biochemical activity of flavone compounds is confirmed by crystal structures solved at 2.06 Å and 2.32 Å resolution, deciphering the way inhibitors bind in the ATP-binding pocket of the kinase, which is driven by the arrangement of hydroxyl moieties. We also demonstrate antidiabetic properties of these biomolecules and prove that they could be further improved by therapy combined with TGF-β inhibitors. Our data will allow future structure-based optimization of the presented scaffolds toward potent, bioavailable and selective anti-diabetic drugs.

Development of novel PEX5-PEX14 protein-protein interaction (PPI) inhibitors based on an oxopiperazine template

Protein-protein interactions (PPIs) constitute an important but challenging class of molecular targets for small molecules. The PEX5-PEX14 PPI has been shown to play a critical role in glycosome biogenesis and its disruption impairs the metabolism in Trpanosoma parasites, eventually leading to their death. Therefore, this PPI is a potential molecular target for new drugs against diseases caused by Trypanosoma infections. Here, we report a new class of peptidomimetic scaffolds to target the PEX5-PEX14 PPI. The molecular design was based on an oxopiperazine template for the α-helical mimetics. A structural simplification along with modifications of the central oxopiperazine scaffold and addressing the lipophilic interactions led to the development of peptidomimetics that inhibit PEX5-TbPEX14 PPI and display cellular activity against T. b. brucei. This approach provides an alternative approach towards the development of trypanocidal agents and may be generally useful for the design of helical mimetics as PPI inhibitors.

Extrinsic stabilization of antiviral ACE2-Fc fusion proteins targeting SARS-CoV-2

The angiotensin-converting enzyme 2 (ACE2) is a viral receptor used by sarbecoviruses to infect cells. Fusion proteins comprising extracellular ACE2 domains and the Fc part of immunoglobulins exhibit high virus neutralization efficiency, but the structure and stability of these molecules are poorly understood. We show that although the hinge between the ACE2 and the IgG4-Fc is highly flexible, the conformational dynamics of the two ACE2 domains is restricted by their association. Interestingly, the conformational stability of the ACE2 moiety is much lower than that of the Fc part. We found that chemical compounds binding to ACE2, such as DX600 and MLN4760, can be used to strongly increase the thermal stability of the ACE2 by different mechanisms. Together, our findings reveal a general concept for stabilizing the labile receptor segments of therapeutic antiviral fusion proteins by chemical compounds.

Silmitasertib (CX-4945), a Clinically Used CK2-Kinase Inhibitor with Additional Effects on GSK3β and DYRK1A Kinases: A Structural Perspective

A clinical casein kinase 2 inhibitor, CX-4945 (silmitasertib), shows significant affinity toward the DYRK1A and GSK3β kinases, involved in down syndrome phenotypes, Alzheimer's disease, circadian clock regulation, and diabetes. This off-target activity offers an opportunity for studying the effect of the DYRK1A/GSK3β kinase system in disease biology and possible line extension. Motivated by the dual inhibition of these kinases, we solved and analyzed the crystal structures of DYRK1A and GSK3β with CX-4945. We built a quantum-chemistry-based model to rationalize the compound affinity for CK2α, DYRK1A, and GSK3β kinases. Our calculations identified a key element for CK2α's subnanomolar affinity to CX-4945. The methodology is expandable to other kinase selectivity modeling. We show that the inhibitor limits DYRK1A- and GSK3β-mediated cyclin D1 phosphorylation and reduces kinase-mediated NFAT signaling in the cell. Given the CX-4945's clinical and pharmacological profile, this inhibitory activity makes it an interesting candidate with potential for application in additional disease areas.

A Buckycatcher in Solution-A Computational Perspective

In this work, we study the buckycatcher (C₆₀H₂₈) in solution using quantum chemical models. We investigate the conformational equilibria in several media and the effects that molecules of solvent might have in interconversion barriers between the different conformers. These are studied in a hypothetical gas phase, in the dielectric of a solvent, as well as with hybrid solvation. In the latter case, due to a disruption of π-stacking interactions, the transition states are destabilized. We also evaluate the complexation of the buckycatcher with solvent-like molecules. In most cases studied, there should be no adducts formed because the enthalpy driving force cannot overcome entropic penalties.

Development of Noncovalent Small-Molecule Keap1-Nrf2 Inhibitors by Fragment-Based Drug Discovery

Targeting the protein-protein interaction (PPI) between the transcription factor nuclear factor erythroid 2-related factor 2 (Nrf2) and its repressor, Kelch-like ECH-associated protein 1 (Keap1), constitutes a promising strategy for treating diseases involving oxidative stress and inflammation. Here, a fragment-based drug discovery (FBDD) campaign resulted in novel, high-affinity (K_i = 280 nM), and cell-active noncovalent small-molecule Keap1-Nrf2 PPI inhibitors. We screened 2500 fragments using orthogonal assays─fluorescence polarization (FP), thermal shift assay (TSA), and surface plasmon resonance (SPR)─and validated the hits by saturation transfer difference (STD) NMR, leading to 28 high-priority hits. Thirteen co-structures showed fragments binding mainly in the P4 and P5 subpockets of Keap1's Kelch domain, and three fluorenone-based fragments featuring a novel binding mode were optimized by structure-based drug discovery. We thereby disclose several fragment hits, including their binding modes, and show how FBDD can be performed to find new small-molecule Keap1-Nrf2 PPI inhibitors.

Exploring the Surface of the Ectodomain of the PD-L1 Immune Checkpoint with Small-Molecule Fragments

Development of small molecules targeting the PD-L1/PD-1 interface is advancing both in industry and academia, but only a few have reached early-stage clinical trials. Here, we take a closer look at the general druggability of PD-L1 using in silico hot spot mapping and nuclear magnetic resonance (NMR)-based characterization. We found that the conformational elasticity of the PD-L1 surface strongly influences the formation of hot spots. We deconstructed several generations of known inhibitors into fragments and examined their binding properties using differential scanning fluorimetry (DSF) and protein-based nuclear magnetic resonance (NMR). These biophysical analyses showed that not all fragments bind to the PD-L1 ectodomain despite having the biphenyl scaffold. Although most of the binding fragments induced PD-L1 oligomerization, two compounds, TAH35 and TAH36, retain the monomeric state of proteins upon binding. Additionally, the presence of the entire ectodomain did not affect the binding of the hit compounds and dimerization of PD-L1. The data demonstrated here provide important information on the PD-L1 druggability and the structure-activity relationship of the biphenyl core moiety and therefore may aid in the design of novel inhibitors and focused fragment libraries for PD-L1.

ULYSSES: An Efficient and Easy to Use Semiempirical Library for C+

We introduce ULYSSES, a user-friendly and robust C++ library for semiempirical quantum chemical calculations. In its current version, ULYSSES is equipped with a large set of different semiempirical models, most of which are based on the Neglect of Diatomic Differential Overlap (NDDO) approximation. Empirical corrections for dispersion and hydrogen bonding are available for most methods, so that higher quality is achieved in the calculation of energies of nonbonded complexes. The library is furthermore equipped with geometry optimization, as well as modules for calculating molecular properties of general interest. Ideal gas thermodynamics is available and allows single structure as well as conformer (multistructure) averaged properties to be calculated. We offer the possibility to use several vibrational partition functions according to the nature of interactions being studied: for covalent systems, the traditional harmonic oscillator approximation is available; for nonbonded complexes, we systematically extended the partition function proposed by Grimme for all thermodynamic functions. The library is also capable of running Born-Oppenheimer molecular dynamics.

Characterization of SARS-CoV-2 replication complex elongation and proofreading activity

The replication complex (RC) of SARS-CoV-2 was recently shown to be one of the fastest RNA-dependent RNA polymerases of any known coronavirus. With this rapid elongation, the RC is more prone to incorporate mismatches during elongation, resulting in a highly variable genomic sequence. Such mutations render the design of viral protein targets difficult, as drugs optimized for a given viral protein sequence can quickly become inefficient as the genomic sequence evolves. Here, we use biochemical experiments to characterize features of RNA template recognition and elongation fidelity of the SARS-CoV-2 RdRp, and the role of the exonuclease, nsp14. Our study highlights the 2'OH group of the RNA ribose as a critical component for RdRp template recognition and elongation. We show that RdRp fidelity is reduced in the presence of the 3' deoxy-terminator nucleotide 3'dATP, which promotes the incorporation of mismatched nucleotides (leading to U:C, U:G, U:U, C:U, and A:C base pairs). We find that the nsp10-nsp14 heterodimer is unable to degrade RNA products lacking free 2'OH or 3'OH ribose groups. Our results suggest the potential use of 3' deoxy-terminator nucleotides in RNA-derived oligonucleotide inhibitors as antivirals against SARS-CoV-2.

Novel Trypanocidal Inhibitors that Block Glycosome Biogenesis by Targeting PEX3-PEX19 Interaction

Human pathogenic trypanosomatid parasites harbor a unique form of peroxisomes termed glycosomes that are essential for parasite viability. We and others previously identified and characterized the essential Trypanosoma brucei ortholog TbPEX3, which is the membrane-docking factor for the cytosolic receptor PEX19 bound to the glycosomal membrane proteins. Knockdown of TbPEX3 expression leads to mislocalization of glycosomal membrane and matrix proteins, and subsequent cell death. As an early step in glycosome biogenesis, the PEX3–PEX19 interaction is an attractive drug target. We established a high-throughput assay for TbPEX3–TbPEX19 interaction and screened a compound library for small-molecule inhibitors. Hits from the screen were further validated using an in vitro ELISA assay. We identified three compounds, which exhibit significant trypanocidal activity but show no apparent toxicity to human cells. Furthermore, we show that these compounds lead to mislocalization of glycosomal proteins, which is toxic to the trypanosomes. Moreover, NMR-based experiments indicate that the inhibitors bind to PEX3. The inhibitors interfering with glycosomal biogenesis by targeting the TbPEX3–TbPEX19 interaction serve as starting points for further optimization and anti-trypanosomal drug development.

Bis-choline tetrathiomolybdate prevents copper-induced blood-brain barrier damage

The blood–brain barrier endothelial cell monolayer becomes permeable to elevated copper loosely bound to albumin, which can be avoided by a high-affinity copper chelator but not by D-penicillamine.

In Wilson disease, excessive copper accumulates in patients’ livers and may, upon serum leakage, severely affect the brain according to current viewpoints. Present remedies aim at avoiding copper toxicity by chelation, for example, by D-penicillamine (DPA) or bis-choline tetrathiomolybdate (ALXN1840), the latter with a very high copper affinity. Hence, ALXN1840 may potentially avoid neurological deterioration that frequently occurs upon DPA treatment. As the etiology of such worsening is unclear, we reasoned that copper loosely bound to albumin, that is, mimicking a potential liver copper leakage into blood, may damage cells that constitute the blood-brain barrier, which was found to be the case in an in vitro model using primary porcine brain capillary endothelial cells. Such blood–brain barrier damage was avoided by ALXN1840, plausibly due to firm protein embedding of the chelator bound copper, but not by DPA. Mitochondrial protection was observed, a prerequisite for blood–brain barrier integrity. Thus, high-affinity copper chelators may minimize such deterioration in the treatment of neurologic Wilson disease.

Multifaceted N-Degron Recognition and Ubiquitylation by GID/CTLH E3 Ligases

N-degron E3 ubiquitin ligases recognize specific residues at the N-termini of substrates. Although molecular details of N-degron recognition are known for several E3 ligases, the range of N-terminal motifs that can bind a given E3 substrate binding domain remains unclear. Here, we discovered capacity of Gid4 and Gid10 substrate receptor subunits of yeast "GID"/human "CTLH" multiprotein E3 ligases to tightly bind a wide range of N-terminal residues whose recognition is determined in part by the downstream sequence context. Screening of phage displaying peptide libraries with exposed N-termini identified novel consensus motifs with non-Pro N-terminal residues binding Gid4 or Gid10 with high affinity. Structural data reveal that conformations of flexible loops in Gid4 and Gid10 complement sequences and folds of interacting peptides. Together with analysis of endogenous substrate degrons, the data show that degron identity, substrate domains harboring targeted lysines, and varying E3 ligase higher-order assemblies combinatorially determine efficiency of ubiquitylation and degradation.

Computer-Aided Design and Synthesis of a New Class of PEX14 Inhibitors: Substituted 2,3,4,5-Tetrahydrobenzo[F][1,4]oxazepines as Potential New Trypanocidal Agents

African and American trypanosomiases are estimated to affect several million people across the world, with effective treatments distinctly lacking. New, ideally oral, treatments with higher efficacy against these diseases are desperately needed. Peroxisomal import matrix (PEX) proteins represent a very interesting target for structure- and ligand-based drug design. The PEX5-PEX14 protein-protein interface in particular has been highlighted as a target, with inhibitors shown to disrupt essential cell processes in trypanosomes, leading to cell death. In this work, we present a drug development campaign that utilizes the synergy between structural biology, computer-aided drug design, and medicinal chemistry in the quest to discover and develop new potential compounds to treat trypanosomiasis by targeting the PEX14-PEX5 interaction. Using the structure of the known lead compounds discovered by Dawidowski et al. as the template for a chemically advanced template search (CATS) algorithm, we performed scaffold-hopping to obtain a new class of compounds with trypanocidal activity, based on 2,3,4,5-tetrahydrobenzo[f][1,4]oxazepines chemistry. The initial compounds obtained were taken forward to a first round of hit-to-lead optimization by synthesis of derivatives, which show activities in the range of low- to high-digit micromolar IC₅₀ in the in vitro tests. The NMR measurements confirm binding to PEX14 in solution, while immunofluorescent microscopy indicates disruption of protein import into the glycosomes, indicating that the PEX14-PEX5 protein-protein interface was successfully disrupted. These studies result in development of a novel scaffold for future lead optimization, while ADME testing gives an indication of further areas of improvement in the path from lead molecules toward a new drug active against trypanosomes.

A fragment-based approach identifies an allosteric pocket that impacts malate dehydrogenase activity

Malate dehydrogenases (MDHs) sustain tumor growth and carbon metabolism by pathogens including Plasmodium falciparum. However, clinical success of MDH inhibitors is absent, as current small molecule approaches targeting the active site are unselective. The presence of an allosteric binding site at oligomeric interface allows the development of more specific inhibitors. To this end we performed a differential NMR-based screening of 1500 fragments to identify fragments that bind at the oligomeric interface. Subsequent biophysical and biochemical experiments of an identified fragment indicate an allosteric mechanism of 4-(3,4-difluorophenyl) thiazol-2-amine (4DT) inhibition by impacting the formation of the active site loop, located >30 Å from the 4DT binding site. Further characterization of the more tractable homolog 4-phenylthiazol-2-amine (4PA) and 16 other derivatives are also reported. These data pave the way for downstream development of more selective molecules by utilizing the oligomeric interfaces showing higher species sequence divergence than the MDH active site.

Romero et al. perform NMR-based screening of 1500 fragments to identify fragments that bind at the oligomeric interface of malate dehydrogenase (MDH). Their study indicates an allosteric mechanism impacting enzymatic activity, paving the way for development of more selective molecules and a starting point for the future development of specific MDH inhibitors.

Latency, thermal stability, and identification of an inhibitory compound of mirolysin, a secretory protease of the human periodontopathogen Tannerella forsythia

Mirolysin is a secretory protease of Tannerella forsythia, a member of the dysbiotic oral microbiota responsible for periodontitis. In this study, we show that mirolysin latency is achieved by a “cysteine-switch” mechanism exerted by Cys23 in the N-terminal profragment. Mutation of Cys23 shortened the time needed for activation of the zymogen from several days to 5 min. The mutation also decreased the thermal stability and autoproteolysis resistance of promirolysin. Mature mirolysin is a thermophilic enzyme and shows optimal activity at 65 °C. Through NMR-based fragment screening, we identified a small molecule (compound (cpd) 9) that blocks promirolysin maturation and functions as a competitive inhibitor (K_i = 3.2 µM), binding to the S1′ subsite of the substrate-binding pocket. Cpd 9 shows superior specificity and does not interact with other T. forsythia proteases or Lys/Arg-specific proteases.

Deconstructing Noncovalent Kelch-like ECH-Associated Protein 1 (Keap1) Inhibitors into Fragments to Reconstruct New Potent Compounds

Targeting the protein-protein interaction (PPI) between nuclear factor erythroid 2-related factor 2 (Nrf2) and Kelch-like ECH-associated protein 1 (Keap1) is a potential therapeutic strategy to control diseases involving oxidative stress. Here, six classes of known small-molecule Keap1-Nrf2 PPI inhibitors were dissected into 77 fragments in a fragment-based deconstruction reconstruction (FBDR) study and tested in four orthogonal assays. This gave 17 fragment hits of which six were shown by X-ray crystallography to bind in the Keap1 Kelch binding pocket. Two hits were merged into compound 8 with a 220-380-fold stronger affinity (K_i = 16 μM) relative to the parent fragments. Systematic optimization resulted in several novel analogues with K_i values of 0.04-0.5 μM, binding modes determined by X-ray crystallography, and enhanced microsomal stability. This demonstrates how FBDR can be used to find new fragment hits, elucidate important ligand-protein interactions, and identify new potent inhibitors of the Keap1-Nrf2 PPI.

Deep learning model predicts water interaction sites on the surface of proteins using limited-resolution data

We develop a residual deep learning model, hotWater (https://pypi.org/project/hotWater/), to identify key water interaction sites on proteins for binding models and drug discovery. This is tested on new crystal structures, as well as cryo-EM and NMR structures from the PDB and in crystallographic refinement with promising results.

Paramagnetic NMR in drug discovery

The presence of an unpaired electron in paramagnetic molecules generates significant effects in NMR spectra, which can be exploited to provide restraints complementary to those used in standard structure-calculation protocols. NMR already occupies a central position in drug discovery for its use in fragment screening, structural biology and validation of ligand-target interactions. Paramagnetic restraints provide unique opportunities, for example, for more sensitive screening to identify weaker-binding fragments. A key application of paramagnetic NMR in drug discovery, however, is to provide new structural restraints in cases where crystallography proves intractable. This is particularly important at early stages in drug-discovery programs where crystal structures of weakly-binding fragments are difficult to obtain and crystallization artefacts are probable, but structural information about ligand poses is crucial to guide medicinal chemistry. Numerous applications show the value of paramagnetic restraints to filter computational docking poses and to generate interaction models. Paramagnetic relaxation enhancements (PREs) generate a distance-dependent effect, while pseudo-contact shift (PCS) restraints provide both distance and angular information. Here, we review strategies for introducing paramagnetic centers and discuss examples that illustrate the utility of paramagnetic restraints in drug discovery. Combined with standard approaches, such as chemical shift perturbation and NOE-derived distance information, paramagnetic NMR promises a valuable source of information for many challenging drug-discovery programs.

A real-time cell-binding assay reveals dynamic features of STxB-Gb3 cointernalization and STxB-mediated cargo delivery into cancer cells

The interaction between the Shiga toxin B-subunit (STxB) and its globotriaosylceramide receptor (Gb3) has a high potential for being exploited for targeted cancer therapy. The primary goal of this study was to evaluate the capacity of STxB to carry small molecules and proteins as cargo into cells. For this purpose, an assay was designed to provide real-time information about the StxB-Gb3 interaction as well as the dynamics and mechanism of the internalization process. The assay revealed the ability to distinguish the process of binding to the cell surface from internalization and presented the importance of receptor and STxB clustering for internalization. The overall setup demonstrated that the binding mechanism is complex, and the concept of affinity is difficult to apply. Hence, time-resolved methods, providing detailed information about the interaction of STxB with cells, are critical for the optimization of intracellular delivery.

Structure-activity relationship and cardiac safety of 2-aryl-2-(pyridin-2-yl)acetamides as a new class of broad-spectrum anticonvulsants derived from Disopyramide

A series of 2-aryl-2-(pyridin-2-yl)acetamides were synthesized and screened for their anticonvulsant activity in animal models of epilepsy. The compounds were broadly active in the 'classical' maximal electroshock seizure (MES) and subcutaneous Metrazol (scMET) tests as well as in the 6 Hz and kindling models of pharmacoresistant seizures. Furthermore, the compounds showed good therapeutic indices between anticonvulsant activity and motor impairment. Structure-activity relationship (SAR) trends clearly showed the highest activity resides in unsubstituted phenyl derivatives or compounds having ortho- and meta- substituents on the phenyl ring. The 2-aryl-2-(pyridin-2-yl)acetamides were derived by redesign of the cardiotoxic sodium channel blocker Disopyramide (DISO). Our results show that the compounds preserve the capability of the parent compound to inhibit voltage gated sodium currents in patch-clamp experiments; however, in contrast to DISO, a representative compound from the series 1 displays high levels of cardiac safety in a panel of in vitro and in vivo experiments.

Introducing the CSP Analyzer: A novel Machine Learning-based application for automated analysis of two-dimensional NMR spectra in NMR fragment-based screening

NMR-based screening, especially fragment-based drug discovery is a valuable approach in early-stage drug discovery. Monitoring fragment-binding in protein-detected 2D NMR experiments requires analysis of hundreds of spectra to detect chemical shift perturbations (CSPs) in the presence of ligands screened. Computational tools are available that simplify the tracking of CSPs in 2D NMR spectra. However, to the best of our knowledge, an efficient automated tool for the assessment and binning of multiple spectra for ligand binding has not yet been described. We present a novel and fast approach for analysis of multiple 2D HSQC spectra based on machine-learning-driven statistical discrimination. The CSP Analyzer features a C# frontend interfaced to a Python ML classifier. The software allows rapid evaluation of 2D screening data from large number of spectra, reducing user-introduced bias in the evaluation. The CSP Analyzer software package is available on GitHub https://github.com/rubbs14/CSP-Analyzer/releases/tag/v1.0 under the GPL license 3.0 and is free to use for academic and commercial uses.

Structure-Activity Relationship in Pyrazolo[4,3-c]pyridines, First Inhibitors of PEX14-PEX5 Protein-Protein Interaction with Trypanocidal Activity

Trypanosoma protists are pathogens leading to a spectrum of devastating infectious diseases. The range of available chemotherapeutics against Trypanosoma is limited, and the existing therapies are partially ineffective and cause serious adverse effects. Formation of the PEX14-PEX5 complex is essential for protein import into the parasites' glycosomes. This transport is critical for parasite metabolism and failure leads to mislocalization of glycosomal enzymes, with fatal consequences for the parasite. Hence, inhibiting the PEX14-PEX5 protein-protein interaction (PPI) is an attractive way to affect multiple metabolic pathways. Herein, we have used structure-guided computational screening and optimization to develop the first line of compounds that inhibit PEX14-PEX5 PPI. The optimization was driven by several X-ray structures, NMR binding data, and molecular dynamics simulations. Importantly, the developed compounds show significant cellular activity against Trypanosoma, including the human pathogen Trypanosoma brucei gambiense and Trypanosoma cruzi parasites.

FSP1 is a glutathione-independent ferroptosis suppressor

Ferroptosis is an iron-dependent form of necrotic cell death marked by oxidative damage to phospholipids1,2. To date, ferroptosis has been thought to be controlled only by the phospholipid hydroperoxide-reducing enzyme glutathione peroxidase 4 (GPX4)3,4 and radical-trapping antioxidants5,6. However, elucidation of the factors that underlie the sensitivity of a given cell type to ferroptosis7 is crucial to understand the pathophysiological role of ferroptosis and how it may be exploited for the treatment of cancer. Although metabolic constraints8 and phospholipid composition9,10 contribute to ferroptosis sensitivity, no cell-autonomous mechanisms have been identified that account for the resistance of cells to ferroptosis. Here we used an expression cloning approach to identify genes in human cancer cells that are able to complement the loss of GPX4. We found that the flavoprotein apoptosis-inducing factor mitochondria-associated 2 (AIFM2) is a previously unrecognized anti-ferroptotic gene. AIFM2, which we renamed ferroptosis suppressor protein 1 (FSP1) and which was initially described as a pro-apoptotic gene11, confers protection against ferroptosis elicited by GPX4 deletion. We further demonstrate that the suppression of ferroptosis by FSP1 is mediated by ubiquinone (also known as coenzyme Q10, CoQ10): the reduced form, ubiquinol, traps lipid peroxyl radicals that mediate lipid peroxidation, whereas FSP1 catalyses the regeneration of CoQ10 using NAD(P)H. Pharmacological targeting of FSP1 strongly synergizes with GPX4 inhibitors to trigger ferroptosis in a number of cancer entities. In conclusion, the FSP1-CoQ10-NAD(P)H pathway exists as a stand-alone parallel system, which co-operates with GPX4 and glutathione to suppress phospholipid peroxidation and ferroptosis.

Crystal Structure of Kluyveromyces lactis Glucokinase (KlGlk1)

Glucose phosphorylating enzymes are crucial in the regulation of basic cellular processes, including metabolism and gene expression. Glucokinases and hexokinases provide a pool of phosphorylated glucose in an adenosine diphosphate (ADP)- and ATP-dependent manner to shape the cell metabolism. The glucose processing enzymes from Kluyveromyces lactis are poorly characterized despite the emerging contribution of this yeast strain to industrial and laboratory scale biotechnology. The first reports on K. lactis glucokinase (KlGlk1) positioned the enzyme as an essential component required for glucose signaling. Nevertheless, no biochemical and structural information was available until now. Here, we present the first crystal structure of KlGlk1 together with biochemical characterization, including substrate specificity and enzyme kinetics. Additionally, comparative analysis of the presented structure and the prior structures of lactis hexokinase (KlHxk1) demonstrates the potential transitions between open and closed enzyme conformations upon ligand binding.

HuR biological function involves RRM3-mediated dimerization and RNA binding by all three RRMs

HuR/ELAVL1 is an RNA-binding protein involved in differentiation and stress response that acts primarily by stabilizing messenger RNA (mRNA) targets. HuR comprises three RNA recognition motifs (RRMs) where the structure and RNA binding of RRM3 and of full-length HuR remain poorly understood. Here, we report crystal structures of RRM3 free and bound to cognate RNAs. Our structural, NMR and biochemical data show that RRM3 mediates canonical RNA interactions and reveal molecular details of a dimerization interface localized on the α-helical face of RRM3. NMR and SAXS analyses indicate that the three RRMs in full-length HuR are flexibly connected in the absence of RNA, while they adopt a more compact arrangement when bound to RNA. Based on these data and crystal structures of tandem RRM1,2-RNA and our RRM3-RNA complexes, we present a structural model of RNA recognition involving all three RRM domains of full-length HuR. Mutational analysis demonstrates that RRM3 dimerization and RNA binding is required for functional activity of full-length HuR in vitro and to regulate target mRNAs levels in human cells, thus providing a fine-tuning for HuR activity in vivo.

Come, sweet death: targeting glycosomal protein import for antitrypanosomal drug development

Glycosomes evolved as specialized system for glycolysis in trypanosomatids. These organelle rely on protein import to maintain function. A machinery of peroxin (PEX) proteins is responsible for recognition and transport of glycosomal proteins to the organelle. Disruption of PEX-based import system was expected to be a strategy against trypanosomatids. Recently, a proof of this hypothesis has been presented. Here, we review current information about trypanosomatids' glycosomal transport components as targets for new trypanocidal therapies.

NMR fragment-based screening for development of the CD44-binding small molecules

The cell-surface protein CD44, a primary receptor for hyaluronic acid (HA), is one of the most promising targets for cancer therapies. It is prominently involved in the process of tumor growth and metastasis. The possibility of modulating the CD44-HA interaction with a pharmacological inhibitor is therefore of great importance, yet until now there are only few small molecules reported to bind to CD44. Here, we describe the results of the NMR fragment-based screening conducted against CD44 by which we found eight new hit compounds that bind to the receptor with the affinity in milimolar range. The NMR-based characterization revealed that there are two possible binding modes for these compounds, and for some of them the binding is no longer possible in the presence of hyaluronic acid. This could provide an interesting starting point for the development of new high-affinity ligands targeting the CD44-HA axis.

Targeting TRAF6 E3 ligase activity with a small-molecule inhibitor combats autoimmunity

Constitutive NF-κB signaling represents a hallmark of chronic inflammation and autoimmune diseases. The E3 ligase TNF receptor-associated factor 6 (TRAF6) acts as a key regulator bridging innate immunity, pro-inflammatory cytokines, and antigen receptors to the canonical NF-κB pathway. Structural analysis and point mutations have unraveled the essential role of TRAF6 binding to the E2-conjugating enzyme ubiquitin-conjugating enzyme E2 N (Ubc13 or UBE2N) to generate Lys⁶³-linked ubiquitin chains for inflammatory and immune signal propagation. Genetic mutations disrupting TRAF6-Ubc13 binding have been shown to reduce TRAF6 activity and, consequently, NF-κB activation. However, to date, no small-molecule modulator is available to inhibit the TRAF6-Ubc13 interaction and thereby counteract NF-κB signaling and associated diseases. Here, using a high-throughput small-molecule screening approach, we discovered an inhibitor of the TRAF6-Ubc13 interaction that reduces TRAF6-Ubc13 activity both in vitro and in cells. We found that this compound, C25-140, impedes NF-κB activation in various immune and inflammatory signaling pathways also in primary human and murine cells. Importantly, C25-140 ameliorated inflammation and improved disease outcomes of autoimmune psoriasis and rheumatoid arthritis in preclinical in vivo mouse models. Hence, the first-in-class TRAF6-Ubc13 inhibitor C25-140 expands the toolbox for studying the impact of the ubiquitin system on immune signaling and underscores the importance of TRAF6 E3 ligase activity in psoriasis and rheumatoid arthritis. We propose that inhibition of TRAF6 activity by small molecules represents a promising novel strategy for targeting autoimmune and chronic inflammatory diseases.

Identification of small-molecule inhibitors of USP2a

USP2a is a deubiquitinating protease that rescues its target proteins from destruction by the proteasome by reversing the process of protein ubiquitination. USP2a shows oncogenic properties in vivo and has been found to be a specific activator of cyclin D1. Many types of cancers are addicted to cyclin D1 expression. Targeting USP2a is a promising strategy for cancer therapy but little progress has been made in the field of inhibition of USP2a. Using NMR-based fragment screening and biophysical binding assays, we have discovered small molecules that bind to USP2a. Iterations of fragment combination and structure-driven design identified two 5-(2-thienyl)-3-isoxazoles as the inhibitors of the USP2a-ubiquitin protein-protein interaction. The affinity of these molecules for the catalytic domain of USP2a parallels their ability to interfere with USP2a binding to ubiquitin in vitro. Altogether, our results establish the 5-(2-thienyl)-3-isoxazole pharmacophore as an attractive starting point for lead optimization.

Inhibitors of PEX14 disrupt protein import into glycosomes and kill Trypanosoma parasites

The parasitic protists of the Trypanosoma genus infect humans and domestic mammals, causing severe mortality and huge economic losses. The most threatening trypanosomiasis is Chagas disease, affecting up to 12 million people in the Americas. We report a way to selectively kill Trypanosoma by blocking glycosomal/peroxisomal import that depends on the PEX14-PEX5 protein-protein interaction. We developed small molecules that efficiently disrupt the PEX14-PEX5 interaction. This results in mislocalization of glycosomal enzymes, causing metabolic catastrophe, and it kills the parasite. High-resolution x-ray structures and nuclear magnetic resonance data enabled the efficient design of inhibitors with trypanocidal activities comparable to approved medications. These results identify PEX14 as an "Achilles' heel" of the Trypanosoma suitable for the development of new therapies against trypanosomiases and provide the structural basis for their development.

Inhibitors of glycosomal protein import provide new leads against trypanosomiasis

Vector-borne trypanosomatid parasite infections in tropical and sub-tropical countries constitute a major threat to humans and livestock. Trypanosoma brucei parasites are transmitted by tsetse fly and lead to African sleeping sickness in humans and Nagana in cattle. In Latin American countries, Trypanosoma cruzi infections spread by triatomine kissing bugs lead to Chagas disease. Various species of Leishmania transmitted to humans by phlebotomine sandflies manifest in a spectrum of diseases termed Leishmaniasis. 20 million people are currently infected with trypanosomatid parasites, leading to over 30,000 deaths annually and half billion people at risk of the infection. It is estimated that 300,000 Chagas infected people reside in the United States and 100,000 in Europe. Glycosomes are peroxisome-like organelles found only in trypanosomatids. Glycolysis occurs in the cytosol in all other organisms, but glycolytic enzymes and other metabolic pathways are compartmentalized inside glycosomes in trypanosomatids. Glycosomes are essential for the parasite survival and hence thought to be an attractive drug target. Our recent study [Dawidowski et al. Science (2017)] is the first to report small molecule inhibitors of glycosomal protein import. Using structure-based drug design, we developed small molecule inhibitors of the Trypanosoma PEX5-PEX14 protein-protein interaction that disrupt glycosomal protein import and kill the parasites. Oral treatment of T. brucei infected mice with PEX14 inhibitor significantly reduced the parasite levels with no adverse effect on mice. The study provides the grounds for further development of the glycosome inhibitors into clinical candidates and validates the parasite protein-protein interactions as drug targets.

The conserved ubiquitin-like protein Hub1 plays a critical role in splicing in human cells

Different from canonical ubiquitin-like proteins, Hub1 does not form covalent conjugates with substrates but binds proteins non-covalently. In Saccharomyces cerevisiae, Hub1 associates with spliceosomes and mediates alternative splicing of SRC1, without affecting pre-mRNA splicing generally. Human Hub1 is highly similar to its yeast homolog, but its cellular function remains largely unexplored. Here, we show that human Hub1 binds to the spliceosomal protein Snu66 as in yeast; however, unlike its S. cerevisiae homolog, human Hub1 is essential for viability. Prolonged in vivo depletion of human Hub1 leads to various cellular defects, including splicing speckle abnormalities, partial nuclear retention of mRNAs, mitotic catastrophe, and consequently cell death by apoptosis. Early consequences of Hub1 depletion are severe splicing defects, however, only for specific splice sites leading to exon skipping and intron retention. Thus, the ubiquitin-like protein Hub1 is not a canonical spliceosomal factor needed generally for splicing, but rather a modulator of spliceosome performance and facilitator of alternative splicing.

Discovery of highly potent p53-MDM2 antagonists and structural basis for anti-acute myeloid leukemia activities

The inhibition of p53-MDM2 interaction is a promising new approach to non-genotoxic cancer treatment. A potential application for drugs blocking the p53-MDM2 interaction is acute myeloid leukemia (AML) due to the occurrence of wild type p53 (wt p53) in the majority of patients. Although there are very promising preclinical results of several p53-MDM2 antagonists in early development, none of the compounds have yet proven the utility as a next generation anticancer agent. Herein we report the design, synthesis and optimization of YH239-EE (ethyl ester of the free carboxylic acid compound YH239), a potent p53-MDM2 antagonizing and apoptosis-inducing agent characterized by a number of leukemia cell lines as well as patient-derived AML blast samples. The structural basis of the interaction between MDM2 (the p53 receptor) and YH239 is elucidated by a co-crystal structure. YH239-EE acts as a prodrug and is the most potent compound that induces apoptosis in AML cells and patient samples. The observed superior activity compared to reference compounds provides the preclinical basis for further investigation and progression of YH239-EE.

Transient protein states in designing inhibitors of the MDM2-p53 interaction

Reactivation of p53 by release of the functional protein from its inhibition by MDM2 provides an efficient, nongenotoxic approach to a wide variety of cancers. We present the cocrystal structures of two complexes of MDM2 with inhibitors based on 6-chloroindole scaffolds. Both molecules bound to a distinct conformational state of MDM2 with nM-μM affinities. In contrast to other structurally characterized antagonists, which mimic three amino acids of p53 (Phe19, Trp23, and Leu26), the compounds induced an additional hydrophobic pocket on the MDM2 surface and unveiled a four-point binding mode. The enlarged interaction interface of the inhibitors resulted in extension of small molecules binding toward the "lid" segment of MDM2 (residues 19-23)--a nascent element that interferes with p53 binding. As supported by protein engineering and molecular dynamics studies, employing these unstable elements of MDM2 provides an efficient and yet unexplored alternative in development of MDM2-p53 association inhibitors.

Structural and functional analysis of the natural JNK1 inhibitor quercetagetin

c-Jun NH2-terminal kinases (JNKs) and phosphatidylinositol 3-kinase (PI3-K) play critical roles in chronic diseases such as cancer, type II diabetes, and obesity. We describe here the binding of quercetagetin (3,3',4',5,6,7-hydroxyflavone), related flavonoids, and SP600125 to JNK1 and PI3-K by ATP-competitive and immobilized metal ion affinity-based fluorescence polarization assays and measure the effect of quercetagetin on JNK1 and PI3-K activities. Quercetagetin attenuated the phosphorylation of c-Jun and AKT, suppressed AP-1 and NF-κB promoter activities, and also reduced cell transformation. It attenuated tumor incidence and reduced tumor volumes in a two-stage skin carcinogenesis mouse model. Our crystallographic structure determination data show that quercetagetin binds to the ATP-binding site of JNK1. Notably, the interaction between Lys55, Asp169, and Glu73 of JNK1 and the catechol moiety of quercetagetin reorients the N-terminal lobe of JNK1, thereby improving compatibility of the ligand with its binding site. The results of a theoretical docking study suggest a binding mode of PI3-K with the hydroxyl groups of the catechol moiety forming hydrogen bonds with the side chains of Asp964 and Asp841 in the p110γ catalytic subunit. These interactions could contribute to the high inhibitory activity of quercetagetin against PI3-K. Our study suggests the potential use of quercetagetin in the prevention or therapy of cancer and other chronic diseases.

Structure of the stapled p53 peptide bound to Mdm2

Mdm2 is a major negative regulator of the tumor suppressor p53 protein, a protein that plays a crucial role in maintaining genome integrity. Inactivation of p53 is the most prevalent defect in human cancers. Inhibitors of the Mdm2-p53 interaction that restore the functional p53 constitute potential nongenotoxic anticancer agents with a novel mode of action. We present here a 2.0 Å resolution structure of the Mdm2 protein with a bound stapled p53 peptide. Such peptides, which are conformationally and proteolytically stabilized with all-hydrocarbon staples, are an emerging class of biologics that are capable of disrupting protein-protein interactions and thus have broad therapeutic potential. The structure represents the first crystal structure of an i, i + 7 stapled peptide bound to its target and reveals that rather than acting solely as a passive conformational brace, a staple can intimately interact with the surface of a protein and augment the binding interface.

Role of the ubiquitin-like protein Hub1 in splice-site usage and alternative splicing

Alternative splicing of pre-messenger RNAs diversifies gene products in eukaryotes and is guided by factors that enable spliceosomes to recognize particular splice sites. Here we report that alternative splicing of Saccharomyces cerevisiae SRC1 pre-mRNA is promoted by the conserved ubiquitin-like protein Hub1. Structural and biochemical data show that Hub1 binds non-covalently to a conserved element termed HIND, which is present in the spliceosomal protein Snu66 in yeast and mammals, and Prp38 in plants. Hub1 binding mildly alters spliceosomal protein interactions and barely affects general splicing in S. cerevisiae. However, spliceosomes that lack Hub1, or are defective in Hub1-HIND interaction, cannot use certain non-canonical 5' splice sites and are defective in alternative SRC1 splicing. Hub1 confers alternative splicing not only when bound to HIND, but also when experimentally fused to Snu66, Prp38, or even the core splicing factor Prp8. Our study indicates a novel mechanism for splice site utilization that is guided by non-covalent modification of the spliceosome by an unconventional ubiquitin-like modifier.

The structure-based design of Mdm2/Mdmx-p53 inhibitors gets serious

The p53 protein is the cell's principal bastion of defense against tumor-associated DNA damage. Commonly referred as a "guardian of the genome", p53 is responsible for determining the fate of the cell when the integrity of its genome is damaged. The development of tumors requires breaching this defense line. All known tumor cells either mutate the p53 gene, or in a similar number of cases, use internal cell p53 modulators, Mdm2 and Mdmx proteins, to disable its function. The release of functional p53 from the inhibition by Mdm2 and Mdmx should in principle provide an efficient, nongenotoxic means of cancer therapy. In recent years substantial progress has been made in developing novel p53-activating molecules thanks to several reported crystal structures of Mdm2/x in complex with p53-mimicking peptides and nonpeptidic drug candidates. Understanding the structural attributes of ligand binding holds the key to developing novel, highly effective, and selective drug candidates. Two low-molecular-weight compounds have just recently progressed into early clinical studies.