Library Design Strategies To Accelerate Fragment-Based Drug Discovery

Fragment-based drug discovery (FBDD) has become an established approach for the generation of early lead candidates. However, despite its success and inherent advantages, hit-to-candidate progression for FBDD is not necessarily faster than that of traditional high-throughput screening. Thus, new technology-driven library design strategies have emerged as a means to facilitate more efficient fragment screening and/or subsequent fragment-to-hit chemistry. This minireview discusses such strategies, which cover the use of labeled fragments for NMR spectroscopy, X-ray crystallographic screening of specialized fragments, covalent linkage for mass spectrometry, dynamic combinatorial chemistry, and fragments optimized for easy elaboration.     

The 3F Library: Fluorinated Fsp3‐Rich Fragments for Expeditious 19F NMR Based Screening

Fragment‐based drug discovery (FBDD) is a popular method in academia and the pharmaceutical industry for the discovery of early lead candidates. Despite its wide‐spread use, the approach still suffers from laborious screening workflows and a limited diversity in the fragments applied. Presented here is the design, synthesis, and biological evaluation of the first fragment library specifically tailored to tackle both these challenges. The 3F library of 115 fluorinated, Fsp³‐rich fragments is shape diverse and natural‐product‐like with desirable physicochemical properties. The library is perfectly suited for rapid and efficient screening by NMR spectroscopy in a two‐stage workflow of ¹⁹F NMR and subsequent ¹H NMR methods. Hits against four diverse protein targets are widely distributed among the fragment scaffolds in the 3F library and a 67 % validation rate was achieved using secondary assays. This collection is the first synthetic fragment library tailor‐made for ¹⁹F NMR screening and the results demonstrate that the approach should find broad application in the FBDD community.     

Rationalizing fragment based drug discovery for BACE1: insights from FB-QSAR,FB-QSSR,multi objective (MO-QSPR) and MIF studies

The ability to identify fragments that interact with a biological target is a key step in FBDD. To date, the concept of fragment based drug design (FBDD) is increasingly driven by bio-physical methods. To expand the boundaries of QSAR paradigm, and to rationalize FBDD using In silico approach, we propose a fragment based QSAR methodology referred here in as FB-QSAR. The FB-QSAR methodology was validated on a dataset consisting of 52 Hydroxy ethylamine (HEA) inhibitors, disclosed by GlaxoSmithKline Pharmaceuticals as potential anti-Alzheimer agents. To address the issue of target selectivity, a major confounding factor in the development of selective BACE1 inhibitors, FB-QSSR models were developed using the reported off target activity values. A heat map constructed, based on the activity and selectivity profile of the individual R-group fragments, and was in turn used to identify superior R-group fragments. Further, simultaneous optimization of multiple properties, an issue encountered in real-world drug discovery scenario, and often overlooked in QSAR approaches, was addressed using a Multi Objective (MO-QSPR) method that balances properties, based on the defined objectives. MO-QSPR was implemented using Derringer and Suich desirability algorithm to identify the optimal level of independent variables (X) that could confer a trade-off between selectivity and activity. The results obtained from FB-QSAR were further substantiated using MIF (Molecular Interaction Fields) studies. To exemplify the potentials of FB-QSAR and MO-QSPR in a pragmatic fashion, the insights gleaned from the MO-QSPR study was reverse engineered using Inverse-QSAR in a combinatorial fashion to enumerate some prospective novel, potent and selective BACE1 inhibitors.     

Virtual Screening with Gnina 1.0

Virtual screening—predicting which compounds within a specified compound library bind to a target molecule, typically a protein—is a fundamental task in the field of drug discovery. Doing virtual screening well provides tangible practical benefits, including reduced drug development costs, faster time to therapeutic viability, and fewer unforeseen side effects. As with most applied computational tasks, the algorithms currently used to perform virtual screening feature inherent tradeoffs between speed and accuracy. Furthermore, even theoretically rigorous, computationally intensive methods may fail to account for important effects relevant to whether a given compound will ultimately be usable as a drug. Here we investigate the virtual screening performance of the recently released Gnina molecular docking software, which uses deep convolutional networks to score protein-ligand structures. We find, on average, that Gnina outperforms conventional empirical scoring. The default scoring in Gnina outperforms the empirical AutoDock Vina scoring function on 89 of the 117 targets of the DUD-E and LIT-PCBA virtual screening benchmarks with a median 1% early enrichment factor that is more than twice that of Vina. However, we also find that issues of bias linger in these sets, even when not used directly to train models, and this bias obfuscates to what extent machine learning models are achieving their performance through a sophisticated interpretation of molecular interactions versus fitting to non-informative simplistic property distributions.     

The multiple roles of computational chemistry in fragment-based drug design

Fragment-based drug discovery (FBDD) represents a change in strategy from the screening of molecules with higher molecular weights and physical properties more akin to fully drug-like compounds, to the screening of smaller, less complex molecules. This is because it has been recognised that fragment hit molecules can be efficiently grown and optimised into leads, particularly after the binding mode to the target protein has been first determined by 3D structural elucidation, e.g. by NMR or X-ray crystallography. Several studies have shown that medicinal chemistry optimisation of an already drug-like hit or lead compound can result in a final compound with too high molecular weight and lipophilicity. The evolution of a lower molecular weight fragment hit therefore represents an attractive alternative approach to optimisation as it allows better control of compound properties. Computational chemistry can play an important role both prior to a fragment screen, in producing a target focussed fragment library, and post-screening in the evolution of a drug-like molecule from a fragment hit, both with and without the available fragment-target co-complex structure. We will review many of the current developments in the area and illustrate with some recent examples from successful FBDD discovery projects that we have conducted.     

Novel strategy for three-dimensional fragment-based lead discovery

Fragment-based drug design (FBDD) is considered a promising approach in lead discovery. However, for a practical application of this approach, problems remain to be solved. Hence, a novel practical strategy for three-dimensional lead discovery is presented in this work. Diverse fragments with spatial positions and orientations retained in separately adjacent regions were generated by deconstructing well-aligned known inhibitors in the same target active site. These three-dimensional fragments retained their original binding modes in the process of new molecule construction by fragment linking and merging. Root-mean-square deviation (rmsd) values were used to evaluate the conformational changes of the component fragments in the final compounds and to identify the potential leads as the main criteria. Furthermore, the successful validation of our strategy is presented on the basis of two relevant tumor targets (CDK2 and c-Met), demonstrating the potential of our strategy to facilitate lead discovery against some drug targets.     

Identification,synthesis and evaluation of CSF1R inhibitors using fragment based drug design

Colony-stimulating factor 1 receptor is a type III receptor protein tyrosine kinase belonging to PDGFR family. CSF1R signaling is essential for differentiation, proliferation and survival of macrophages. Aberrant expression of CSF1R appears to be an attractive target in several cancer types. Higher expression of CSF1R ligands correlates to tumor progression. CSF1R inhibitors have been shown to suppress cancers. We have attempted an in silico fragment derived drug discovery approach by screening ˜25,000 in-house compounds as potential CSF1R inhibitors. Using FBDD approach we have identified six diverse fragments that exhibit affinity towards hinge region of CSF1R. Some of the fragments 5-nitroindole and 7-azaindole and their derivatives were synthesized for further evaluation. The in silico and in vitro enzyme activity studies reveal moderate inhibition of CSF1R kinase activity by 5-nitroindole and good inhibition by 7-azaindole fragments. Bio and chemiinformatics studies have shown that 7-azaindole compounds have better membrane permeability and enzyme inhibition properties. Molecular docking studies show that the amino acid residues 664–666 in the hinge region of the cytosolic domain of CSF1R to be the preferred region of binding for nitroindole and azaindole derivatives. Further optimization and biological analysis would identify these fragments as potential and promising leads as CSF1R inhibitors.     

Metabolic stability studies of lead compounds supported by separation techniques and chemometrics analysis

With metabolism being one of the main routes of drug elimination from the body (accounting for removal of around 75% of known drugs), it is crucial to understand and study metabolic stability of drug candidates. Metabolically unstable compounds are uncomfortable to administer (requiring repetitive dosage during therapy), while overly stable drugs increase risk of adverse drug reactions. Additionally, biotransformation reactions can lead to formation of toxic or pharmacologically active metabolites (either less‐active than parent drug, or even with different action). There were numerous approaches in estimating metabolic stability, including in vitro, in vivo, in silico, and high‐throughput screening to name a few. This review aims at describing separation techniques used in in vitro metabolic stability estimation, as well as chemometric techniques allowing for creation of predictive models which enable high‐throughput screening approach for estimation of metabolic stability. With a very low rate of drug approval, it is important to understand in silico methods that aim at supporting classical in vitro approach. Predictive models that allow assessment of certain biological properties of drug candidates allow for cutting not only cost, but also time required to synthesize compounds predicted to be unstable or inactive by in silico models.     

Identification of New Potential Inhibitors of Quorum Sensing through a Specialized Multi-Level Computational Approach

Biofilms are aggregates of microorganisms anchored to a surface and embedded in a self-produced matrix of extracellular polymeric substances and have been associated with 80% of all bacterial infections in humans. Because bacteria in biofilms are less amenable to antibiotic treatment, biofilms have been associated with developing antibiotic resistance, a problem that urges developing new therapeutic options and approaches. Interfering with quorum-sensing (QS), an important process of cell-to-cell communication by bacteria in biofilms is a promising strategy to inhibit biofilm formation and development. Here we describe and apply an in silico computational protocol for identifying novel potential inhibitors of quorum-sensing, using CviR—the quorum-sensing receptor from Chromobacterium violaceum—as a model target. This in silico approach combines protein-ligand docking (with 7 different docking programs/scoring functions), receptor-based virtual screening, molecular dynamic simulations, and free energy calculations. Particular emphasis was dedicated to optimizing the discrimination ability between active/inactive molecules in virtual screening tests using a target-specific training set. Overall, the optimized protocol was used to evaluate 66,461 molecules, including those on the ZINC/FDA-Approved database and to the Mu.Ta.Lig Virtual Chemotheca. Multiple promising compounds were identified, yielding good prospects for future experimental validation and for drug repurposing towards QS inhibition.     

Drug Design—Past,Present, Future

Drug design is a complex pharmaceutical science with a long history. Many achievements have been made in the field of drug design since the end of 19th century, when Emil Fisher suggested that the drug–receptor interaction resembles the key and lock interplay. Gradually, drug design has been transformed into a coherent and well-organized science with a solid theoretical background and practical applications. Now, drug design is the most advanced approach for drug discovery. It utilizes the innovations in science and technology and includes them in its wide-ranging arsenal of methods and tools in order to achieve the main goal: discovery of effective, specific, non-toxic, safe and well-tolerated drugs. Drug design is one of the most intensively developing modern sciences and its progress is accelerated by the implication of artificial intelligence. The present review aims to capture some of the most important milestones in the development of drug design, to outline some of the most used current methods and to sketch the future perspective according to the author’s point of view. Without pretending to cover fully the wide range of drug design topics, the review introduces the reader to the content of Molecules’ Special Issue “Drug Design—Science and Practice”.     

To Affinity and Beyond: A Personal Reflection on the Design and Discovery of Drugs

Faced with new and as yet unmet medical need, the stark underperformance of the pharmaceutical discovery process is well described if not perfectly understood. Driven primarily by profit rather than societal need, the search for new pharmaceutical products—small molecule drugs, biologicals, and vaccines—is neither properly funded nor sufficiently systematic. Many innovative approaches remain significantly underused and severely underappreciated, while dominant methodologies are replete with problems and limitations. Design is a component of drug discovery that is much discussed but seldom realised. In and of itself, technical innovation alone is unlikely to fulfil all the possibilities of drug discovery if the necessary underlying infrastructure remains unaltered. A fundamental revision in attitudes, with greater reliance on design powered by computational approaches, as well as a move away from the commercial imperative, is thus essential to capitalise fully on the potential of pharmaceutical intervention in healthcare.     

Open-Source Browser-Based Tools for Structure-Based Computer-Aided Drug Discovery

We here outline the importance of open-source, accessible tools for computer-aided drug discovery (CADD). We begin with a discussion of drug discovery in general to provide context for a subsequent discussion of structure-based CADD applied to small-molecule ligand discovery. Next, we identify usability challenges common to many open-source CADD tools. To address these challenges, we propose a browser-based approach to CADD tool deployment in which CADD calculations run in modern web browsers on users’ local computers. The browser app approach eliminates the need for user-initiated download and installation, ensures broad operating system compatibility, enables easy updates, and provides a user-friendly graphical user interface. Unlike server apps—which run calculations “in the cloud” rather than on users’ local computers—browser apps do not require users to upload proprietary information to a third-party (remote) server. They also eliminate the need for the difficult-to-maintain computer infrastructure required to run user-initiated calculations remotely. We conclude by describing some CADD browser apps developed in our lab, which illustrate the utility of this approach. Aside from introducing readers to these specific tools, we are hopeful that this review highlights the need for additional browser-compatible, user-friendly CADD software.     

Improving the Cellular Selectivity of a Membrane-Disrupting Antimicrobial Agent by Monomer Control and by Taming

Antimicrobial resistance represents a significant world-wide health threat that is looming. To meet this challenge, new classes of antimicrobial agents and the redesign of existing ones will be required. This review summarizes some of the studies that have been carried out in my own laboratories involving membrane-disrupting agents. A major discovery that we made, using a Triton X-100 as a prototypical membrane-disrupting molecule and cholesterol-rich liposomes as model systems, was that membrane disruption can occur by two distinct processes, depending on the state of aggregation of the attacking agent. Specifically, we found that monomers induced leakage, while attack by aggregates resulted in a catastrophic rupture of the membrane. This discovery led us to design of a series of derivatives of the clinically important antifungal agent, Amphotericin B, where we demonstrated the feasibility of separating antifungal from hemolytic activity by decreasing the molecule’s tendency to aggregate, i.e., by controlling its monomer concentration. Using an entirely different approach (i.e., a “taming” strategy), we found that by covalently attaching one or more facial amphiphiles (“floats”) to Amphotericin B, its aggregate forms were much less active in lysing red blood cells while maintaining high antifungal activity. The possibility of applying such “monomer control” and “taming” strategies to other membrane-disrupting antimicrobial agents is briefly discussed.     

Expanding medicinal chemistry into 3D space: metallofragments as 3D scaffolds for fragment-based drug discovery

Fragment-based drug discovery (FBDD) is a powerful strategy for the identification of new bioactive molecules. FBDD relies on fragment libraries, generally of modest size, but of high chemical diversity. Although good chemical diversity in FBDD libraries has been achieved in many respects, achieving shape diversity – particularly fragments with three-dimensional (3D) structures – has remained challenging. A recent analysis revealed that >75% of all conventional, organic fragments are predominantly 1D or 2D in shape. However, 3D fragments are desired because molecular shape is one of the most important factors in molecular recognition by a biomolecule. To address this challenge, the use of inert metal complexes, so-called ‘metallofragments’ (mFs), to construct a 3D fragment library is introduced. A modest library of 71 compounds has been prepared with rich shape diversity as gauged by normalized principle moment of inertia (PMI) analysis. PMI analysis shows that these metallofragments occupy an area of fragment space that is unique and highly underrepresented when compared to conventional organic fragment libraries that are comprised of orders of magnitude more molecules. The potential value of this metallofragment library is demonstrated by screening against several different types of proteins, including an antiviral, an antibacterial, and an anticancer target. The suitability of the metallofragments for future hit-to-lead development was validated through the determination of IC₅₀ and thermal shift values for select fragments against several proteins. These findings demonstrate the utility of metallofragment libraries as a means of accessing underutilized 3D fragment space for FBDD against a variety of protein targets.

Fragment-based drug discovery (FBDD) using 3-dimensional metallofragments is a new strategy for the identification of bioactive molecules.     

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.     

Pharmacophore Mapping Combined with dbCICA Reveal New Structural Features for the Development of Novel Ligands Targeting α4β2 and α7 Nicotinic Acetylcholine Receptors

The neuronal nicotinic acetylcholine receptors (nAChRs) belong to the ligand-gated ion channel (GLIC) group, presenting a crucial role in several biological processes and neuronal disorders. The α4β2 and α7 nAChRs are the most abundant in the central nervous system (CNS), being involved in challenging diseases such as epilepsy, Alzheimer’s disease, schizophrenia, and anxiety disorder, as well as alcohol and nicotine dependencies. In addition, in silico-based strategies may contribute to revealing new insights into drug design and virtual screening to find new drug candidates to treat CNS disorders. In this context, the pharmacophore maps were constructed and validated for the orthosteric sites of α4β2 and α7 nAChRs, through a docking-based Comparative Intermolecular Contacts Analysis (dbCICA). In this sense, bioactive ligands were retrieved from the literature for each receptor. A molecular docking protocol was developed for all ligands in both receptors by using GOLD software, considering GoldScore, ChemScore, ASP, and ChemPLP scoring functions. Output GOLD results were post-processed through dbCICA to identify critical contacts involved in protein-ligand interactions. Moreover, Crossminer software was used to construct a pharmacophoric map based on the most well-behaved ligands and negative contacts from the dbCICA model for each receptor. Both pharmacophore maps were validated by using a ROC curve. The results revealed important features for the ligands, such as the presence of hydrophobic regions, a planar ring, and hydrogen bond donor and acceptor atoms for α4β2. Parallelly, a non-planar ring region was identified for α7. These results can enable fragment-based drug design (FBDD) strategies, such as fragment growing, linking, and merging, allowing an increase in the activity of known fragments. Thus, our results can contribute to a further understanding of structural subunits presenting the potential for key ligand-receptor interactions, favoring the search in molecular databases and the design of novel ligands.     

The 3F Library: Fluorinated Fsp3-Rich Fragments for Expeditious 19F NMR Based Screening

Fragment-based drug discovery (FBDD) is a popular method in academia and the pharmaceutical industry for the discovery of early lead candidates. Despite its wide-spread use, the approach still suffers from laborious screening workflows and a limited diversity in the fragments applied. Presented here is the design, synthesis, and biological evaluation of the first fragment library specifically tailored to tackle both these challenges. The 3F library of 115 fluorinated, Fsp³-rich fragments is shape diverse and natural-product-like with desirable physicochemical properties. The library is perfectly suited for rapid and efficient screening by NMR spectroscopy in a two-stage workflow of ¹⁹F NMR and subsequent ¹H NMR methods. Hits against four diverse protein targets are widely distributed among the fragment scaffolds in the 3F library and a 67 % validation rate was achieved using secondary assays. This collection is the first synthetic fragment library tailor-made for ¹⁹F NMR screening and the results demonstrate that the approach should find broad application in the FBDD community.     

Effective Use of Empirical Data for Virtual Screening against APJR GPCR Receptor

Alzheimer’s disease is a neurodegenerative disorder incompatible with normal daily activity, affecting one in nine people. One of its potential targets is the apelin receptor (APJR), a G-protein coupled receptor, which presents considerably high expression levels in the central nervous system. In silico studies of APJR drug-like molecule binding are in small numbers while high throughput screenings (HTS) are already sufficiently many to devise efficient drug design strategies. This presents itself as an opportunity to optimize different steps in future large scale virtual screening endeavours. Here, we ran a first stage docking simulation against a library of 95 known binders and 3829 generated decoys in an effort to improve the rescoring stage. We then analyzed receptor binding site structure and ligands binding poses to describe their interactions. As a result, we devised a simple and straightforward virtual screening Stage II filtering score based on search space extension followed by a geometric estimation of the ligand—binding site fitness. Having this score, we used an ensemble of receptors generated by Hamiltonian Monte Carlo simulation and reported the results. The improvements shown herein prove that our ensemble docking protocol is suited for APJR and can be easily extrapolated to other GPCRs.     

In Silico Study towards Repositioning of FDA-Approved Drug Candidates for Anticoronaviral Therapy: Molecular Docking,Molecular Dynamics and Binding Free Energy Calculations

The SARS-CoV-2 targets were evaluated for a set of FDA-approved drugs using a combination of drug repositioning and rigorous computational modeling methodologies such as molecular docking and molecular dynamics (MD) simulations followed by binding free energy calculations. Six FDA-approved drugs including, Ouabain, Digitoxin, Digoxin, Proscillaridin, Salinomycin and Niclosamide with promising anti-SARS-CoV-2 activity were screened in silico against four SARS-CoV-2 proteins—papain-like protease (PLpro), RNA-dependent RNA polymerase (RdRp), SARS-CoV-2 main protease (Mpro), and adaptor-associated kinase 1 (AAK1)—in an attempt to define their promising targets. The applied computational techniques suggest that all the tested drugs exhibited excellent binding patterns with higher scores and stable complexes compared to the native protein cocrystallized inhibitors. Ouabain was suggested to act as a dual inhibitor for both PLpro and Mpro enzymes, while Digitoxin bonded perfectly to RdRp. In addition, Salinomycin targeted PLpro. Particularly, Niclosamide was found to target AAK1 with greater affinity compared to the reference drug. Our study provides comprehensive molecular-level insights for identifying or designing novel anti-COVID-19 drugs.     

Fragment Linking and Optimization of Inhibitors of the Aspartic Protease Endothiapepsin: Fragment‐Based Drug Design Facilitated by Dynamic Combinatorial Chemistry

Fragment‐based drug design (FBDD) affords active compounds for biological targets. While there are numerous reports on FBDD by fragment growing/optimization, fragment linking has rarely been reported. Dynamic combinatorial chemistry (DCC) has become a powerful hit‐identification strategy for biological targets. We report the synergistic combination of fragment linking and DCC to identify inhibitors of the aspartic protease endothiapepsin. Based on X‐ray crystal structures of endothiapepsin in complex with fragments, we designed a library of bis‐acylhydrazones and used DCC to identify potent inhibitors. The most potent inhibitor exhibits an IC₅₀ value of 54 nm , which represents a 240‐fold improvement in potency compared to the parent hits. Subsequent X‐ray crystallography validated the predicted binding mode, thus demonstrating the efficiency of the combination of fragment linking and DCC as a hit‐identification strategy. This approach could be applied to a range of biological targets, and holds the potential to facilitate hit‐to‐lead optimization.