Structure-based druggability assessment--identifying suitable targets for small molecule therapeutics

The IKK-binding domain of NEMO is an irregular coiled coil with a dynamic binding interface

NEMO is an essential component in the activation of the canonical NF-κB pathway and exerts its function by recruiting the IκB kinases (IKK) to the IKK complex. Inhibition of the NEMO/IKKs interaction is an attractive therapeutic paradigm for diseases related to NF-κB mis-regulation, but a difficult endeavor because of the extensive protein-protein interface. Here we report the high-resolution structure of the unbound IKKβ-binding domain of NEMO that will greatly facilitate the design of NEMO/IKK inhibitors. The structures of unbound NEMO show a closed conformation that partially occludes the three binding hot-spots and suggest a facile transition to an open state that can accommodate ligand binding. By fusing coiled-coil adaptors to the IKKβ-binding domain of NEMO, we succeeded in creating a protein with improved solution behavior, IKKβ-binding affinity and crystallization compatibility, which will enable the structural characterization of new NEMO/inhibitor complexes.

In Silico Target Druggability Assessment: From Structural to Systemic Approaches

This chapter will focus on today's in silico direct and indirect approaches to assess therapeutic target druggability. The direct approach tries to infer from the 3D structure the capacity of the target protein to bind small molecule in order to modulate its biological function. Algorithms to recognize and characterize the quality of the ligand interaction sites whether within buried protein cavities or within large protein-protein interface will be reviewed in the first part of the paper. In the case a ligand-binding site is already identified, indirect aspects of target druggability can be assessed. These indirect approaches focus first on target promiscuity and the potential difficulties in developing specific drugs. It is based on large-scale comparison of protein-binding sites. The second aspect concerns the capacity of the target to induce resistant pathway once it is inhibited or activated by a drug. The emergence of drug-resistant pathways can be assessed through systemic analysis of biological networks implementing metabolism and/or cell regulation signaling.

P2Rank: machine learning based tool for rapid and accurate prediction of ligand binding sites from protein structure

Background

Ligand binding site prediction from protein structure has many applications related to elucidation of protein function and structure based drug discovery. It often represents only one step of many in complex computational drug design efforts. Although many methods have been published to date, only few of them are suitable for use in automated pipelines or for processing large datasets. These use cases require stability and speed, which disqualifies many of the recently introduced tools that are either template based or available only as web servers.

Results

We present P2Rank, a stand-alone template-free tool for prediction of ligand binding sites based on machine learning. It is based on prediction of ligandability of local chemical neighbourhoods that are centered on points placed on the solvent accessible surface of a protein. We show that P2Rank outperforms several existing tools, which include two widely used stand-alone tools (Fpocket, SiteHound), a comprehensive consensus based tool (MetaPocket 2.0), and a recent deep learning based method (DeepSite). P2Rank belongs to the fastest available tools (requires under 1 s for prediction on one protein), with additional advantage of multi-threaded implementation.

Conclusions

P2Rank is a new open source software package for ligand binding site prediction from protein structure. It is available as a user-friendly stand-alone command line program and a Java library. P2Rank has a lightweight installation and does not depend on other bioinformatics tools or large structural or sequence databases. Thanks to its speed and ability to make fully automated predictions, it is particularly well suited for processing large datasets or as a component of scalable structural bioinformatics pipelines.

Electronic supplementary material

The online version of this article (10.1186/s13321-018-0285-8) contains supplementary material, which is available to authorized users.

Principles for targeting RNA with drug-like small molecules

Recent studies have indicated the potential to develop small-molecule drugs that act on RNA targets, leading to burgeoning interest in the field. This article discusses general principles for discovering small-molecule drugs that target RNA and argues that the overarching challenge is to identify appropriate target structures in disease-causing RNAs that have high information content and, consequently, appropriate ligand-binding pockets.

RNA molecules are essential for cellular information transfer and gene regulation, and RNAs have been implicated in many human diseases. Messenger and non-coding RNAs contain highly structured elements, and evidence suggests that many of these structures are important for function. Targeting these RNAs with small molecules offers opportunities to therapeutically modulate numerous cellular processes, including those linked to 'undruggable' protein targets. Despite this promise, there is currently only a single class of human-designed small molecules that target RNA used clinically — the linezolid antibiotics. However, a growing number of small-molecule RNA ligands are being identified, leading to burgeoning interest in the field. Here, we discuss principles for discovering small-molecule drugs that target RNA and argue that the overarching challenge is to identify appropriate target structures — namely, in disease-causing RNAs that have high information content and, consequently, appropriate ligand-binding pockets. If focus is placed on such druggable binding sites in RNA, extensive knowledge of the typical physicochemical properties of drug-like small molecules could then enable small-molecule drug discovery for RNA targets to become (only) roughly as difficult as for protein targets.

Subtractive Genomics, Molecular Docking and Molecular Dynamics Simulation Revealed LpxC as a Potential Drug Target Against Multi-Drug Resistant Klebsiella pneumoniae

The emergence and dissemination of pan drug resistant clones of Klebsiella pneumoniae are great threat to public health. In this regard new therapeutic targets must be highlighted to pave the path for novel drug discovery and development. Subtractive proteomic pipeline brought forth UDP-3-O-[3-hydroxymyristoyl] N-acetylglucosamine deacetylase (LpxC), a Zn⁺² dependent cytoplasmic metalloprotein and catalyze the rate limiting deacetylation step of lipid A biosynthesis pathway. Primary sequence analysis followed by 3-dimensional (3-D) structure elucidation of the protein led to the detection of K. pneumoniae LpxC (KpLpxC) topology distinct from its orthologous counterparts in other bacterial species. Molecular docking study of the protein recognized receptor antagonist compound 106, a uridine-based LpxC inhibitory compound, as a ligand best able to fit the binding pocket with a Gold Score of 67.53. Molecular dynamics simulation of docked KpLpxC revealed an alternate binding pattern of ligand in the active site. The ligand tail exhibited preferred binding to the domain I residues as opposed to the substrate binding hydrophobic channel of subdomain II, usually targeted by inhibitory compounds. Comparison with the undocked KpLpxC system demonstrated ligand induced high conformational changes in the hydrophobic channel of subdomain II in KpLpxC. Hence, ligand exerted its inhibitory potential by rendering the channel unstable for substrate binding.

Systematic interrogation of diverse Omic data reveals interpretable, robust, and generalizable transcriptomic features of clinically successful therapeutic targets

Target selection is the first and pivotal step in drug discovery. An incorrect choice may not manifest itself for many years after hundreds of millions of research dollars have been spent. We collected a set of 332 targets that succeeded or failed in phase III clinical trials, and explored whether Omic features describing the target genes could predict clinical success. We obtained features from the recently published comprehensive resource: Harmonizome. Nineteen features appeared to be significantly correlated with phase III clinical trial outcomes, but only 4 passed validation schemes that used bootstrapping or modified permutation tests to assess feature robustness and generalizability while accounting for target class selection bias. We also used classifiers to perform multivariate feature selection and found that classifiers with a single feature performed as well in cross-validation as classifiers with more features (AUROC = 0.57 and AUPR = 0.81). The two predominantly selected features were mean mRNA expression across tissues and standard deviation of expression across tissues, where successful targets tended to have lower mean expression and higher expression variance than failed targets. This finding supports the conventional wisdom that it is favorable for a target to be present in the tissue(s) affected by a disease and absent from other tissues. Overall, our results suggest that it is feasible to construct a model integrating interpretable target features to inform target selection. We anticipate deeper insights and better models in the future, as researchers can reuse the data we have provided to improve methods for handling sample biases and learn more informative features. Code, documentation, and data for this study have been deposited on GitHub at https://github.com/arouillard/omic-features-successful-targets.

Tapping the RNA world for therapeutics

A recent revolution in RNA biology has led to the identification of new RNA classes with unanticipated functions, new types of RNA modifications, an unexpected multiplicity of alternative transcripts and widespread transcription of extragenic regions. This development in basic RNA biology has spawned a corresponding revolution in RNA-based strategies to generate new types of therapeutics. Here, I review RNA-based drug design and discuss barriers to broader applications and possible ways to overcome them. Because they target nucleic acids rather than proteins, RNA-based drugs promise to greatly extend the domain of 'druggable' targets beyond what can be achieved with small molecules and biologics.

Computer-Aided Drug Design in Epigenetics

Epigenetic dysfunction has been widely implicated in several diseases especially cancers thus highlights the therapeutic potential for chemical interventions in this field. With rapid development of computational methodologies and high-performance computational resources, computer-aided drug design has emerged as a promising strategy to speed up epigenetic drug discovery. Herein, we make a brief overview of major computational methods reported in the literature including druggability prediction, virtual screening, homology modeling, scaffold hopping, pharmacophore modeling, molecular dynamics simulations, quantum chemistry calculation, and 3D quantitative structure activity relationship that have been successfully applied in the design and discovery of epi-drugs and epi-probes. Finally, we discuss about major limitations of current virtual drug design strategies in epigenetics drug discovery and future directions in this field.

Big Data in Drug Discovery

Interpretation of Big Data in the drug discovery community should enhance project timelines and reduce clinical attrition through improved early decision making. The issues we encounter start with the sheer volume of data and how we first ingest it before building an infrastructure to house it to make use of the data in an efficient and productive way. There are many problems associated with the data itself including general reproducibility, but often, it is the context surrounding an experiment that is critical to success. Help, in the form of artificial intelligence (AI), is required to understand and translate the context. On the back of natural language processing pipelines, AI is also used to prospectively generate new hypotheses by linking data together. We explain Big Data from the context of biology, chemistry and clinical trials, showcasing some of the impressive public domain sources and initiatives now available for interrogation.

Screening of candidate substrates and coupling ions of transporters by thermostability shift assays

Substrates of most transport proteins have not been identified, limiting our understanding of their role in physiology and disease. Traditional identification methods use transport assays with radioactive compounds, but they are technically challenging and many compounds are unavailable in radioactive form or are prohibitively expensive, precluding large-scale trials. Here, we present a high-throughput screening method that can identify candidate substrates from libraries of unlabeled compounds. The assay is based on the principle that transport proteins recognize substrates through specific interactions, which lead to enhanced stabilization of the transporter population in thermostability shift assays. Representatives of three different transporter (super)families were tested, which differ in structure as well as transport and ion coupling mechanisms. In each case, the substrates were identified correctly from a large set of chemically related compounds, including stereo-isoforms. In some cases, stabilization by substrate binding was enhanced further by ions, providing testable hypotheses on energy coupling mechanisms.

Molecular Dynamics Simulation and Prediction of Druggable Binding Sites

Binding site identification and druggability evaluation are two essential steps in structure-based drug design. A druggable binding site tends to have high binding affinity to drug-like molecules. Predicting such sites can have a significant impact on a drug design campaign. This chapter focuses on summarizing the different methods that are used to predict druggable binding sites. The chapter also discusses the importance of including protein flexibility in the search process and the use of molecular dynamics simulations to address this aspect. Case studies from the literature are also summarized and discussed. We hope that this chapter would provide an overview on the different methods employed in binding site identification evaluation.

Bioinformatics in translational drug discovery

Bioinformatics approaches are becoming ever more essential in translational drug discovery both in academia and within the pharmaceutical industry. Computational exploitation of the increasing volumes of data generated during all phases of drug discovery is enabling key challenges of the process to be addressed. Here, we highlight some of the areas in which bioinformatics resources and methods are being developed to support the drug discovery pipeline. These include the creation of large data warehouses, bioinformatics algorithms to analyse 'big data' that identify novel drug targets and/or biomarkers, programs to assess the tractability of targets, and prediction of repositioning opportunities that use licensed drugs to treat additional indications.

Nanocaged platforms: modification, drug delivery and nanotoxicity. Opening synthetic cages to release the tiger

Nanocages (NCs) have emerged as a new class of drug-carriers, with a wide range of possibilities in multi-modality medical treatments and theranostics. Nanocages can overcome such limitations as high toxicity caused by anti-cancer chemotherapy or by the nanocarrier itself, due to their unique characteristics. These properties consist of: (1) a high loading-capacity (spacious interior); (2) a porous structure (analogous to openings between the bars of the cage); (3) enabling smart release (a key to unlock the cage); and (4) a low likelihood of unfavorable immune responses (the outside of the cage is safe). In this review, we cover different classes of NC structures such as virus-like particles (VLPs), protein NCs, DNA NCs, supramolecular nanosystems, hybrid metal-organic NCs, gold NCs, carbon-based NCs and silica NCs. Moreover, NC-assisted drug delivery including modification methods, drug immobilization, active targeting, and stimulus-responsive release mechanisms are discussed, highlighting the advantages, disadvantages and challenges. Finally, translation of NCs into clinical applications, and an up-to-date assessment of the nanotoxicology considerations of NCs are presented.

Mapping the 3D structures of small molecule binding sites

Background

Analysis of the 3D structures of protein–ligand binding sites can provide valuable insights for drug discovery. Binding site comparison (BSC) studies can be employed to elucidate the function of orphan proteins or to predict the potential for polypharmacology. Many previous binding site analyses only consider binding sites surrounding an experimentally observed bound ligand.

Results

To encompass potential protein–ligand binding sites that do not have ligands known to bind, we have incorporated fpocket cavity detection software and assessed the impact of this inclusion on BSC performance. Using fpocket, we generated a database of ligand-independent potential binding sites and applied the BSC tool, SiteHopper, to analyze similarity relationships between protein binding sites. We developed a method for clustering potential binding sites using a curated dataset of structures for six therapeutically relevant proteins from diverse protein classes in the protein data bank. Two clustering methods were explored; hierarchical clustering and a density-based method adept at excluding noise and outliers from a dataset. We introduce circular plots to visualize binding site structure space. From the datasets analyzed in this study, we highlight a structural relationship between binding sites of cationic trypsin and prothrombin, protein targets known to bind structurally similar small molecules, exemplifying the potential utility of objectively and holistically mapping binding site space from the structural proteome.

Conclusions

We present a workflow for the objective mapping of potential protein–ligand binding sites derived from the currently available structural proteome. We show that ligand-independent binding site detection tools can be introduced without excessive penalty on BSC performance. Clustering combined with intuitive visualization tools can be applied to map relationships between the 3D structures of protein binding sites.Graphical abstract

The Intersection of Structural and Chemical Biology - An Essential Synergy

The continual improvement in our ability to generate high resolution structural models of biological molecules has stimulated and supported innovative chemical biology projects that target increasingly challenging ligand interaction sites. In this review we outline some of the recent developments in chemical biology and rational ligand design and show selected examples that illustrate the synergy between these research areas.

Can We Rely on Computational Predictions To Correctly Identify Ligand Binding Sites on Novel Protein Drug Targets? Assessment of Binding Site Prediction Methods and a Protocol for Validation of Predicted Binding Sites

In the field of medicinal chemistry there is increasing focus on identifying key proteins whose biochemical functions can firmly be linked to serious diseases. Such proteins become targets for drug or inhibitor molecules that could treat or halt the disease through therapeutic action or by blocking the protein function respectively. The protein must be targeted at the relevant biologically active site for drug or inhibitor binding to be effective. As insufficient experimental data is available to confirm the biologically active binding site for novel protein targets, researchers often rely on computational prediction methods to identify binding sites. Presented herein is a short review on structure-based computational methods that (i) predict putative binding sites and (ii) assess the druggability of predicted binding sites on protein targets. This review briefly covers the principles upon which these methods are based, where they can be accessed and their reliability in identifying the correct binding site on a protein target. Based on this review, we believe that these methods are useful in predicting putative binding sites, but as they do not account for the dynamic nature of protein-ligand binding interactions, they cannot definitively identify the correct site from a ranked list of putative sites. To overcome this shortcoming, we strongly recommend using molecular docking to predict the most likely protein-ligand binding site(s) and mode(s), followed by molecular dynamics simulations and binding thermodynamics calculations to validate the docking results. This protocol provides a valuable platform for experimental and computational efforts to design novel drugs and inhibitors that target disease-related proteins.

Exploring the role of water in molecular recognition: predicting protein ligandability using a combinatorial search of surface hydration sites

The interaction between any two biological molecules must compete with their interaction with water molecules. This makes water the most important molecule in medicine, as it controls the interactions of every therapeutic with its target. A small molecule binding to a protein is able to recognize a unique binding site on a protein by displacing bound water molecules from specific hydration sites. Quantifying the interactions of these water molecules allows us to estimate the potential of the protein to bind a small molecule. This is referred to as ligandability. In the study, we describe a method to predict ligandability by performing a search of all possible combinations of hydration sites on protein surfaces. We predict ligandability as the summed binding free energy for each of the constituent hydration sites, computed using inhomogeneous fluid solvation theory. We compared the predicted ligandability with the maximum observed binding affinity for 20 proteins in the human bromodomain family. Based on this comparison, it was determined that effective inhibitors have been developed for the majority of bromodomains, in the range from 10 to 100 nM. However, we predict that more potent inhibitors can be developed for the bromodomains BPTF and BRD7 with relative ease, but that further efforts to develop inhibitors for ATAD2 will be extremely challenging. We have also made predictions for the 14 bromodomains with no reported small molecule K d values by isothermal titration calorimetry. The calculations predict that PBRM1(1) will be a challenging target, while others such as TAF1L(2), PBRM1(4) and TAF1(2), should be highly ligandable. As an outcome of this work, we assembled a database of experimental maximal K d that can serve as a community resource assisting medicinal chemistry efforts focused on BRDs. Effective prediction of ligandability would be a very useful tool in the drug discovery process.

Benzene Probes in Molecular Dynamics Simulations Reveal Novel Binding Sites for Ligand Design

Protein flexibility poses a major challenge in binding site identification. Several computational pocket detection methods that utilize small-molecule probes in molecular dynamics (MD) simulations have been developed to address this issue. Although they have proven hugely successful at reproducing experimental structural data, their ability to predict new binding sites that are yet to be identified and characterized has not been demonstrated. Here, we report the use of benzenes as probe molecules in ligand-mapping MD (LMMD) simulations to predict the existence of two novel binding sites on the surface of the oncoprotein MDM2. One of them was serendipitously confirmed by biophysical assays and X-ray crystallography to be important for the binding of a new family of hydrocarbon stapled peptides that were specifically designed to target the other putative site. These results highlight the predictive power of LMMD and suggest that predictions derived from LMMD simulations can serve as a reliable basis for the identification of novel ligand binding sites in structure-based drug design.

Computational prediction of formulation strategies for beyond-rule-of-5 compounds

The physicochemical properties of some contemporary drug candidates are moving towards higher molecular weight, and coincidentally also higher lipophilicity in the quest for biological selectivity and specificity. These physicochemical properties move the compounds towards beyond rule-of-5 (B-r-o-5) chemical space and often result in lower water solubility. For such B-r-o-5 compounds non-traditional delivery strategies (i.e. those other than conventional tablet and capsule formulations) typically are required to achieve adequate exposure after oral administration. In this review, we present the current status of computational tools for prediction of intestinal drug absorption, models for prediction of the most suitable formulation strategies for B-r-o-5 compounds and models to obtain an enhanced understanding of the interplay between drug, formulation and physiological environment. In silico models are able to identify the likely molecular basis for low solubility in physiologically relevant fluids such as gastric and intestinal fluids. With this baseline information, a formulation scientist can, at an early stage, evaluate different orally administered, enabling formulation strategies. Recent computational models have emerged that predict glass-forming ability and crystallisation tendency and therefore the potential utility of amorphous solid dispersion formulations. Further, computational models of loading capacity in lipids, and therefore the potential for formulation as a lipid-based formulation, are now available. Whilst such tools are useful for rapid identification of suitable formulation strategies, they do not reveal drug localisation and molecular interaction patterns between drug and excipients. For the latter, Molecular Dynamics simulations provide an insight into the interplay between drug, formulation and intestinal fluid. These different computational approaches are reviewed. Additionally, we analyse the molecular requirements of different targets, since these can provide an early signal that enabling formulation strategies will be required. Based on the analysis we conclude that computational biopharmaceutical profiling can be used to identify where non-conventional gateways, such as prediction of 'formulate-ability' during lead optimisation and early development stages, are important and may ultimately increase the number of orally tractable contemporary targets.

Translating Neurogenomics Into New Medicines

Brain disorders remain one of the defining challenges of modern medicine and among the most poorly served with new therapeutics. Advances in human neurogenetics have begun to shed light on the genomic architecture of complex diseases of mood, cognition, brain development, and neurodegeneration. From genome-wide association studies to rare variants, these findings hold promise for defining the pathogenesis of brain disorders that have resisted simple molecular description. However, the path from genetics to new medicines is far from clear and can take decades, even for the most well-understood genetic disorders. In this review, we define three challenges for the field of neurogenetics that we believe must be addressed to translate human genetics efficiently into new therapeutics for brain disorders.

Ultra-High-Throughput Structure-Based Virtual Screening for Small-Molecule Inhibitors of Protein-Protein Interactions

Protein-protein interactions play important roles in virtually all cellular processes, making them enticing targets for modulation by small-molecule therapeutics: specific examples have been well validated in diseases ranging from cancer and autoimmune disorders, to bacterial and viral infections. Despite several notable successes, however, overall these remain a very challenging target class. Protein interaction sites are especially challenging for computational approaches, because the target protein surface often undergoes a conformational change to enable ligand binding: this confounds traditional approaches for virtual screening. Through previous studies, we demonstrated that biased "pocket optimization" simulations could be used to build collections of low-energy pocket-containing conformations, starting from an unbound protein structure. Here, we demonstrate that these pockets can further be used to identify ligands that complement the protein surface. To do so, we first build from a given pocket its "exemplar": a perfect, but nonphysical, pseudoligand that would optimally match the shape and chemical features of the pocket. In our previous studies, we used these exemplars to quantitatively compare protein surface pockets to one another. Here, we now introduce this exemplar as a template for pharmacophore-based screening of chemical libraries. Through a series of benchmark experiments, we demonstrate that this approach exhibits comparable performance as traditional docking methods for identifying known inhibitors acting at protein interaction sites. However, because this approach is predicated on ligand/exemplar overlays, and thus does not require explicit calculation of protein-ligand interactions, exemplar screening provides a tremendous speed advantage over docking: 6 million compounds can be screened in about 15 min on a single 16-core, dual-GPU computer. The extreme speed at which large compound libraries can be traversed easily enables screening against a "pocket-optimized" ensemble of protein conformations, which in turn facilitates identification of more diverse classes of active compounds for a given protein target.

Prediction of Druggable Proteins Using Machine Learning and Systems Biology: A Mini-Review

The emergence of -omics technologies has allowed the collection of vast amounts of data on biological systems. Although, the pace of such collection has been exponential, the impact of these data remains small on many critical biomedical applications such as drug development. Limited resources, high costs, and low hit-to-lead ratio have led researchers to search for more cost effective methodologies. A possible alternative is to incorporate computational methods of potential drug target prediction early during drug discovery workflow. Computational methods based on systems approaches have the advantage of taking into account the global properties of a molecule not limited to its sequence, structure or function. Machine learning techniques are powerful tools that can extract relevant information from massive and noisy data sets. In recent years the scientific community has explored the combined power of these fields to propose increasingly accurate and low cost methods to propose interesting drug targets. In this mini-review, we describe promising approaches based on the simultaneous use of systems biology and machine learning to access gene and protein druggability. Moreover, we discuss the state-of-the-art of this emerging and interdisciplinary field, discussing data sources, algorithms and the performance of the different methodologies. Finally, we indicate interesting avenues of research and some remaining open challenges.

Lessons from Hot Spot Analysis for Fragment-Based Drug Discovery

Analysis of binding energy hot spots at protein surfaces can provide crucial insights into the prospects for successful application of fragment-based drug discovery (FBDD), and whether a fragment hit can be advanced into a high-affinity, drug-like ligand. The key factor is the strength of the top ranking hot spot, and how well a given fragment complements it. We show that published data are sufficient to provide a sophisticated and quantitative understanding of how hot spots derive from a protein 3D structure, and how their strength, number, and spatial arrangement govern the potential for a surface site to bind to fragment-sized and larger ligands. This improved understanding provides important guidance for the effective application of FBDD in drug discovery.

A Call for Systematic Research on Solute Carriers

Solute carrier (SLC) membrane transport proteins control essential physiological functions, including nutrient uptake, ion transport, and waste removal. SLCs interact with several important drugs, and a quarter of the more than 400 SLC genes are associated with human diseases. Yet, compared to other gene families of similar stature, SLCs are relatively understudied. The time is right for a systematic attack on SLC structure, specificity, and function, taking into account kinship and expression, as well as the dependencies that arise from the common metabolic space.

To Hit or Not to Hit, That Is the Question - Genome-wide Structure-Based Druggability Predictions for Pseudomonas aeruginosa Proteins

Pseudomonas aeruginosa is a Gram-negative bacterium known to cause opportunistic infections in immune-compromised or immunosuppressed individuals that often prove fatal. New drugs to combat this organism are therefore sought after. To this end, we subjected the gene products of predicted perturbative genes to structure-based druggability predictions using DrugPred. Making this approach suitable for large-scale predictions required the introduction of new methods for calculation of descriptors, development of a workflow to identify suitable pockets in homologous proteins and establishment of criteria to obtain valid druggability predictions based on homologs. We were able to identify 29 perturbative proteins of P. aeruginosa that may contain druggable pockets, including some of them with no or no drug-like inhibitors deposited in ChEMBL. These proteins form promising novel targets for drug discovery against P. aeruginosa.

New Frontiers in Druggability

A powerful early approach to evaluating the druggability of proteins involved determining the hit rate in NMR-based screening of a library of small compounds. Here, we show that a computational analog of this method, based on mapping proteins using small molecules as probes, can reliably reproduce druggability results from NMR-based screening and can provide a more meaningful assessment in cases where the two approaches disagree. We apply the method to a large set of proteins. The results show that, because the method is based on the biophysics of binding rather than on empirical parametrization, meaningful information can be gained about classes of proteins and classes of compounds beyond those resembling validated targets and conventionally druglike ligands. In particular, the method identifies targets that, while not druggable by druglike compounds, may become druggable using compound classes such as macrocycles or other large molecules beyond the rule-of-five limit.

DARC: Mapping Surface Topography by Ray-Casting for Effective Virtual Screening at Protein Interaction Sites

Protein-protein interactions represent an exciting and challenging target class for therapeutic intervention using small molecules. Protein interaction sites are often devoid of the deep surface pockets presented by "traditional" drug targets, and crystal structures reveal that inhibitors typically engage these sites using very shallow binding modes. As a consequence, modern virtual screening tools developed to identify inhibitors of traditional drug targets do not perform as well when they are instead deployed at protein interaction sites. To address the need for novel inhibitors of important protein interactions, here we introduce an alternate docking strategy specifically designed for this regime. Our method, termed DARC (Docking Approach using Ray-Casting), matches the topography of a surface pocket "observed" from within the protein to the topography "observed" when viewing a potential ligand from the same vantage point. We applied DARC to carry out a virtual screen against the protein interaction site of human antiapoptotic protein Mcl-1 and found that four of the top-scoring 21 compounds showed clear inhibition in a biochemical assay. The Ki values for these compounds ranged from 1.2 to 21 μM, and each had ligand efficiency comparable to promising small-molecule inhibitors of other protein-protein interactions. These hit compounds do not resemble the natural (protein) binding partner of Mcl-1, nor do they resemble any known inhibitors of Mcl-1. Our results thus demonstrate the utility of DARC for identifying novel inhibitors of protein-protein interactions.

PockDrug: A Model for Predicting Pocket Druggability That Overcomes Pocket Estimation Uncertainties

Predicting protein druggability is a key interest in the target identification phase of drug discovery. Here, we assess the pocket estimation methods' influence on druggability predictions by comparing statistical models constructed from pockets estimated using different pocket estimation methods: a proximity of either 4 or 5.5 Å to a cocrystallized ligand or DoGSite and fpocket estimation methods. We developed PockDrug, a robust pocket druggability model that copes with uncertainties in pocket boundaries. It is based on a linear discriminant analysis from a pool of 52 descriptors combined with a selection of the most stable and efficient models using different pocket estimation methods. PockDrug retains the best combinations of three pocket properties which impact druggability: geometry, hydrophobicity, and aromaticity. It results in an average accuracy of 87.9% ± 4.7% using a test set and exhibits higher accuracy (∼5-10%) than previous studies that used an identical apo set. In conclusion, this study confirms the influence of pocket estimation on pocket druggability prediction and proposes PockDrug as a new model that overcomes pocket estimation variability.

In silico design of low molecular weight protein-protein interaction inhibitors: Overall concept and recent advances

Protein-protein interactions (PPIs) are carrying out diverse functions in living systems and are playing a major role in the health and disease states. Low molecular weight (LMW) "drug-like" inhibitors of PPIs would be very valuable not only to enhance our understanding over physiological processes but also for drug discovery endeavors. However, PPIs were deemed intractable by LMW chemicals during many years. But today, with the new experimental and in silico technologies that have been developed, about 50 PPIs have already been inhibited by LMW molecules. Here, we first focus on general concepts about protein-protein interactions, present a consensual view about ligandable pockets at the protein interfaces and the possibilities of using fast and cost effective structure-based virtual screening methods to identify PPI hits. We then discuss the design of compound collections dedicated to PPIs. Recent financial analyses of the field suggest that LMW PPI modulators could be gaining momentum over biologics in the coming years supporting further research in this area.

Therapeutic opportunities within the DNA damage response

The DNA damage response (DDR) is essential for maintaining the genomic integrity of the cell, and its disruption is one of the hallmarks of cancer. Classically, defects in the DDR have been exploited therapeutically in the treatment of cancer with radiation therapies or genotoxic chemotherapies. More recently, protein components of the DDR systems have been identified as promising avenues for targeted cancer therapeutics. Here, we present an in-depth analysis of the function, role in cancer and therapeutic potential of 450 expert-curated human DDR genes. We discuss the DDR drugs that have been approved by the US Food and Drug Administration (FDA) or that are under clinical investigation. We examine large-scale genomic and expression data for 15 cancers to identify deregulated components of the DDR, and we apply systematic computational analysis to identify DDR proteins that are amenable to modulation by small molecules, highlighting potential novel therapeutic targets.

Genetic Approaches to Facilitate Antibacterial Drug Development

Very few chemically novel agents have been approved for antibacterial chemotherapies during the last 50 yr. Yet new antibacterial drugs are needed to reduce the impact on global health of an increasing number of drug-resistant infections, including highly drug-resistant forms of tuberculosis. This review discusses how genetic approaches can be used to study the mechanism of action of whole-cell screening hits and facilitate target-driven strategies for antimicrobial drug development.

Synthesis and tau RNA binding evaluation of ametantrone-containing ligands

We describe the synthesis and characterization of ametantrone-containing RNA ligands based on the derivatization of this intercalator with two neamine moieties (Amt-Nea,Nea) or with one azaquinolone heterocycle and one neamine (Amt-Nea,Azq) as well as its combination with guanidinoneamine (Amt-NeaG4). Biophysical studies revealed that guanidinylation of the parent ligand (Amt-Nea) had a positive effect on the binding of the resulting compound for Tau pre-mRNA target as well as on the stabilization upon complexation of some of the mutated RNA sequences associated with the development of tauopathies. Further studies by NMR revealed the existence of a preferred binding site in the stem-loop structure, in which ametantrone intercalates in the characteristic bulged region. Regarding doubly-functionalized ligands, binding affinity and stabilizing ability of Amt-Nea,Nea were similar to those of the guanidinylated ligand, but the two aminoglycoside fragments seem to interfere with its accommodation in a single binding site. However, Amt-Nea,Azq binds at the bulged region in a similar way than Amt-NeaG4. Overall, these results provide new insights on fine-tuning RNA binding properties of ametantrone by single or double derivatization with other RNA recognition motifs, which could help in the future design of new ligands with improved selectivity for disease-causing RNA molecules.

Binding site druggability assessment in fragment-based drug design

Target druggability refers to the propensity that a particular target is amenable to bind high-affinity drug-like molecules. A robust yet accurate computational assessment of target druggability would greatly benefit the fields of chemical genomics and drug discovery. Here, we illustrate a structure-based computational protocol to quantitatively assess the target binding-site druggability via in silico screening a fragment-like compound library. In particular, we provide guidelines, suggestions, and critical thoughts on different aspects of this computational protocol, including: construction of fragment library, preparation of target structure, in silico fragment screening, and analysis of druggability.

Disulfide-mediated stabilization of the IκB kinase binding domain of NF-κB essential modulator (NEMO)

Human NEMO (NF-κB essential modulator) is a 419 residue scaffolding protein that, together with catalytic subunits IKKα and IKKβ, forms the IκB kinase (IKK) complex, a key regulator of NF-κB pathway signaling. NEMO is an elongated homodimer comprising mostly α-helix. It has been shown that a NEMO fragment spanning residues 44–111, which contains the IKKα/β binding site, is structurally disordered in the absence of bound IKKβ. Herein we show that enforcing dimerization of NEMO_1–120 or NEMO_44–111 constructs through introduction of one or two interchain disulfide bonds, through oxidation of the native Cys54 residue and/or at position 107 through a Leu107Cys mutation, induces a stable α-helical coiled-coil structure that is preorganized to bind IKKβ with high affinity. Chemical and thermal denaturation studies showed that, in the context of a covalent dimer, the ordered structure was stabilized relative to the denatured state by up to 3 kcal/mol. A full-length NEMO-L107C protein formed covalent dimers upon treatment of mammalian cells with H₂O₂. Furthermore, NEMO-L107C bound endogenous IKKβ in A293T cells, reconstituted TNF-induced NF-κB signaling in NEMO-deficient cells, and interacted with TRAF6. Our results indicate that the IKKβ binding domain of NEMO possesses an ordered structure in the unbound state, provided that it is constrained within a dimer as is the case in the constitutively dimeric full-length NEMO protein. The stability of the NEMO coiled coil is maintained by strong interhelix interactions in the region centered on residue 54. The disulfide-linked constructs we describe herein may be useful for crystallization of NEMO’s IKKβ binding domain in the absence of bound IKKβ, thereby facilitating the structural characterization of small-molecule inhibitors.

Using chiral molecules as an approach to address low-druggability recognition sites

The content of chiral carbon atoms or structural complexity, which is known to correlate well with relevant physicochemical properties of small molecules, represents a promising descriptor that could fill the gap in existing drug discovery between ligand library filtering rules and the corresponding properties of the target's recognition site. Herein, we present an in silico study on the yet unclear underlying correlations between molecular complexity and other more sophisticated physicochemical and biological properties. By analyzing thousands of protein-ligand complexes from DrugBank, we show that increasing molecular complexity of drugs is an approach to addressing particularly low-druggability and polar recognition sites. We also show that biologically relevant protein classes characteristically bind molecules with a certain degree of structural complexity. Three distinct behaviors toward drug recognition are described. The reported results set the basis for a better understanding of protein-drug recognition, and open the possibility of including target information in the filtering of large ligand libraries for screening.

Fractal Dimensions of Macromolecular Structures

Quantifying the properties of macromolecules is a prerequisite for understanding their roles in biochemical processes. One of the less-explored geometric features of macromolecules is molecular surface irregularity, or 'roughness', which can be measured in terms of fractal dimension (D). In this study, we demonstrate that surface roughness correlates with ligand binding potential. We quantified the surface roughnesses of biological macromolecules in a large-scale survey that revealed D values between 2.0 and 2.4. The results of our study imply that surface patches involved in molecular interactions, such as ligand-binding pockets and protein-protein interfaces, exhibit greater local fluctuations in their fractal dimensions than 'inert' surface areas. We expect approximately 22 % of a protein's surface outside of the crystallographically known ligand binding sites to be ligandable. These findings provide a fresh perspective on macromolecular structure and have considerable implications for drug design as well as chemical and systems biology.

SPILLO-PBSS: detecting hidden binding sites within protein 3D-structures through a flexible structure-based approach

The study reports a flexible structure-based approach aimed at identifying binding sites within target proteins starting from a well-defined reference binding site. The method, named SPILLO potential binding sites searcher (SPILLO-PBSS), includes a suitably designed tolerance which allows an efficient recognition of the potential binding sites regardless of both involved residues and protein conformation. Hence, the proposed method overcomes the rigidity which affects the available approaches and which prevents a proper analysis of distorted binding sites. We apply SPILLO-PBSS to several test cases, including the search for the guanosine diphosphate binding site in distorted H-Ras proteins and the identification of acetylcholine binding proteins from among a library of heterogeneous resolved proteins. Tests are also performed to compare SPILLO-PBSS with other related and available methods. The encouraging results confirm the notable potentialities of this approach and lay the foundation for its use to analyze and predict target proteins on a proteome-wide scale.

Bioengineered vaults: self-assembling protein shell-lipophilic core nanoparticles for drug delivery

We report a novel approach to a new class of bioengineered, monodispersed, self-assembling vault nanoparticles consisting of a protein shell exterior with a lipophilic core interior designed for drug and probe delivery. Recombinant vaults were engineered to contain a small amphipathic α-helix derived from the nonstructural protein 5A of hepatitis C virus, thereby creating within the vault lumen a lipophilic microenvironment into which lipophilic compounds could be reversibly encapsulated. Multiple types of electron microscopy showed that attachment of this peptide resulted in larger than expected additional mass internalized within the vault lumen attributable to incorporation of host lipid membrane constituents spanning the vault waist (>35 nm). These bioengineered lipophilic vaults reversibly associate with a sample set of therapeutic compounds, including all-trans retinoic acid, amphotericin B, and bryostatin 1, incorporating hundreds to thousands of drug molecules per vault nanoparticle. Bryostatin 1 is of particular therapeutic interest because of its ability to potently induce expression of latent HIV, thus representing a preclinical lead in efforts to eradicate HIV/AIDS. Vaults loaded with bryostatin 1 released free drug, resulting in activation of HIV from provirus latency in vitro and induction of CD69 biomarker expression following intravenous injection into mice. The ability to preferentially and reversibly encapsulate lipophilic compounds into these novel bioengineered vault nanoparticles greatly advances their potential use as drug delivery systems.

Structure-based druggability assessment of the mammalian structural proteome with inclusion of light protein flexibility

Advances reported over the last few years and the increasing availability of protein crystal structure data have greatly improved structure-based druggability approaches. However, in practice, nearly all druggability estimation methods are applied to protein crystal structures as rigid proteins, with protein flexibility often not directly addressed. The inclusion of protein flexibility is important in correctly identifying the druggability of pockets that would be missed by methods based solely on the rigid crystal structure. These include cryptic pockets and flexible pockets often found at protein-protein interaction interfaces. Here, we apply an approach that uses protein modeling in concert with druggability estimation to account for light protein backbone movement and protein side-chain flexibility in protein binding sites. We assess the advantages and limitations of this approach on widely-used protein druggability sets. Applying the approach to all mammalian protein crystal structures in the PDB results in identification of 69 proteins with potential druggable cryptic pockets.

Modelling three-dimensional protein structures for applications in drug design

A structural perspective of drug target and anti-target proteins, and their molecular interactions with biologically active molecules, largely advances many areas of drug discovery, including target validation, hit and lead finding and lead optimisation. In the absence of experimental 3D structures, protein structure prediction often offers a suitable alternative to facilitate structure-based studies. This review outlines recent methodical advances in homology modelling, with a focus on those techniques that necessitate consideration of ligand binding. In this context, model quality estimation deserves special attention because the accuracy and reliability of different structure prediction techniques vary considerably, and the quality of a model ultimately determines its usefulness for structure-based drug discovery. Examples of G-protein-coupled receptors (GPCRs) and ADMET-related proteins were selected to illustrate recent progress and current limitations of protein structure prediction. Basic guidelines for good modelling practice are also provided.

Drug-Like Protein-Protein Interaction Modulators: Challenges and Opportunities for Drug Discovery and Chemical Biology

[Formula: see text] Fundamental processes in living cells are largely controlled by macromolecular interactions and among them, protein-protein interactions (PPIs) have a critical role while their dysregulations can contribute to the pathogenesis of numerous diseases. Although PPIs were considered as attractive pharmaceutical targets already some years ago, they have been thus far largely unexploited for therapeutic interventions with low molecular weight compounds. Several limiting factors, from technological hurdles to conceptual barriers, are known, which, taken together, explain why research in this area has been relatively slow. However, this last decade, the scientific community has challenged the dogma and became more enthusiastic about the modulation of PPIs with small drug-like molecules. In fact, several success stories were reported both, at the preclinical and clinical stages. In this review article, written for the 2014 International Summer School in Chemoinformatics (Strasbourg, France), we discuss in silico tools (essentially post 2012) and databases that can assist the design of low molecular weight PPI modulators (these tools can be found at www.vls3d.com). We first introduce the field of protein-protein interaction research, discuss key challenges and comment recently reported in silico packages, protocols and databases dedicated to PPIs. Then, we illustrate how in silico methods can be used and combined with experimental work to identify PPI modulators.

High-resolution structure of the M14-type cytosolic carboxypeptidase from Burkholderia cenocepacia refined exploiting PDB_REDO strategies

A potential cytosolic metallocarboxypeptidase from Burkholderia cenocepacia has been crystallized and a synchrotron-radiation microfocus beamline allowed the acquisition of diffraction data to 1.9 Å resolution. The asymmetric unit comprises a tetramer containing over 1500 amino acids, and the high-throughput automated protocols embedded in PDB_REDO were coupled with model-map inspections in refinement. This approach has highlighted the value of such protocols for efficient analyses. The subunit is constructed from two domains. The N-terminal domain has previously only been observed in cytosolic carboxypeptidase (CCP) proteins. The C-terminal domain, which carries the Zn2+-containing active site, serves to classify this protein as a member of the M14D subfamily of carboxypeptidases. Although eukaryotic CCPs possess deglutamylase activity and are implicated in processing modified tubulin, the function and substrates of the bacterial family members remain unknown. The B. cenocepacia protein did not display deglutamylase activity towards a furylacryloyl glutamate derivative, a potential substrate. Residues previously shown to coordinate the divalent cation and that contribute to peptide-bond cleavage in related enzymes such as bovine carboxypeptidase are conserved. The location of a conserved basic patch in the active site adjacent to the catalytic Zn2+, where an acetate ion is identified, suggests recognition of the carboxy-terminus in a similar fashion to other carboxypeptidases. However, there are significant differences that indicate the recognition of substrates with different properties. Of note is the presence of a lysine in the S1' recognition subsite that suggests specificity towards an acidic substrate.

Identifying druggable targets by protein microenvironments matching: application to transcription factors

Druggability of a protein is its potential to be modulated by drug-like molecules. It is important in the target selection phase. We hypothesize that: (i) known drug-binding sites contain advantageous physicochemical properties for drug binding, or “druggable microenvironments” and (ii) given a target, the presence of multiple druggable microenvironments similar to those seen previously is associated with a high likelihood of druggability. We developed DrugFEATURE to quantify druggability by assessing the microenvironments in potential small-molecule binding sites. We benchmarked DrugFEATURE using two data sets. One data set measures druggability using NMR-based screening. DrugFEATURE correlates well with this metric. The second data set is based on historical drug discovery outcomes. Using the DrugFEATURE cutoffs derived from the first, we accurately discriminated druggable and difficult targets in the second. We further identified novel druggable transcription factors with implications for cancer therapy. DrugFEATURE provides useful insight for drug discovery, by evaluating druggability and suggesting specific regions for interacting with drug-like molecules.

Computational methods in drug discovery

Computer-aided drug discovery/design methods have played a major role in the development of therapeutically important small molecules for over three decades. These methods are broadly classified as either structure-based or ligand-based methods. Structure-based methods are in principle analogous to high-throughput screening in that both target and ligand structure information is imperative. Structure-based approaches include ligand docking, pharmacophore, and ligand design methods. The article discusses theory behind the most important methods and recent successful applications. Ligand-based methods use only ligand information for predicting activity depending on its similarity/dissimilarity to previously known active ligands. We review widely used ligand-based methods such as ligand-based pharmacophores, molecular descriptors, and quantitative structure-activity relationships. In addition, important tools such as target/ligand data bases, homology modeling, ligand fingerprint methods, etc., necessary for successful implementation of various computer-aided drug discovery/design methods in a drug discovery campaign are discussed. Finally, computational methods for toxicity prediction and optimization for favorable physiologic properties are discussed with successful examples from literature.

Structural properties of non-traditional drug targets present new challenges for virtual screening

Traditional drug targets have historically included signaling proteins that respond to small molecules and enzymes that use small molecules as substrates. Increasing attention is now being directed toward other types of protein targets, in particular those that exert their function by interacting with nucleic acids or other proteins rather than small-molecule ligands. Here, we systematically compare existing examples of inhibitors of protein-protein interactions to inhibitors of traditional drug targets. While both sets of inhibitors bind with similar potency, we find that the inhibitors of protein-protein interactions typically bury a smaller fraction of their surface area upon binding to their protein targets. The fact that an average atom is less buried suggests that more atoms are needed to achieve a given potency, explaining the observation that ligand efficiency is typically poor for inhibitors of protein-protein interactions. We then carried out a series of docking experiments and found a further consequence of these relatively exposed binding modes is that structure-based virtual screening may be more difficult: such binding modes do not provide sufficient clues to pick out active compounds from decoy compounds. Collectively, these results suggest that the challenges associated with such non-traditional drug targets may not lie with identifying compounds that potently bind to the target protein surface, but rather with identifying compounds that bind in a sufficiently buried manner to achieve good ligand efficiency and, thus, good oral bioavailability. While the number of available crystal structures of distinct protein interaction sites bound to small-molecule inhibitors is relatively small at present (only 21 such complexes were included in this study), these are sufficient to draw conclusions based on the current state of the field; as additional data accumulate it will be exciting to refine the viewpoint presented here. Even with this limited perspective however, we anticipate that these insights, together with new methods for exploring protein conformational fluctuations, may prove useful for identifying the "low-hanging fruit" among non-traditional targets for therapeutic intervention.