Efficient virtual screening using multiple protein conformations described as negative images of the ligand-binding site

Deciphering the Lexicon of Protein Targets: A Review on Multifaceted Drug Discovery in the Era of Artificial Intelligence

Understanding protein sequence and structure is essential for understanding protein-protein interactions (PPIs), which are essential for many biological processes and diseases. Targeting protein binding hot spots, which regulate signaling and growth, with rational drug design is promising. Rational drug design uses structural data and computational tools to study protein binding sites and protein interfaces to design inhibitors that can change these interactions, thereby potentially leading to therapeutic approaches. Artificial intelligence (AI), such as machine learning (ML) and deep learning (DL), has advanced drug discovery and design by providing computational resources and methods. Quantum chemistry is essential for drug reactivity, toxicology, drug screening, and quantitative structure-activity relationship (QSAR) properties. This review discusses the methodologies and challenges of identifying and characterizing hot spots and binding sites. It also explores the strategies and applications of artificial-intelligence-based rational drug design technologies that target proteins and protein-protein interaction (PPI) binding hot spots. It provides valuable insights for drug design with therapeutic implications. We have also demonstrated the pathological conditions of heat shock protein 27 (HSP27) and matrix metallopoproteinases (MMP2 and MMP9) and designed inhibitors of these proteins using the drug discovery paradigm in a case study on the discovery of drug molecules for cancer treatment. Additionally, the implications of benzothiazole derivatives for anticancer drug design and discovery are deliberated.

Virtual Screening Strategy to Identify Retinoic Acid-Related Orphan Receptor γt Modulators

Molecular docking is a key method used in virtual screening (VS) campaigns to identify small-molecule ligands for drug discovery targets. While docking provides a tangible way to understand and predict the protein-ligand complex formation, the docking algorithms are often unable to separate active ligands from inactive molecules in practical VS usage. Here, a novel docking and shape-focused pharmacophore VS protocol is demonstrated for facilitating effective hit discovery using retinoic acid receptor-related orphan receptor gamma t (RORγt) as a case study. RORγt is a prospective target for treating inflammatory diseases such as psoriasis and multiple sclerosis. First, a commercial molecular database was flexibly docked. Second, the alternative docking poses were rescored against the shape/electrostatic potential of negative image-based (NIB) models that mirror the target's binding cavity. The compositions of the NIB models were optimized via iterative trimming and benchmarking using a greedy search-driven algorithm or brute force NIB optimization. Third, a pharmacophore point-based filtering was performed to focus the hit identification on the known RORγt activity hotspots. Fourth, free energy binding affinity evaluation was performed on the remaining molecules. Finally, twenty-eight compounds were selected for in vitro testing and eight compounds were determined to be low μM range RORγt inhibitors, thereby showing that the introduced VS protocol generated an effective hit rate of ~29%.

Optimization of Cavity-Based Negative Images to Boost Docking Enrichment in Virtual Screening

Molecular docking is a key in silico method used routinely in modern drug discovery projects. Although docking provides high-quality ligand binding predictions, it regularly fails to separate the active compounds from the inactive ones. In negative image-based rescoring (R-NiB), the shape/electrostatic potential (ESP) of docking poses is compared to the negative image of the protein's ligand binding cavity. While R-NiB often improves the docking yield considerably, the cavity-based models do not reach their full potential without expert editing. Accordingly, a greedy search-driven methodology, brute force negative image-based optimization (BR-NiB), is presented for optimizing the models via iterative editing and benchmarking. Thorough and unbiased training, testing and stringent validation with a multitude of drug targets, and alternative docking software show that BR-NiB ensures excellent docking efficacy. BR-NiB can be considered as a new type of shape-focused pharmacophore modeling, where the optimized models contain only the most vital cavity information needed for effectively filtering docked actives from the inactive or decoy compounds. Finally, the BR-NiB code for performing the automated optimization is provided free-of-charge under MIT license via GitHub (https://github.com/jvlehtonen/brutenib) for boosting the success rates of docking-based virtual screening campaigns.

Negative Image-Based Rescoring: Using Cavity Information to Improve Docking Screening

Molecular docking produces often lackluster results in real-life virtual screening assays that aim to discover novel drug candidates or hit compounds. The problem lies in the inability of the default docking scoring to properly estimate the Gibbs free energy of binding, which impairs the recognition of the best binding poses and the separation of active ligands from inactive compounds. Negative image-based rescoring (R-NiB) provides both effective and efficient way for re-ranking the outputted flexible docking poses to improve the virtual screening yield. Importantly, R-NiB has been shown to work with multiple genuine drug targets and six popular docking algorithms using demanding benchmark test sets. The effectiveness of the R-NiB methodology relies on the shape/electrostatics similarity between the target protein's ligand-binding cavity and the docked ligand poses. In this chapter, the R-NiB method is described with practical usability in mind.

Negative Image-Based Screening: Rigid Docking Using Cavity Information

Rational drug discovery relies heavily on molecular docking-based virtual screening, which samples flexibly the ligand binding poses against the target protein's structure. The upside of flexible docking is that the geometries of the generated docking poses are adjusted to match the residue alignment inside the target protein's ligand-binding pocket. The downside is that the flexible docking requires plenty of computing resources and, regardless, acquiring a decent level of enrichment typically demands further rescoring or post-processing. Negative image-based screening is a rigid docking technique that is ultrafast and computationally light but also effective as proven by vast benchmarking and screening experiments. In the NIB screening, the target protein cavity's shape/electrostatics is aligned and compared against ab initio-generated ligand 3D conformers. In this chapter, the NIB methodology is explained at the practical level and both its weaknesses and strengths are discussed candidly.

A Perspective: Active Role of Lipids in Neurotransmitter Dynamics

Synaptic neurotransmission is generally considered as a function of membrane-embedded receptors and ion channels in response to the neurotransmitter (NT) release and binding. This perspective aims to widen the protein-centric view by including another vital component—the synaptic membrane—in the discussion. A vast set of atomistic molecular dynamics simulations and biophysical experiments indicate that NTs are divided into membrane-binding and membrane-nonbinding categories. The binary choice takes place at the water-membrane interface and follows closely the positioning of the receptors’ binding sites in relation to the membrane. Accordingly, when a lipophilic NT is on route to a membrane-buried binding site, it adheres on the membrane and, then, travels along its plane towards the receptor. In contrast, lipophobic NTs, which are destined to bind into receptors with extracellular binding sites, prefer the water phase. This membrane-based sorting splits the neurotransmission into membrane-independent and membrane-dependent mechanisms and should make the NT binding into the receptors more efficient than random diffusion would allow. The potential implications and notable exceptions to the mechanisms are discussed here. Importantly, maintaining specific membrane lipid compositions (MLCs) at the synapses, especially regarding anionic lipids, affect the level of NT-membrane association. These effects provide a plausible link between the MLC imbalances and neurological diseases such as depression or Parkinson’s disease. Moreover, the membrane plays a vital role in other phases of the NT life cycle, including storage and release from the synaptic vesicles, transport from the synaptic cleft, as well as their synthesis and degradation.

Getting Docking into Shape Using Negative Image-Based Rescoring

The failure of default scoring functions to ensure virtual screening enrichment is a persistent problem for the molecular docking algorithms used in structure-based drug discovery. To remedy this problem, elaborate rescoring and postprocessing schemes have been developed with a varying degree of success, specificity, and cost. The negative image-based rescoring (R-NiB) has been shown to improve the flexible docking performance markedly with a variety of drug targets. The yield improvement is achieved by comparing the alternative docking poses against the negative image of the target protein's ligand-binding cavity. In other words, the shape and electrostatics of the binding pocket is directly used in the similarity comparison to rank the explicit docking poses. Here, the PANTHER/ShaEP-based R-NiB methodology is tested with six popular docking softwares, including GLIDE, PLANTS, GOLD, DOCK, AUTODOCK, and AUTODOCK VINA, using five validated benchmark sets. Overall, the results indicate that R-NiB outperforms the default docking scoring consistently and inexpensively, demonstrating that the methodology is ready for wide-scale virtual screening usage.

Fragment- and negative image-based screening of phosphodiesterase 10A inhibitors

A novel virtual screening methodology called fragment- and negative image-based (F-NiB) screening is introduced and tested experimentally using phosphodiesterase 10A (PDE10A) as a case study. Potent PDE10A-specific small-molecule inhibitors are actively sought after for their antipsychotic and neuroprotective effects. The F-NiB combines features from both fragment-based drug discovery and negative image-based (NIB) screening methodologies to facilitate rational drug discovery. The selected structural parts of protein-bound ligand(s) are seamlessly combined with the negative image of the target's ligand-binding cavity. This cavity- and fragment-based hybrid model, namely its shape and electrostatics, is used directly in the rigid docking of ab initio generated ligand 3D conformers. In total, 14 compounds were acquired using the F-NiB methodology, 3D quantitative structure-activity relationship modeling, and pharmacophore modeling. Three of the small molecules inhibited PDE10A at ~27 to ~67 μM range in a radiometric assay. In a larger context, the study shows that the F-NiB provides a flexible way to incorporate small-molecule fragments into the drug discovery.

A Practical Perspective: The Effect of Ligand Conformers on the Negative Image-Based Screening

Negative image-based (NIB) screening is a rigid molecular docking methodology that can also be employed in docking rescoring. During the NIB screening, a negative image is generated based on the target protein’s ligand-binding cavity by inverting its shape and electrostatics. The resulting NIB model is a drug-like entity or pseudo-ligand that is compared directly against ligand 3D conformers, as is done with a template compound in the ligand-based screening. This cavity-based rigid docking has been demonstrated to work with genuine drug targets in both benchmark testing and drug candidate/lead discovery. Firstly, the study explores in-depth the applicability of different ligand 3D conformer generation software for acquiring the best NIB screening results using cyclooxygenase-2 (COX-2) as the example system. Secondly, the entire NIB workflow from the protein structure preparation, model build-up, and ligand conformer generation to the similarity comparison is performed for COX-2. Accordingly, hands-on instructions are provided on how to employ the NIB methodology from start to finish, both with the rigid docking and docking rescoring using noncommercial software. The practical aspects of the NIB methodology, especially the effect of ligand conformers, are discussed thoroughly, thus, making the methodology accessible for new users.

Discovery of Retinoic Acid-Related Orphan Receptor γt Inverse Agonists via Docking and Negative Image-Based Screening

Retinoic acid-related orphan receptor γt (RORγt) has a vital role in the differentiation of T-helper 17 (TH17) cells. Potent and specific RORγt inverse agonists are sought for treating TH17-related diseases such as psoriasis, rheumatoid arthritis, and type 1 diabetes. Here, the aim was to discover novel RORγt ligands using both standard molecular docking and negative image-based screening. Interestingly, both of these in silico techniques put forward mostly the same compounds for experimental testing. In total, 11 of the 34 molecules purchased for testing were verified as RORγt inverse agonists, thus making the effective hit rate 32%. The pIC₅₀ values for the compounds varied from 4.9 (11 μM) to 6.2 (590 nM). Importantly, the fact that the verified hits represent four different cores highlights the structural diversity of the RORγt inverse agonism and the ability of the applied screening methodologies to facilitate much-desired scaffold hopping for drug design.

Blocking oestradiol synthesis pathways with potent and selective coumarin derivatives

A comprehensive set of 3-phenylcoumarin analogues with polar substituents was synthesised for blocking oestradiol synthesis by 17-β-hydroxysteroid dehydrogenase 1 (HSD1) in the latter part of the sulphatase pathway. Five analogues produced ≥62% HSD1 inhibition at 5 µM and, furthermore, three of them produced ≥68% inhibition at 1 µM. A docking-based structure-activity relationship analysis was done to determine the molecular basis of the inhibition and the cross-reactivity of the analogues was tested against oestrogen receptor, aromatase, cytochrome P450 1A2, and monoamine oxidases. Most of the analogues are only modestly active with 17-β-hydroxysteroid dehydrogenase 2 – a requirement for lowering effective oestradiol levels in vivo. Moreover, the analysis led to the synthesis and discovery of 3-imidazolecoumarin as a potent aromatase inhibitor. In short, coumarin core can be tailored with specific ring and polar moiety substitutions to block either the sulphatase pathway or the aromatase pathway for treating breast cancer and endometriosis.

Improving Docking Performance Using Negative Image-Based Rescoring

Despite the large computational costs of molecular docking, the default scoring functions are often unable to recognize the active hits from the inactive molecules in large-scale virtual screening experiments. Thus, even though a correct binding pose might be sampled during the docking, the active compound or its biologically relevant pose is not necessarily given high enough score to arouse the attention. Various rescoring and post-processing approaches have emerged for improving the docking performance. Here, it is shown that the very early enrichment (number of actives scored higher than 1% of the highest ranked decoys) can be improved on average 2.5-fold or even 8.7-fold by comparing the docking-based ligand conformers directly against the target protein's cavity shape and electrostatics. The similarity comparison of the conformers is performed without geometry optimization against the negative image of the target protein's ligand-binding cavity using the negative image-based (NIB) screening protocol. The viability of the NIB rescoring or the R-NiB, pioneered in this study, was tested with 11 target proteins using benchmark libraries. By focusing on the shape/electrostatics complementarity of the ligand-receptor association, the R-NiB is able to improve the early enrichment of docking essentially without adding to the computing cost. By implementing consensus scoring, in which the R-NiB and the original docking scoring are weighted for optimal outcome, the early enrichment is improved to a level that facilitates effective drug discovery. Moreover, the use of equal weight from the original docking scoring and the R-NiB scoring improves the yield in most cases.

Identification of estrogen receptor α ligands with virtual screening techniques

Utilization of computer-aided molecular discovery methods in virtual screening (VS) is a cost-effective approach to identify novel bioactive small molecules. Unfortunately, no universal VS strategy can guarantee high hit rates for all biological targets, but each target requires distinct, fine-tuned solutions. Here, we have studied in retrospective manner the effectiveness and usefulness of common pharmacophore hypothesis, molecular docking and negative image-based screening as potential VS tools for a widely applied drug discovery target, estrogen receptor α (ERα). The comparison of the methods helps to demonstrate the differences in their ability to identify active molecules. For example, structure-based methods identified an already known active ligand from the widely-used bechmarking decoy molecule set. Although prospective VS against one commercially available database with around 100,000 drug-like molecules did not retrieve many testworthy hits, one novel hit molecule with pIC50 value of 6.6, was identified. Furthermore, our small in-house compound collection of easy-to-synthesize molecules was virtually screened against ERα, yielding to five hit candidates, which were found to be active in vitro having pIC50 values from 5.5 to 6.5.

Ultrafast protein structure-based virtual screening with Panther

Molecular docking is by far the most common method used in protein structure-based virtual screening. This paper presents Panther, a novel ultrafast multipurpose docking tool. In Panther, a simple shape-electrostatic model of the ligand-binding area of the protein is created by utilizing the protein crystal structure. The features of the possible ligands are then compared to the model by using a similarity search algorithm. On average, one ligand can be processed in a few minutes by using classical docking methods, whereas using Panther processing takes <1 s. The presented Panther protocol can be used in several applications, such as speeding up the early phases of drug discovery projects, reducing the number of failures in the clinical phase of the drug development process, and estimating the environmental toxicity of chemicals. Panther-code is available in our web pages (http://www.jyu.fi/panther) free of charge after registration.

Pocket-space maps to identify novel binding-site conformations in proteins

The identification of novel binding-site conformations can greatly assist the progress of structure-based ligand design projects. Diverse pocket shapes drive medicinal chemistry to explore a broader chemical space and thus present additional opportunities to overcome key drug discovery issues such as potency, selectivity, toxicity, and pharmacokinetics. We report a new automated approach to diverse pocket selection, PocketAnalyzer(PCA), which applies principal component analysis and clustering to the output of a grid-based pocket detection algorithm. Since the approach works directly with pocket shape descriptors, it is free from some of the problems hampering methods that are based on proxy shape descriptors, e.g. a set of atomic positional coordinates. The approach is technically straightforward and allows simultaneous analysis of mutants, isoforms, and protein structures derived from multiple sources with different residue numbering schemes. The PocketAnalyzer(PCA) approach is illustrated by the compilation of diverse sets of pocket shapes for aldose reductase and viral neuraminidase. In both cases this allows identification of novel computationally derived binding-site conformations that are yet to be observed crystallographically. Indeed, known inhibitors capable of exploiting these novel binding-site conformations are subsequently identified, thereby demonstrating the utility of PocketAnalyzer(PCA) for rationalizing and improving the understanding of the molecular basis of protein-ligand interaction and bioactivity. A Python program implementing the PocketAnalyzer(PCA) approach is available for download under an open-source license ( http://sourceforge.net/projects/papca/ or http://cpclab.uni-duesseldorf.de/downloads ).

Comparison of virtual high-throughput screening methods for the identification of phosphodiesterase-5 inhibitors

Reliable and effective virtual high-throughput screening (vHTS) methods are desperately needed to minimize the expenses involved in drug discovery projects. Here, we present an improvement to the negative image-based (NIB) screening: the shape, the electrostatics, and the solvation state of the target protein's ligand-binding site are included into the vHTS. Additionally, the initial vHTS results are postprocessed with molecular mechanics/generalized Born surface area (MMGBSA) calculations to estimate the favorability of ligand-protein interactions. The results show that docking produces very good early enrichment for phosphodiesterase-5 (PDE-5); however, in general, the NIB and the ligand-based screening performed better with or without the added electrostatics. Furthermore, the postprocessing of the NIB screening results using MMGBSA calculations improved the early enrichment for the PDE-5 considerably, thus, making hit discovery affordable.