Successful generation of structural information for fragment-based drug discovery

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.

Fluorine and chlorine substituted adamantyl-urea as molecular tools for inhibition of human soluble epoxide hydrolase with picomolar efficacy

Series of 1,3-disubstituted ureas and diadamantyl disubstituted diureas with fluorinated and chlorinated adamantane residues were shown to inhibit human soluble epoxide hydrolase (sEH) with inhibition potency ranging from 40 pM to 9.2 nM. The measured IC50 values for some molecules were below the accuracy limit of the existing in vitro assays. Such an increase in activity was achieved by minimal structural modifications to the molecules of known inhibitors, including 4-[trans-4-(1-adamantylcarbamoylamino)cyclohexyl]oxybenzoic acid. For the chlorinated homologue of the latter the sharp jump in inhibitory activity can be (according to molecular dynamics data) the result of interactions - Cl-π interaction. Considering the extremely high inhibitory activity, acceptable solubility and partial blockage of metabolically sensitive centres in their structures, some compounds are of interest for further in vivo biotesting.

Screening and Biological Evaluation of Soluble Epoxide Hydrolase Inhibitors: Assessing the Role of Hydrophobicity in the Pharmacophore-Guided Search of Novel Hits

The human soluble epoxide hydrolase (sEH) is a bifunctional enzyme that modulates the levels of regulatory epoxy lipids. The hydrolase activity is carried out by a catalytic triad located at the center of a wide L-shaped binding site, which contains two hydrophobic subpockets at both sides. On the basis of these structural features, it can be assumed that desolvation is a major factor in determining the maximal achievable affinity that can be attained for this pocket. Accordingly, hydrophobic descriptors may be better suited to the search of novel hits targeting this enzyme. This study examines the suitability of quantum mechanically derived hydrophobic descriptors in the discovery of novel sEH inhibitors. To this end, three-dimensional quantitative structure-activity relationship (3D-QSAR) pharmacophores were generated by combining electrostatic and steric or alternatively hydrophobic and hydrogen-bond parameters in conjunction with a tailored list of 76 known sEH inhibitors. The pharmacophore models were then validated by using two external sets chosen (i) to rank the potency of four distinct series of compounds and (ii) to discriminate actives from decoys, using in both cases datasets taken from the literature. Finally, a prospective study was performed including a virtual screening of two chemical libraries to identify new potential hits, which were subsequently experimentally tested for their inhibitory activity on human, rat, and mouse sEH. The use of hydrophobic-based descriptors led to the identification of six compounds as inhibitors of the human enzyme with IC50 < 20 nM, including two with IC50 values of 0.4 and 0.7 nM. The results support the use of hydrophobic descriptors as a valuable tool in the search of novel scaffolds that encode a proper hydrophilic/hydrophobic distribution complementary to the target's binding site.

Benzoxazolone-5-Urea Derivatives as Human Soluble Epoxide Hydrolase (sEH) Inhibitors

Inhibition of soluble epoxide hydrolase (sEH) is indicated as a new therapeutic modality against a variety of inflammatory diseases, including metabolic, renal, and cardiovascular disorders. In our ongoing research on sEH inhibitors, we synthesized novel benzoxazolone-5-urea analogues with highly potent sEH inhibitory properties inspired by the crystallographic fragment scaffolds incorporating a single H-bond donor/acceptor pair. The tractable SAR results indicated that the aryl or benzyl fragments flanking the benzoxazolone-urea scaffold conferred potent sEH inhibition, and compounds 31-39 inhibited the sEH activity with IC₅₀ values in the range of 0.39-570 nM. Docking studies and molecular dynamics simulations with the most potent analogue 33 provided valuable insights into potential binding interactions of the inhibitor in the sEH binding region. In conclusion, benzoxazolone-5-ureas furnished with benzyl groups on the urea function can be regarded as novel lead structures, which allow the development of advanced analogues with enhanced properties against sEH.

Synthesis, In Vitro Profiling, and In Vivo Evaluation of Benzohomoadamantane-Based Ureas for Visceral Pain: A New Indication for Soluble Epoxide Hydrolase Inhibitors

The soluble epoxide hydrolase (sEH) has been suggested as a pharmacological target for the treatment of several diseases, including pain-related disorders. Herein, we report further medicinal chemistry around new benzohomoadamantane-based sEH inhibitors (sEHI) in order to improve the drug metabolism and pharmacokinetics properties of a previous hit. After an extensive in vitro screening cascade, molecular modeling, and in vivo pharmacokinetics studies, two candidates were evaluated in vivo in a murine model of capsaicin-induced allodynia. The two compounds showed an anti-allodynic effect in a dose-dependent manner. Moreover, the most potent compound presented robust analgesic efficacy in the cyclophosphamide-induced murine model of cystitis, a well-established model of visceral pain. Overall, these results suggest painful bladder syndrome as a new possible indication for sEHI, opening a new range of applications for them in the visceral pain field.

Repositioning of Quinazolinedione-Based Compounds on Soluble Epoxide Hydrolase (sEH) through 3D Structure-Based Pharmacophore Model-Driven Investigation

The development of new bioactive compounds represents one of the main purposes of the drug discovery process. Various tools can be employed to identify new drug candidates against pharmacologically relevant biological targets, and the search for new approaches and methodologies often represents a critical issue. In this context, in silico drug repositioning procedures are required even more in order to re-evaluate compounds that already showed poor biological results against a specific biological target. 3D structure-based pharmacophoric models, usually built for specific targets to accelerate the identification of new promising compounds, can be employed for drug repositioning campaigns as well. In this work, an in-house library of 190 synthesized compounds was re-evaluated using a 3D structure-based pharmacophoric model developed on soluble epoxide hydrolase (sEH). Among the analyzed compounds, a small set of quinazolinedione-based molecules, originally selected from a virtual combinatorial library and showing poor results when preliminarily investigated against heat shock protein 90 (Hsp90), was successfully repositioned against sEH, accounting the related built 3D structure-based pharmacophoric model. The promising results here obtained highlight the reliability of this computational workflow for accelerating the drug discovery/repositioning processes.

Allosteric regulation of the soluble epoxide hydrolase by nitro fatty acids using a combined experimental and computational approach

The human soluble epoxide hydrolase (hsEH) is a key regulator of epoxy fatty acid (EpFA) metabolism. Inhibition of sEH can maintain endogenous levels of beneficial EpFAs and reduce the levels of their corresponding diol products, thus ameliorating a variety of pathological conditions including cardiovascular, central nervous system and metabolic diseases. The quest for orthosteric drugs that bind directly to the catalytic crevice of hsEH has been prolonged and sustained over the past decades, but the disappointing outcome of clinical trials to date warrants alternative pharmacological approaches. Previously, we have shown that hsEH can be allosterically inhibited by the endogenous electrophilic lipid 15-deoxy-Δ^12,14-Prostaglandin-J₂, via covalent adduction to two cysteines, C423 and C522. In this study, we explore the properties and behaviour of three electrophilic lipids belonging to the class of the nitro fatty acids, namely 9- and 10-nitrooleate and 10-nitrolinoleate. Biochemical and biophysical investigations revealed that, in addition to C423 and C522, nitro fatty acids can covalently bind to additional nucleophilic residues in hsEH C-terminal domain (CTD), two of which predicted in this study to be latent allosteric sites. Systematic mapping of the protein mutational space and evaluation of possible propagation pathways delineated selected residues, both in the allosteric patches and in other regions of the enzyme, envisaged to play a role on allosteric signalling. The responses elicited by the ligands on the covalent adduction sites supports future fragment-based design studies of new allosteric effectors for hsEH with increased efficacy and selectivity.

The biosynthetic origin of ribofuranose in bacterial polysaccharides

Bacterial surface polysaccharides are assembled by glycosyltransferase enzymes that typically use sugar nucleotide or polyprenyl-monophosphosugar activated donors. Characterized representatives exist for many monosaccharides but neither the donor nor the corresponding glycosyltransferases have been definitively identified for ribofuranose residues found in some polysaccharides. Klebsiella pneumoniae O-antigen polysaccharides provided prototypes to identify dual-domain ribofuranosyltransferase proteins catalyzing a two-step reaction sequence. Phosphoribosyl-5-phospho-D-ribosyl-α-1-diphosphate serves as the donor for a glycan acceptor-specific phosphoribosyl transferase (gPRT), and a more promiscuous phosphoribosyl-phosphatase (PRP) then removes the residual 5'-phosphate. The 2.5-Å resolution crystal structure of a dual-domain ribofuranosyltransferase ortholog from Thermobacillus composti revealed a PRP domain that conserves many features of the phosphatase members of the haloacid dehalogenase family, and a gPRT domain that diverges substantially from all previously characterized phosphoribosyl transferases. The gPRT represents a new glycosyltransferase fold conserved in the most abundant ribofuranosyltransferase family.

Computational insights into the known inhibitors of human soluble epoxide hydrolase

Human soluble epoxide hydrolase (hsEH) is involved in the hydrolysis of epoxyeicosatrienoic acids (EETs), which have potent anti-inflammatory properties. Given that EET conversion generates nonbioactive molecules, inhibition of this enzyme would be beneficial. Past decades of work on hsEH inhibitors resulted in numerous potential compounds, of which a hundred hsEH-ligand complexes were crystallized and deposited in the Protein Data Bank (PDB). We analyzed all deposited hsEH-ligand complexes to gain insight into the binding of inhibitors and to provide feedback on the future drug design processes. We also reviewed computationally driven strategies that were used to propose novel hsEH inhibitors.

From the Design to the In Vivo Evaluation of Benzohomoadamantane-Derived Soluble Epoxide Hydrolase Inhibitors for the Treatment of Acute Pancreatitis

The pharmacological inhibition of soluble epoxide hydrolase (sEH) is efficient for the treatment of inflammatory and pain-related diseases. Numerous potent sEH inhibitors (sEHIs) present adamantyl or phenyl moieties, such as the clinical candidates AR9281 or EC5026. Herein, in a new series of sEHIs, these hydrophobic moieties have been merged in a benzohomoadamantane scaffold. Most of the new sEHIs have excellent inhibitory activities against sEH. Molecular dynamics simulations suggested that the addition of an aromatic ring into the adamantane scaffold produced conformational rearrangements in the enzyme to stabilize the aromatic ring of the benzohomoadamantane core. A screening cascade permitted us to select a candidate for an in vivo efficacy study in a murine model of cerulein-induced acute pancreatitis. The administration of 22 improved the health status of the animals and reduced pancreatic damage, demonstrating that the benzohomoadamantane unit is a promising scaffold for the design of novel sEHIs.

Molecular Interactions of Pyrazine-Based Compounds to Proteins

Pyrazine-based compounds are of great importance in medicinal chemistry. Due to their heteroaromatic nature, they uniquely combine properties of heteroatoms (polar interactions) with the properties of aromatic moieties (nonpolar interactions). This review summarizes results of a systematic analysis of RCSB PDB database focused on important binding interactions of pyrazine-based ligands cocrystallized in protein targets. The most frequent interaction of pyrazine was hydrogen bond to pyrazine nitrogen atom as an acceptor, followed by weak hydrogen bond with pyrazine hydrogen as donor. We also identified intramolecular hydrogen bonds within pyrazine ligands, π-interactions, coordination to metal ions, and few halogen bonds in chloropyrazines. In many cases the binding mode of the pyrazine fragment was complex, involving a combination of several interactions. We conclude that pyrazine as a molecular fragment should not be perceived as a simple aromatic isostere but rather as a readily interacting moiety of drug-like molecules with high potential for interactions to proteins.

2-Oxaadamant-1-yl Ureas as Soluble Epoxide Hydrolase Inhibitors: In Vivo Evaluation in a Murine Model of Acute Pancreatitis

In vivo pharmacological inhibition of soluble epoxide hydrolase (sEH) reduces inflammatory diseases, including acute pancreatitis (AP). Adamantyl ureas are very potent sEH inhibitors, but the lipophilicity and metabolism of the adamantane group compromise their overall usefulness. Herein, we report that the replacement of a methylene unit of the adamantane group by an oxygen atom increases the solubility, permeability, and stability of three series of urea-based sEH inhibitors. Most of these oxa-analogues are nanomolar inhibitors of both the human and murine sEH. Molecular dynamics simulations rationalize the molecular basis for their activity and suggest that the presence of the oxygen atom on the adamantane scaffold results in active site rearrangements to establish a weak hydrogen bond. The 2-oxaadamantane 22, which has a good solubility, microsomal stability, and selectivity for sEH, was selected for further in vitro and in vivo studies in models of cerulein-induced AP. Both in prophylactic and treatment studies, 22 diminished the overexpression of inflammatory and endoplasmic reticulum stress markers induced by cerulein and reduced the pancreatic damage.

FragMAX: the fragment-screening platform at the MAX IV Laboratory

Advances in synchrotron storage rings and beamline automation have pushed data-collection rates to thousands of data sets per week. With this increase in throughput, massive projects such as in-crystal fragment screening have become accessible to a larger number of research groups. The quality of support offered at large-scale facilities allows medicinal chemistry-focused or biochemistry-focused groups to supplement their research with structural biology. Preparing the experiment, analysing multiple data sets and prospecting for interesting complexes of protein and fragments require, for both newcomers and experienced users, efficient management of the project and extensive computational power for data processing and structure refinement. Here, FragMAX, a new complete platform for fragment screening at the BioMAX beamline of the MAX IV Laboratory, is described. The ways in which users are assisted in X-ray-based fragment screenings and in which the fourth-generation storage ring available at the facility is best exploited are also described.

Combining Machine Learning and Enhanced Sampling Techniques for Efficient and Accurate Calculation of Absolute Binding Free Energies

Calculating absolute binding free energies is challenging and important. In this paper, we test some recently developed metadynamics-based methods and develop a new combination with a Hamiltonian replica-exchange approach. The methods were tested on 18 chemically diverse ligands with a wide range of different binding affinities to a complex target; namely, human soluble epoxide hydrolase. The results suggest that metadynamics with a funnel-shaped restraint can be used to calculate, in a computationally affordable and relatively accurate way, the absolute binding free energy for small fragments. When used in combination with an optimal pathlike variable obtained using machine learning or with the Hamiltonian replica-exchange algorithm SWISH, this method can achieve reasonably accurate results for increasingly complex ligands, with a good balance of computational cost and speed. An additional benefit of using the combination of metadynamics and SWISH is that it also provides useful information about the role of water in the binding mechanism.

Fragment Pose Prediction Using Non-equilibrium Candidate Monte Carlo and Molecular Dynamics Simulations

Part of early stage drug discovery involves determining how molecules may bind to the target protein. Through understanding where and how molecules bind, chemists can begin to build ideas on how to design improvements to increase binding affinities. In this retrospective study, we compare how computational approaches like docking, molecular dynamics (MD) simulations, and a non-equilibrium candidate Monte Carlo (NCMC)-based method (NCMC + MD) perform in predicting binding modes for a set of 12 fragment-like molecules, which bind to soluble epoxide hydrolase. We evaluate each method's effectiveness in identifying the dominant binding mode and finding additional binding modes (if any). Then, we compare our predicted binding modes to experimentally obtained X-ray crystal structures. We dock each of the 12 small molecules into the apo-protein crystal structure and then run simulations up to 1 μs each. Small and fragment-like molecules likely have smaller energy barriers separating different binding modes by virtue of relatively fewer and weaker interactions relative to drug-like molecules and thus likely undergo more rapid binding mode transitions. We expect, thus, to see more rapid transitions between binding modes in our study. Following this, we build Markov State Models to define our stable ligand binding modes. We investigate if adequate sampling of ligand binding modes and transitions between them can occur at the microsecond timescale using traditional MD or a hybrid NCMC+MD simulation approach. Our findings suggest that even with small fragment-like molecules, we fail to sample all the crystallographic binding modes using microsecond MD simulations, but using NCMC+MD, we have better success in sampling the crystal structure while obtaining the correct populations.

Protostane-type triterpenoids as natural soluble epoxide hydrolase inhibitors: Inhibition potentials and molecular dynamics

The inhibition of soluble epoxide hydrolase (sEH) is a promising therapeutic approach to treat inflammation and other disorders. In our present investigation on searching for sEH inhibitors from traditional Chinese medicines, we found that Alisma orientale displayed inhibition of sEH. We constructed a small library of protostane-type triterpenoids (1-25) isolated from A. orientale, and screened their inhibitory activities. Alismanin B (1), 11-deoxy-25-anhydro alisol E (4), 11-deoxy alisol B (5), and 25-O-ethyl alisol A (15) displayed concentration-dependently inhibitory activities against sEH with IC₅₀ values from 3.40 ± 0.57 μM to 9.57 ± 0.88 μM. 11-Deoxy-25-anhydro alisol E (4) and 11-deoxy alisol B (5) were defined as mixed-type competitive inhibitors with Ki values of 12.6 and 3.48 μM, respectively, based on the result of inhibition kinetics. The potential interaction mechanism of 11-deoxy alisol B (5) with sEH was analyzed by molecular docking and molecular dynamics, revealing that amino acid residues Trp336 and Tyr466 were vital for its inhibitory activity.

Discovery of the First in Vivo Active Inhibitors of the Soluble Epoxide Hydrolase Phosphatase Domain

The emerging pharmacological target soluble epoxide hydrolase (sEH) is a bifunctional enzyme exhibiting two different catalytic activities that are located in two distinct domains. Although the physiological role of the C-terminal hydrolase domain is well-investigated, little is known about its phosphatase activity, located in the N-terminal phosphatase domain of sEH (sEH-P). Herein we report the discovery and optimization of the first inhibitor of human and rat sEH-P that is applicable in vivo. X-ray structure analysis of the sEH phosphatase domain complexed with an inhibitor provides insights in the molecular basis of small-molecule sEH-P inhibition and helps to rationalize the structure-activity relationships. 4-(4-(3,4-Dichlorophenyl)-5-phenyloxazol-2-yl)butanoic acid (22b, SWE101) has an excellent pharmacokinetic and pharmacodynamic profile in rats and enables the investigation of the physiological and pathophysiological role of sEH-P in vivo.

15-deoxy-Δ12,14-Prostaglandin J2 inhibits human soluble epoxide hydrolase by a dual orthosteric and allosteric mechanism

Human soluble epoxide hydrolase (hsEH) is an enzyme responsible for the inactivation of bioactive epoxy fatty acids, and its inhibition is emerging as a promising therapeutical strategy to target hypertension, cardiovascular disease, pain and insulin sensitivity. Here, we uncover the molecular bases of hsEH inhibition mediated by the endogenous 15-deoxy-Δ12,14-Prostaglandin J2 (15d-PGJ2). Our data reveal a dual inhibitory mechanism, whereby hsEH can be inhibited by reversible docking of 15d-PGJ2 in the catalytic pocket, as well as by covalent locking of the same compound onto cysteine residues C423 and C522, remote to the active site. Biophysical characterisations allied with in silico investigations indicate that the covalent modification of the reactive cysteines may be part of a hitherto undiscovered allosteric regulatory mechanism of the enzyme. This study provides insights into the molecular modes of inhibition of hsEH epoxy-hydrolytic activity and paves the way for the development of new allosteric inhibitors.

Elucidating the druggability of the human proteome with eFindSite

Identifying the viability of protein targets is one of the preliminary steps of drug discovery. Determining the ability of a protein to bind drugs in order to modulate its function, termed the druggability, requires a non-trivial amount of time and resources. Inability to properly measure druggability has accounted for a significant portion of failures in drug discovery. This problem is only further exacerbated by the large sample space of proteins involved in human diseases. With these barriers, the druggability space within the human proteome remains unexplored and has made it difficult to develop drugs for numerous diseases. Hence, we present a new feature developed in eFindSite that employs supervised machine learning to predict the druggability of a given protein. Benchmarking calculations against the Non-Redundant data set of Druggable and Less Druggable binding sites demonstrate that an AUC for druggability prediction with eFindSite is as high as 0.88. With eFindSite, we elucidated the human druggability space to be 10,191 proteins. Considering the disease space from the Open Targets Platform and excluding already known targets from the predicted data set reveal 2731 potentially novel therapeutic targets. eFindSite is freely available as a stand-alone software at https://github.com/michal-brylinski/efindsite .

Structure-based virtual screening for novel chymase inhibitors by in silico fragment mapping

The term chymase refers to a family of chymotrypsin-like serine proteases stored within the secretory granules of mast cells. Recently, a variety of small molecule inhibitors for chymase have been developed with a primary focus on the treatment of cardiovascular diseases. Despite the expected therapeutic benefit of these chymase inhibitors, they have not been used clinically. Here, we attempted to identify new chymase inhibitors using a multistep structure-based virtual screening protocol combined with our knowledge-based in silico fragment mapping technique. The mapping procedure identified fragments with novel modes of interaction at the oxyanion hole of chymase. Next, we constructed a three-dimensional (3D) pharmacophore model and retrieved eight candidate chymase inhibitors from a commercial database that included approximately five million compounds. This selection was achieved using a multistep virtual screening protocol, which combined a 3D pharmacophore-based search, docking calculations, and analyses of binding free energy. One of the eight compounds exhibited concentration-dependent chymase inhibitory activity, which could be further optimized to develop more potent chymase inhibitors.

Precipitant-ligand exchange technique reveals the ADP binding mode in Mycobacterium tuberculosis dethiobiotin synthetase

Dethiobiotin synthetase from Mycobacterium tuberculosis (MtDTBS) is a promising antituberculosis drug target. Small-molecule inhibitors that target MtDTBS provide a route towards new therapeutics for the treatment of antibiotic-resistant tuberculosis. Adenosine diphosphate (ADP) is an inhibitor of MtDTBS; however, structural studies into its mechanism of inhibition have been unsuccessful owing to competitive binding to the enzyme by crystallographic precipitants such as citrate and sulfate. Here, a crystallographic technique termed precipitant-ligand exchange has been developed to exchange protein-bound precipitants with ligands of interest. Proof of concept for the exchange method was demonstrated using cytidine triphosphate (CTP), which adopted the same binding mechanism as that obtained with traditional crystal-soaking techniques. Precipitant-ligand exchange also yielded the previously intractable structure of MtDTBS in complex with ADP solved to 2.4 Å resolution. This result demonstrates the utility of precipitant-ligand exchange, which may be widely applicable to protein crystallography.

Expression, purification, and characterisation of human soluble Epoxide Hydrolase (hsEH) and of its functional C-terminal domain

The human soluble Epoxide Hydrolase (hsEH) is an enzyme involved in the hydrolysis of endogenous anti-inflammatory and cardio-protective signalling mediators known as epoxyeicosatrienoic acids (EETs). EETs’ conversion into the corresponding diols by hsEH generates non-bioactive molecules, thereby the enzyme inhibition would be expected to enhance the EETs bioavailability, and their beneficial properties. Numerous inhibitors have been developed to target the enzyme, some of which are showing promising antihypertensive and anti-inflammatory properties in vivo. Thus far, the preparation of the recombinant enzyme for enzymatic and structural in vitro studies has been performed mainly using a baculovirus expression system. More recently, it was reported that the enzyme could be exogenously expressed and isolated from E. coli, although limited amounts of active protein were obtained. We herein describe two novel methods to yield pure recombinant enzyme. The first describes the expression and purification of the full-length enzyme from eukaryotic cells HEK293-F, whilst the second concerns the C-terminal domain of hsEH obtained from the cost-effective and rapid E. coli prokaryotic system. The two methods successfully generated satisfactory amounts of functional enzyme, with virtually identical enzymatic activity. Overall, the protocols described in this paper can be employed for the recombinant expression and purification of active hsEH, to be used in future biomedical investigations and for high-throughput screening of inhibitors for potential use in the treatment of cardiovascular disease.

•

hsEH is a key regulator of cardiovascular homeostasis.

•

A HEK293-F mammalian expression system for hsEH full-length (FL) was developed.

•

An E. coli expression system for the hsEH C-terminal Domain (CTD) was established.

•

Both proteins exhibited the same enzymatic specific activity in vitro.

•

The CTD preparation provides benefits of easy operation, and high yield and purity.

In silico fragment-mapping method: a new tool for fragment-based/structure-based drug discovery

Here, we propose an in silico fragment-mapping method as a potential tool for fragment-based/structure-based drug discovery (FBDD/SBDD). For this method, we created a database named Canonical Subsite-Fragment DataBase (CSFDB) and developed a knowledge-based fragment-mapping program, Fsubsite. CSFDB consists of various pairs of subsite-fragments derived from X-ray crystal structures of known protein-ligand complexes. Using three-dimensional similarity-matching between subsites on one protein and another, Fsubsite compares the surface of a target protein with all subsites in CSFDB. When a local topography similar to the subsite is found on the surface, Fsubsite places a fragment combined with the subsite in CSFDB on the target protein. For validation purposes, we applied the method to the apo-structure of cyclin-dependent kinase 2 (CDK2) and identified four compounds containing three mapped fragments that existed in the list of known inhibitors of CDK2. Next, the utility of our fragment-mapping method for fragment-growing was examined on the complex structure of tRNA-guanine transglycosylase with a small ligand. Fsubsite mapped appropriate fragments on the same position as the binding ligand or in the vicinity of the ligand. Finally, a 3D-pharmacophore model was constructed from the fragments mapped on the apo-structure of heat shock protein 90-α (HSP90α). Then, 3D pharmacophore-based virtual screening was carried out using a commercially available compound database. The resultant hit compounds were very similar to a known ligand of HSP90α. As a result of these findings, this in silico fragment-mapping method seems to be a useful tool for computational FBDD and SBDD.

In silico fragment-based drug design with SEED

We report on two fragment-based drug design protocols, SEED2XR and ALTA, which start by high-throughput docking. SEED2XR is a two-stage protocol for fragment-based drug design. The first stage is in silico and consists of the automatic docking of 10³-10⁴ fragments using SEED, which requires about 1 s per fragment. SEED is a docking software developed specifically for fragment docking and binding energy evaluation by a force field with implicit solvent. In the second stage of SEED2XR, the 10-10² fragments with the most favorable predicted binding energies are validated by protein X-ray crystallography. The recent applications of SEED2XR to bromodomains demonstrate that the whole SEED2XR protocol can be carried out in about a week of working time, with hit rates ranging from 10% to 40%. Information on fragment-target interactions generated by the SEED2XR protocol or directly from SEED docking has been used for the discovery of hundreds of hits. ALTA is a computational protocol for screening which identifies candidate ligands that preserve the interactions between the optimal SEED fragments and the protein target. Medicinal chemistry optimization of ligands predicted by ALTA has resulted in pre-clinical candidates for protein kinases and bromodomains. The high-throughput, very low cost, sustainability, and high hit rate of the SEED-based protocols, unreachable by purely experimental techniques, make them perfectly suitable for both academic and industrial drug discovery research.

Structural biology data archiving - where we are and what lies ahead

For almost 50 years, structural biology has endeavoured to conserve and share its experimental data and their interpretations (usually, atomistic models) through global public archives such as the Protein Data Bank, Electron Microscopy Data Bank and Biological Magnetic Resonance Data Bank (BMRB). These archives are treasure troves of freely accessible data that document our quest for molecular or atomic understanding of biological function and processes in health and disease. They have prepared the field to tackle new archiving challenges as more and more (combinations of) techniques are being utilized to elucidate structure at ever increasing length scales. Furthermore, the field has made substantial efforts to develop validation methods that help users to assess the reliability of structures and to identify the most appropriate data for their needs. In this Review, we present an overview of public data archives in structural biology and discuss the importance of validation for users and producers of structural data. Finally, we sketch our efforts to integrate structural data with bioimaging data and with other sources of biological data. This will make relevant structural information available and more easily discoverable for a wide range of scientists.

Phosphatase activity of soluble epoxide hydrolase

Soluble epoxide hydrolase (sEH) is a bifunctional enzyme that exhibits lipid epoxide hydrolase (sEH-H) and lipid phosphatase activity (sEH-P), with each being located in its own distinct domain. While the epoxide hydrolase activity is well-investigated, the role of the phosphatase domain remains unclear. This article briefly summarizes the evolution, structure and SNPs of the human sEH, with a special focus on the function and postulated role of the N-terminal domain of sEH. Furthermore, the article provides an overview of tools to study sEH-P.

Guidelines for the successful generation of protein-ligand complex crystals

With continuous technical improvements at synchrotron facilities, data-collection rates have increased dramatically. This makes it possible to collect diffraction data for hundreds of protein-ligand complexes within a day, provided that a suitable crystal system is at hand. However, developing a suitable crystal system can prove challenging, exceeding the timescale of data collection by several orders of magnitude. Firstly, a useful crystallization construct of the protein of interest needs to be chosen and its expression and purification optimized, before screening for suitable crystallization and soaking conditions can start. This article reviews recent publications analysing large data sets of crystallization trials, with the aim of identifying factors that do or do not make a good crystallization construct, and gives guidance in the design of an expression construct. It provides an overview of common protein-expression systems, addresses how ligand binding can be both help and hindrance for protein purification, and describes ligand co-crystallization and soaking, with an emphasis on troubleshooting.

Probing the orientation of inhibitor and epoxy-eicosatrienoic acid binding in the active site of soluble epoxide hydrolase

Soluble epoxide hydrolase (sEH) is an important therapeutic target of many diseases, such as chronic obstructive pulmonary disease (COPD) and diabetic neuropathic pain. It acts by hydrolyzing and thus regulating specific bioactive long chain polyunsaturated fatty acid epoxides (lcPUFA), like epoxyeicosatrienoic acids (EETs). To better predict which epoxides could be hydrolyzed by sEH, one needs to dissect the important factors and structural requirements that govern the binding of the substrates to sEH. This knowledge allows further exploration of the physiological role played by sEH. Unfortunately, a crystal structure of sEH with a substrate bound has not yet been reported. In this report, new photoaffinity mimics of a sEH inhibitor and EET regioisomers were prepared and used in combination with peptide sequencing and computational modeling, to identify the binding orientation of different regioisomers and enantiomers of EETs into the catalytic cavity of sEH. Results indicate that the stereochemistry of the epoxide plays a crucial role in dictating the binding orientation of the substrate.

Challenges of docking in large, flexible and promiscuous binding sites

After decades of work, the correct determination of the binding mode of a small molecule into a target protein is still a challenging problem, whose difficulty depends on: (i) the sizes of the binding site and the ligand; (ii) the flexibility of both interacting partners, and (iii) the differential solvation of bound and unbound partners. We have evaluated the performance of standard rigid(receptor)/flexible(ligand) docking approaches with respect to last-generation fully flexible docking methods to obtain reasonable poses in a very challenging case: soluble Epoxide Hydrolase (sEH), a flexible protein showing different binding sites. We found that full description of the flexibility of both protein and ligand and accurate description of solvation leads to significant improvement in the ability of docking to reproduce well known binding modes, and at the same time capture the intrinsic binding promiscuity of the protein.

High-Throughput Crystallography: Reliable and Efficient Identification of Fragment Hits

Today the identification of lead structures for drug development often starts from small fragment-like molecules raising the chances to find compounds that successfully pass clinical trials. At the heart of the screening for fragments binding to a specific target, crystallography delivers structural information essential for subsequent drug design. While it is common to search for bound ligands in electron densities calculated directly after an initial refinement cycle, we raise the important question whether this strategy is viable for fragments characterized by low affinities. Here, we describe and provide a collection of high-quality diffraction data obtained from 364 protein crystals treated with diverse fragments. Subsequent data analysis showed that ∼25% of all hits would have been missed without further refining the resulting structures. To enable fast and reliable hit identification, we have designed an automated refinement pipeline that will inspire the development of optimized tools facilitating the successful application of fragment-based methods.

Applications of (19)F-NMR in Fragment-Based Drug Discovery

(19)F-NMR has proved to be a valuable tool in fragment-based drug discovery. Its applications include screening libraries of fluorinated fragments, assessing competition among elaborated fragments and identifying the binding poses of promising hits. By observing fluorine in both the ligand and the target protein, useful information can be obtained on not only the binding pose but also the dynamics of ligand-protein interactions. These applications of (19)F-NMR will be illustrated in this review with studies from our fragment-based drug discovery campaigns against protein targets in parasitic and infectious diseases.

Six Biophysical Screening Methods Miss a Large Proportion of Crystallographically Discovered Fragment Hits: A Case Study

Fragment-based lead discovery (FBLD) has become a pillar in drug development. Typical applications of this method comprise at least two biophysical screens as prefilter and a follow-up crystallographic experiment on a subset of fragments. Clearly, structural information is pivotal in FBLD, but a key question is whether such a screening cascade strategy will retrieve the majority of fragment-bound structures. We therefore set out to screen 361 fragments for binding to endothiapepsin, a representative of the challenging group of aspartic proteases, employing six screening techniques and crystallography in parallel. Crystallography resulted in the very high number of 71 structures. Yet alarmingly, 44% of these hits were not detected by any biophysical screening approach. Moreover, any screening cascade, building on the results from two or more screening methods, would have failed to predict at least 73% of these hits. We thus conclude that, at least in the present case, the frequently applied biophysical prescreening filters deteriorate the number of possible X-ray hits while only the immediate use of crystallography enables exhaustive retrieval of a maximum of fragment structures, which represent a rich source guiding hit-to-lead-to-drug evolution.

Stacking with No Planarity?

This viewpoint describes the results obtained from matched molecular pair analyses and quantum mechanics calculations that show unsaturated rings found in drug-like molecules may be replaced with their saturated counterparts without losing potency even if they are engaged in stacking interactions with the side chains of aromatic residues.

Compound high-quality criteria: a new vision to guide the development of drugs, current situation

For several decades, the pharmaceutical industry has suffered due to major issues such as reductions of the number of FDA approved drugs and biologics. Several analyses have been highlighted that the 'druglikeness' is one of the strategies to improve succeed rates of screening such as, for instance, high-throughput screening (HTS), and then hits (as starting point), leads and clinical candidates. It is clear that the improvement of compound quality accelerates the drug discovery projects. The monitoring of several indices to avoid 'molecular obesity' (ADMET problems) of final drugs from good-quality 'low-fat' starting points represents today a powerful strategy of optimization process. The development of the new guides to find drugs highlighting attempts at improving the attrition rate from hits to final medicines by focusing on how to improve the druggability of hits, leads and drugs during the drug discovery process represents a key approach to design next better generation of medicines.

Fragment Screening of Soluble Epoxide Hydrolase for Lead Generation-Structure-Based Hit Evaluation and Chemistry Exploration

Soluble epoxide hydrolase (sEH) is involved in the regulation of many biological processes by metabolizing the key bioactive lipid mediator, epoxyeicosatrienoic acids. For the development of sEH inhibitors with improved physicochemical properties, we performed both a fragment screening and a high-throughput screening aiming at an integrated hit evaluation and lead generation. Followed by a joint dose-response analysis to confirm the hits, the identified actives were then effectively triaged by a structure-based hit-classification approach to three prioritized series. Two distinct scaffolds were identified as tractable starting points for potential lead chemistry work. The oxoindoline series bind at the right-hand side of the active-site pocket with hydrogen bonds to the protein. The 2-phenylbenzimidazole-4-sulfonamide series bind at the central channel with significant induced fit, which has not been previously reported. On the basis of the encouraging initial results, we envision that a new lead series with improved properties could be generated if a vector is found that could merge the cyclohexyl functionality of the oxoindoline series with the trifluoromethyl moiety of the 2-phenylbenzimidazole-4-sulfonamide series.

MeshAndCollect: an automated multi-crystal data-collection workflow for synchrotron macromolecular crystallography beamlines

Here, an automated procedure is described to identify the positions of many cryocooled crystals mounted on the same sample holder, to rapidly predict and rank their relative diffraction strengths and to collect partial X-ray diffraction data sets from as many of the crystals as desired. Subsequent hierarchical cluster analysis then allows the best combination of partial data sets, optimizing the quality of the final data set obtained. The results of applying the method developed to various systems and scenarios including the compilation of a complete data set from tiny crystals of the membrane protein bacteriorhodopsin and the collection of data sets for successful structure determination using the single-wavelength anomalous dispersion technique are also presented.

Ion mobility coupled to native mass spectrometry as a relevant tool to investigate extremely small ligand-induced conformational changes

Native and ion-mobility mass spectrometry reveal the conformational evolution over time of a peptide deformylase binding different ligands, which is consistent with slow-tight inhibition of the enzyme.