The binding database: overview and user's guide

Machine learning-assisted search for novel coagulants: When machine learning can be efficient even if data availability is low

Design of new drugs is a challenging process: a candidate molecule should satisfy multiple conditions to act properly and make the least side-effect-perfect candidates selectively attach to and influence only targets, leaving off-targets intact. The amount of experimental data about various properties of molecules constantly grows, promoting data-driven approaches. However, the applicability of typical predictive machine learning techniques can be substantially limited by a lack of experimental data about a particular target. For example, there are many known Thrombin inhibitors (acting as anticoagulants), but a very limited number of known Protein C inhibitors (coagulants). In this study, we present our approach to suggest new inhibitor candidates by building an effective representation of chemical space. For this aim, we developed a deep learning model-autoencoder, trained on a large set of molecules in the SMILES format to map the chemical space. Further, we applied different sampling strategies to generate novel coagulant candidates. Symmetrically, we tested our approach on anticoagulant candidates, where we were able to predict their inhibition towards Thrombin. We also compare our approach with MegaMolBART-another deep learning generative model, but exploiting similar principles of navigation in a chemical space.

Integrating single-cell and multi-omic approaches reveals Euphorbiae Humifusae Herba-dependent mitochondrial dysfunction in non-small-cell lung cancer

Euphorbiae Humifusae Herba (EHH) is a pivotal therapeutic agent with diverse pharmacological effects. However, a substantial gap exists in understanding its pharmacological properties and anti-tumour mechanisms. This study aimed to address this gap by exploring EHH's pharmacological properties, identifying NSCLC therapy-associated protein targets, and elucidating how EHH induces mitochondrial disruption in NSCLC cells, offering insights into novel NSCLC treatment strategies. String database was utilized to explore protein-protein interactions. Subsequently, single-cell analysis and multi-omics further unveiled the impact of EHH-targeted genes on the immune microenvironment of NSCLC, as well as their influence on immunotherapeutic responses. Finally, both in vivo and in vitro experiments elucidated the anti-tumour mechanisms of EHH, specifically through the assessment of mitochondrial ROS levels and alterations in mitochondrial membrane potential. EHH exerts its influence through engagement with a cluster of 10 genes, including the apoptotic gene CASP3. This regulatory impact on the immune milieu within NSCLC holds promise as an indicator for predicting responses to immunotherapy. Besides, EHH demonstrated the capability to induce mitochondrial ROS generation and perturbations in mitochondrial membrane potential in NSCLC cells, ultimately leading to mitochondrial dysfunction and consequent apoptosis of tumour cells. EHH induces mitochondrial disruption in NSCLC cells, leading to cell apoptosis to inhibit the progress of NSCLC.

Artificial intelligence methods in kinase target profiling: Advances and challenges

Kinases have a crucial role in regulating almost the full range of cellular processes, making them essential targets for therapeutic interventions against various diseases. Accurate kinase-profiling prediction is vital for addressing the selectivity/specificity challenges in kinase drug discovery, which is closely related to lead optimization, drug repurposing, and the understanding of potential drug side effects. In this review, we provide an overview of the latest advancements in machine learning (ML)-based and deep learning (DL)-based quantitative structure-activity relationship (QSAR) models for kinase profiling. We highlight current trends in this rapidly evolving field and discuss the existing challenges and future directions regarding experimental data set construction and model architecture design. Our aim is to offer practical insights and guidance for the development and utilization of these approaches.

ALPACA: A machine Learning Platform for Affinity and selectivity profiling of CAnnabinoids receptors modulators

The development of small molecules that selectively target the cannabinoid receptor subtype 2 (CB2R) is emerging as an intriguing therapeutic strategy to treat neurodegeneration, as well as to contrast the onset and progression of cancer. In this context, in-silico tools able to predict CB2R affinity and selectivity with respect to the subtype 1 (CB1R), whose modulation is responsible for undesired psychotropic effects, are highly desirable. In this work, we developed a series of machine learning classifiers trained on high-quality bioactivity data of small molecules acting on CB2R and/or CB1R extracted from ChEMBL v30. Our classifiers showed strong predictive power in accurately determining CB2R affinity, CB1R affinity, and CB2R/CB1R selectivity. Among the built models, those obtained using random forest as algorithm proved to be the top-performing ones (AUC in validation ≥0.96) and were made freely accessible through a user-friendly web platform developed ad hoc and called ALPACA (https://www.ba.ic.cnr.it/softwareic/alpaca/). Due to its user-friendly interface and robust predictive power, ALPACA can be a valuable tool in saving both time and resources involved in the design of selective CB2R modulators.

Can molecular dynamics simulations improve predictions of protein-ligand binding affinity with machine learning?

Binding affinity prediction largely determines the discovery efficiency of lead compounds in drug discovery. Recently, machine learning (ML)-based approaches have attracted much attention in hopes of enhancing the predictive performance of traditional physics-based approaches. In this study, we evaluated the impact of structural dynamic information on the binding affinity prediction by comparing the models trained on different dimensional descriptors, using three targets (i.e. JAK1, TAF1-BD2 and DDR1) and their corresponding ligands as the examples. Here, 2D descriptors are traditional ECFP4 fingerprints, 3D descriptors are the energy terms of the Smina and NNscore scoring functions and 4D descriptors contain the structural dynamic information derived from the trajectories based on molecular dynamics (MD) simulations. We systematically investigate the MD-refined binding affinity prediction performance of three classical ML algorithms (i.e. RF, SVR and XGB) as well as two common virtual screening methods, namely Glide docking and MM/PBSA. The outcomes of the ML models built using various dimensional descriptors and their combinations reveal that the MD refinement with the optimized protocol can improve the predictive performance on the TAF1-BD2 target with considerable structural flexibility, but not for the less flexible JAK1 and DDR1 targets, when taking docking poses as the initial structure instead of the crystal structures. The results highlight the importance of the initial structures to the final performance of the model through conformational analysis on the three targets with different flexibility.

Sunsetting Binding MOAD with its last data update and the addition of 3D-ligand polypharmacology tools

Binding MOAD is a database of protein-ligand complexes and their affinities with many structured relationships across the dataset. The project has been in development for over 20 years, but now, the time has come to bring it to a close. Currently, the database contains 41,409 structures with affinity coverage for 15,223 (37%) complexes. The website BindingMOAD.org provides numerous tools for polypharmacology exploration. Current relationships include links for structures with sequence similarity, 2D ligand similarity, and binding-site similarity. In this last update, we have added 3D ligand similarity using ROCS to identify ligands which may not necessarily be similar in two dimensions but can occupy the same three-dimensional space. For the 20,387 different ligands present in the database, a total of 1,320,511 3D-shape matches between the ligands were added. Examples of the utility of 3D-shape matching in polypharmacology are presented. Finally, plans for future access to the project data are outlined.

Identification of 11-Hydroxytephrosin and Torosaflavone A as Potential Inhibitors of 3-Phosphoinositide-Dependent Protein Kinase 1 (PDPK1): Toward Anticancer Drug Discovery

The 3-phosphoinositide-dependent protein kinase 1 (PDPK1) has a significant role in cancer progression and metastasis as well as other inflammatory disorders, and has been proposed as a promising therapeutic target for several malignancies. In this work, we conducted a systematic virtual screening of natural compounds from the IMPPAT database to identify possible PDPK1 inhibitors. Primarily, the Lipinski rules, ADMET, and PAINS filter were applied and then the binding affinities, docking scores, and selectivity were carried out to find effective hits against PDPK1. Finally, we identified two natural compounds, 11-Hydroxytephrosin and Torosaflavone A, bearing substantial affinity with PDPK1. Both compounds showed drug-likeness as predicted by the ADMET analysis and their physicochemical parameters. These compounds preferentially bind to the ATP-binding pocket of PDPK1 and interact with functionally significant residues. The conformational dynamics and complex stability of PDPK1 with the selected compounds were then studied using interaction analysis and molecular dynamics (MD) simulations for 100 ns. The simulation results revealed that PDPK1 forms stable docked complexes with the elucidated compounds. The findings show that the newly discovered 11-Hydroxytephrosin and Torosaflavone A bind to PDPK1 in an ATP-competitive manner, suggesting that they could one day be used as therapeutic scaffolds against PDPK1-associated diseases including cancer.

Complex machine learning model needs complex testing: Examining predictability of molecular binding affinity by a graph neural network

Drug discovery pipelines typically involve high-throughput screening of large amounts of compounds in a search of potential drugs candidates. As a chemical space of small organic molecules is huge, a "navigation" over it urges for fast and lightweight computational methods, thus promoting machine-learning approaches for processing huge pools of candidates. In this contribution, we present a graph-based deep neural network for prediction of protein-drug binding affinity and assess its predictive power under thorough testing conditions. Within the suggested approach, both protein and drug molecules are represented as graphs and passed to separate graph sub-networks, then concatenated and regressed towards a binding affinity. The neural network is trained on two binding affinity datasets-PDBbind and data imported from RCSB Protein Data Bank. In order to explore the generalization capabilities of the model we go beyond traditional random or leave-cluster-out techniques and demonstrate the need for more elaborate model performance assessment - six different strategies for test/train data partitioning (random, time- and property-arranged, protein- and ligand-clustered) with a k-fold cross-validation are engaged. Finally, we discuss the model performance in terms of a set of metrics for different split strategies and fold arrangement. Our code is available at https://github.com/SoftServeInc/affinity-by-GNN.

Inhibition of mutant RAS-RAF interaction by mimicking structural and dynamic properties of phosphorylated RAS

Undruggability of RAS proteins has necessitated alternative strategies for the development of effective inhibitors. In this respect, phosphorylation has recently come into prominence as this reversible post-translational modification attenuates sensitivity of RAS towards RAF. As such, in this study, we set out to unveil the impact of phosphorylation on dynamics of HRAS^WT and aim to invoke similar behavior in HRAS^G12D mutant by means of small therapeutic molecules. To this end, we performed molecular dynamics (MD) simulations using phosphorylated HRAS and showed that phosphorylation of Y32 distorted Switch I, hence the RAS/RAF interface. Consequently, we targeted Switch I in HRAS^G12D by means of approved therapeutic molecules and showed that the ligands enabled detachment of Switch I from the nucleotide-binding pocket. Moreover, we demonstrated that displacement of Switch I from the nucleotide-binding pocket was energetically more favorable in the presence of the ligand. Importantly, we verified computational findings in vitro where HRAS^G12D/RAF interaction was prevented by the ligand in HEK293T cells that expressed HRAS^G12D mutant protein. Therefore, these findings suggest that targeting Switch I, hence making Y32 accessible might open up new avenues in future drug discovery strategies that target mutant RAS proteins.

A review on network pharmacology based phytotherapy in treating diabetes- An environmental perspective

Diabetes has become common lifestyle disorder associated with obesity and cardiovascular diseases. Environmental factors like physical inactivity, polluted surroundings and unhealthy dieting also plays a vital role in diabetes pathogenesis. As the current anti-diabetic drugs possess unprecedented side effects, traditional herbal medicine can be used an alternative therapy. The paramount challenge with the herbal formulation usage is the lack of standardized procedure, entangled with little knowledge on drug safety and mechanism of drug action. Heavy metal contamination is a major environmental hazard where plants tend to accumulate toxic metals like nickel, chromium and lead through industrial and agricultural activities. It becomes inappropriate to use these plants for phytotherapy as it may affect the human health on long term consumption. This review discuss about the environmental risk factors related to diabetes and better implication of medicinal plants in anti-diabetic therapy using network pharmacology. It is an in silico analytical tool that helps to unravel the multi-targeted action of herbal formulations rich in secondary metabolites. Also, a special focus is attempted to pool the databases regarding the medicinal plants for diabetes and associated diseases, their bioactive compounds, possible diabetic targets, drug-target interaction and toxicology reports that may open an aisle in safer, effective and toxicity-free drug discovery.

Assessing the Anti-inflammatory Mechanism of Reduning Injection by Network Pharmacology

Reduning Injection (RDNI) is a traditional Chinese medicine formula indicated for the treatment of inflammatory diseases. However, the molecular mechanism of RDNI is unclear. The information of RDNI ingredients was collected from previous studies. Targets of them were obtained by data mining and molecular docking. The information of targets and related pathways was collected in UniProt and KEGG. Networks were constructed and analyzed by Cytoscape to identify key compounds, targets, and pathways. Data mining and molecular docking identified 11 compounds, 84 targets, and 201 pathways that are related to the anti-inflammatory activity of RDNI. Network analysis identified two key compounds (caffeic acid and ferulic acid), five key targets (Bcl-2, eNOS, PTGS2, PPARA, and MMPs), and four key pathways (estrogen signaling pathway, PI3K-AKT signaling pathway, cGMP-PKG signaling pathway, and calcium signaling pathway) which would play critical roles in the treatment of inflammatory diseases by RDNI. The cross-talks among pathways provided a deeper understanding of anti-inflammatory effect of RDNI. RDNI is capable of regulating multiple biological processes and treating inflammation at a systems level. Network pharmacology is a practical approach to explore the therapeutic mechanism of TCM for complex disease.

A Computational Approach with Biological Evaluation: Combinatorial Treatment of Curcumin and Exemestane Synergistically Regulates DDX3 Expression in Cancer Cell Lines

DDX3 belongs to RNA helicase family that demonstrates oncogenic properties and has gained wider attention due to its role in cancer progression, proliferation and transformation. Mounting reports have evidenced the role of DDX3 in cancers making it a promising target to abrogate DDX3 triggered cancers. Dual pharmacophore models were generated and were subsequently validated. They were used as 3D queries to screen the InterBioScreen database, resulting in the selection of curcumin that was escalated to molecular dynamics simulation studies. In vitro anti-cancer analysis was conducted on three cell lines such as MCF-7, MDA-MB-231 and HeLa, which were evaluated along with exemestane. Curcumin was docked into the active site of the protein target (PDB code 2I4I) to estimate the binding affinity. The compound has interacted with two key residues and has displayed stable molecular dynamics simulation results. In vitro analysis has demonstrated that both the candidate compounds have reduced the expression of DDX3 in three cell lines. However, upon combinatorial treatment of curcumin (10 and 20 μM) and exemestane (50 μM) a synergism was exhibited, strikingly downregulating the DDX3 expression and has enhanced apoptosis in three cell lines. The obtained results illuminate the use of curcumin as an alternative DDX3 inhibitor and can serve as a chemical scaffold to design new small molecules.

Review and comparative assessment of similarity-based methods for prediction of drug–protein interactions in the druggable human proteome

Drug–protein interactions (DPIs) underlie the desired therapeutic actions and the adverse side effects of a significant majority of drugs. Computational prediction of DPIs facilitates research in drug discovery, characterization and repurposing. Similarity-based methods that do not require knowledge of protein structures are particularly suitable for druggable genome-wide predictions of DPIs. We review 35 high-impact similarity-based predictors that were published in the past decade. We group them based on three types of similarities and their combinations that they use. We discuss and compare key aspects of these methods including source databases, internal databases and their predictive models. Using our novel benchmark database, we perform comparative empirical analysis of predictive performance of seven types of representative predictors that utilize each type of similarity individually and all possible combinations of similarities. We assess predictive quality at the database-wide DPI level and we are the first to also include evaluation over individual drugs. Our comprehensive analysis shows that predictors that use more similarity types outperform methods that employ fewer similarities, and that the model combining all three types of similarities secures area under the receiver operating characteristic curve of 0.93. We offer a comprehensive analysis of sensitivity of predictive performance to intrinsic and extrinsic characteristics of the considered predictors. We find that predictive performance is sensitive to low levels of similarities between sequences of the drug targets and several extrinsic properties of the input drug structures, drug profiles and drug targets. The benchmark database and a webserver for the seven predictors are freely available at http://biomine.cs.vcu.edu/servers/CONNECTOR/.

Small molecule purine and pseudopurine derivatives: synthesis, cytostatic evaluations and investigation of growth inhibitory effect in non-small cell lung cancer A549

Novel halogenated purines and pseudopurines with diverse aryl-substituted 1,2,3-triazoles were prepared. While p-(trifluoromethyl)-substituted 1,2,3-triazole in N-9 alkylated purine and 3-deazapurine was critical for strong albeit unselective activity on pancreatic adenocarcinoma cells CFPAC-1,1-(p-fluorophenyl)-1,2,3-triazole derivative of 7-deazapurine showed selective cytostatic effect on metastatic colon cancer cells SW620. Importantly, 1-(p-chlorophenyl)-1,2,3-triazole-tagged benzimidazole displayed the most pronounced and highly selective inhibitory effect in nM range on non-small cell lung cancer A549. This compound revealed to target molecular processes at the extracellular side and inside the plasma membrane regulated by GPLD1 and growth factor receptors PDGFR and IGF-1R leading to the inhibition of cell proliferation and induction of apoptosis mediated by p38 MAP kinase and NF-κB, respectively. Further optimisation of this compound as to reduce its toxicity in normal cells may lead to the development of novel agent effective against lung cancer.

Design, synthesis and biological evaluation of novel benzimidazole amidines as potent multi-target inhibitors for the treatment of non-small cell lung cancer

A series of novel amidino 2-substituted benzimidazoles linked to 1,4-disubstituted 1,2,3-triazoles were synthesized by implementation of microwave and ultrasound irradiation in click reaction and subsequent condensation of thus obtained 4-(1,2,3-triazol-1-yl)benzaldehyde with o-phenylenediamines. In vitro antiproliferative screening of compounds performed on human cancer cell lines revealed that p-chlorophenyl-substituted 1,2,3-triazolyl N-isopropylamidine 10c and benzyl-substituted 1,2,3-triazolyl imidazoline 11f benzimidazoles had selective and potent cytostatic activities in the low nM range against non-small cell lung cancer cell line A549, which could be attributed to induction of apoptosis and primary necrosis. Additional Western blot analyses showed different mechanisms of cytostatic activity between compounds 10c and 11f that could be associated with the nature of aromatic substituent at 1-(1,2,3-triazolyl) and amidino moiety at C-5 position of benzimidazole ring. Specifically, compound 11f abrogated the activity of several protein kinases including TGM2, CDK9, SK1 and p38 MAPK, whereas compound 10c did not have profound effect on the activities of CDK9 and TGM2, but instead showed moderate downregulation of SK1 activity concomitant with a significant reduction in p38 MAPK. Further in silico structural analysis demonstrated that compound 11f bound slightly better to the ATP binding site of p38 MAPK compared to 10c, which correlated well with observed stronger decrement in the expression level of phospho-p38 MAPK elicited by 11f in comparison with 10c.

An integrated anti-arrhythmic target network of a Chinese medicine compound, Wenxin Keli, revealed by combined machine learning and molecular pathway analysis

Deciphering the compound Wenxin Keli's anti-arrhythmic pharmacological mechanism by integrating network pharmacology and experimental verification methods.

Lead identification for the K-Ras protein: virtual screening and combinatorial fragment-based approaches

OBJECTIVE

Kirsten rat sarcoma (K-Ras) protein is a member of Ras family belonging to the small guanosine triphosphatases superfamily. The members of this family share a conserved structure and biochemical properties, acting as binary molecular switches. The guanosine triphosphate-bound active K-Ras interacts with a range of effectors, resulting in the stimulation of downstream signaling pathways regulating cell proliferation, differentiation, and apoptosis. Efforts to target K-Ras have been unsuccessful until now, placing it among high-value molecules against which developing a therapy would have an enormous impact. K-Ras transduces signals when it binds to guanosine triphosphate by directly binding to downstream effector proteins, but in case of guanosine diphosphate-bound conformation, these interactions get disrupted.

METHODS

In the present study, we targeted the nucleotide-binding site in the "on" and "off" state conformations of the K-Ras protein to find out suitable lead compounds. A structure-based virtual screening approach has been used to screen compounds from different databases, followed by a combinatorial fragment-based approach to design the apposite lead for the K-Ras protein.

RESULTS

Interestingly, the designed compounds exhibit a binding preference for the "off" state over "on" state conformation of K-Ras protein. Moreover, the designed compounds' interactions are similar to guanosine diphosphate and, thus, could presumably act as a potential lead for K-Ras. The predicted drug-likeness properties of these compounds suggest that these compounds follow the Lipinski's rule of five and have tolerable absorption, distribution, metabolism, excretion and toxicity values.

CONCLUSION

Thus, through the current study, we propose targeting only "off" state conformations as a promising strategy for the design of reversible inhibitors to pharmacologically inhibit distinct conformations of K-Ras protein.

Connecting proteins with drug-like compounds: Open source drug discovery workflows with BindingDB and KNIME

Today's large, public databases of protein-small molecule interaction data are creating important new opportunities for data mining and integration. At the same time, new graphical user interface-based workflow tools offer facile alternatives to custom scripting for informatics and data analysis. Here, we illustrate how the large protein-ligand database BindingDB may be incorporated into KNIME workflows as a step toward the integration of pharmacological data with broader biomolecular analyses. Thus, we describe a collection of KNIME workflows that access BindingDB data via RESTful webservices and, for more intensive queries, via a local distillation of the full BindingDB dataset. We focus in particular on the KNIME implementation of knowledge-based tools to generate informed hypotheses regarding protein targets of bioactive compounds, based on notions of chemical similarity. A number of variants of this basic approach are tested for seven existing drugs with relatively ill-defined therapeutic targets, leading to replication of some previously confirmed results and discovery of new, high-quality hits. Implications for future development are discussed. Database URL: www.bindingdb.org.

Activity cliffs and activity cliff generators based on chemotype-related activity landscapes

Activity cliffs have large impact in drug discovery; therefore, their detection and quantification are of major importance. This work introduces the metric activity cliff enrichment factor and expands the previously reported activity cliff generator concept by adding chemotype information to representations of the activity landscape. To exemplify these concepts, three molecular databases with multiple biological activities were characterized. Compounds in each database were grouped into chemotype classes. Then, pairwise comparisons of structure similarities and activity differences were calculated for each compound and used to construct chemotype-based structure-activity similarity (SAS) maps. Different landscape distributions among four major regions of the SAS maps were observed for different subsets of molecules grouped in chemotypes. Based on this observation, the activity cliff enrichment factor was calculated to numerically detect chemotypes enriched in activity cliffs. Several chemotype classes were detected having major proportion of activity cliffs than the entire database. In addition, some chemotype classes comprising compounds with smooth structure activity relationships (SAR) were detected. Finally, the activity cliff generator concept was applied to compounds grouped in chemotypes to extract valuable SAR information.

A structural chemogenomics analysis of aminergic GPCRs: lessons for histamine receptor ligand design

BACKGROUND AND PURPOSE

Chemogenomics focuses on the discovery of new connections between chemical and biological space leading to the discovery of new protein targets and biologically active molecules. G-protein coupled receptors (GPCRs) are a particularly interesting protein family for chemogenomics studies because there is an overwhelming amount of ligand binding affinity data available. The increasing number of aminergic GPCR crystal structures now for the first time allows the integration of chemogenomics studies with high-resolution structural analyses of GPCR-ligand complexes.

EXPERIMENTAL APPROACH

In this study, we have combined ligand affinity data, receptor mutagenesis studies, and amino acid sequence analyses to high-resolution structural analyses of (hist)aminergic GPCR-ligand interactions. This integrated structural chemogenomics analysis is used to more accurately describe the molecular and structural determinants of ligand affinity and selectivity in different key binding regions of the crystallized aminergic GPCRs, and histamine receptors in particular.

KEY RESULTS

Our investigations highlight interesting correlations and differences between ligand similarity and ligand binding site similarity of different aminergic receptors. Apparent discrepancies can be explained by combining detailed analysis of crystallized or predicted protein-ligand binding modes, receptor mutation studies, and ligand structure-selectivity relationships that identify local differences in essential pharmacophore features in the ligand binding sites of different receptors.

CONCLUSIONS AND IMPLICATIONS

We have performed structural chemogenomics studies that identify links between (hist)aminergic receptor ligands and their binding sites and binding modes. This knowledge can be used to identify structure-selectivity relationships that increase our understanding of ligand binding to (hist)aminergic receptors and hence can be used in future GPCR ligand discovery and design.

CASE plots for the chemotype-based activity and selectivity analysis: a CASE study of cyclooxygenase inhibitors

Structure-activity characterization of molecular databases plays a central role in drug discovery. However, the characterization of large databases containing structurally diverse molecules with several end-points represents a major challenge. For this purpose, the use of chemoinformatic methods plays an important role to elucidate structure-activity relationships. Herein, a general methodology, namely Chemotype Activity and Selectivity Enrichment plots, is presented. Chemotype Activity and Selectivity Enrichment plots provide graphical information concerning the activity and selectivity patterns of particular chemotypes contained in structurally diverse databases. As a case study, we analyzed a set of 658 compounds screened against cyclooxygenase-1 and cyclooxygenase-2. Chemotype Activity and Selectivity Enrichment plots analysis highlighted chemotypes enriched with active and selective molecules against cyclooxygenase-2; all this in a simple 2D graphical representation. Additionally, the most active and selective chemotypes detected in Chemotype Activity and Selectivity Enrichment plots were analyzed separately using the previously reported dual activity-difference maps. These findings indicate that Chemotype Activity and Selectivity Enrichment plots and dual activity-difference maps are complementary chemoinformatic tools to explore the structure-activity relationships of structurally diverse databases screened against two biological end-points.

A prospective cross-screening study on G-protein-coupled receptors: lessons learned in virtual compound library design

We present the systematic prospective evaluation of a protein-based and a ligand-based virtual screening platform against a set of three G-protein-coupled receptors (GPCRs): the β-2 adrenoreceptor (ADRB2), the adenosine A(2A) receptor (AA2AR), and the sphingosine 1-phosphate receptor (S1PR1). Novel bioactive compounds were identified using a consensus scoring procedure combining ligand-based (frequent substructure ranking) and structure-based (Snooker) tools, and all 900 selected compounds were screened against all three receptors. A striking number of ligands showed affinity/activity for GPCRs other than the intended target, which could be partly attributed to the fuzziness and overlap of protein-based pharmacophore models. Surprisingly, the phosphodiesterase 5 (PDE5) inhibitor sildenafil was found to possess submicromolar affinity for AA2AR. Overall, this is one of the first published prospective chemogenomics studies that demonstrate the identification of novel cross-pharmacology between unrelated protein targets. The lessons learned from this study can be used to guide future virtual ligand design efforts.

2D QSAR Studies of Several Potent Aminopyridine, Anilinopyrimidine and Pyridine Carboxamide-based JNK Inhibitors

The c-Jan N-terminal kinases are members of the mitogen activated protein kinase family of signaling proteins. Amino pyridine based compounds, 4-anilino pyrimidine derivatives, and 2-pyridine carboxamide derivatives have been identified as potent JNK inhibitors with good cellular activity. In this study we calculated molecular topological and quantum chemical descriptors of 15 training compounds and three quantitative structure activity relationships models have been constructed. The significance of three models is judged on the basis of correlation, Fischer F test and quality factor (Q). This study is helpful for screening potent inhibitors of protein kinases.

Snooker: a structure-based pharmacophore generation tool applied to class A GPCRs

G-protein coupled receptors (GPCRs) are important drug targets for various diseases and of major interest to pharmaceutical companies. The function of individual members of this protein family can be modulated by the binding of small molecules at the extracellular side of the structurally conserved transmembrane (TM) domain. Here, we present Snooker, a structure-based approach to generate pharmacophore hypotheses for compounds binding to this extracellular side of the TM domain. Snooker does not require knowledge of ligands, is therefore suitable for apo-proteins, and can be applied to all receptors of the GPCR protein family. The method comprises the construction of a homology model of the TM domains and prioritization of residues on the probability of being ligand binding. Subsequently, protein properties are converted to ligand space, and pharmacophore features are generated at positions where protein ligand interactions are likely. Using this semiautomated knowledge-driven bioinformatics approach we have created pharmacophore hypotheses for 15 different GPCRs from several different subfamilies. For the beta-2-adrenergic receptor we show that ligand poses predicted by Snooker pharmacophore hypotheses reproduce literature supported binding modes for ∼75% of compounds fulfilling pharmacophore constraints. All 15 pharmacophore hypotheses represent interactions with essential residues for ligand binding as observed in mutagenesis experiments and compound selections based on these hypotheses are shown to be target specific. For 8 out of 15 targets enrichment factors above 10-fold are observed in the top 0.5% ranked compounds in a virtual screen. Additionally, prospectively predicted ligand binding poses in the human dopamine D3 receptor based on Snooker pharmacophores were ranked among the best models in the community wide GPCR dock 2010.

Computational systems chemical biology

There is a critical need for improving the level of chemistry awareness in systems biology. The data and information related to modulation of genes and proteins by small molecules continue to accumulate at the same time as simulation tools in systems biology and whole body physiologically based pharmacokinetics (PBPK) continue to evolve. We called this emerging area at the interface between chemical biology and systems biology systems chemical biology (SCB) (Nat Chem Biol 3: 447-450, 2007).The overarching goal of computational SCB is to develop tools for integrated chemical-biological data acquisition, filtering and processing, by taking into account relevant information related to interactions between proteins and small molecules, possible metabolic transformations of small molecules, as well as associated information related to genes, networks, small molecules, and, where applicable, mutants and variants of those proteins. There is yet an unmet need to develop an integrated in silico pharmacology/systems biology continuum that embeds drug-target-clinical outcome (DTCO) triplets, a capability that is vital to the future of chemical biology, pharmacology, and systems biology. Through the development of the SCB approach, scientists will be able to start addressing, in an integrated simulation environment, questions that make the best use of our ever-growing chemical and biological data repositories at the system-wide level. This chapter reviews some of the major research concepts and describes key components that constitute the emerging area of computational systems chemical biology.

Ligand efficiency indices for an effective mapping of chemico-biological space: the concept of an atlas-like representation

We propose a numerical framework that permits an effective atlas-like representation of chemico-biological space based on a series of Cartesian planes mapping the ligands with the corresponding targets connected by an affinity parameter (K(i) or related). The numerical framework is derived from the concept of ligand efficiency indices, which provide a natural coordinate system combining the potency toward the target (biological space) with the physicochemical properties of the ligand (chemical space). This framework facilitates navigation in the multidimensional drug discovery space using map-like representations based on pairs of combined variables related to the efficiency of the ligands per Dalton (molecular weight or number of non-hydrogen atoms) and per unit of polar surface area (or number of polar atoms).

Binding mechanisms in supramolecular complexes

Supramolecular chemistry has expanded dramatically in recent years both in terms of potential applications and in its relevance to analogous biological systems. The formation and function of supramolecular complexes occur through a multiplicity of often difficult to differentiate noncovalent forces. The aim of this Review is to describe the crucial interaction mechanisms in context, and thus classify the entire subject. In most cases, organic host-guest complexes have been selected as examples, but biologically relevant problems are also considered. An understanding and quantification of intermolecular interactions is of importance both for the rational planning of new supramolecular systems, including intelligent materials, as well as for developing new biologically active agents.

Host-guest complexes with protein-ligand-like affinities: computational analysis and design

It has recently been discovered that guests combining a nonpolar core with cationic substituents bind cucurbit[7]uril (CB[7]) in water with ultrahigh affinities. The present study uses the Mining Minima algorithm to study the physics of these extraordinary associations and to computationally test a new series of CB[7] ligands designed to bind with similarly high affinity. The calculations reproduce key experimental observations regarding the affinities of ferrocene-based guests with CB[7] and beta-cyclodextrin and provide a coherent view of the roles of electrostatics and configurational entropy as determinants of affinity in these systems. The newly designed series of compounds is based on a bicyclo[2.2.2]octane core, which is similar in size and polarity to the ferrocene core of the existing series. Mining Minima predicts that these new compounds will, like the ferrocenes, bind CB[7] with extremely high affinities.

Protein structure analysis online

Computational biology/chemistry tools are used in most areas of life/health science research. These methods are continually being developed and their use can present difficulties for both experienced and novice investigators. To facilitate the use of these applications, many packages have been implemented online during these last 5 years. This unit focuses on online computational methods with a special emphasis on structural refinement/atomic simulations, protein electrostatic calculations, searches for functional sites, searches for druggable pockets, protein docking and small molecule docking, and prediction of potential impact of amino acid variations on the structure and function of the protein molecules.

BindingDB: a web-accessible database of experimentally determined protein-ligand binding affinities

BindingDB () is a publicly accessible database currently containing ∼20 000 experimentally determined binding affinities of protein–ligand complexes, for 110 protein targets including isoforms and mutational variants, and ∼11 000 small molecule ligands. The data are extracted from the scientific literature, data collection focusing on proteins that are drug-targets or candidate drug-targets and for which structural data are present in the Protein Data Bank. The BindingDB website supports a range of query types, including searches by chemical structure, substructure and similarity; protein sequence; ligand and protein names; affinity ranges and molecular weight. Data sets generated by BindingDB queries can be downloaded in the form of annotated SDfiles for further analysis, or used as the basis for virtual screening of a compound database uploaded by the user. The data in BindingDB are linked both to structural data in the PDB via PDB IDs and chemical and sequence searches, and to the literature in PubMed via PubMed IDs.

Internet resources related to drug action and human response: a review

It has been demonstrated that numerous proteins interact with drugs or their metabolites. Knowledge of these proteins is necessary to understand the mechanisms of drug action and human response. Progress in modern genetics, molecular biology, biochemistry and pharmacology is generating a comprehensive mechanistic understanding of drug-target interaction on the molecular level. This is valuable for researchers and pharmaceutical companies in their efforts to improve the efficacy of existing drugs and to discover new ones. Most recently, the integration of a systems biology approach into drug discovery processes calls for more holistic knowledge and easily accessible resources of the proteins that are important in drug action and human response. We have reviewed many publicly accessible internet resources of these proteins, according to their roles in drug action and human response, such as therapeutic effect, adverse reaction, absorption, distribution, metabolism and excretion.

Cataloging the relationships between proteins: a review of interaction databases

By organizing and making widely accessible the increasing amounts of data from high-throughput analyses, protein interaction databases have become an integral resource for the biological community in relating sequence data with higher-order function. To provide a sense of the use and applicability of these databases, we describe each of the major comprehensive interaction databases as well as some of the more specialized ones. Content description, search/browse functionalities, and data presentation are discussed. A succinct explanation of database contents helps the user quickly identify whether the database contains applicable information to their research interest. Broad levels of search/browse functions as well as descriptions/examples allow users to quickly find and access pertinent data. At this point, clear presentation of search results as well as the primary content is necessary. Many databases display information graphically or divided into smaller digestible parts over a number of tabbed/linked pages. In addition, cross-linking between the databases promotes interconnectivity of the data and is an added layer of relational data for the user. Overall, although these protein interaction databases are under continual improvement, their current state shows that much time and effort has gone into organizing and presenting these large sets of data-describing protein interactions.

How good is your screening library?

Efficient library design is an ongoing challenge for investigators seeking novel ligands for proteins, whether for drug discovery or chemical biology. Strategies that add neglected chemistry or exclude unproductive compounds are two dominant recent themes, as is a growing awareness of molecular complexity and its implications. The choice of how complex molecules in screening libraries should be often amounts to how big they should be. Small, simple molecules have lower affinities and must be screened at high concentration, but they will also have higher hit rates. Larger compounds, on the other hand, will often more closely resemble final drugs, but because they are more highly functionalized and specific, they will have much lower hit rates. The best general-purpose screening libraries may well be those of intermediate complexity that are free of artifact-causing nuisance compounds.

Binding MOAD (Mother Of All Databases)

Binding MOAD (Mother of All Databases) is the largest collection of high-quality, protein-ligand complexes available from the Protein Data Bank. At this time, Binding MOAD contains 5331 protein-ligand complexes comprised of 1780 unique protein families and 2630 unique ligands. We have searched the crystallography papers for all 5000+ structures and compiled binding data for 1375 (26%) of the protein-ligand complexes. The binding-affinity data ranges 13 orders of magnitude. This is the largest collection of binding data reported to date in the literature. We have also addressed the issue of redundancy in the data. To create a nonredundant dataset, one protein from each of the 1780 protein families was chosen as a representative. Representatives were chosen by tightest binding, best resolution, etc. For the 1780 "best" complexes that comprise the nonredundant version of Binding MOAD, 475 (27%) have binding data. This significant collection of protein-ligand complexes will be very useful in elucidating the biophysical patterns of molecular recognition and enzymatic regulation. The complexes with binding-affinity data will help in the development of improved scoring functions and structure-based drug discovery techniques. The dataset can be accessed at http://www.BindingMOAD.org.

A survey of the year 2002 literature on applications of isothermal titration calorimetry

Isothermal titration calorimetry (ITC) is becoming widely accepted as a key instrument in any laboratory in which quantification of biomolecular interactions is a requisite. The method has matured with respect to general acceptance and application development over recent years. The number of publications on ITC has grown exponentially over the last 10 years, reflecting the general utility of the method. Here all the published works of the year 2002 in this area have been surveyed. We review the broad range of systems to which ITC is being directed and classify these into general areas highlighting key publications of interest. This provides an overview of what can be achieved using this method and what developments are likely to occur in the near future.

Enzyme classification by ligand binding

The problem of assigning a biochemical function to newly discovered proteins has been traditionally approached by expert enzymological analysis, sequence analysis, and structural modeling. In recent years, the appearance of databases containing protein-ligand interaction data for large numbers of protein classes and chemical compounds have provided new ways of investigating proteins for which the biochemical function is not completely understood. In this work, we introduce a method that utilizes ligand-binding data for functional classification of enzymes. The method makes use of the existing Enzyme Commission (EC) classification scheme and the data on interactions of small molecules with enzymes from the BRENDA database. A set of ligands that binds to an enzyme with unknown biochemical function serves as a query to search a protein-ligand interaction database for enzyme classes that are known to interact with a similar set of ligands. These classes provide hypotheses of the query enzyme's function and complement other computational annotations that take advantage of sequence and structural information. Similarity between sets of ligands is computed using point set similarity measures based upon similarity between individual compounds. We present the statistics of classification of the enzymes in the database by a cross-validation procedure and illustrate the application of the method on several examples.

The evolution of structural databases

Starting with the Protein Data Bank (PDB) as a common ancestor, the evolution of structural databases has been driven by the rapprochement of the structural world and the practical applications. The result is an impressive number of secondary structural databases that is welcomed by structural biologists and bioinformaticians but runs the risk of producing an embarrassment of riches among non-specialist users. Given that any profit depends on the number of customers, efficient interfaces between many structural data banks must be available to make their contents easily accessible. Increasing the information content of central structural repositories might be the best way to guide users through the many, sometimes overlapping databases.