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.

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.

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.

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.

Steric hindrance and structural flexibility shape the functional properties of a guanine-rich oligonucleotide.

Ligand/protein molecular recognition involves a dynamic process, whereby both partners require a degree of structural plasticity to regulate the binding/unbinding event. Here, we present the characterization of the interaction between a highly dynamic G-rich oligonucleotide, M08s-1, and its target protein, human α-thrombin. M08s-1 is the most active anticoagulant aptamer selected thus far. Circular dichroism and gel electrophoresis analyses indicate that both intramolecular and intermolecular G-quadruplex structures are populated in solution. The presence of thrombin stabilises the antiparallel intramolecular chair-like G-quadruplex conformation, that provides by far the main contribution to the biological activity of the aptamer. The crystal structure of the thrombin-oligonucleotide complex reveals that M08s-1 adopts a kinked structural organization formed by a G-quadruplex domain and a long duplex module, linked by a stretch of five purine bases. The quadruplex motif hooks the exosite I region of thrombin and the duplex region is folded towards the surface of the protein. This structural feature, which has never been observed in other anti-exosite I aptamers with a shorter duplex motif, hinders the approach of a protein substrate to the active site region and may well explain the significant increase in the anticoagulant activity of M08s-1 compared to the other anti-exosite I aptamers.

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.

A Buckycatcher in Solution-A Computational Perspective.

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.

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.

Structure-based design, synthesis and evaluation of a novel family of PEX5-PEX14 interaction inhibitors against Trypanosoma.

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.

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

Small molecule mediated inhibition of protein cargo recognition by peroxisomal transport receptor PEX5 is toxic to Trypanosoma.

Trypanosomiases are life-threatening infections of humans and livestock, and novel effective therapeutic approaches are needed. Trypanosoma compartmentalize glycolysis into specialized organelles termed glycosomes. Most of the trypanosomal glycolytic enzymes harbor a peroxisomal targeting signal-1 (PTS1) which is recognized by the soluble receptor PEX5 to facilitate docking and translocation of the cargo into the glycosomal lumen. Given its pivotal role in the glycosomal protein import, the PEX5-PTS1 interaction represents a potential target to inhibit import of glycolytic enzymes and thus kill the parasite. We developed a fluorescence polarization (FP)-based assay for monitoring the PEX5-PTS1 interaction and performed a High Throughput Screening (HTS) campaign to identify small molecule inhibitors of the interaction. Six of the identified hits passed orthogonal selection criteria and were found to inhibit parasite growth in cell culture. Our results validate PEX5 as a target for small molecule inhibitors and provide scaffolds suitable for further pre-clinical development of novel trypanocidal compounds.

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.

NudC guides client transfer between the Hsp40/70 and Hsp90 chaperone systems.

In the eukaryotic cytosol, the Hsp70 and the Hsp90 chaperone machines work in tandem with the maturation of a diverse array of client proteins. The transfer of nonnative clients between these systems is essential to the chaperoning process, but how it is regulated is still not clear. We discovered that NudC is an essential transfer factor with an unprecedented mode of action: NudC interacts with Hsp40 in Hsp40-Hsp70-client complexes and displaces Hsp70. Then, the interaction of NudC with Hsp90 allows the direct transfer of Hsp40-bound clients to Hsp90 for further processing. Consistent with this mechanism, NudC increases client activation in vitro as well as in cells and is essential for cellular viability. Together, our results show the complexity of the cooperation between the major chaperone machineries in the eukaryotic cytosol.

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.