A pose prediction approach based on ligand 3D shape similarity

Building shape-focused pharmacophore models for effective docking screening

The performance of molecular docking can be improved by comparing the shape similarity of the flexibly sampled poses against the target proteins' inverted binding cavities. The effectiveness of these pseudo-ligands or negative image-based models in docking rescoring is boosted further by performing enrichment-driven optimization. Here, we introduce a novel shape-focused pharmacophore modeling algorithm O-LAP that generates a new class of cavity-filling models by clumping together overlapping atomic content via pairwise distance graph clustering. Top-ranked poses of flexibly docked active ligands were used as the modeling input and multiple alternative clustering settings were benchmark-tested thoroughly with five demanding drug targets using random training/test divisions. In docking rescoring, the O-LAP modeling typically improved massively on the default docking enrichment; furthermore, the results indicate that the clustered models work well in rigid docking. The C+ +/Qt5-based algorithm O-LAP is released under the GNU General Public License v3.0 via GitHub ( https://github.com/jvlehtonen/overlap-toolkit ). SCIENTIFIC CONTRIBUTION: This study introduces O-LAP, a C++/Qt5-based graph clustering software for generating new type of shape-focused pharmacophore models. In the O-LAP modeling, the target protein cavity is filled with flexibly docked active ligands, the overlapping ligand atoms are clustered, and the shape/electrostatic potential of the resulting model is compared against the flexibly sampled molecular docking poses. The O-LAP modeling is shown to ensure high enrichment in both docking rescoring and rigid docking based on comprehensive benchmark-testing.

Biased Docking for Protein-Ligand Pose Prediction

The interaction between a protein and its ligands is one of the basic and most important processes in biological chemistry. Docking methods aim to predict the molecular 3D structure of protein-ligand complexes starting from coordinates of the protein and the ligand separately. They are widely used in both industry and academia, especially in the context of drug development projects. AutoDock4 is one of the most popular docking tools and, as for any docking method, its performance is highly system dependent. Knowledge about specific protein-ligand interactions on a particular target can be used to successfully overcome this limitation. Here, we describe how to apply the AutoDock Bias protocol, a simple and elegant strategy that allows users to incorporate target-specific information through a modified scoring function that biases the ligand structure towards those poses (or conformations) that establish selected interactions. We discuss two examples using different bias sources. In the first, we show how to steer dockings towards interactions derived from crystal structures of the receptor with different ligands; in the second example, we define and apply hydrophobic biases derived from Molecular Dynamics simulations in mixed solvents. Finally, we discuss general concepts of biased docking, its performance in pose prediction, and virtual screening campaigns as well as other potential applications.

Scope of β-Secretase (BACE1)-Targeted Therapy in Alzheimer's Disease: Emphasizing the Flavonoid Based Natural Scaffold for BACE1 Inhibition

Alzheimer's disease (AD) is a progressive neurodegenerative disease and the most common form of dementia in the world. Studies report the presence of extracellular amyloid plaques consisting of β-amyloid peptide and intracellular tangles consisting of hyperphosphorylated tau proteins as the histopathological indicators of AD. The process of β-amyloid peptide generation by sequential cleavage of amyloid precursor protein by β-secretase (BACE1) and γ-secretase, followed by its aggregation to form amyloid plaques, is the mechanistic basis of the amyloid hypothesis. Other popular hypotheses related to the pathogenesis of AD include the tau hypothesis and the oxidative stress hypothesis. Various targets of the amyloid cascade are now in prime focus to develop drugs for AD. Many BACE1 inhibitors, β-amyloid aggregation inhibitors, and Aβ clearance strategies using monoclonal antibodies are in various stages of clinical trials. This review provides an in-depth evaluation of the role of BACE1 in disease pathogenesis and also highlights the therapeutic approaches developed to find more potent but less toxic inhibitors for BACE1, particularly emphasizing the natural scaffold as a nontoxic lead for BACE1 inhibition. Cellular targets and signaling cascades involving BACE1 have been highlighted to understand the physiological role of BACE1. This knowledge is extremely crucial to understand the toxicity evaluations for BACE1-targeted therapy. We have particularly highlighted the scope of flavonoids as a new generation of nontoxic BACE1 inhibitory scaffolds. The structure-activity relationship of BACE1 inhibition for this group of compounds has been highlighted to provide a guideline to design more selective highly potent inhibitors. The review aims to provide a holistic overview of BACE1-targeted therapy for AD that paves the way for future drug development.

Merging Ligand-Based and Structure-Based Methods in Drug Discovery: An Overview of Combined Virtual Screening Approaches

Virtual screening (VS) is an outstanding cornerstone in the drug discovery pipeline. A variety of computational approaches, which are generally classified as ligand-based (LB) and structure-based (SB) techniques, exploit key structural and physicochemical properties of ligands and targets to enable the screening of virtual libraries in the search of active compounds. Though LB and SB methods have found widespread application in the discovery of novel drug-like candidates, their complementary natures have stimulated continued efforts toward the development of hybrid strategies that combine LB and SB techniques, integrating them in a holistic computational framework that exploits the available information of both ligand and target to enhance the success of drug discovery projects. In this review, we analyze the main strategies and concepts that have emerged in the last years for defining hybrid LB + SB computational schemes in VS studies. Particularly, attention is focused on the combination of molecular similarity and docking, illustrating them with selected applications taken from the literature.

Insights Into the Explication of Potent Tyrosinase Inhibitors with Reference to Computational Studies

Background:

Pigment melanin has primarily a photo defensive role in human skin, its unnecessary production and irregular distribution can cause uneven skin tone ultimately results in hyper pigmentation. Melanin biosynthesis is initiated by tyrosine oxidation through tyrosinase, the key enzyme for melanogenesis. Not only in humans, tyrosinase is also widely distributed in plants and liable for browning of vegetables and fruits. Search for the inhibitors of tyrosinase have been an important target to facilitate development of therapies for the prevention of hyperpigmentary disorders and an undesired browning of vegetables and fruits.

Methods:

Different natural and synthetic chemical compounds have been tested as potential tyrosinase inhibitors, but the mechanism of inhibition is not known, and the quest for information regarding interaction between tyrosinase and its inhibitors is one of the recent areas of research. Computer based methods hence are useful to overcome such issues. Successful utilization of in silico tools like molecular docking simulations make it possible to interpret the tyrosinase and its inhibitor’s intermolecular interactions and helps in identification and development of new and potent tyrosinase inhibitors.

Results:

The present review has pointed out the prominent role of computer aided approaches for the explication of promising tyrosinase inhibitors with a focus on molecular docking approach. Highlighting certain examples of natural compounds whose antityrosinase effects has been evaluated using computational simulations.

Conclusion:

The investigation of new and potent inhibitors of tyrosinase using computational chemistry and bioinformatics will ultimately help millions of peoples to get rid of hyperpigmentary disorders as well as browning of fruits and vegetables.

Improving ligand 3D shape similarity-based pose prediction with a continuum solvent model

In order to improve the pose prediction performance of docking methods, we have previously developed the pose prediction using shape similarity (PoPSS) method. It identifies a ligand conformation of the highest shape similarity with target protein crystal ligands. The identified ligand conformation is then placed into the target protein binding pocket and refined using side-chain repacking and Monte Carlo energy minimization. Subsequently, we have reported a modification to PoPSS, named as PoPSS-Lite, using a simple grid-based energy minimization for side-chain repacking and Tversky correlation coefficient as the similarity metric. This modification has improved the pose prediction performance and PoPSS-Lite was one of the top performers in D3R GC3. Here we report a further modification to PoPSS that utilizes a continuum solvent model to account for water mediated protein ligand interactions. In this approach, named as PoPSS-PB, the ligand conformation of the highest shape similarity with crystal ligands is refined along with the target protein binding site by incorporating the Poisson-Boltzmann electrostatics. The performance of PoPSS-PB along with PoPSS and PoPSS-Lite was prospectively evaluated in D3R GC4. PoPSS-PB not only demonstrated excellent performance with mean and median RMSDs of 1.20 and 1.13 Å but also achieved improved performance over PoPSS and PoPSS-Lite. Furthermore, the comparison with other D3R GC4 pose prediction submissions revealed admirable performance. Our results showed that the binding poses of ligands with unknown binding modes can be successfully predicted by utilizing ligand 3D shape similarity with known crystallographic ligands and that taking the solvation into consideration improves pose prediction.

Binding mode information improves fragment docking

Docking is commonly used in drug discovery to predict how ligand binds to protein target. Best programs are generally able to generate a correct solution, yet often fail to identify it. In the case of drug-like molecules, the correct and incorrect poses can be sorted by similarity to the crystallographic structure of the protein in complex with reference ligands. Fragments are particularly sensitive to scoring problems because they are weak ligands which form few interactions with protein. In the present study, we assessed the utility of binding mode information in fragment pose prediction. We compared three approaches: interaction fingerprints, 3D-matching of interaction patterns and 3D-matching of shapes. We prepared a test set composed of high-quality structures of the Protein Data Bank. We generated and evaluated the docking poses of 586 fragment/protein complexes. We observed that the best approach is twice as accurate as the native scoring function, and that post-processing is less effective for smaller fragments. Interestingly, fragments and drug-like molecules both proved to be useful references. In the discussion, we suggest the best conditions for a successful pose prediction with the three approaches.

D3R Grand Challenge 3: blind prediction of protein-ligand poses and affinity rankings

The Drug Design Data Resource aims to test and advance the state of the art in protein-ligand modeling by holding community-wide blinded, prediction challenges. Here, we report on our third major round, Grand Challenge 3 (GC3). Held 2017-2018, GC3 centered on the protein Cathepsin S and the kinases VEGFR2, JAK2, p38-α, TIE2, and ABL1, and included both pose-prediction and affinity-ranking components. GC3 was structured much like the prior challenges GC2015 and GC2. First, Stage 1 tested pose prediction and affinity ranking methods; then all available crystal structures were released, and Stage 2 tested only affinity rankings, now in the context of the available structures. Unique to GC3 was the addition of a Stage 1b self-docking subchallenge, in which the protein coordinates from all of the cocrystal structures used in the cross-docking challenge were released, and participants were asked to predict the pose of CatS ligands using these newly released structures. We provide an overview of the outcomes and discuss insights into trends and best-practices.

Shape similarity guided pose prediction: lessons from D3R Grand Challenge 3

To extend the utility of ligand 3D shape similarity into pose prediction and virtual screening, we have previously developed CDVS and PoPSS methods. Both of them utilize ligand 3D shape similarity with the crystallographic ligands to improve pose prediction. While CDVS utilizes shape similarity to select suitable receptor structures for molecular docking, PoPSS places a ligand conformation of the highest shape similarity with crystal ligands into the target protein binding pocket which is then refined by side-chain repacking and Monte Carlo energy minimization. Analyses of PoPSS revealed some drawbacks in ligand conformation generation and the scoring scheme used. Moreover, as PoPSS does not sample the ligand conformation after placing it in the binding pocket, it relies solely on conformation generation methods to produce native like conformations. To address these limitations of PoPSS method, we report here a modified approach named as PoPSS-Lite, where side-chain repacking was replaced by a simple grid-based energy minimization. This modification also allowed the sampling of terminal functional groups while keeping the core scaffold fixed. Furthermore, shape similarity calculations were improved by increasing the number of ligand conformations and using a different similarity metric. The performance of PoPSS-Lite was prospectively evaluated in D3R GC3. Comparison of PoPSS-Lite demonstrated superior performance over PoPSS and CDVS with lower mean and median RMSDs. Furthermore, comparison with other D3R GC3 pose prediction submissions revealed top performance for PoPSS-Lite. Our D3R GC3 result extends our perspective that ligand 3D shape similarity with known crystallographic information can be successfully used to predict the binding pose of ligands with unknown binding modes. Our D3R GC3 results further highlight the necessity for improvement in conformer generation methods in order to improve shape similarity guided pose prediction.

Advances in the Development of Shape Similarity Methods and Their Application in Drug Discovery

Molecular similarity is a key concept in drug discovery. It is based on the assumption that structurally similar molecules frequently have similar properties. Assessment of similarity between small molecules has been highly effective in the discovery and development of various drugs. Especially, two-dimensional (2D) similarity approaches have been quite popular due to their simplicity, accuracy and efficiency. Recently, the focus has been shifted toward the development of methods involving the representation and comparison of three-dimensional (3D) conformation of small molecules. Among the 3D similarity methods, evaluation of shape similarity is now gaining attention for its application not only in virtual screening but also in molecular target prediction, drug repurposing and scaffold hopping. A wide range of methods have been developed to describe molecular shape and to determine the shape similarity between small molecules. The most widely used methods include atom distance-based methods, surface-based approaches such as spherical harmonics and 3D Zernike descriptors, atom-centered Gaussian overlay based representations. Several of these methods demonstrated excellent virtual screening performance not only retrospectively but also prospectively. In addition to methods assessing the similarity between small molecules, shape similarity approaches have been developed to compare shapes of protein structures and binding pockets. Additionally, shape comparisons between atomic models and 3D density maps allowed the fitting of atomic models into cryo-electron microscopy maps. This review aims to summarize the methodological advances in shape similarity assessment highlighting advantages, disadvantages and their application in drug discovery.

Discovery of new erbB4 inhibitors: Repositioning an orphan chemical library by inverse virtual screening

Inverse Virtual Screening (IVS) is a docking based approach aimed to the evaluation of the virtual ability of a single compound to interact with a library of proteins. For the first time, we applied this methodology to a library of synthetic compounds, which proved to be inactive towards the target they were initially designed for. Trifluoromethyl-benzenesulfonamides 3-21 were repositioned by means of IVS identifying new lead compounds (14-16, 19 and 20) for the inhibition of erbB4 in the low micromolar range. Among these, compound 20 exhibited an interesting value of IC50 on MCF7 cell lines, thus validating IVS in lead repurposing.

A cross docking pipeline for improving pose prediction and virtual screening performance

Pose prediction and virtual screening performance of a molecular docking method depend on the choice of protein structures used for docking. Multiple structures for a target protein are often used to take into account the receptor flexibility and problems associated with a single receptor structure. However, the use of multiple receptor structures is computationally expensive when docking a large library of small molecules. Here, we propose a new cross-docking pipeline suitable to dock a large library of molecules while taking advantage of multiple target protein structures. Our method involves the selection of a suitable receptor for each ligand in a screening library utilizing ligand 3D shape similarity with crystallographic ligands. We have prospectively evaluated our method in D3R Grand Challenge 2 and demonstrated that our cross-docking pipeline can achieve similar or better performance than using either single or multiple-receptor structures. Moreover, our method displayed not only decent pose prediction performance but also better virtual screening performance over several other methods.

D3R grand challenge 2015: Evaluation of protein-ligand pose and affinity predictions

The Drug Design Data Resource (D3R) ran Grand Challenge 2015 between September 2015 and February 2016. Two targets served as the framework to test community docking and scoring methods: (1) HSP90, donated by AbbVie and the Community Structure Activity Resource (CSAR), and (2) MAP4K4, donated by Genentech. The challenges for both target datasets were conducted in two stages, with the first stage testing pose predictions and the capacity to rank compounds by affinity with minimal structural data; and the second stage testing methods for ranking compounds with knowledge of at least a subset of the ligand-protein poses. An additional sub-challenge provided small groups of chemically similar HSP90 compounds amenable to alchemical calculations of relative binding free energy. Unlike previous blinded Challenges, we did not provide cognate receptors or receptors prepared with hydrogens and likewise did not require a specified crystal structure to be used for pose or affinity prediction in Stage 1. Given the freedom to select from over 200 crystal structures of HSP90 in the PDB, participants employed workflows that tested not only core docking and scoring technologies, but also methods for addressing water-mediated ligand-protein interactions, binding pocket flexibility, and the optimal selection of protein structures for use in docking calculations. Nearly 40 participating groups submitted over 350 prediction sets for Grand Challenge 2015. This overview describes the datasets and the organization of the challenge components, summarizes the results across all submitted predictions, and considers broad conclusions that may be drawn from this collaborative community endeavor.