Exploring the potential of protein-based pharmacophore models in ligand pose prediction and ranking

Revolutionizing Drug Targeting Strategies: Integrating Artificial Intelligence and Structure-Based Methods in PROTAC Development

PROteolysis TArgeting Chimera (PROTAC) is an emerging technology in chemical biology and drug discovery. This technique facilitates the complete removal of the target proteins that are "undruggable" or challenging to target through chemical molecules via the Ubiquitin-Proteasome System (UPS). PROTACs have been widely explored and outperformed not only in cancer but also in other diseases. During the past few decades, several academic institutes and pharma companies have poured more efforts into PROTAC-related technologies, setting the stage for several major degrader trial readouts in clinical phases. Despite their promising results, the formation of robust ternary orientation, off-target activity, poor permeability, and binding affinity are some of the limitations that hinder their development. Recent advancements in computational technologies have facilitated progress in the development of PROTACs. Researchers have been able to utilize these technologies to explore a wider range of E3 ligases and optimize linkers, thereby gaining a better understanding of the effectiveness and safety of PROTACs in clinical settings. In this review, we briefly explore the computational strategies reported to date for the formation of PROTAC components and discuss the key challenges and opportunities for further research in this area.

Recent Advances in Computer-Aided Structure-Based Drug Design on Ion Channels

Ion channels play important roles in fundamental biological processes, such as electric signaling in cells, muscle contraction, hormone secretion, and regulation of the immune response. Targeting ion channels with drugs represents a treatment option for neurological and cardiovascular diseases, muscular degradation disorders, and pathologies related to disturbed pain sensation. While there are more than 300 different ion channels in the human organism, drugs have been developed only for some of them and currently available drugs lack selectivity. Computational approaches are an indispensable tool for drug discovery and can speed up, especially, the early development stages of lead identification and optimization. The number of molecular structures of ion channels has considerably increased over the last ten years, providing new opportunities for structure-based drug development. This review summarizes important knowledge about ion channel classification, structure, mechanisms, and pathology with the main focus on recent developments in the field of computer-aided, structure-based drug design on ion channels. We highlight studies that link structural data with modeling and chemoinformatic approaches for the identification and characterization of new molecules targeting ion channels. These approaches hold great potential to advance research on ion channel drugs in the future.

Apo2ph4: A Versatile Workflow for the Generation of Receptor-based Pharmacophore Models for Virtual Screening

Pharmacophore models are widely used as efficient virtual screening (VS) filters for the target-directed enrichment of large compound libraries. However, the generation of pharmacophore models that have the power to discriminate between active and inactive molecules traditionally requires structural information about ligand-target complexes or at the very least knowledge of one active ligand. The fact that the discovery of the first known active ligand of a newly investigated target represents a major hurdle at the beginning of every drug discovery project underscores the need for methods that are able to derive high-quality pharmacophore models even without the prior knowledge of any active ligand structures. In this work, we introduce a novel workflow, called apo2ph4, that enables the rapid derivation of pharmacophore models solely from the three-dimensional structure of the target receptor. The utility of this workflow is demonstrated retrospectively for the generation of a pharmacophore model for the M2 muscarinic acetylcholine receptor. Furthermore, in order to show the general applicability of apo2ph4, the workflow was employed for all 15 targets of the recently published LIT-PCBA dataset. Pharmacophore-based VS runs using the apo2ph4-derived models achieved a significant enrichment of actives for 13 targets. In the last presented example, a pharmacophore model derived from the etomidate site of the α1β2γ2 GABAA receptor was used in VS campaigns. Subsequent in vitro testing of selected hits revealed that 19 out of 20 (95%) tested compounds were able to significantly enhance GABA currents, which impressively demonstrates the applicability of apo2ph4 for real-world drug design projects.

FitDock: protein-ligand docking by template fitting

Protein-ligand docking is an essential method in computer-aided drug design and structural bioinformatics. It can be used to identify active compounds and reveal molecular mechanisms of biological processes. A successful docking usually requires thorough conformation sampling and scoring, which are computationally expensive and difficult. Recent studies demonstrated that it can be beneficial to docking with the guidance of existing similar co-crystal structures. In this work, we developed a protein-ligand docking method, named FitDock, which fits initial conformation to the given template using a hierarchical multi-feature alignment approach, subsequently explores the possible conformations and finally outputs refined docking poses. In our comprehensive benchmark tests, FitDock showed 40%-60% improvement in terms of docking success rate and an order of magnitude faster over popular docking methods, if template structures exist (> 0.5 ligand similarity). FitDock has been implemented in a user-friendly program, which could serve as a convenient tool for drug design and molecular mechanism exploration. It is now freely available for academic users at http://cao.labshare.cn/fitdock/.

Molecular Modeling Techniques Applied to the Design of Multitarget Drugs: Methods and Applications

Multifactorial diseases, such as cancer and diabetes present a challenge for the traditional "one-target, one disease" paradigm due to their complex pathogenic mechanisms. Although a combination of drugs can be used, a multitarget drug may be a better choice face of its efficacy, lower adverse effects and lower chance of resistance development. The computer-based design of these multitarget drugs can explore the same techniques used for single-target drug design, but the difficulties associated to the obtention of drugs that are capable of modulating two or more targets with similar efficacy impose new challenges, whose solutions involve the adaptation of known techniques and also to the development of new ones, including machine-learning approaches. In this review, some SBDD and LBDD techniques for the multitarget drug design are discussed, together with some cases where the application of such techniques led to effective multitarget ligands.

AncPhore: A versatile tool for anchor pharmacophore steered drug discovery with applications in discovery of new inhibitors targeting metallo-β-lactamases and indoleamine/tryptophan 2,3-dioxygenases

We herein describe AncPhore, a versatile tool for drug discovery, which is characterized by pharmacophore feature analysis and anchor pharmacophore (i.e., most important pharmacophore features) steered molecular fitting and virtual screening. Comparative analyses of numerous protein–ligand complexes using AncPhore revealed that anchor pharmacophore features are biologically important, commonly associated with protein conservative characteristics, and have significant contributions to the binding affinity. Performance evaluation of AncPhore showed that it had substantially improved prediction ability on different types of target proteins including metalloenzymes by considering the specific contributions and diversity of anchor pharmacophore features. To demonstrate the practicability of AncPhore, we screened commercially available chemical compounds and discovered a set of structurally diverse inhibitors for clinically relevant metallo-β-lactamases (MBLs); of them, 4 and 6 manifested potent inhibitory activity to VIM-2, NDM-1 and IMP-1 MBLs. Crystallographic analyses of VIM-2:4 complex revealed the precise inhibition mode of 4 with VIM-2, highly consistent with the defined anchor pharmacophore features. Besides, we also identified new hit compounds by using AncPhore for indoleamine/tryptophan 2,3-dioxygenases (IDO/TDO), another class of clinically relevant metalloenzymes. This work reveals anchor pharmacophore as a valuable concept for target-centered drug discovery and illustrates the potential of AncPhore to efficiently identify new inhibitors for different types of protein targets.

AutoPH4: An Automated Method for Generating Pharmacophore Models from Protein Binding Pockets

Pharmacophore models are widely used in computational drug discovery (e.g., in the virtual screening of drug molecules) to capture essential information about interactions between ligands and a target protein. Generating pharmacophore models from protein structures is typically a manual process, but there has been growing interest in automated pharmacophore generation methods. Automation makes feasible the processing of large numbers of protein conformations, such as those generated by molecular dynamics (MD) simulations, and thus may help achieve the longstanding goal of incorporating protein flexibility into virtual screening workflows. Here, we present AutoPH4, a new automated method for generating pharmacophore models based on protein structures; we show that a virtual screening workflow incorporating AutoPH4 ranks compounds more accurately than any other pharmacophore-based virtual screening workflow for which results on a public benchmark have been reported. The strong performance of the virtual screening workflow indicates that the AutoPH4 component of the workflow generates high-quality pharmacophores, making AutoPH4 promising for use in future virtual screening workflows as well, such as ones that use conformations generated by MD simulations.

Conservation and Relevance of Pharmacophore Point Types

Pharmacophore models in general use a variety of features for distinct chemical characteristics, such as hydrogen-bond properties, lipohilicity, and ionizability. Usually, features have to match onto their identical type. To clarify if this stringent one-to-one assignment is justified, we investigated a set of 581 unique ligands from the BindingDB with known orientation inside the respective binding pockets and conducted a statistical analysis of the likelihood of observed exchanges in between the pharmacophore features, respectively their degree of conservation. To find out if certain features are obsolete, we derived a ranking to determine the most relevant ones. We found that the most conserved one-to-one feature is the negative ionizable (acids), followed by hydrogen-bond donor, positive ionizable (basic nitrogens), hydrogen-bond acceptor, aromatic, nonaromatic π-systems, and other lipophilic characteristics. The most likely exchanges were found between carboxylate groups and hydrogen-bond acceptors and likewise between basic nitrogens and hydrogen-bond donors, which reflects the characteristics of Lewis acids and bases. Exchanges between hydrogen-bond donors and hydrogen-bond acceptors are hardly more likely than by chance. The kind of target (e.g., kinase, phosphatase, protease, phosphodiesterase, nuclear receptor, metal-containing, or transmembrane protein) did not show substantial influence on the degree of conservation. The relevance of the actual pharmacophore features was found to be strongly dependent on the applied ranking scheme. Mutual information ranks all hydrophobic features as least important, whereas the aromatic feature is put into second place by using a geometric series. Both ranking schemes see the negative ionizable feature of higher significance than the positively ionizable feature.

All in One: Cavity Detection, Druggability Estimate, Cavity-Based Pharmacophore Perception, and Virtual Screening

Discovering the very first ligands of pharmacologically important targets in a fast and cost-efficient manner is an important issue in drug discovery. In the absence of structural information on either endogenous or synthetic ligands, computational chemists classically identify the very first hits by docking compound libraries to a binding site of interest, with well-known biases arising from the usage of scoring functions. We herewith propose a novel computational method tailored to ligand-free protein structures and consisting in the generation of simple cavity-based pharmacophores to which potential ligands could be aligned by the use of a smooth Gaussian function. The method, embedded in the IChem toolkit, automatically detects ligand-binding cavities, then predicts their structural druggability, and last creates a structure-based pharmacophore for predicted druggable binding sites. A companion tool (Shaper2) was designed to align ligands to cavity-derived pharmacophoric features. The proposed method is as efficient as state-of-the-art virtual screening methods (ROCS, Surflex-Dock) in both posing and virtual screening challenges. Interestingly, IChem-Shaper2 is clearly orthogonal to these latter methods in retrieving unique chemotypes from high-throughput virtual screening data.

Truly Target-Focused Pharmacophore Modeling: A Novel Tool for Mapping Intermolecular Surfaces

Pharmacophore models are an accurate and minimal tridimensional abstraction of intermolecular interactions between chemical structures, usually derived from a group of molecules or from a ligand-target complex. Only a limited amount of solutions exists to model comprehensive pharmacophores using the information of a particular target structure without knowledge of any binding ligand. In this work, an automated and customable tool for truly target-focused (T²F) pharmacophore modeling is introduced. Key molecular interaction fields of a macromolecular structure are calculated using the AutoGRID energy functions. The most relevant points are selected by a newly developed filtering cascade and clustered to pharmacophore features with a density-based algorithm. Using five different protein classes, the ability of this method to identify essential pharmacophore features was compared to structure-based pharmacophores derived from ligand-target interactions. This method represents an extremely valuable instrument for drug design in a situation of scarce ligand information available, but also in the case of underexplored therapeutic targets, as well as to investigate protein allosteric pockets and protein-protein interactions.

Development of pharmacophore-based classification model for activators of constitutive androstane receptor

Constitutive androstane receptor (CAR) is predominantly expressed in the liver and is important for regulating drug metabolism and transport. Despite its biological importance, there have been few attempts to develop in silico models to predict the activity of CAR modulated by chemical compounds. The number of in silico studies of CAR may be limited because of CAR's constitutive activity under normal conditions, which makes it difficult to elucidate the key structural features of the interaction between CAR and its ligands. In this study, to address these limitations, we introduced 3D pharmacophore-based descriptors with an integrated ligand and structure-based pharmacophore features, which represent the receptor-ligand interaction. Machine learning methods (support vector machine and artificial neural network) were applied to develop an in silico model with the descriptors containing significant information regarding the ligand binding positions. The best classification model built with a solvent accessibility volume-based filter and the support vector machine showed good predictabilities of 87%, and 85.4% for the training set and validation set, respectively. This demonstrates that our model can be used to accurately predict CAR activators and offers structural information regarding ligand/protein interactions.

Multi-level structure-based pharmacophore modelling of caspase-3-non-peptide complexes: Extracting essential pharmacophore features and its application to virtual screening

Enormous caspase-3-non-peptide crystal structures have been developed to study the structural basis of caspase-3 enzyme inhibition using active site directed small molecular design. These complexes have not been explored thoroughly to decipher the essential non-covalent interactions made by crystal ligands. We present here a multi-level analysis of these caspase-3 complexes using structure-based pharmacophore approach wherein numerous candidate pharmacophore hypotheses were assessed for its ability to cover available caspase-3 small molecular inhibitor dataset. The reliability of the resultant pharmacophores was evaluated using three different validation sets comprising focussed caspase-3 inhibitors, focussed + random decoys, and focussed + structurally similar random decoys and its performance was measured by the Güner-Henry (GH) scoring and enrichment statistics. Furthermore, the effect on excluded volumes toward caspase-3 inhibitors mapping was investigated by an iterative deletion in the structure-based models and created optimal structure-based pharmacophore models to enable effective design of caspase-3 small molecular inhibitor design.

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.

Selectivity by small-molecule inhibitors of protein interactions can be driven by protein surface fluctuations

Small-molecules that inhibit interactions between specific pairs of proteins have long represented a promising avenue for therapeutic intervention in a variety of settings. Structural studies have shown that in many cases, the inhibitor-bound protein adopts a conformation that is distinct from its unbound and its protein-bound conformations. This plasticity of the protein surface presents a major challenge in predicting which members of a protein family will be inhibited by a given ligand. Here, we use biased simulations of Bcl-2-family proteins to generate ensembles of low-energy conformations that contain surface pockets suitable for small molecule binding. We find that the resulting conformational ensembles include surface pockets that mimic those observed in inhibitor-bound crystal structures. Next, we find that the ensembles generated using different members of this protein family are overlapping but distinct, and that the activity of a given compound against a particular family member (ligand selectivity) can be predicted from whether the corresponding ensemble samples a complementary surface pocket. Finally, we find that each ensemble includes certain surface pockets that are not shared by any other family member: while no inhibitors have yet been identified to take advantage of these pockets, we expect that chemical scaffolds complementing these “distinct” pockets will prove highly selective for their targets. The opportunity to achieve target selectivity within a protein family by exploiting differences in surface fluctuations represents a new paradigm that may facilitate design of family-selective small-molecule inhibitors of protein-protein interactions.

Despite intense interest and considerable effort, there are few examples of small molecules that directly inhibit protein-protein interactions. Crystal structures of early successes have highlighted the plasticity of the protein surface, as some inhibitor-bound proteins are captured in conformations that are distinct from both their unbound and their protein-bound conformations. The lack of a single well-defined structure presents a challenge for predicting the members of a protein family to which a given compound will show activity (ligand selectivity). Here we generate ensembles of conformations from simulation, and show that we can predict ligand selectivity based on which family members sample conformations complementary to the ligand. This approach may present a new avenue for designing highly-selective inhibitors of protein-protein interactions.

Pharmacophore-based similarity scoring for DOCK

Pharmacophore modeling incorporates geometric and chemical features of known inhibitors and/or targeted binding sites to rationally identify and design new drug leads. In this study, we have encoded a three-dimensional pharmacophore matching similarity (FMS) scoring function into the structure-based design program DOCK. Validation and characterization of the method are presented through pose reproduction, crossdocking, and enrichment studies. When used alone, FMS scoring dramatically improves pose reproduction success to 93.5% (∼20% increase) and reduces sampling failures to 3.7% (∼6% drop) compared to the standard energy score (SGE) across 1043 protein-ligand complexes. The combined FMS+SGE function further improves success to 98.3%. Crossdocking experiments using FMS and FMS+SGE scoring, for six diverse protein families, similarly showed improvements in success, provided proper pharmacophore references are employed. For enrichment, incorporating pharmacophores during sampling and scoring, in most cases, also yield improved outcomes when docking and rank-ordering libraries of known actives and decoys to 15 systems. Retrospective analyses of virtual screenings to three clinical drug targets (EGFR, IGF-1R, and HIVgp41) using X-ray structures of known inhibitors as pharmacophore references are also reported, including a customized FMS scoring protocol to bias on selected regions in the reference. Overall, the results and fundamental insights gained from this study should benefit the docking community in general, particularly researchers using the new FMS method to guide computational drug discovery with DOCK.

PL-PatchSurfer: a novel molecular local surface-based method for exploring protein-ligand interactions

Structure-based computational methods have been widely used in exploring protein-ligand interactions, including predicting the binding ligands of a given protein based on their structural complementarity. Compared to other protein and ligand representations, the advantages of a surface representation include reduced sensitivity to subtle changes in the pocket and ligand conformation and fast search speed. Here we developed a novel method named PL-PatchSurfer (Protein-Ligand PatchSurfer). PL-PatchSurfer represents the protein binding pocket and the ligand molecular surface as a combination of segmented surface patches. Each patch is characterized by its geometrical shape and the electrostatic potential, which are represented using the 3D Zernike descriptor (3DZD). We first tested PL-PatchSurfer on binding ligand prediction and found it outperformed the pocket-similarity based ligand prediction program. We then optimized the search algorithm of PL-PatchSurfer using the PDBbind dataset. Finally, we explored the utility of applying PL-PatchSurfer to a larger and more diverse dataset and showed that PL-PatchSurfer was able to provide a high early enrichment for most of the targets. To the best of our knowledge, PL-PatchSurfer is the first surface patch-based method that treats ligand complementarity at protein binding sites. We believe that using a surface patch approach to better understand protein-ligand interactions has the potential to significantly enhance the design of new ligands for a wide array of drug-targets.

PharmDock: a pharmacophore-based docking program

Background

Protein-based pharmacophore models are enriched with the information of potential interactions between ligands and the protein target. We have shown in a previous study that protein-based pharmacophore models can be applied for ligand pose prediction and pose ranking. In this publication, we present a new pharmacophore-based docking program PharmDock that combines pose sampling and ranking based on optimized protein-based pharmacophore models with local optimization using an empirical scoring function.

Results

Tests of PharmDock on ligand pose prediction, binding affinity estimation, compound ranking and virtual screening yielded comparable or better performance to existing and widely used docking programs. The docking program comes with an easy-to-use GUI within PyMOL. Two features have been incorporated in the program suite that allow for user-defined guidance of the docking process based on previous experimental data. Docking with those features demonstrated superior performance compared to unbiased docking.

Conclusion

A protein pharmacophore-based docking program, PharmDock, has been made available with a PyMOL plugin. PharmDock and the PyMOL plugin are freely available from http://people.pharmacy.purdue.edu/~mlill/software/pharmdock.

Development of purely structure-based pharmacophores for the topoisomerase I-DNA-ligand binding pocket

Purely structure-based pharmacophores (SBPs) are an alternative method to ligand-based approaches and have the advantage of describing the entire interaction capability of a binding pocket. Here, we present the development of SBPs for topoisomerase I, an anticancer target with an unusual ligand binding pocket consisting of protein and DNA atoms. Different approaches to cluster and select pharmacophore features are investigated, including hierarchical clustering and energy calculations. In addition, the performance of SBPs is evaluated retrospectively and compared to the performance of ligand- and complex-based pharmacophores. SBPs emerge as a valid method in virtual screening and a complementary approach to ligand-focussed methods. The study further reveals that the choice of pharmacophore feature clustering and selection methods has a large impact on the virtual screening hit lists. A prospective application of the SBPs in virtual screening reveals that they can be used successfully to identify novel topoisomerase inhibitors.