Advantages of crystallographic fragment screening: functional and mechanistic insights from a powerful platform for efficient drug discovery

ProteinsPlus: a publicly available resource for protein structure mining

The openly accessible ProteinsPlus web server at https://proteins.plus is a unique resource enabling protein structure mining and modeling, focussing on protein-ligand interactions. Since its launch in 2017, the number of tools is steadily increasing. Currently, the server comprises six methods for protein structure analysis, four tools for mining the Protein Data Bank (PDB), and five prediction approaches regarding protein-ligand complex modeling. Users can use experimental structures from the PDB or computationally predicted structures from the AlphaFold Protein Structure Database as starting points. Alternatively, they can upload individual structure files. Recent updates include novel methods for detecting binding sites and predicting solvent channels in crystallographic structures, as well as updates of tools for protein-ligand interaction depiction in 2D and binding site mining. Given these updates, we present a real-life application scenario that underpins the novelties and applicability of the web server's tools for modern structure-based design projects. It also highlights the next steps for the web server, which will be redesigned using a different technology stack to improve the inter-usability of the tools, ease maintainability, and make it future-proof.

Status and perspective of protein crystallography at the first multi-bend achromat based synchrotron MAX IV

The first multi-bend achromat based synchrotron MAX IV operates two protein crystallography beamlines, BioMAX and MicroMAX. BioMAX is designed as a versatile, stable, high-throughput beamline catering for most protein crystallography experiments. MicroMAX is a more ambitious beamline dedicated to serial crystallography including time-resolved experiments. Both beamlines exploit the special characteristics of fourth-generation beamlines provided by the 3 GeV ring of MAX IV. In addition, the fragment-based drug discovery platform, FragMAX, is hosted and, at the FemtoMAX beamline, protein diffraction experiments exploring ultrafast time resolution can be performed. A technical and operational overview of the different beamlines and the platform is given as well as an outlook for protein crystallography embedded in the wider possibilities that MAX IV offers to users in the life sciences.

NMR2-Based Drug Discovery Pipeline Presented on the Oncogenic Protein KRAS

Fragment-based drug discovery has emerged as a powerful approach for developing therapeutics against challenging targets, including the GTPase KRAS. Here, we report an NMR-based screening campaign employing state-of-the-art techniques to evaluate a library of 890 fragments against the oncogenic KRAS G12V mutant bound to GMP-PNP. Further HSQC titration experiments identified hits with low millimolar affinities binding within the SI/SII switch region, which forms the binding interface for the effector proteins. To elucidate the binding modes, we applied NMR molecular replacement (NMR2) structure calculations, bypassing the need for a conventional protein resonance assignment. Traditionally, NMR2 relies on isotope-filtered nuclear Overhauser effect spectroscopy experiments requiring double-labeled [13C,15N]-protein. We introduce a cost-efficient alternative using a relaxation-based filter that eliminates isotope labeling while preserving structural accuracy. Validation against standard isotopically labeled workflows confirmed the equivalence of the derived protein-ligand structures. This approach enabled the determination of 12 NMR2 KRAS-fragment complex structures, providing critical insights into structure-activity relationships to guide ligand optimization. These results demonstrate the streamlined integration of NMR2 into a fragment-based drug discovery pipeline composed of screening, binding characterization, and rapid structural elucidation with or without isotopic labeling.

Deciphering the enigmas of non-nucleoside reverse transcriptase inhibitors (NNRTIs): A medicinal chemistry expedition towards combating HIV drug resistance

The pivotal involvement of reverse transcriptase activity in the pathogenesis of the progressive HIV virus has stimulated gradual advancements in drug discovery initiatives spanning three decades. Consequently, nonnucleoside reverse transcriptase inhibitors (NNRTIs) have emerged as a preeminent category of therapeutic agents for HIV management. Academic institutions and pharmaceutical companies have developed numerous NNRTIs, an essential component of antiretroviral therapy. Six NNRTIs have received Food and Drug Administration approval and are widely used in clinical practice, significantly improving the quality of HIV patients. However, the rapid emergence of drug resistance has limited the effectiveness of these medications, underscoring the necessity for perpetual research and development of novel therapeutic alternatives. To supplement the existing literatures on NNRTIs, a comprehensive review has been compiled to synthesize this extensive dataset into a comprehensible format for the medicinal chemistry community. In this review, a thorough investigation and meticulous analysis were conducted on the progressions achieved in NNRTIs within the past 8 years (2016-2023), and the experiences and insights gained in the development of inhibitors with varying chemical structures were also summarized. The provision of a crucial point of reference for the development of wide-ranging anti-HIV medications is anticipated.

Combining MicroED and native mass spectrometry for structural discovery of enzyme-biosynthetic inhibitor complexes

With the goal of accelerating the discovery of small molecule-protein complexes, we leverage fast, low-dose, event based electron counting microcrystal electron diffraction (MicroED) data collection and native mass spectrometry. This approach resolves structures of the epoxide-based cysteine protease inhibitor, and natural product, E-64, and its biosynthetic analogs bound to the model cysteine protease, papain. The combined structural power of MicroED and the analytical capabilities of native mass spectrometry (ED-MS) allows assignment of papain structures bound to E-64-like ligands with data obtained from crystal slurries soaked with mixtures of known inhibitors, and crude biosynthetic reactions. ED-MS further discriminates the highest-affinity ligand soaked into microcrystals from a broad inhibitor cocktail, and identifies multiple similarly high-affinity ligands soaked into microcrystals simultaneously. This extends to libraries of printed ligands dispensed directly onto TEM grids and later soaked with papain microcrystal slurries. ED-MS identifies papain binding to its preferred natural products, by showing that two analogues of E-64 outcompete others in binding to papain crystals, and by detecting papain bound to E-64 and an analogue from crude biosynthetic reactions, without purification. This illustrates the utility of ED-MS for natural product ligand discovery and for structure-based screening of small molecule binders to macromolecular targets.

Ligand Screening and Discovery using Cocktail Soaking and Automated MicroED

Cocktail soaking using single-crystal X-ray diffraction (SC-XRD) has previously allowed high-throughput crystallographic screening of ligands against protein targets. However, protein microcrystals are not amenable to this approach if they are too small to yield strong diffraction patterns. In this study, we developed a workflow integrating cocktail soaking with automated microcrystal electron diffraction (MicroED) to allow rapid ligand screening, structure determination, and binding analysis directly from microcrystals. This can improve the successful hit rate, because binding is often more efficient when smaller crystals are soaked in the ligand. The approach was validated with known ligands of thermolysin and identified novel binding interactions for ligands of proteinase K. The structures of multiple protein-ligand complexes, including ligands with weak binding affinities, could be solved quickly. Their estimated relative binding affinities are in good agreement with previous work and independent microscale thermophoresis (MST) measurements.

Fragment Screening Identifies Novel Allosteric Binders and Binding Sites in the VHR (DUSP3) Phosphatase

The human Vaccinia H1-related phosphatase (VHR; DUSP3) is a critical positive regulator of the innate immune response. Recent studies suggest that inhibiting VHR could be beneficial in treating sepsis and septic shock. VHR belongs to the superfamily of protein tyrosine phosphatases (PTPs), a large class of enzymes that are notoriously difficult to target with small molecules. Fragment-based drug discovery (FBDD) has emerged as an effective strategy for generating potent ligands, even for challenging drug targets. Here, we present a fluorine NMR-based discovery platform for identifying fragments that bind to VHR. This platform encompasses automated library assembly, mixture formation, quantitative material transfer, fluorine NMR screening, and biophysical hit confirmation. We demonstrate that this streamlined, integrated screening workflow produces validated hits with diverse chemical matter and tangible structure-activity relationships (SAR). Crystal structures yielded detailed information on the fragment-protein interactions and provide a basis for future structurally enabled ligand optimization. Notably, we discovered novel ligand binding sites on VHR, distant from the conserved active site, facilitating the generation of selective VHR modulators. This fragment discovery platform can be applied to other PTPs and holds significant potential for identifying potent and selective ligands.

Structural Biology for Target Identification and Validation

Structural biology is catalyzing a paradigm shift in drug discovery towards rational approaches in target identification and validation. Leveraging structural insights obtained through cryo-EM or X-ray crystallography not only enhances the efficiency of drug discovery projects in terms of time and cost, but also significantly improves the likelihood of achieving market approval.Initiating a successful project necessitates more than just a robust package for target credentialing; it demands a comprehensive strategy for the identification and optimization of potential drugs. The critical evaluation of target druggability is markedly enhanced when supported by experimentally derived structural information. This nuanced approach ensures a more thorough understanding of the technical feasibility of drug development from the project's inception.

TDP1 represents a promising therapeutic target for overcoming tumor resistance to chemotherapeutic agents: progress and potential

Tyrosyl-DNA phosphodiesterase 1 (TDP1) is an enzyme that plays a crucial role in repairing DNA lesions caused by the entrapment of DNA topoisomerase IB (TOP1)-DNA break-associated crosslinks. TDP1 inhibitors exhibit synergistic effects with TOP1 inhibitors in cancer cells, effectively overcoming resistance to TOP1 inhibitors. Therefore, this approach presents a promising strategy for reversing tumor resistance to TOP1 inhibitors. This review comprehensively outlines the structural and biological features of TDP1, the substrates involved in its catalytic hydrolysis, and its potential as a therapeutic target in oncology. Additionally, we summarize the various screening methods used to identify TDP1 inhibitors, alongside the latest advancements in TDP1 inhibitor research.

Identification of Potent CHK2 Inhibitors‐Modulators for Therapeutic Application in Cancer: A Machine Learning Integrated Fragment‐Based Drug Design Approach

The CHK2 protein regulates the cell division cycle and responds to DNA damage. Additionally, it facilitates the repair of DNA damage and maintains the integrity of its biological processes. Dysregulation of the CHK2 protein is associated with a predisposition to harmful diseases. The current research protocol was designed to identify novel hit molecules as CHK2 inhibitors and disrupt the normal biological function of the CHK2 protein via a fragment‐based drug discovery approach. The protocol involved generating fragments using the MacFrag tool, followed by a chemical similarity search utilizing RDKit to identify fragment molecules analogous to previously established CHK2 inhibitor scaffolds. The bioactive molecules were constructed using the Fragmenstein tool, followed by molecular docking simulations to investigate their binding affinity. In addition, pharmacokinetic properties were analyzed, and a molecular dynamics simulation study was conducted to assess the stability of selected compounds with CHK2 protein. Finally, five novel compounds were identified as excellent CHK2 inhibitors through the FBDD and show good binding interactions at active sites of CHK2 with beneficial ADMET properties. This research work presents novel CHK2 inhibitor molecules that have the potential to be utilized in drug discovery, serving as key leads for future advancements in healthcare industries and sectors.

Identifying novel chemical matter against the Chikungunya virus nsP3 macrodomain through crystallographic fragment screening

Chikungunya virus (CHIKV) causes severe fever, rash and debilitating joint pain that can last for months1,2or even years. Millions of people have been infected with CHIKV, mostly in low and middle-income countries, and the virus continues to spread into new areas due to the geographical expansion of its mosquito hosts. Its genome encodes a macrodomain, which functions as an ADP-ribosyl hydrolase, removing ADPr from viral and host-cell proteins interfering with the innate immune response. Mutational studies have shown that the CHIKV nsP3 macrodomain is necessary for viral replication, making it a potential target for the development of antiviral therapeutics. We, therefore, performed a high-throughput crystallographic fragment screen against the CHIKV nsP3 macrodomain, yielding 109 fragment hits covering the ADPr-binding site and two adjacent subsites that are absent in the homologous macrodomain of SARS-CoV-2 but may be present in other alphaviruses, such as Venezuelan equine encephalitis virus (VEEV) and eastern equine encephalitis virus (EEEV). Finally, a subset of overlapping fragments was used to manually design three fragment merges covering the adenine and oxyanion subsites. The rich dataset of chemical matter and structural information discovered from this fragment screen is publicly available and can be used as a starting point for developing a CHIKV nsP3 macrodomain inhibitor.

Development of a Crystallographic Screening to Identify Sudan Virus VP40 Ligands

The matrix protein VP40 of the highly pathogenic Sudan virus (genus Orthoebolavirus) is a multifunctional protein responsible for the recruitment of viral nucleocapsids to the plasma membrane and the budding of infectious virions. In addition to its role in assembly, VP40 also downregulates viral genome replication and transcription. VP40's existence in various homo-oligomeric states is presumed to underpin its diverse functional capabilities during the viral life cycle. Given the absence of licensed therapeutics targeting the Sudan virus, our study focused on inhibiting VP40 dimers, the structural precursors to critical higher-order oligomers, as a novel antiviral strategy. We have established a crystallographic screening pipeline for the identification of small-molecule fragments capable of binding to VP40. Dimeric VP40 of the Sudan virus was recombinantly expressed in bacteria, purified, crystallized, and soaked in a solution of 96 different preselected fragments. Salicylic acid was identified as a crystallographic hit with clear electron density in the pocket between the N- and the C-termini of the VP40 dimer. The binding interaction is predominantly coordinated by amino acid residues leucine 158 (L158) and arginine 214 (R214), which are key in stabilizing salicylic acid within the binding pocket. While salicylic acid displayed minimal impact on the functional aspects of VP40, we delved deeper into characterizing the druggability of the identified binding pocket. We analyzed the influence of residues L158 and R214 on the formation of virus-like particles and viral RNA synthesis. Site-directed mutagenesis of these residues to alanine markedly affected both VP40's budding activity and its effect on viral RNA synthesis, underscoring the potential of the salicylic acid binding pocket as a drug target. In summary, our findings lay the foundation for structure-guided drug design to provide lead compounds against Sudan virus VP40.

Fragment-based drug discovery campaigns guided by native mass spectrometry

Native mass spectrometry (nMS) is well established as a biophysical technique for characterising biomolecules and their interactions with endogenous or investigational small molecule ligands. The high sensitivity mass measurements make nMS particularly well suited for applications in fragment-based drug discovery (FBDD) screening campaigns where the detection of weakly binding ligands to a target biomolecule is crucial. We first reviewed the contributions of nMS to guiding FBDD hit identification in 2013, providing a comprehensive perspective on the early adoption of nMS for fragment screening. Here we update this initial progress with a focus on contributions of nMS that have guided FBDD for the period 2014 until end of 2023. We highlight the development of nMS adoption in FBDD in the context of other biophysical fragment screening techniques. We also discuss the roadmap for increased adoption of nMS for fragment screening beyond soluble proteins, including for guiding the discovery of fragments supporting advances in PROTAC discovery, RNA-binding small molecules and covalent therapeutic drug discovery.

Novel Fragment Inhibitors of PYCR1 from Docking-Guided X-ray Crystallography

The proline biosynthetic enzyme Δ1-pyrroline-5-carboxylate (P5C) reductase 1 (PYCR1) is one of the most consistently upregulated enzymes across multiple cancer types and central to the metabolic rewiring of cancer cells. Herein, we describe a fragment-based, structure-first approach to the discovery of PYCR1 inhibitors. Thirty-seven fragment-like carboxylic acids in the molecular weight range of 143-289 Da were selected from docking and then screened using X-ray crystallography as the primary assay. Strong electron density was observed for eight compounds, corresponding to a crystallographic hit rate of 22%. The fragments are novel compared to existing proline analog inhibitors in that they block both the P5C substrate pocket and the NAD(P)H binding site. Four hits showed inhibition of PYCR1 in kinetic assays, and one has lower apparent IC50 than the current best proline analog inhibitor. These results show proof-of-concept for our inhibitor discovery approach and provide a basis for fragment-to-lead optimization.

Fragment library screening by X-ray crystallography and binding site analysis on thioredoxin glutathione reductase of Schistosoma mansoni

Schistosomiasis is caused by parasites of the genus Schistosoma, which infect more than 200 million people. Praziquantel (PZQ) has been the main drug for controlling schistosomiasis for over four decades, but despite that it is ineffective against juvenile worms and size and taste issues with its pharmaceutical forms impose challenges for treating school-aged children. It is also important to note that PZQ resistant strains can be generated in laboratory conditions and observed in the field, hence its extensive use in mass drug administration programs raises concerns about resistance, highlighting the need to search for new schistosomicidal drugs. Schistosomes survival relies on the redox enzyme thioredoxin glutathione reductase (TGR), a validated target for the development of new anti-schistosomal drugs. Here we report a high-throughput fragment screening campaign of 768 compounds against S. mansoni TGR (SmTGR) using X-ray crystallography. We observed 49 binding events involving 35 distinct molecular fragments which were found to be distributed across 16 binding sites. Most sites are described for the first time within SmTGR, a noteworthy exception being the "doorstop pocket" near the NADPH binding site. We have compared results from hotspots and pocket druggability analysis of SmTGR with the experimental binding sites found in this work, with our results indicating only limited coincidence between experimental and computational results. Finally, we discuss that binding sites at the doorstop/NADPH binding site and in the SmTGR dimer interface, should be prioritized for developing SmTGR inhibitors as new antischistosomal drugs.

Geometric Deep Learning for Structure-Based Ligand Design

A pervasive challenge in drug design is determining how to expand a ligand-a small molecule that binds to a target biomolecule-in order to improve various properties of the ligand. Adding single chemical groups, known as fragments, is important for lead optimization tasks, and adding multiple fragments is critical for fragment-based drug design. We have developed a comprehensive framework that uses machine learning and three-dimensional protein-ligand structures to address this challenge. Our method, FRAME, iteratively determines where on a ligand to add fragments, selects fragments to add, and predicts the geometry of the added fragments. On a comprehensive benchmark, FRAME consistently improves predicted affinity and selectivity relative to the initial ligand, while generating molecules with more drug-like chemical properties than docking-based methods currently in widespread use. FRAME learns to accurately describe molecular interactions despite being given no prior information on such interactions. The resulting framework for quality molecular hypothesis generation can be easily incorporated into the workflows of medicinal chemists for diverse tasks, including lead optimization, fragment-based drug discovery, and de novo drug design.

Covalent fragment libraries in drug discovery-Design, synthesis, and screening methods

As the development of drugs with a covalent mode of action is becoming increasingly popular, well-validated covalent fragment-based drug discovery (FBDD) methods have been comparatively slow to keep up with the demand. In this chapter the principles of covalent fragment reactivity, library design, synthesis, and screening methods are explored in depth, focussing on literature examples with direct applications to practical covalent fragment library design and screening. Further, questions about the future of the field are explored and potential useful advances are proposed.

A Benchmark Study of Protein-Fragment Complex Structure Calculations with NMR2

Protein-fragment complex structures are particularly sought after in medicinal chemistry to rationally design lead molecules. These structures are usually derived using X-ray crystallography, but the failure rate is non-neglectable. NMR is a possible alternative for the calculation of weakly interacting complexes. Nevertheless, the time-consuming protein signal assignment step remains a barrier to its routine application. NMR Molecular Replacement (NMR2) is a versatile and rapid method that enables the elucidation of a protein-ligand complex structure. It has been successfully applied to peptides, drug-like molecules, and more recently to fragments. Due to the small size of the fragments, ca < 300 Da, solving the structures of the protein-fragment complexes is particularly challenging. Here, we present the expected performances of NMR2 when applied to protein-fragment complexes. The NMR2 approach has been benchmarked with the SERAPhic fragment library to identify the technical challenges in protein-fragment NMR structure calculation. A straightforward strategy is proposed to increase the method's success rate further. The presented work confirms that NMR2 is an alternative method to X-ray crystallography for solving protein-fragment complex structures.

Fast fragment and compound screening pipeline at the Swiss Light Source

Crystallography-based fragment screening is a highly effective technique employed in structure-based drug discovery to expand the range of lead development opportunities. It allows screening and sorting of weakly binding, low molecular mass fragments, which can be developed into larger high-affinity lead compounds. Technical improvements at synchrotron beamlines, design of innovative libraries mapping chemical space efficiently, effective soaking methods and enhanced data analysis have enabled the implementation of high-throughput fragment screening pipelines at multiple synchrotron facilities. This widened access to CBFS beyond the pharma industry has allowed academic users to rapidly screen large quantities of fragment-soaked protein crystals. The positive outcome of a CBFS campaign is a set of structures that present the three-dimensional arrangement of fragment-protein complexes in detail, thereby providing information on the location and the mode of interaction of bound fragments. Through this review, we provide users with a comprehensive guide that sets clear expectations before embarking on a crystallography-based fragment screening campaign. We present a list of essential pre-requirements that must be assessed, including the suitability of your current crystal system for a fragment screening campaign. Furthermore, we extensively discuss the available methodological options, addressing their limitations and providing strategies to overcome them. Additionally, we provide a brief perspective on how to proceed once hits are obtained. Notably, we emphasize the solutions we have implemented for instrumentation and software development within our Fast Fragment and Compound Screening pipeline. We also highlight third-party software options that can be utilized for rapid refinement and hit assessment.

LifeSoaks: a tool for analyzing solvent channels in protein crystals and obstacles for soaking experiments

Due to the structural complexity of proteins, their corresponding crystal arrangements generally contain a significant amount of solvent-occupied space. These areas allow a certain degree of intracrystalline protein flexibility and mobility of solutes. Therefore, knowledge of the geometry of solvent-filled channels and cavities is essential whenever the dynamics inside a crystal are of interest. Especially in soaking experiments for structure-based drug design, ligands must be able to traverse the crystal solvent channels and reach the corresponding binding pockets. Unsuccessful screenings are sometimes attributed to the geometry of the crystal packing, but the underlying causes are often difficult to understand. This work presents LifeSoaks, a novel tool for analyzing and visualizing solvent channels in protein crystals. LifeSoaks uses a Voronoi diagram-based periodic channel representation which can be efficiently computed. The size and location of channel bottlenecks, which might hinder molecular diffusion, can be directly derived from this representation. This work presents the calculated bottleneck radii for all crystal structures in the PDB and the analysis of a new, hand-curated data set of structures obtained by soaking experiments. The results indicate that the consideration of bottleneck radii and the visual inspection of channels are beneficial for planning soaking experiments.

Halo Library, a Tool for Rapid Identification of Ligand Binding Sites on Proteins Using Crystallographic Fragment Screening

X-ray crystallographic fragment screening (XCFS) uses fragment-sized molecules (∼60 to 300 Da) to access binding sites on proteins that may be inaccessible to larger drug-like molecules (>300 Da). Previous studies have shown that fragments containing halogen atoms bind more often to proteins than non-halogenated fragments. Here, we designed the Halo Library containing 46 halogenated fragments (including the "universal fragment" 4-bromopyrazole), a majority of which have been reported to bind to or inhibit one or more targets. The library was screened against the crystals of HIV-1 reverse transcriptase with the drug rilpivirine, yielding an overall hit rate of 26%. Two new binding sites were discovered, and several hot spots were identified. This small library may thus provide a convenient tool for rapidly assessing the feasibility of a target for XCFS, mapping hot spots and cryptic sites, as well as finding fragment binders that can be useful for developing drug leads.

Molecular-evaluated and explainable drug repurposing for COVID-19 using ensemble knowledge graph embedding

The search for an effective drug is still urgent for COVID-19 as no drug with proven clinical efficacy is available. Finding the new purpose of an approved or investigational drug, known as drug repurposing, has become increasingly popular in recent years. We propose here a new drug repurposing approach for COVID-19, based on knowledge graph (KG) embeddings. Our approach learns "ensemble embeddings" of entities and relations in a COVID-19 centric KG, in order to get a better latent representation of the graph elements. Ensemble KG-embeddings are subsequently used in a deep neural network trained for discovering potential drugs for COVID-19. Compared to related works, we retrieve more in-trial drugs among our top-ranked predictions, thus giving greater confidence in our prediction for out-of-trial drugs. For the first time to our knowledge, molecular docking is then used to evaluate the predictions obtained from drug repurposing using KG embedding. We show that Fosinopril is a potential ligand for the SARS-CoV-2 nsp13 target. We also provide explanations of our predictions thanks to rules extracted from the KG and instanciated by KG-derived explanatory paths. Molecular evaluation and explanatory paths bring reliability to our results and constitute new complementary and reusable methods for assessing KG-based drug repurposing.

Exiting the tunnel of uncertainty: crystal soak to validated hit

Crystallographic fragment screens provide an efficient and effective way to identify small-molecule ligands of a crystallized protein. Due to their low molecular weight, such hits tend to have low, often unquantifiable, affinity for their target, complicating the twin challenges of validating the hits as authentic solution-phase ligands of the target and identifying the `best' hit(s) for further elaboration. In this article, approaches that address these challenges are assessed. Using retrospective analysis of a recent ATAD2 hit-identification campaign, alongside other examples of successful fragment-screening campaigns, it is suggested that hit validation and prioritization are best achieved by a `triangulation' approach in which the results of multiple available biochemical and biophysical techniques are correlated to develop qualitative structure-activity relationships (SARs). Such qualitative SARs may indeed be the only means by which to navigate a project through the tunnel of uncertainty that prevails before on-scale biophysical, biochemical and/or biological measurements become possible.

Development of New Potential Inhibitors of β1 Integrins through In Silico Methods-Screening and Computational Validation

Integrins are transmembrane receptors that play a critical role in many biological processes which can be therapeutically modulated using integrin blockers, such as peptidomimetic ligands. This work aimed to develop new potential β1 integrin antagonists using modeled receptors based on the aligned crystallographic structures and docked with three lead compounds (BIO1211, BIO5192, and TCS2314), widely known as α4β1 antagonists. Lead-compound complex optimization was performed by keeping intact the carboxylate moiety of the ligand, adding substituents in two other regions of the molecule to increase the affinity with the target. Additionally, pharmacokinetic predictions were performed for the ten best ligands generated, with the lowest docking interaction energy obtained for α4β1 and BIO5192. Results revealed an essential salt bridge between the BIO5192 carboxylate group and the Mg²⁺ MIDAS ion of the integrin. We then generated more than 200 new BIO5192 derivatives, some with a greater predicted affinity to α4β1. Furthermore, the significance of retaining the pyrrolidine core of the ligand and increasing the therapeutic potential of the new compounds is emphasized. Finally, one novel molecule (1592) was identified as a potential drug candidate, with appropriate pharmacokinetic profiles, similar dynamic behavior at the integrin interaction site compared with BIO5192, and a higher predicted affinity to VLA-4.

SAMPL7 protein-ligand challenge: A community-wide evaluation of computational methods against fragment screening and pose-prediction

A novel crystallographic fragment screening data set was generated and used in the SAMPL7 challenge for protein-ligands. The SAMPL challenges prospectively assess the predictive power of methods involved in computer-aided drug design. Application of various methods to fragment molecules are now widely used in the search for new drugs. However, there is little in the way of systematic validation specifically for fragment-based approaches. We have performed a large crystallographic high-throughput fragment screen against the therapeutically relevant second bromodomain of the Pleckstrin-homology domain interacting protein (PHIP2) that revealed 52 different fragments bound across 4 distinct sites, 47 of which were bound to the pharmacologically relevant acetylated lysine (Kac) binding site. These data were used to assess computational screening, binding pose prediction and follow-up enumeration. All submissions performed randomly for screening. Pose prediction success rates (defined as less than 2 Å root mean squared deviation against heavy atom crystal positions) ranged between 0 and 25% and only a very few follow-up compounds were deemed viable candidates from a medicinal-chemistry perspective based on a common molecular descriptors analysis. The tight deadlines imposed during the challenge led to a small number of submissions suggesting that the accuracy of rapidly responsive workflows remains limited. In addition, the application of these methods to reproduce crystallographic fragment data still appears to be very challenging. The results show that there is room for improvement in the development of computational tools particularly when applied to fragment-based drug design.

Using Structure-guided Fragment-Based Drug Discovery to Target Pseudomonas aeruginosa Infections in Cystic Fibrosis

Cystic fibrosis (CF) is progressive genetic disease that predisposes lungs and other organs to multiple long-lasting microbial infections. Pseudomonas aeruginosa is the most prevalent and deadly pathogen among these microbes. Lung function of CF patients worsens following chronic infections with P. aeruginosa and is associated with increased mortality and morbidity. Emergence of multidrug-resistant, extensively drug-resistant and pandrug-resistant strains of P. aeruginosa due to intrinsic and adaptive antibiotic resistance mechanisms has failed the current anti-pseudomonal antibiotics. Hence new antibacterials are urgently needed to treat P. aeruginosa infections. Structure-guided fragment-based drug discovery (FBDD) is a powerful approach in the field of drug development that has succeeded in delivering six FDA approved drugs over the past 20 years targeting a variety of biological molecules. However, FBDD has not been widely used in the development of anti-pseudomonal molecules. In this review, we first give a brief overview of our structure-guided FBDD pipeline and then give a detailed account of FBDD campaigns to combat P. aeruginosa infections by developing small molecules having either bactericidal or anti-virulence properties. We conclude with a brief overview of the FBDD efforts in our lab at the University of Cambridge towards targeting P. aeruginosa infections.

Fragment Hotspot Mapping to Identify Selectivity-Determining Regions between Related Proteins

Selectivity is a crucial property in small molecule development. Binding site comparisons within a protein family are a key piece of information when aiming to modulate the selectivity profile of a compound. Binding site differences can be exploited to confer selectivity for a specific target, while shared areas can provide insights into polypharmacology. As the quantity of structural data grows, automated methods are needed to process, summarize, and present these data to users. We present a computational method that provides quantitative and data-driven summaries of the available binding site information from an ensemble of structures of the same protein. The resulting ensemble maps identify the key interactions important for ligand binding in the ensemble. The comparison of ensemble maps of related proteins enables the identification of selectivity-determining regions within a protein family. We applied the method to three examples from the well-researched human bromodomain and kinase families, demonstrating that the method is able to identify selectivity-determining regions that have been used to introduce selectivity in past drug discovery campaigns. We then illustrate how the resulting maps can be used to automate comparisons across a target protein family.

Fragment-Based Drug Discovery for RNA Targets

Rapid development within the fields of both fragment-based drug discovery (FBDD) and medicinal targeting of RNA provides possibilities for combining technologies and methods in novel ways. This review provides an overview of fragment-based screening (FBS) against RNA targets, including a discussion of the most recently used screening and hit validation methods such as NMR spectroscopy, X-ray crystallography, and virtual screening methods. A discussion of fragment library design based on research from small-molecule RNA binders provides an overview on both the currently limited guidelines within RNA-targeting fragment library design, and future possibilities. Finally, future perspectives are provided on screening and hit validation methods not yet used in combination with both fragment screening and RNA targets.

Probing the Surface of a Parasite Drug Target Thioredoxin Glutathione Reductase Using Small Molecule Fragments

Fragment screening is a powerful drug discovery approach particularly useful for enzymes difficult to inhibit selectively, such as the thiol/selenol-dependent thioredoxin reductases (TrxRs), which are essential and druggable in several infectious diseases. Several known inhibitors are reactive electrophiles targeting the selenocysteine-containing C-terminus and thus often suffering from off-target reactivity in vivo. The lack of structural information on the interaction modalities of the C-terminus-targeting inhibitors, due to the high mobility of this domain and the lack of alternative druggable sites, prevents the development of selective inhibitors for TrxRs. In this work, fragments selected from actives identified in a large screen carried out against Thioredoxin Glutathione Reductase from Schistosoma mansoni (SmTGR) were probed by X-ray crystallography. SmTGR is one of the most promising drug targets for schistosomiasis, a devastating, neglected disease. Utilizing a multicrystal method to analyze electron density maps, structural analysis, and functional studies, three binding sites were characterized in SmTGR: two sites are close to or partially superposable with the NADPH binding site, while the third one is found between two symmetry related SmTGR subunits of the crystal lattice. Surprisingly, one compound bound to this latter site stabilizes, through allosteric effects mediated by the so-called guiding bar residues, the crucial redox active C-terminus of SmTGR, making it finally visible at high resolution. These results further promote fragments as small molecule probes for investigating functional aspects of the target protein, exemplified by the allosteric effect on the C-terminus, and providing fundamental chemical information exploitable in drug discovery.

Crystal structures of Val58Ile tryptophan repressor in a domain-swapped array in the presence and absence of L-tryptophan

The crystal structures of domain-swapped tryptophan repressor (TrpR) variant Val58Ile before and after soaking with the physiological ligand L-tryptophan (L-Trp) indicate that L-Trp occupies the same location in the domain-swapped form as in native dimeric TrpR and makes equivalent residue contacts. This result is unexpected because the ligand binding-site residues arise from three separate polypeptide chains in the domain-swapped form. This work represents the first published structure of a domain-swapped form of TrpR with L-Trp bound. The presented structures also show that the protein amino-terminus, whether or not it bears a disordered extension of about 20 residues, is accessible in the large solvent channels of the domain-swapped crystal form, as in the structures reported previously in this form for TrpR without N-terminal extensions. These findings inspire the exploration of L-Trp analogs and N-terminal modifications as labels to orient guest proteins that cannot otherwise be crystallized in the solvent channels of crystalline domain-swapped TrpR hosts for potential diffraction analysis.

Comprehensive Analysis of Binding Sites in Tubulin

Tubulin plays essential roles in vital cellular activities and is the target of a wide range of proteins and ligands. Here, using a combined computational and crystallographic fragment screening approach, we addressed the question of how many binding sites exist in tubulin. We identified 27 distinct sites, of which 11 have not been described previously, and analyzed their relationship to known tubulin-protein and tubulin-ligand interactions. We further observed an intricate pocket communication network and identified 56 chemically diverse fragments that bound to 10 distinct tubulin sites. Our results offer a unique structural basis for the development of novel small molecules for use as tubulin modulators in basic research applications or as drugs. Furthermore, our method lays down a framework that may help to discover new pockets in other pharmaceutically important targets and characterize them in terms of chemical tractability and allosteric modulation.

Exploring protein hotspots by optimized fragment pharmacophores

Fragment-based drug design has introduced a bottom-up process for drug development, with improved sampling of chemical space and increased effectiveness in early drug discovery. Here, we combine the use of pharmacophores, the most general concept of representing drug-target interactions with the theory of protein hotspots, to develop a design protocol for fragment libraries. The SpotXplorer approach compiles small fragment libraries that maximize the coverage of experimentally confirmed binding pharmacophores at the most preferred hotspots. The efficiency of this approach is demonstrated with a pilot library of 96 fragment-sized compounds (SpotXplorer0) that is validated on popular target classes and emerging drug targets. Biochemical screening against a set of GPCRs and proteases retrieves compounds containing an average of 70% of known pharmacophores for these targets. More importantly, SpotXplorer0 screening identifies confirmed hits against recently established challenging targets such as the histone methyltransferase SETD2, the main protease (3CLPro) and the NSP3 macrodomain of SARS-CoV-2.

Fragment Screening Reveals Starting Points for Rational Design of Galactokinase 1 Inhibitors to Treat Classic Galactosemia

Classic galactosemia is caused by loss-of-function mutations in galactose-1-phosphate uridylyltransferase (GALT) that lead to toxic accumulation of its substrate, galactose-1-phosphate. One proposed therapy is to inhibit the biosynthesis of galactose-1-phosphate, catalyzed by galactokinase 1 (GALK1). Existing inhibitors of human GALK1 (hGALK1) are primarily ATP-competitive with limited clinical utility to date. Here, we determined crystal structures of hGALK1 bound with reported ATP-competitive inhibitors of the spiro-benzoxazole series, to reveal their binding mode in the active site. Spurred by the need for additional chemotypes of hGALK1 inhibitors, desirably targeting a nonorthosteric site, we also performed crystallography-based screening by soaking hundreds of hGALK1 crystals, already containing active site ligands, with fragments from a custom library. Two fragments were found to bind close to the ATP binding site, and a further eight were found in a hotspot distal from the active site, highlighting the strength of this method in identifying previously uncharacterized allosteric sites. To generate inhibitors of improved potency and selectivity targeting the newly identified binding hotspot, new compounds were designed by merging overlapping fragments. This yielded two micromolar inhibitors of hGALK1 that were not competitive with respect to either substrate (ATP or galactose) and demonstrated good selectivity over hGALK1 homologues, galactokinase 2 and mevalonate kinase. Our findings are therefore the first to demonstrate inhibition of hGALK1 from an allosteric site, with potential for further development of potent and selective inhibitors to provide novel therapeutics for classic galactosemia.

Avoiding Drug Resistance in HIV Reverse Transcriptase

HIV reverse transcriptase (RT) is an enzyme that plays a major role in the replication cycle of HIV and has been a key target of anti-HIV drug development efforts. Because of the high genetic diversity of the virus, mutations in RT can impart resistance to various RT inhibitors. As the prevalence of drug resistance mutations is on the rise, it is necessary to design strategies that will lead to drugs less susceptible to resistance. Here we provide an in-depth review of HIV reverse transcriptase, current RT inhibitors, novel RT inhibitors, and mechanisms of drug resistance. We also present novel strategies that can be useful to overcome RT's ability to escape therapies through drug resistance. While resistance may not be completely avoidable, designing drugs based on the strategies and principles discussed in this review could decrease the prevalence of drug resistance.

Use of Fluorescent Nucleotides to Map RNA-Binding Sites on Protein Surface

Currently, studies of RNA/protein interactions occupy a prominent place in molecular biology and medicine. The structures of RNA-protein complexes may be determined by X-ray crystallography or NMR for further analyses. These methods are time-consuming and difficult due to the versatility and dynamics of the RNA structure. Furthermore, due to the need to solve the "phase problem" for each dataset in crystallography, crystallographic structures of RNA are still underrepresented. Structure determination of single ribonucleotide-protein complexes is a useful tool to identify the position of single-stranded RNA-binding sites in proteins. We describe here a structural approach that incorporates affinity measurement of a protein for various single ribonucleotides, ranking the RNA/protein complexes according to their stability. This chapter describes how to perform these measurements, including a perspective for the analysis of RNA-binding sites in protein and single-nucleotide crystal soaking.

Fragment-based drug discovery: opportunities for organic synthesis

This Review describes the increasing demand for organic synthesis to facilitate fragment-based drug discovery (FBDD), focusing on polar, unprotected fragments. In FBDD, X-ray crystal structures are used to design target molecules for synthesis with new groups added onto a fragment via specific growth vectors. This requires challenging synthesis which slows down drug discovery, and some fragments are not progressed into optimisation due to synthetic intractability. We have evaluated the output from Astex's fragment screenings for a number of programs, including urokinase-type plasminogen activator, hematopoietic prostaglandin D2 synthase, and hepatitis C virus NS3 protease-helicase, and identified fragments that were not elaborated due, in part, to a lack of commercially available analogues and/or suitable synthetic methodology. This represents an opportunity for the development of new synthetic research to enable rapid access to novel chemical space and fragment optimisation.

Structure driven compound optimization in targeted protein degradation

Small molecule induced protein degradation has created tremendous excitement in drug discovery within recent years. Not being confined to target inhibition and being able to remove disease-causing protein targets via engagement and subsequent ubiquitination has provided scientists with a powerful tool to expand the druggable space. At the center of this approach sits the ternary complex formed between an E3 ubiquitin ligase, the small molecule degrader, and the target protein. A productive ternary complex is pivotal for a ubiquitin to be transferred to a surface lysine of the target protein resulting in poly-ubiquitination which enables recognition and finally degradation by the proteasome. As understanding the ternary complex means understanding the degradation process, many efforts are put into obtaining structural information of the ternary complex and getting a snapshot of the underlying conformations and molecular contacts. Locking this transient trimeric intermediate in a crystalline state has proven to be very demanding but the obtained results have tremendously improved our understanding of small molecule degraders. This review discusses target protein degradation from a structural perspective and highlights the evolution of certain degraders based on the obtained structural insights.

BANΔIT: B'-Factor Analysis for Drug Design and Structural Biology

The analysis of B‐factor profiles from X‐ray protein structures can be utilized for structure‐based drug design since protein mobility changes have been associated with the quality of protein‐ligand interactions. With the BANΔIT (B’‐factor analysis and ΔB’ interpretation toolkit), we have developed a JavaScript‐based browser application that provides a graphical user interface for the normalization and analysis of B’‐factor profiles. To emphasize the usability for rational drug design applications, we have analyzed a selection of crystallographic protein‐ligand complexes and have given exemplary conclusions for further drug optimization including the development of a B’‐factor‐supported pharmacophore model for SARS CoV‐2 main protease inhibitors. BANΔIT is available online at https://bandit.uni‐mainz.de. The source code can be downloaded from https://github.com/FBarthels/BANDIT.

Demonstration of the utility of DOS-derived fragment libraries for rapid hit derivatisation in a multidirectional fashion

Organic synthesis underpins the evolution of weak fragment hits into potent lead compounds. Deficiencies within current screening collections often result in the requirement of significant synthetic investment to enable multidirectional fragment growth, limiting the efficiency of the hit evolution process. Diversity-oriented synthesis (DOS)-derived fragment libraries are constructed in an efficient and modular fashion and thus are well-suited to address this challenge. To demonstrate the effective nature of such libraries within fragment-based drug discovery, we herein describe the screening of a 40-member DOS library against three functionally distinct biological targets using X-Ray crystallography. Firstly, we demonstrate the importance for diversity in aiding hit identification with four fragment binders resulting from these efforts. Moreover, we also exemplify the ability to readily access a library of analogues from cheap commercially available materials, which ultimately enabled the exploration of a minimum of four synthetic vectors from each molecule. In total, 10-14 analogues of each hit were rapidly accessed in three to six synthetic steps. Thus, we showcase how DOS-derived fragment libraries enable efficient hit derivatisation and can be utilised to remove the synthetic limitations encountered in early stage fragment-based drug discovery.

Protein-fragment complex structures derived by NMR molecular replacement

Recently we have established an NMR molecular replacement method, which is capable of solving the structure of the interaction site of protein-ligand complexes in a fully automated manner. While the method was successfully applied for ligands with strong and weak binding affinities, including small molecules and peptides, its applicability on ligand fragments remains to be shown. Structures of fragment-protein complexes are more challenging for the method since fragments contain only few protons. Here we show a successful application of the NMR molecular replacement method in solving structures of complexes between three derivatives of a ligand fragment and the protein receptor PIN1. We anticipate that this approach will find a broad application in fragment-based lead discovery.

Structural characterization of melatonin as an inhibitor of the Wnt deacylase Notum

The hormone melatonin, secreted from the pineal gland, mediates multiple physiological effects including modulation of Wnt/β-catenin signalling. The Wnt palmitoleate lipid modification is essential for its signalling activity, while the carboxylesterase Notum can remove the lipid from Wnt and inactivate it. Notum enzyme inhibition can therefore upregulate Wnt signalling. While searching for Notum inhibitors by crystallographic fragment screening, a hit compound N-[2-(5-fluoro-1H-indol-3-yl)ethyl]acetamide that is structurally similar to melatonin came to our attention. We then soaked melatonin and its precursor N-acetylserotonin into Notum crystals and obtained high-resolution structures (≤1.5 Å) of their complexes. In each of the structures, two compound molecules bind with Notum: one at the enzyme's catalytic pocket, overlapping the space occupied by the acyl tail of the Wnt palmitoleate lipid, and the other at the edge of the pocket opposite the substrate entrance. Although the inhibitory activity of melatonin shown by in vitro enzyme assays is low (IC50 75 µmol/L), the structural information reported here provides a basis for the design of potent and brain accessible drugs for neurodegenerative diseases such as Alzheimer's disease, in which upregulation of Wnt signalling may be beneficial.

In silico Strategies to Support Fragment-to-Lead Optimization in Drug Discovery

Fragment-based drug (or lead) discovery (FBDD or FBLD) has developed in the last two decades to become a successful key technology in the pharmaceutical industry for early stage drug discovery and development. The FBDD strategy consists of screening low molecular weight compounds against macromolecular targets (usually proteins) of clinical relevance. These small molecular fragments can bind at one or more sites on the target and act as starting points for the development of lead compounds. In developing the fragments attractive features that can translate into compounds with favorable physical, pharmacokinetics and toxicity (ADMET-absorption, distribution, metabolism, excretion, and toxicity) properties can be integrated. Structure-enabled fragment screening campaigns use a combination of screening by a range of biophysical techniques, such as differential scanning fluorimetry, surface plasmon resonance, and thermophoresis, followed by structural characterization of fragment binding using NMR or X-ray crystallography. Structural characterization is also used in subsequent analysis for growing fragments of selected screening hits. The latest iteration of the FBDD workflow employs a high-throughput methodology of massively parallel screening by X-ray crystallography of individually soaked fragments. In this review we will outline the FBDD strategies and explore a variety of in silico approaches to support the follow-up fragment-to-lead optimization of either: growing, linking, and merging. These fragment expansion strategies include hot spot analysis, druggability prediction, SAR (structure-activity relationships) by catalog methods, application of machine learning/deep learning models for virtual screening and several de novo design methods for proposing synthesizable new compounds. Finally, we will highlight recent case studies in fragment-based drug discovery where in silico methods have successfully contributed to the development of lead compounds.

Theory and applications of differential scanning fluorimetry in early-stage drug discovery

Differential scanning fluorimetry (DSF) is an accessible, rapid, and economical biophysical technique that has seen many applications over the years, ranging from protein folding state detection to the identification of ligands that bind to the target protein. In this review, we discuss the theory, applications, and limitations of DSF, including the latest applications of DSF by ourselves and other researchers. We show that DSF is a powerful high-throughput tool in early drug discovery efforts. We place DSF in the context of other biophysical methods frequently used in drug discovery and highlight their benefits and downsides. We illustrate the uses of DSF in protein buffer optimization for stability, refolding, and crystallization purposes and provide several examples of each. We also show the use of DSF in a more downstream application, where it is used as an in vivo validation tool of ligand-target interaction in cell assays. Although DSF is a potent tool in buffer optimization and large chemical library screens when it comes to ligand-binding validation and optimization, orthogonal techniques are recommended as DSF is prone to false positives and negatives.

Identification of a ligand binding hot spot and structural motifs replicating aspects of tyrosyl-DNA phosphodiesterase I (TDP1) phosphoryl recognition by crystallographic fragment cocktail screening

Tyrosyl DNA-phosphodiesterase I (TDP1) repairs type IB topoisomerase (TOP1) cleavage complexes generated by TOP1 inhibitors commonly used as anticancer agents. TDP1 also removes DNA 3' end blocking lesions generated by chain-terminating nucleosides and alkylating agents, and base oxidation both in the nuclear and mitochondrial genomes. Combination therapy with TDP1 inhibitors is proposed to synergize with topoisomerase targeting drugs to enhance selectivity against cancer cells exhibiting deficiencies in parallel DNA repair pathways. A crystallographic fragment screening campaign against the catalytic domain of TDP1 was conducted to identify new lead compounds. Crystal structures revealed two fragments that bind to the TDP1 active site and exhibit inhibitory activity against TDP1. These fragments occupy a similar position in the TDP1 active site as seen in prior crystal structures of TDP1 with bound vanadate, a transition state mimic. Using structural insights into fragment binding, several fragment derivatives have been prepared and evaluated in biochemical assays. These results demonstrate that fragment-based methods can be a highly feasible approach toward the discovery of small-molecule chemical scaffolds to target TDP1, and for the first time, we provide co-crystal structures of small molecule inhibitors bound to TDP1, which could serve for the rational development of medicinal TDP1 inhibitors.

Target 2035: probing the human proteome

Biomedical scientists tend to focus on only a small fraction of the proteins encoded by the human genome despite overwhelming genetic evidence that many understudied proteins are important for human disease. One of the best ways to interrogate the function of a protein and to determine its relevance as a drug target is by using a pharmacological modulator, such as a chemical probe or an antibody. If these tools were available for most human proteins, it should be possible to translate the tremendous advances in genomics into a greater understanding of human health and disease, and catalyze the creation of innovative new medicines. Target 2035 is a global federation for developing and applying new technologies with the goal of creating chemogenomic libraries, chemical probes, and/or functional antibodies for the entire proteome.

Relative Binding Energies Predict Crystallographic Binding Modes of Ethionamide Booster Lead Compounds

Transcriptional repressor EthR from Mycobacterium tuberculosis is a valuable target for antibiotic booster drugs. We previously reported a virtual screening campaign to identify EthR inhibitors for development. Two ligand binding orientations were often proposed, though only the top scoring pose was utilized for filtering of the large data set. We obtained biophysically validated hits, some of which yielded complex crystal structures. In some cases, the crystallized binding mode and top scoring mode agree, while for others an alternate ligand binding orientation was found. In this contribution, we combine rigid docking, molecular dynamics simulations, and the linear interaction energy method to calculate binding free energies and derive relative binding energies for a number of EthR inhibitors in both modes. This strategy allowed us to correctly predict the most favorable orientation. Therefore, this widely applicable approach will be suitable to triage multiple binding modes within EthR and other potential drug targets with similar characteristics.