Search

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.

Large language models (LLMs) have greatly enhanced our ability to understand biology and chemistry. Yet, relatively few robust methods have been reported for structure-based drug discovery. Highly precise biomolecule-ligand interaction datasets are urgently needed in particular for LLMs, that require extensive training data. We present MISATO, the first dataset that combines quantum mechanics properties of small molecules and associated molecular dynamics simulations of about 20000 experimental protein-ligand complexes. Starting from the PDBbind dataset, semi-empirical quantum mechanics was used to systematically refine these structures. The largest collection to date of molecular dynamics traces of protein-ligand complexes in explicit water are included, accumulating to 170 μs. We give ML baseline models and simple Python data loaders, and aim to foster a thriving community around MISATO (https://github.com/t7morgen/misato-dataset). An easy entry point for ML experts is provided without the need of deep domain expertise to enable the next generation of drug discovery AI models.

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.

Targeting and import of peroxisomal proteins depends on PEX5, PEX14 and PEX13. We present a biochemical and structural characterization of the PEX13 C-terminal region. By combining NMR spectroscopy, X-ray crystallography and biochemical methods, we show that the PEX13 SH3 domain mediates intramolecular interactions with a newly identified proximal FxxxF motif and also binds to WxxxF peptide motifs from the PEX5 NTD, demonstrating evolutionary conservation of this interaction from yeast to human. Strikingly, the C-terminal FxxxF motif autoinhibits the WxxxF/Y binding surface on the PEX13 SH3 domain. This is supported by high-resolution crystal structures, which show FxxxF or WxxxF/Y binding to the same, non-canonical surface on the SH3 domain. The FxxxF motif also binds the PEX14 NTD with micromolar affinity. Surprisingly, the canonical binding surface for PxxP motifs on the human PEX13 SH3 fold does not recognize PxxP motifs in PEX14, distinct from the yeast ortholog. The dynamic network of PEX13, PEX14 and PEX5 interactions mediated by diaromatic peptide motifs fine-tunes and modulates peroxisomal matrix import in cells.

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.

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.

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.

Trypanosomiases are neglected tropical diseases caused by Trypanosoma (sub)species. Available treatments are limited and have considerable adverse effects and questionable efficacy in the chronic stage of the disease, urgently calling for the identification of new targets and drug candidates. Recently, we have shown that impairment of glycosomal protein import by the inhibition of the PEX5-PEX14 protein-protein interaction (PPI) is lethal to Trypanosoma. Here, we report the development of a novel dibenzo[b,f][1,4]oxazepin-11(10H)-one scaffold for small molecule inhibitors of PEX5-PEX14 PPI. The initial hit was identified by a high throughput screening (HTS) of a library of compounds. A bioisosteric replacement approach allowed to replace the metabolically unstable sulphur atom from the initial dibenzo[b,f][1,4]thiazepin-11(10H)-one HTS hit with oxygen. A crystal structure of the hit compound bound to PEX14 surface facilitated the rational design of the compound series accessible by a straightforward chemistry for the initial structure-activity relationship (SAR) analysis. This guided the design of compounds with trypanocidal activity in cell-based assays providing a promising starting point for the development of new drug candidates to tackle trypanosomiases.

In this work, we study the buckycatcher (C60H28) 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.

Coronaviruses protect their single-stranded RNA genome with a methylated cap during replication. 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 where nsp10 acts as a co-factor to both. Aditionally, 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 should be impossible due to steric clashes. Here, we show that nsp14, nsp10, and nsp16 can form a heterotrimer complex. This interaction is expected to encourage formation of mature capped viral mRNA, modulating the 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.

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. This research has been supported by Grants Maestro 2017/26/A/ST5/00572 (to T.A.H.) , Sonata UMO-2020/39/D/ST4/01344 (to E.S.) , Preludium UMO-2021/41/N/ST4/03485 (to M.Z.) , and Preludium UMO-2020/37/N/ST4/02691 (to D.M.) from the National Science Centre, Poland. X.d.C. thanks the Basque Country Government for the predoctoral and EGONLABUR grants.

We show that a water envelope network plays a critical role in protein-protein interactions (PPI). The potency of a PPI inhibitor is modulated by orders of magnitude on manipulation of the solvent envelope alone. The structure-activity relationship of PEX14 inhibitors was analyzed as an example using in silico and X-ray data.

BET proteins such as BRD3 are oncogenic transcriptional coactivators. SPOP binding triggers their proteasomal degradation. In both endometrial and prostate cancers, SPOP mutations occur in the MATH domain, but with opposed influence on drug susceptibility. In prostate cancer, SPOP mutations presumably cause increased BET levels, decreasing BET inhibitor drug susceptibility. As opposed, in endometrial cancer, decreased BET levels concomitant with higher BET inhibitor drug susceptibility were observed. Here, we present the to our knowledge first co-crystal structure of SPOP and a bromodomain containing protein (BRD3). Our structural and biophysical data confirm the suggested loss-of-function in prostate cancer-associated SPOP mutants and provide mechanistic explanation. As opposed to previous literature, our data on endometrial cancer-associated SPOP mutants do not show altered binding behavior compared to wild-type SPOP, indicating a more complex regulatory mechanism. SPOP mutation screening may thus be considered a valuable personalized medicine tool for effective antitumor therapy. (C) 2019 Elsevier Ltd. All rights reserved.

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. (c) 2021 The Author(s). Published by Elsevier Ltd. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).

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 (Ki = 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.

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 (Ki = 16 μM) relative to the parent fragments. Systematic optimization resulted in several novel analogues with Ki 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. 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 (Ki = 16 μM) relative to the parent fragments. Systematic optimization resulted in several novel analogues with Ki 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.

SummaryThe COVID-19 pandemic caused by SARS-CoV-2 has been socially and economically devastating. Despite an unprecedented research effort, effective therapeutics are still missing to limit severe disease and mortality. Using high-throughput screening, we identified acriflavine as a potent papain-like protease (PLpro) inhibitor. NMR titrations and a co-crystal structure confirm that acriflavine blocks the PLprocatalytic pocket in an unexpected binding mode. We show that the drug inhibits viral replication at nanomolar concentration in cellular models,in vivoin mice andex vivoin human airway epithelia, with broad range activity against SARS-CoV-2 and other betacoronaviruses. Considering that acriflavine is an inexpensive drug approved in some countries, it may be immediately tested in clinical trials and play an important role during the current pandemic and future outbreaks.

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 (Ki = 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.

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.

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.

Multidrug resistance among Gram-negative bacteria is a major global public health threat. Metallo-β-lactamases (MBLs) target the most widely used antibiotic class, the β-lactams, including the most recent generation of carbapenems. Interspecies spread renders these enzymes a serious clinical threat, and there are no clinically available inhibitors. We present the crystal structures of IMP-13, a structurally uncharacterized MBL from the Gram-negative bacterium Pseudomonas aeruginosa found in clinical outbreaks globally, and characterize the binding using solution nuclear magnetic resonance spectroscopy and molecular dynamics simulations., Multidrug resistance among Gram-negative bacteria is a major global public health threat. Metallo-β-lactamases (MBLs) target the most widely used antibiotic class, the β-lactams, including the most recent generation of carbapenems. Interspecies spread renders these enzymes a serious clinical threat, and there are no clinically available inhibitors. We present the crystal structures of IMP-13, a structurally uncharacterized MBL from the Gram-negative bacterium Pseudomonas aeruginosa found in clinical outbreaks globally, and characterize the binding using solution nuclear magnetic resonance spectroscopy and molecular dynamics simulations. The crystal structures of apo IMP-13 and IMP-13 bound to four clinically relevant carbapenem antibiotics (doripenem, ertapenem, imipenem, and meropenem) are presented. Active-site plasticity and the active-site loop, where a tryptophan residue stabilizes the antibiotic core scaffold, are essential to the substrate-binding mechanism. The conserved carbapenem scaffold plays the most significant role in IMP-13 binding, explaining the broad substrate specificity. The observed plasticity and substrate-locking mechanism provide opportunities for rational drug design of novel metallo-β-lactamase inhibitors, essential in the fight against antibiotic resistance.

The surfaces of proteins are generally hydrophilic but there have been reports of sites that exhibit an exceptionally high affinity for individual water molecules. Not only do such molecules often fulfil critical biological functions, but also, they may alter the binding of newly designed drugs. In crystal structures, sites consistently occupied in each unit cell yield electron density clouds that represent water molecule presence. These are recorded in virtually all high-resolution structures obtained through X-ray diffraction. In this work, we utilized the wealth of data from the RCSB Protein Data Bank to train a residual deep learning model named ‘hotWater’ to identify sites on the surface of proteins that are most likely to bind water, the so-called water hot spots. The model can be used to score existing water molecules from a PDB file to provide their ranking according to the predicted binding strength or to scan the surface of a protein to determine the most likely water hot-spots de novo. This is computationally much more efficient than currently used molecular dynamics simulations. Based on testing the model on three example proteins, which have been resolved using both high-resolution X-ray crystallography (providing accurate positions of trapped waters) as well as low-resolution X-ray diffraction, NMR or CryoEM (where structure refinement does not yield water positions), we were able to show that the hotWater method is able to recover in the “water-free” structures many water binding sites known from the high-resolution structures. A blind test on a newly solved protein structure with waters removed from the PDB also showed good prediction of the crystal water positions. This was compared to two known algorithms that use electron density and was shown to have higher recall at resolutions >2.6 Å. We also show that the algorithm can be applied to novel proteins such as the RNA polymerase complex from SARS-CoV-2, which could be of use in drug discovery. The hotWater model is freely available at (https://pypi.org/project/hotWater/).

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 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.

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., Computational and Structural Biotechnology Journal, 18, ISSN:2001-0370

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.

Lipases show high stability in lipophilic solvents and catalyze reactions at the water-oil interfaces, which are of great industrial interest. One promising application of lipases is the production of epoxides from alkenes and hydrogen peroxide. So far, little attention has been given to uncover the reaction mechanism for this in detail at the molecular level. Here, we present structural and mutational analysis of a lipase from Penicillium camembertii that indicates a two-stage synergistic mechanism for this reaction. Surprisingly, a mutant devoid of the catalytic serine retains a fraction of activity while histidine from the catalytic triad is absolutely critical to maintain the enzyme activity. Histidine appears to perform a dual-activation role acting both towards hydrogen peroxide and the catalytic serine. These results thus allow a better understanding of enzymatic epoxidation and engineering of more potent, stable and selective enzymes.

Fungal lipases are efficient and environment-friendly biocatalysts for many industrially relevant processes. One of the most widely applied lipases in the manufacturing industry is Candida antarctica lipase B (CALB). Here, we report the biochemical and structural characterization of a novel CALB-like lipase from an important human pathogen-Aspergillus fumigatus (AFLB), which has high sn-1,3-specificity toward triolein. AFLB crystal structure displays a CALB-like catalytic domain and hosts a unique tightly closed 'lid' domain that contains a disulfide bridge, as well as an extra N-terminal subdomain composed of residues 1-128 (including the helix α1-α5 located above the active site). To gain insight into the function of this novel lid and N-terminal subdomain, we constructed and characterized a series of mutants in these two domains. Deleting the protruding bulk lid's residues, replacing the bulk and tight lid with a small and loose lid from CALB, or breaking the disulfide bridge increased the affinity of CALB for glyceride substrates and improved its catalytic activity, along with the loss of enzyme fold stability and thermostability. N-terminal truncation mutants revealed that the N-terminal peptide (residues 1-59) is a strong inhibitor of AFLB binding to lipid films. This peptide thus limits AFLB's penetration power and specific activity, revealing a unique enzyme activity regulatory mechanism. Our findings on the functional and structural properties of AFLB provide a better understanding of the functions of the CALB-like lipases and pave the way for its future protein engineering. DATABASE: Structural data are available in the Protein Data Bank under the accession numbers 6IDY.

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.

Summary Pdx1 is a transcription factor crucial for development and maintenance of a functional pancreas. It regulates insulin expression and glucose homeostasis. SPOP is an E3-ubiquitin ligase adaptor protein that binds Pdx1, thus triggering its ubiquitination and proteasomal degradation. However, the underlying mechanisms are not well understood. Here, we present the crystal structure of the SPOP-Pdx1 complex. We show that Pdx1 residues 223–233 bind to SPOP MATH domain with low micromolar affinity. The interface is extended compared to other SPOP-client proteins. Previously, Pdx1 phosphorylation has been proposed to have a regulatory function. In this respect we show that phosphorylation lowers the affinity of Pdx1 to SPOP by isothermal titration calorimetry and nuclear magnetic resonance data. Our data provide insights into a critical protein-protein interaction that regulates cellular Pdx1 levels by SPOP-mediated decay. A reduction of Pdx1 levels in β cells is linked to apoptosis and considered a hallmark of type 2 diabetes.

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 fulllength HuR in vitro and to regulate target mRNAs levels in human cells, thus providing a fine-tuning for HuR activity in vivo.

Staphylococcal SplB protease belongs to the chymotrypsin family. Chymotrypsin zymogen is activated by proteolytic processing at the N terminus, resulting in significant structural rearrangement at the active site. Here, we demonstrate that the molecular mechanism of SplB protease activation differs significantly and we characterize the novel mechanism in detail. Using peptide and protein substrates we show that the native signal peptide, or any N-terminal extension, has an inhibitory effect on SplB. Only precise N-terminal processing releases the full proteolytic activity of the wild type analogously to chymotrypsin. However, comparison of the crystal structures of mature SplB and a zymogen mimic show no rearrangement at the active site whatsoever. Instead, only the formation of a unique hydrogen bond network, distant form the active site, by the new N-terminal glutamic acid of mature SplB is observed. The importance of this network and influence of particular hydrogen bond interactions at the N terminus on the catalytic process is demonstrated by evaluating the kinetics of a series of mutants. The results allow us to propose a consistent model where changes in the overall protein dynamics rather than structural rearrangement of the active site are involved in the activation process.

Staphylococcus aureus is responsible for a variety of human infections, including life-threatening, systemic conditions. Secreted proteome, including a range of proteases, constitutes the major virulence factor of the bacterium. However, the functions of individual enzymes, in particular SplA protease, remain poorly characterized. Here, we report development of specific inhibitors of SplA protease. The design, synthesis, and activity of a series of α-aminoalkylphosphonate diaryl esters and their peptidyl derivatives are described. Potent inhibitors of SplA are reported, which may facilitate future investigation of physiological function of the protease. The binding modes of the high-affinity compounds Cbz-PheP-(OC6H4−4-SO2CH3)2 and Suc-Val-Pro-PheP-(OC6H5)2 are revealed by high-resolution crystal structures of complexes with the protease. Surprisingly, the binding mode of both compounds deviates from previously characterized canonical interaction of α-aminoalkylphosphonate peptidyl derivatives and family S1 serine proteases.

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.

Das zellulare p53-Protein ist das Bollwerk im Kampf gegen Tumor-assoziierte DNA-Schaden. Allgemein auch als “Beschutzer des Genoms” bekannt, ist p53 fur das Schicksal der Zelle verantwortlich, wenn eine substanzielle DNA-Schadigung vorliegt. Tumoren konnen sich nur entwickeln, wenn diese Verteidigungslinie durchbrochen wird. Alle bekannten Tumorzellen mutieren das p53-Gen oder nutzen interne p53-Modulatoren, wie die beiden Proteine Mdm2 und Mdmx, um p53 zu desaktivieren. Das Freisetzen von p53 aus den Komplexen mit Mdm2 oder Mdmx sollte daher eine effiziente und nicht gentoxische Krebstherapie ermoglichen. Dank mehrerer Kristallstrukturen von Mdm2/x im Komplex mit p53-imitierenden Molekulen, sowohl Peptiden als auch niedermolekularen Verbindungen, gab es in den letzten Jahren fundamentale Fortschritte bei der Entwicklung neuer p53-aktivierender Molekule. Das Verstandnis der Bindungseigenschaften des Liganden ermoglicht die Entwicklung neuer, hoch aktiver und selektiver Wirkstoff-Kandidaten. Zwei niedermolekulare Verbindungen werden seit Kurzem in klinischen Studien untersucht.

Three classes of proteins are known to nucleate new filaments: the Arp2/3 complex, formins, and the third group of proteins that contain ca. 25 amino acid long actin-binding Wiskott-Aldrich syndrome protein homology 2 domains, called the WH2 repeats. Crystal structures of the complexes between the actin-binding WH2 repeats of the Spire protein and actin were determined for the Spire single WH2 domain D, the double (SpirCD), triple (SpirBCD), quadruple (SpirABCD) domains, and an artificial Spire WH2 construct comprising three identical D repeats (SpirDDD). SpirCD represents the minimal functional core of Spire that can nucleate actin filaments. Packing in the crystals of the actin complexes with SpirCD, SpirBCD, SpirABCD, and SpirDDD shows the presence of two types of assemblies, “side-to-side” and “straight-longitudinal,” which can serve as actin filament nuclei. The principal feature of these structures is their loose, open conformations, in which the sides of actins that normally constitute the inner interface core of a filament are flipped inside out. These Spire structures are distant from those seen in the filamentous nuclei of Arp2/3, formins, and in the F-actin filament.

Intensive anticancer drug discovery efforts have been made to develop small molecule inhibitors of the p53-MDM2 and p53-MDMX interactions. We present here the structures of the most potent inhibitors bound to MDM2 and MDMX that are based on the new imidazo-indole scaffold. In addition, the structure of the recently reported spiro-oxindole inhibitor bound to MDM2 is described. The structures indicate how the substituents of a small molecule that bind to the three subpockets of the MDM2/X-p53 interaction should be optimized for effective binding to MDM2 and/or MDMX. While the spiro-oxindole inhibitor triggers significant ligand-induced changes in MDM2, the imidazo-indoles share similar binding modes for MDMX and MDM2, but cause only minimal induced-fit changes in the structures of both proteins. Our study includes the first structure of the complex between MDMX and a small molecule and should aid in developing efficient scaffolds for binding to MDMX and/or MDM2.

Staphylococcus aureus is a dangerous human pathogen whose antibiotic resistance is steadily increasing and no efficient vaccine is as yet available. This serious threat drives extensive studies on staphylococcal physiology and pathogenicity pathways, especially virulence factors. Spl (serine protease-like) proteins encoded by an operon containing up to six genes are a good example of poorly characterized secreted proteins probably involved in virulence. In the present study, we describe an efficient heterologous expression system for SplA and detailed biochemical and structural characterization of the recombinant SplA protease. The enzyme shares a significant sequence homology to V8 protease and epidermolytic toxins which are well documented staphylococcal virulence factors. SplA has a very narrow substrate specificity apparently imposed by the precise recognition of three amino acid residues positioned N-terminal to the hydrolysed peptide bond. To explain determinants of this extended specificity we resolve the crystal structure of SplA and define the consensus model of substrate binding. Furthermore we demonstrate that artificial N-terminal elongation of mature SplA mimicking a naturally present signal peptide abolishes enzymatic activity. The probable physiological role of the process is discussed. Of interest, even though precise N-terminal trimming is a common regulatory mechanism among S1 family enzymes, the crystal structure of SplA reveals novel significantly different mechanistic details.

Upon exposure to DNA damage the p53 tumor suppressor is accumulated and activated to stall cellular growth. For this to occur, p53 must be relieved from its major inhibitors, Mdm2 (Hdm2 in humans) and Mdmx (Mdm4; Hdmx in humans). A key mechanism controlling this relief is the post-translational modifications of p53 and its inhibitors. We have previously demonstrated that the stress-activated tyrosine kinase, c-Abl, contributes to the relief of p53 from Hdm2. Because Hdmx is the major inhibitor of p53 activity, the additional possibility that c-Abl protects p53 through targeting Hdmx was explored in this study. c-Abl was found to interact with and to phosphorylate Hdmx. This phosphorylation was enhanced in response to DNA damage. Importantly, we mapped the sites of phosphorylation to the p53 binding domain of Hdmx. One of these phosphorylations, on tyrosine 99, inhibited Hdmx interaction with p53. This inhibition is consistent with the predicted role of this residue in the interaction with p53 based on the crystal structure of the interaction site. Our results show that c-Abl not only targets Hdm2, but also Hdmx, which together contribute to p53 activation in response to DNA damage.

NMR-based drug screening methods provide the most reliable characterization of binding propensities of ligands to their target proteins. They are, however, one of the least effective methods in terms of the amount of protein required and the time needed for acquiring an NMR experiment. We show here that the introduction of tryptophan to proteins permits rapid screening by monitoring a simple 1D proton NMR signal of the NH side chain ((N)H(epsilon)) of the tryptophan. The method could also provide quantitative characterization of the antagonist-protein and antagonist-protein-protein interactions in the form of KDs and fractions of the released proteins from their mutual binding. We illustrate the method with the lead compounds that block the Mdm2-p53 interaction and by studying inhibitors that bind to cyclin-dependent kinase 2 (CDK2).

Aims The p53 protein, also called guardian of the genome, has a key role in cell cycle regulation. It is activated under stressful circumstances, such as DNA damage which results in permanent arrest or cell death. The protein is disabled in several types of human cancer due to over-expression of the two regulators, Mdm2 and Mdmx. As a result, inhibiting Mdm subtypes could reactivate p53 and bring about a promising therapeutic strategy in cancers. Main methods Here a structure-based pharmacophore search and docking simulation are presented in order to filter our in-house library which contains 1035 compounds to find novel scaffolds that inhibit Mdm2 and Mdmx concomitantly. Afterwards, fluorescence polarization binding assay was used to obtain inhibition constant of final compounds. Key findings Thirty two ligands were introduced to bioassay as a result of in-silico methods. Twelve of them inhibit both proteins with almost balanced K i value ranging from 18 to 162 μM for Mdm2 and 18 to 233 μM for Mdmx. It was observed that all compounds fill Phe19 and Trp23 pockets of Mdm2/x binding sites and form a hydrogen bond with Trp23 pocket's neighbor amino acids in a manner similar to p53 protein. Additionally, it was concluded that Trp23 pocket of Mdmx has a bigger hydrophobic volume comparing with the one of Mdm2. Significance Three structure–activity relationship patterns are supposed which one of them presents usefulness features and can be used in future studies. This study presents first qualitative SAR for dual inhibitors against Mdm2/x.

This paper describes the identification of compounds, which stimulate the activity of the protein phosphatase PPH-5 and addresses the influence of the identified compounds on the enzymatic properties and the potential mechanism of their action., Protein phosphatase 5 (PP5) is an evolutionary conserved serine/threonine phosphatase. Its dephosphorylation activity modulates a diverse set of cellular factors including protein kinases and the microtubule-associated tau protein involved in neurodegenerative disorders. It is auto-regulated by its heat-shock protein (Hsp90)-interacting tetratricopeptide repeat (TPR) domain and its C-terminal α-helix. In the present study, we report the identification of five specific PP5 activators [PP5 small-molecule activators (P5SAs)] that enhance the phosphatase activity up to 8-fold. The compounds are allosteric modulators accelerating efficiently the turnover rate of PP5, but do barely affect substrate binding or the interaction between PP5 and the chaperone Hsp90. Enzymatic studies imply that the compounds bind to the phosphatase domain of PP5. For the most promising compound crystallographic comparisons of the apo PP5 and the PP5–P5SA-2 complex indicate a relaxation of the auto-inhibited state of PP5. Residual electron density and mutation analyses in PP5 suggest activator binding to a pocket in the phosphatase/TPR domain interface, which may exert regulatory functions. These compounds thus may expose regulatory mechanisms in the PP5 enzyme and serve to develop optimized activators based on these scaffolds.

Despite the fact that most lipases are believed to be active against triacylglycerides, there is a small group of lipases that are active only on mono- and diacylglycerides. The reason for this difference in substrate scope is not clear. We tried to identify the reasons for this in the lipase from Malassezia globosa. By protein engineering, and with only one mutation, we managed to convert this enzyme into a typical triacylglycerol lipase (the wild-type lipase does not accept triacylglycerides). The variant Q282L accepts a broad spectrum of triacylglycerides, although the catalytic behavior is altered to some extent. From in silico analysis it seems that specific hydrophobic interactions are key to the altered substrate specificity.

Summary The III-A intermediate constitutes the major rate-determining step in the oxidative folding of leech carboxypeptidase inhibitor (LCI). In this work, III-A has been directly purified from the folding reaction and structurally characterized by NMR spectroscopy. This species, containing three native disulfides, displays a highly native-like structure; however, it lacks some secondary structure elements, making it more flexible than native LCI. III-A represents a structurally determined example of a disulfide-insecure intermediate; direct oxidation of this species to the fully native protein seems to be restricted by the burial of its two free cysteine residues inside a native-like structure. We also show that theoretical approaches based on topological constraints predict with good accuracy the presence of this folding intermediate. Overall, the derived results suggest that, as it occurs with non-disulfide bonded proteins, native-like interactions between segments of secondary structure rather than the crosslinking of disulfide bonds direct the folding of LCI.

Seven different, but highly conserved 14-3-3 proteins are involved in diverse signaling pathways in human cells. It is unclear how the 14-3-3sigma isoform, a transcriptional target of p53, exerts its inhibitory effect on the cell cycle in the presence of other 14-3-3 isoforms, which are constitutively expressed at high levels. In order to identify structural differences between the 14-3-3 isoforms, we solved the crystal structure of the human 14-3-3sigma protein at a resolution of 2.8 Angstroms and compared it to the known structures of 14-3-3zeta and 14-3-3tau. The global architecture of the 14-3-3sigma fold is similar to the previously determined structures of 14-3-3zeta and 14-3-3t: two 14-3-3sigma molecules form a cup-shaped dimer. Significant differences between these 14-3-3 isoforms were detected adjacent to the amphipathic groove, which mediates the binding to phosphorylated consensus motifs in 14-3-3-ligands. Another specificity determining region is localized between amino-acids 203 to 215. These differences presumably select for the interaction with specific ligands, which may explain the different biological functions of the respective 14-3-3 isoforms. Furthermore, the two 14-3-3sigma molecules forming a dimer differ by the spatial position of the ninth helix, which is shifted to the inside of the ligand interaction surface, thus indicating adaptability of this part of the molecule. In addition, 5 non-conserved residues are located at the interface between two 14-3-3sigma proteins forming a dimer and represent candidate determinants of homo- and hetero-dimerization specificity. The structural differences among the 14-3-3 isoforms described here presumably contribute to isoform-specific interactions and functions.