Similarity searching using 2D structural fingerprints

DeepCompoundNet: enhancing compound-protein interaction prediction with multimodal convolutional neural networks

Virtual screening has emerged as a valuable computational tool for predicting compound-protein interactions, offering a cost-effective and rapid approach to identifying potential candidate drug molecules. Current machine learning-based methods rely on molecular structures and their relationship in the network. The former utilizes information such as amino acid sequences and chemical structures, while the latter leverages interaction network data, such as protein-protein interactions, drug-disease interactions, and protein-disease interactions. However, there has been limited exploration of integrating molecular information with interaction networks. This study presents DeepCompoundNet, a deep learning-based model that integrates protein features, drug properties, and diverse interaction data to predict chemical-protein interactions. DeepCompoundNet outperforms state-of-the-art methods for compound-protein interaction prediction, as demonstrated through performance evaluations. Our findings highlight the complementary nature of multiple interaction data, extending beyond amino acid sequence homology and chemical structure similarity. Moreover, our model's analysis confirms that DeepCompoundNet gets higher performance in predicting interactions between proteins and chemicals not observed in the training samples.Communicated by Ramaswamy H. Sarma.

Accelerating drug discovery and repurposing by combining transcriptional signature connectivity with docking

We present an in silico approach for drug discovery, dubbed connectivity enhanced structure activity relationship (ceSAR). Building on the landmark LINCS library of transcriptional signatures of drug-like molecules and gene knockdowns, ceSAR combines cheminformatic techniques with signature concordance analysis to connect small molecules and their targets and further assess their biophysical compatibility using molecular docking. Candidate compounds are first ranked in a target structure-independent manner, using chemical similarity to LINCS analogs that exhibit transcriptomic concordance with a target gene knockdown. Top candidates are subsequently rescored using docking simulations and machine learning-based consensus of the two approaches. Using extensive benchmarking, we show that ceSAR greatly reduces false-positive rates, while cutting run times by multiple orders of magnitude and further democratizing drug discovery pipelines. We further demonstrate the utility of ceSAR by identifying and experimentally validating inhibitors of BCL2A1, an important antiapoptotic target in melanoma and preterm birth-associated inflammation.

Do Molecular Fingerprints Identify Diverse Active Drugs in Large-Scale Virtual Screening? (No)

Computational approaches for small-molecule drug discovery now regularly scale to the consideration of libraries containing billions of candidate small molecules. One promising approach to increased the speed of evaluating billion-molecule libraries is to develop succinct representations of each molecule that enable the rapid identification of molecules with similar properties. Molecular fingerprints are thought to provide a mechanism for producing such representations. Here, we explore the utility of commonly used fingerprints in the context of predicting similar molecular activity. We show that fingerprint similarity provides little discriminative power between active and inactive molecules for a target protein based on a known active-while they may sometimes provide some enrichment for active molecules in a drug screen, a screened data set will still be dominated by inactive molecules. We also demonstrate that high-similarity actives appear to share a scaffold with the query active, meaning that they could more easily be identified by structural enumeration. Furthermore, even when limited to only active molecules, fingerprint similarity values do not correlate with compound potency. In sum, these results highlight the need for a new wave of molecular representations that will improve the capacity to detect biologically active molecules based on their similarity to other such molecules.

On the relevance of query definition in the performance of 3D ligand-based virtual screening

Ligand-based virtual screening (LBVS) methods are widely used to explore the vast chemical space in the search of novel compounds resorting to a variety of properties encoded in 1D, 2D or 3D descriptors. The success of 3D-LBVS is affected by the overlay of molecular pairs, thus making selection of the template compound, search of accessible conformational space and choice of the query conformation to be potential factors that modulate the successful retrieval of actives. This study examines the impact of adopting different choices for the query conformation of the template, paying also attention to the influence exerted by the structural similarity between templates and actives. The analysis is performed using PharmScreen, a 3D LBVS tool that relies on similarity measurements of the hydrophobic/philic pattern of molecules, and Phase Shape, which is based on the alignment of atom triplets followed by refinement of the volume overlap. The study is performed for the original DUD-E+ database and a Morgan Fingerprint filtered version (denoted DUD-E+-Diverse; available in https://github.com/Pharmacelera/Query-models-to-3DLBVS ), which was prepared to minimize the 2D resemblance between template and actives. Although in most cases the query conformation exhibits a mild influence on the overall performance, a critical analysis is made to disclose factors, such as the content of structural features between template and actives and the induction of conformational strain in the template, that underlie the drastic impact of the query definition in the recovery of actives for certain targets. The findings of this research also provide valuable guidance for assisting the selection of the query definition in 3D LBVS campaigns.

Metabolic Connectome and Its Role in the Prediction, Diagnosis, and Treatment of Complex Diseases

The interconnectivity of advanced biological systems is essential for their proper functioning. In modern connectomics, biological entities such as proteins, genes, RNA, DNA, and metabolites are often represented as nodes, while the physical, biochemical, or functional interactions between them are represented as edges. Among these entities, metabolites are particularly significant as they exhibit a closer relationship to an organism's phenotype compared to genes or proteins. Moreover, the metabolome has the ability to amplify small proteomic and transcriptomic changes, even those from minor genomic changes. Metabolic networks, which consist of complex systems comprising hundreds of metabolites and their interactions, play a critical role in biological research by mediating energy conversion and chemical reactions within cells. This review provides an introduction to common metabolic network models and their construction methods. It also explores the diverse applications of metabolic networks in elucidating disease mechanisms, predicting and diagnosing diseases, and facilitating drug development. Additionally, it discusses potential future directions for research in metabolic networks. Ultimately, this review serves as a valuable reference for researchers interested in metabolic network modeling, analysis, and their applications.

Chemical similarity of molecules with physiological response

Measuring the similarity among molecules is an important task in various chemically oriented problems. This elusive concept is hard to define and quantify. Moreover, the complexity of the problem is elevated by bifurcating the notion of molecular similarity to structural and chemical similarity. While the structural similarity of molecules is being extensively researched, the so-called chemical similarity is being mentioned scarcely. Here, we propose a way of converting the physico-chemical properties into molecular fingerprints. Then, using the apparatus of measuring the structural similarity, the chemical similarity can be assessed. The proof of a concept is demonstrated on a set of molecules that induce diverse physiological responses.

PrePCI: A structure- and chemical similarity-informed database of predicted protein compound interactions

We describe the Predicting Protein-Compound Interactions (PrePCI) database which comprises over 5 billion predicted interactions between 6.8 million chemical compounds and 19,797 human proteins. PrePCI relies on a proteome-wide database of structural models based on both traditional modeling techniques and the AlphaFold Protein Structure Database. Sequence- and structural similarity-based metrics are established between template proteins, T, in the Protein Data Bank that bind compounds, C, and query proteins in the model database, Q. When the metrics exceed threshold values, it is assumed that C also binds to Q with a likelihood ratio (LR) derived from machine learning. If the relationship is based on structural similarity, the LR is based on a scoring function that measures the extent to which C is compatible with the binding site of Q as described in the LT-scanner algorithm. For every predicted complex derived in this way, chemical similarity based on the Tanimoto coefficient identifies other small molecules that may bind to Q. An overall LR for the binding of C to Q is obtained from Naive Bayesian statistics. The PrePCI database can be queried by entering a UniProt ID or gene name for a protein to obtain a list of compounds predicted to bind to it along with associated LRs. Alternatively, entering an identifier for the compound outputs a list of proteins it is predicted to bind. Specific applications of the database to lead discovery, elucidation of drug mechanism of action, and biological function annotation are described.

Exploration of a Large Virtual Chemical Space: Identification of Potent Inhibitors of Lactate Dehydrogenase-A against Pancreatic Cancer

It is imperative to explore the gigantic available chemical space to identify new scaffolds for drug lead discovery. Identifying potent hits from virtual screening of large chemical databases is challenging and computationally demanding. Rather than the traditional two-dimensional (2D)/three-dimensional (3D) approaches on smaller chemical libraries of a few hundred thousand compounds, we screened a ZINC library of 15 million compounds using multiple computational methods. Here, we present the successful application of a virtual screening methodology that identifies several chemotypes as starting hits against lactate dehydrogenase-A (LDHA). From 29 compounds identified from virtual screening, 17 (58%) showed IC₅₀ values < 63 μM, two showed single-digit micromolar inhibition, and the most potent hit compound had IC₅₀ down to 117 nM. We enriched the database and employed an ensemble approach by combining 2D fingerprint similarity searches, pharmacophore modeling, molecular docking, and molecular dynamics. WaterMap calculations were carried out to explore the thermodynamics of surface water molecules and gain insights into the LDHA binding pocket. The present work has led to the discovery of two new chemical classes, including compounds with a succinic acid monoamide moiety or a hydroxy pyrimidinone ring system. Selected hits block lactate production in cells and inhibit pancreatic cancer cell lines with cytotoxicity IC₅₀ down to 12.26 μM against MIAPaCa-2 cells and 14.64 μM against PANC-1, which, under normoxic conditions, is already comparable or more potent than most currently available known LDHA inhibitors.

The Discovery of Potential SARS-CoV-2 Natural Inhibitors among 4924 African Metabolites Targeting the Papain-like Protease: A Multi-Phase In Silico Approach

Four compounds, hippacine, 4,2'-dihydroxy-4'-methoxychalcone, 2',5'-dihydroxy-4-methoxychalcone, and wighteone, were selected from 4924 African natural metabolites as potential inhibitors against SARS-CoV-2 papain-like protease (PLpro, PDB ID: 3E9S). A multi-phased in silico approach was employed to select the most similar metabolites to the co-crystallized ligand (TTT) of the PLpro through molecular fingerprints and structural similarity studies. Followingly, to examine the binding of the selected metabolites with the PLpro (molecular docking. Further, to confirm this binding through molecular dynamics simulations. Finally, in silico ADMET and toxicity studies were carried out to prefer the most convenient compounds and their drug-likeness. The obtained results could be a weapon in the battle against COVID-19 via more in vitro and in vivo studies.

Structure-Based Virtual Screening, Docking, ADMET, Molecular Dynamics, and MM-PBSA Calculations for the Discovery of Potential Natural SARS-CoV-2 Helicase Inhibitors from the Traditional Chinese Medicine

Continuing our antecedent work against COVID-19, a set of 5956 compounds of traditional Chinese medicine have been virtually screened for their potential against SARS-CoV-2 helicase (PDB ID: 5RMM). Initially, a fingerprint study with VXG, the ligand of the target enzyme, disclosed the similarity of 187 compounds. Then, a molecular similarity study declared the most similar 40 compounds. Subsequently, molecular docking studies were carried out to examine the binding modes and energies. Then, the most appropriate 26 compounds were subjected to in silico ADMET and toxicity studies to select the most convenient inhibitors to be: (1R,2S)-ephedrine (57), (1R,2S)-norephedrine (59), 2-(4-(pyrrolidin-1-yl)phenyl)acetic acid (84), 1-phenylpropane-1,2-dione (195), 2-methoxycinnamic acid (246), 2-methoxybenzoic acid (364), (R)-2-((R)-5-oxopyrrolidin-3-yl)-2-phenylacetic acid (405), (Z)-6-(3-hydroxy-4-methoxystyryl)-4-methoxy-2H-pyran-2-one (533), 8-chloro-2-(2-phenylethyl)-5,6,7-trihydroxy-5,6,7,8-tetrahydrochromone (637), 3-((1R,2S)-2-(dimethylamino)-1-hydroxypropyl)phenol (818), (R)-2-ethyl-4-(1-hydroxy-2-(methylamino)ethyl)phenol (5159), and (R)-2-((1S,2S,5S)-2-benzyl-5-hydroxy-4-methylcyclohex-3-en-1-yl)propane-1,2-diol (5168). Among the selected 12 compounds, the metabolites, compound 533 showed the best docking scores. Interestingly, the MD simulation studies for compound 533, the one with the highest docking score, over 100 ns showed its correct binding to SARS-CoV-2 helicase with low energy and optimum dynamics. Finally, MM-PBSA studies showed that 533 bonded favorably to SARS-CoV-2 helicase with a free energy value of −83 kJ/mol. Further, the free energy decomposition study determined the essential amino acid residues that contributed favorably to the binding process. The obtained results give a huge hope to find a cure for COVID-19 through further in vitro and in vivo studies for the selected compounds.

Artificial intelligence and machine-learning approaches in structure and ligand-based discovery of drugs affecting central nervous system

CNS disorders are indications with a very high unmet medical needs, relatively smaller number of available drugs, and a subpar satisfaction level among patients and caregiver. Discovery of CNS drugs is extremely expensive affair with its own unique challenges leading to extremely high attrition rates and low efficiency. With explosion of data in information age, there is hardly any aspect of life that has not been touched by data driven technologies such as artificial intelligence (AI) and machine learning (ML). Drug discovery is no exception, emergence of big data via genomic, proteomic, biological, and chemical technologies has driven pharmaceutical giants to collaborate with AI oriented companies to revolutionise drug discovery, with the goal of increasing the efficiency of the process. In recent years many examples of innovative applications of AI and ML techniques in CNS drug discovery has been reported. Research on therapeutics for diseases such as schizophrenia, Alzheimer's and Parkinsonism has been provided with a new direction and thrust from these developments. AI and ML has been applied to both ligand-based and structure-based drug discovery and design of CNS therapeutics. In this review, we have summarised the general aspects of AI and ML from the perspective of drug discovery followed by a comprehensive coverage of the recent developments in the applications of AI/ML techniques in CNS drug discovery.

Rings in Clinical Trials and Drugs: Present and Future

We present a comprehensive analysis of all ring systems (both heterocyclic and nonheterocyclic) in clinical trial compounds and FDA-approved drugs. We show 67% of small molecules in clinical trials comprise only ring systems found in marketed drugs, which mirrors previously published findings for newly approved drugs. We also show there are approximately 450 000 unique ring systems derived from 2.24 billion molecules currently available in synthesized chemical space, and molecules in clinical trials utilize only 0.1% of this available pool. Moreover, there are fewer ring systems in drugs compared with those in clinical trials, but this is balanced by the drug ring systems being reused more often. Furthermore, systematic changes of up to two atoms on existing drug and clinical trial ring systems give a set of 3902 future clinical trial ring systems, which are predicted to cover approximately 50% of the novel ring systems entering clinical trials.

Computer-Aided Methods for Molecular Classification

The study aims to analyze the degree of similarity of some molecules belonging to two subgroups of Aminoalkylindoles. After extracting the molecules’ characteristics using Cheminformatics methods, and the computation of the Tanimoto coefficients, dendrograms and heatmaps were built to reveal the degree of similarity of the analyzed drugs. Some atom-pair similarities between the molecules in the same group were detected. The clusters determined by the k-means method divided the Benzoylindoles into two subgroups but kept all the Phenylacetylindoles together in the same set. The activity spectrum of the elements in each group was also analyzed, and similarities have been emphasized. The clustering has been validated using the Kruskal–Wallis test on the series of computed probabilities of the main effects.

Similarity-Based Virtual Screen Using Enhanced Siamese Deep Learning Methods

Traditional drug production is a long and complex process that leads to new drug production. The virtual screening technique is a computational method that allows chemical compounds to be screened at an acceptable time and cost. Several databases contain information on various aspects of biologically active substances. Simple statistical tools are difficult to use because of the enormous amount of information and complex data samples of molecules that are structurally heterogeneous recorded in these databases. Many techniques for capturing the biological similarity between a test compound and a known target ligand in LBVS have been established. However, despite the good performances of the above methods compared to their prior, especially when dealing with molecules that have homogeneous active structural elements, they are not satisfied when dealing with molecules that are structurally heterogeneous. Deep learning models have recently achieved considerable success in a variety of disciplines due to their powerful generalization and feature extraction capabilities. Also, the Siamese network has been used in similarity models for more complicated data samples, especially with heterogeneous data samples. The main aim of this study is to enhance the performance of similarity searching, especially with molecules that are structurally heterogeneous. The Siamese architecture will be enhanced using two similarity distance layers with one fusion layer to further improve the similarity measurements between molecules and then adding many layers after the fusion layer for some models to improve the retrieval recall. In this architecture, several methods of deep learning have been used, which are long short-term memory (LSTM), gated recurrent unit (GRU), convolutional neural network-one dimension (CNN1D), and convolutional neural network-two dimensions (CNN2D). A series of experiments have been carried out on real-world data sets, and the results have shown that the proposed methods outperformed the existing methods.

Comparing the performance and coverage of selected in silico (liver) metabolism tools relative to reported studies in the literature to inform analogue selection in read-across: A case study

Changes in the regulatory landscape of chemical safety assessment call for the use of New Approach Methodologies (NAMs) including read-across to fill data gaps. One critical aspect of analogue evaluation is the extent to which target and source analogues are metabolically similar. In this study, a set of 37 structurally diverse chemicals were compiled from the EPA ToxCast inventory to compare and contrast a selection of metabolism in silico tools, in terms of their coverage and performance relative to metabolism information reported in the literature. The aim was to build understanding of the scope and capabilities of these tools and how they could be utilised in a read-across assessment. The tools were Systematic Generation of Metabolites (SyGMa), Meteor Nexus, BioTransformer, Tissue Metabolism Simulator (TIMES), OECD Toolbox, and Chemical Transformation Simulator (CTS). Performance was characterised by sensitivity and precision determined by comparing predictions against literature reported metabolites (from 44 publications). A coverage score was derived to provide a relative quantitative comparison between the tools. Meteor, TIMES, Toolbox, and CTS predictions were run in batch mode, using default settings. SyGMa and BioTransformer were run with user-defined settings, (two passes of phase I and one pass of phase II). Hierarchical clustering revealed high similarity between TIMES and Toolbox. SyGMa had the highest coverage, matching an average of 38.63% of predictions generated by the other tools though was prone to significant overprediction. It generated 5,125 metabolites, which represented 54.67% of all predictions. Precision and sensitivity values ranged from 1.1-29% and 14.7-28.3% respectively. The Toolbox had the highest performance overall. A case study was presented for 3,4-Toluenediamine (3,4-TDA), assessed for the derivation of screening-level Provisional Peer Reviewed Toxicity Values (PPRTVs), was used to demonstrate the practical role in silico metabolism information can play in analogue evaluation as part of a read-across approach.

Computer-Aided Classification of New Psychoactive Substances

The appearance on the free market of synthetic cannabinoids raised the researchers’ interest in establishing their molecular similarity by QSAR analysis. A rigorous criterion for classifying drugs is their chemical structure. Therefore, this article presents the structural similarity of two groups of drugs: benzoylindoles and phenylacetylindoles. Statistical analysis and clustering of the molecules are performed based on their numerical characteristics extracted using Cheminformatics methods. Their similarities/dissimilarities are emphasized using the dendrograms and heat map. The highest discrepancies are found in the phenylacetylindoles group.

Rapid Identification of Potential Drug Candidates from Multi-Million Compounds' Repositories. Combination of 2D Similarity Search with 3D Ligand/Structure Based Methods and In Vitro Screening

Rapid in silico selection of target focused libraries from commercial repositories is an attractive and cost-effective approach in early drug discovery. If structures of active compounds are available, rapid 2D similarity search can be performed on multimillion compounds' databases. This approach can be combined with physico-chemical parameter and diversity filtering, bioisosteric replacements, and fragment-based approaches for performing a first round biological screening. Our objectives were to investigate the combination of 2D similarity search with various 3D ligand and structure-based methods for hit expansion and validation, in order to increase the hit rate and novelty. In the present account, six case studies are described and the efficiency of mixing is evaluated. While sequentially combined 2D/3D similarity approach increases the hit rate significantly, sequential combination of 2D similarity with pharmacophore model or 3D docking enriched the resulting focused library with novel chemotypes. Parallel integrated approaches allowed the comparison of the various 2D and 3D methods and revealed that 2D similarity-based and 3D ligand and structure-based techniques are often complementary, and their combinations represent a powerful synergy. Finally, the lessons we learnt including the advantages and pitfalls of the described approaches are discussed.

Improving Measures of Chemical Structural Similarity Using Machine Learning on Chemical-Genetic Interactions

A common strategy for identifying molecules likely to possess a desired biological activity is to search large databases of compounds for high structural similarity to a query molecule that demonstrates this activity, under the assumption that structural similarity is predictive of similar biological activity. However, efforts to systematically benchmark the diverse array of available molecular fingerprints and similarity coefficients have been limited by a lack of large-scale datasets that reflect biological similarities of compounds. To elucidate the relative performance of these alternatives, we systematically benchmarked 11 different molecular fingerprint encodings, each combined with 13 different similarity coefficients, using a large set of chemical–genetic interaction data from the yeast Saccharomyces cerevisiae as a systematic proxy for biological activity. We found that the performance of different molecular fingerprints and similarity coefficients varied substantially and that the all-shortest path fingerprints paired with the Braun-Blanquet similarity coefficient provided superior performance that was robust across several compound collections. We further proposed a machine learning pipeline based on support vector machines that offered a fivefold improvement relative to the best unsupervised approach. Our results generally suggest that using high-dimensional chemical–genetic data as a basis for refining molecular fingerprints can be a powerful approach for improving prediction of biological functions from chemical structures.

Potential Biological Targets Prediction, ADME Profiling, & Molecular Docking studies of Novel Steroidal Products from Cunninghamella Blakesleana

BACKGROUND

New potential biological targets prediction through inverse molecular docking technique is an another smart strategy to forecast the possibility of compounds being biologically active against various target receptors.

OBJECTIVES

In this case of designed study, we screened our recently obtained novel acetylinic steroidal biotransformed products [(1) 8-β-methyl-14-α-hydroxy∆4tibolone (2) 9-α-Hydroxy∆4 tibolone (3) 8-β-methyl-11-β-hydroxy∆4tibolone (4) 6-β-hydroxy∆4tibolone, (5) 6-β-9-α-dihydroxy∆4tibolone (6) 7-β-hydroxy∆4tibolone) ] from fungi Cunninghemella Blakesleana to predict their possible biological targets and profiling of ADME properties.

METHOD

The prediction of pharmacokinetics properties membrane permeability as well as bioavailability radar properties were carried out by using Swiss target prediction, and Swiss ADME tools, respectively these metabolites were also subjected to predict the possible mechanism of action along with associated biological network pathways by using Reactome data-base.

RESULTS

All the six screened compounds possess excellent drug ability criteria, and exhibited exceptionally excellent non inhibitory potential against all five isozymes of CYP450 enzyme complex, including (CYP1A2, CYP2C19, CYP2C9, CYP2D6, and CYP3A4) respectively. All the screened compounds are lying within the acceptable pink zone of bioavailability radar and showing excellent descriptive properties. Compounds [1-4 & 6] are showing high BBB (Blood Brain Barrier) permeation, while compound 5 is exhibiting high HIA (Human Intestinal Absorption) property of (Egan Egg).

CONCLUSION

In conclusion, the results of this study smartly reveals that in-silico based studies are considered to provide robustness towards a rational drug designing and development approach, therefore in this way it helps to avoid the possibility of failure of drug candidates in the later experimental stages of drug development phases.

Resources and computational strategies to advance small molecule SARS-CoV-2 discovery: Lessons from the pandemic and preparing for future health crises

There is an urgent need to identify new therapies that prevent SARS-CoV-2 infection and improve the outcome of COVID-19 patients. This pandemic has thus spurred intensive research in most scientific areas and in a short period of time, several vaccines have been developed. But, while the race to find vaccines for COVID-19 has dominated the headlines, other types of therapeutic agents are being developed. In this mini-review, we report several databases and online tools that could assist the discovery of anti-SARS-CoV-2 small chemical compounds and peptides. We then give examples of studies that combined in silico and in vitro screening, either for drug repositioning purposes or to search for novel bioactive compounds. Finally, we question the overall lack of discussion and plan observed in academic research in many countries during this crisis and suggest that there is room for improvement.

Integrating 3D structural information into systems biology

Systems biology is a data-heavy field that focuses on systems-wide depictions of biological phenomena necessarily sacrificing a detailed characterization of individual components. As an example, genome-wide protein interaction networks are widely used in systems biology and continuously extended and refined as new sources of evidence become available. Despite the vast amount of information about individual protein structures and protein complexes that has accumulated in the past 50 years in the Protein Data Bank, the data, computational tools, and language of structural biology are not an integral part of systems biology. However, increasing effort has been devoted to this integration, and the related literature is reviewed here. Relationships between proteins that are detected via structural similarity offer a rich source of information not available from sequence similarity, and homology modeling can be used to leverage Protein Data Bank structures to produce 3D models for a significant fraction of many proteomes. A number of structure-informed genomic and cross-species (i.e., virus–host) interactomes will be described, and the unique information they provide will be illustrated with a number of examples. Tissue- and tumor-specific interactomes have also been developed through computational strategies that exploit patient information and through genetic interactions available from increasingly sensitive screens. Strategies to integrate structural information with these alternate data sources will be described. Finally, efforts to link protein structure space with chemical compound space offer novel sources of information in drug design, off-target identification, and the identification of targets for compounds found to be effective in phenotypic screens.

A Quantitative Survey of Bacterial Persistence in the Presence of Antibiotics: Towards Antipersister Antimicrobial Discovery

Background: Bacterial persistence to antibiotics relates to the phenotypic ability to survive lethal concentrations of otherwise bactericidal antibiotics. The quantitative nature of the time-kill assay, which is the sector's standard for the study of antibiotic bacterial persistence, is an invaluable asset for global, unbiased, and cross-species analyses. Methods: We compiled the results of antibiotic persistence from antibiotic-sensitive bacteria during planktonic growth. The data were extracted from a sample of 187 publications over the last 50 years. The antibiotics used in this compilation were also compared in terms of structural similarity to fluorescent molecules known to accumulate in Escherichia coli. Results: We reviewed in detail data from 54 antibiotics and 36 bacterial species. Persistence varies widely as a function of the type of antibiotic (membrane-active antibiotics admit the fewest), the nature of the growth phase and medium (persistence is less common in exponential phase and rich media), and the Gram staining of the target organism (persistence is more common in Gram positives). Some antibiotics bear strong structural similarity to fluorophores known to be taken up by E. coli, potentially allowing competitive assays. Some antibiotics also, paradoxically, seem to allow more persisters at higher antibiotic concentrations. Conclusions: We consolidated an actionable knowledge base to support a rational development of antipersister antimicrobials. Persistence is seen as a step on the pathway to antimicrobial resistance, and we found no organisms that failed to exhibit it. Novel antibiotics need to have antipersister activity. Discovery strategies should include persister-specific approaches that could find antibiotics that preferably target the membrane structure and permeability of slow-growing cells.

VAE-Sim: A Novel Molecular Similarity Measure Based on a Variational Autoencoder

Molecular similarity is an elusive but core "unsupervised" cheminformatics concept, yet different "fingerprint" encodings of molecular structures return very different similarity values, even when using the same similarity metric. Each encoding may be of value when applied to other problems with objective or target functions, implying that a priori none are "better" than the others, nor than encoding-free metrics such as maximum common substructure (MCSS). We here introduce a novel approach to molecular similarity, in the form of a variational autoencoder (VAE). This learns the joint distribution p(z|x) where z is a latent vector and x are the (same) input/output data. It takes the form of a "bowtie"-shaped artificial neural network. In the middle is a "bottleneck layer" or latent vector in which inputs are transformed into, and represented as, a vector of numbers (encoding), with a reverse process (decoding) seeking to return the SMILES string that was the input. We train a VAE on over six million druglike molecules and natural products (including over one million in the final holdout set). The VAE vector distances provide a rapid and novel metric for molecular similarity that is both easily and rapidly calculated. We describe the method and its application to a typical similarity problem in cheminformatics.

Virtual screening web servers: designing chemical probes and drug candidates in the cyberspace

The interplay between life sciences and advancing technology drives a continuous cycle of chemical data growth; these data are most often stored in open or partially open databases. In parallel, many different types of algorithms are being developed to manipulate these chemical objects and associated bioactivity data. Virtual screening methods are among the most popular computational approaches in pharmaceutical research. Today, user-friendly web-based tools are available to help scientists perform virtual screening experiments. This article provides an overview of internet resources enabling and supporting chemical biology and early drug discovery with a main emphasis on web servers dedicated to virtual ligand screening and small-molecule docking. This survey first introduces some key concepts and then presents recent and easily accessible virtual screening and related target-fishing tools as well as briefly discusses case studies enabled by some of these web services. Notwithstanding further improvements, already available web-based tools not only contribute to the design of bioactive molecules and assist drug repositioning but also help to generate new ideas and explore different hypotheses in a timely fashion while contributing to teaching in the field of drug development.

The rcdk and cluster R packages applied to drug candidate selection

The aim of this article is to show how thevpower of statistics and cheminformatics can be combined, in R, using two packages: rcdk and cluster.We describe the role of clustering methods for identifying similar structures in a group of 23 molecules according to their fingerprints. The most commonly used method is to group the molecules using a "score" obtained by measuring the average distance between them. This score reflects the similarity/non-similarity between compounds and helps us identify active or potentially toxic substances through predictive studies.Clustering is the process by which the common characteristics of a particular class of compounds are identified. For clustering applications, we are generally measure the molecular fingerprint similarity with the Tanimoto coefficient. Based on the molecular fingerprints, we calculated the molecular distances between the methotrexate molecule and the other 23 molecules in the group, and organized them into a matrix. According to the molecular distances and Ward 's method, the molecules were grouped into 3 clusters. We can presume structural similarity between the compounds and their locations in the cluster map. Because only 5 molecules were included in the methotrexate cluster, we considered that they might have similar properties and might be further tested as potential drug candidates.

Promising Terpenes as Natural Antagonists of Cancer: An In-Silico Approach

Overexpression of murine double minute 2 (MDM2) results in the inactivation of p53 and causes cancer which is a leading cause of death in recent era. In recent decades, much attention has been paid to discover potential inhibitors against MDM2 in order to cure cancer. Outcomes from studies proposes that the MDM2 is a hot target to screen potent antagonists. Thus, this study aims at discovering natural compounds using several computational approaches to inhibit the MDM2 and to eliminate p53-MDM2 interaction, which would result in the reactivation of p53 activity. A library of 500 terpenes was prepared and several virtual screening approaches were employed to find out the best hits which could serve as p53-MDM2 antagonists. On the basis of the designed protocol, three terpenes were selected. In the present study, for the stability and validation of selected three protein-ligand complexes 20 ns molecular dynamics simulations and principal component analyses (PCA) were performed. Results found that the selected terpenes hits (3-trans-p-coumaroyl maslinic acid, Silvestrol and Betulonic acid) are potential inhibitors of p53–MDM2 interaction and could serve as potent antagonists.

Design and Selection of Novel C1s Inhibitors by In Silico and In Vitro Approaches

The complement system is associated with various diseases such as inflammation or auto-immune diseases. Complement-targeted drugs could provide novel therapeutic intervention against the above diseases. C1s, a serine protease, plays an important role in the CS and could be an attractive target since it blocks the system at an early stage of the complement cascade. Designing C1 inhibitors is particularly challenging since known inhibitors are restricted to a narrow bioactive chemical space in addition selectivity over other serine proteases is an important requirement. The typical architecture of a small molecule inhibitor of C1s contains an amidine (or guanidine) residue, however, the discovery of non-amidine inhibitors might have high value, particularly if novel chemotypes and/or compounds displaying improved selectivity are identified. We applied various virtual screening approaches to identify C1s focused libraries that lack the amidine/guanidine functionalities, then the in silico generated libraries were evaluated by in vitro biological assays. While 3D structure-based methods were not suitable for virtual screening of C1s inhibitors, and a 2D similarity search did not lead to novel chemotypes, pharmacophore model generation allowed us to identify two novel chemotypes with submicromolar activities. In three screening rounds we tested altogether 89 compounds and identified 20 hit compounds (<10 μM activities; overall hit rate: 22.5%). The highest activity determined was 12 nM (1,2,4-triazole), while for the newly identified chemotypes (1,3-benzoxazin-4-one and thieno[2,3-d][1,3]oxazin-4-one) it was 241 nM and 549 nM, respectively.

An in silico Workflow that Yields Experimentally Comparable Inhibitors for Human Dihydroorotate Dehydrogenase

INTRODUCTION

Rheumatoid Arthritis [RA] is an autoimmune disease that can cause chronic inflammation of the joints. Human DiHydroOrotate DeHydrogenase [DHODH] is a clinically validated drug target for the treatment of Rheumatoid Arthritis. DHODH inhibition results in beneficial immunosuppressant and anti-proliferative effects.

MATERIALS AND METHODS

Leflunomide [LEF] and Brequinar Sodium [BREQ], drugs used in the treatment of RA, suppresses the immune cells responsible for inflammation but has several side-effects, most predominant being symptomatic liver damage and toxicity. An existing scaffold based on structural analogies with LEF and BREQ was used to screen out potent inhibitors of DHODH, in ZINC Database using 2D binary fingerprint. 10 structures similar to the scaffold were shortlisted due to their Tanimoto similarity coefficient. Selected structures were docked using the tools AutoDock, Ligand fit and iGEMDOCK with target human DHODH. High scoring compounds having similar interactions as that of scaffold were checked to evaluate their Drug-Likeliness.

RESULTS

The five shortlisted compounds were then subjected to Molecular Dynamics Simulation studies for 50ns using GROMACS. Measures of structural similarity based on 2D Fingerprint Screening and Molecular Dynamics Simulation studies can suggest good leads for drug designing. The novelty of this study is that the workflow used here yields the same results that are at par with the experimental data.

CONCLUSION

This suggests the use of the 2D fingerprint similarity search in various databases, followed by multiple docking algorithms and dynamics as a workflow that will lead to finding novel compounds that a structurally and functionally similar to LEF and BREQ.

Scaffold hopping computational approach for searching novel β-lactamase inhibitors

We present a novel computational ligand-based virtual screening approach with scaffold hopping capabilities for the identification of novel inhibitors of β-lactamases which confer bacterial resistance to β-lactam antibiotics. The structures of known β-lactamase inhibitors were used as query ligands, and a virtual in silico screening a database of 8 million drug-like compounds was performed in order to select the ligands with similar shape and charge distribution. A set of numerical descriptors was used such as chirality, eigen spectrum of matrices of interatomic distances and connectivity together with higher order moment invariants that showed their efficiency in the field of pattern recognition but have not yet been employed in drug discovery. The developed scaffold-hopping approach was applied for the discovery of analogues of four allosteric inhibitors of serine β-lactamases. After a virtual in silico screening, the effect of two selected ligands on the activity of TEM type β-lactamase was studied experimentally. New non-β-lactam inhibitors were found that showed more effective inhibition of β-lactamases compared to query ligands.

A Massively Parallel Selection of Small Molecule-RNA Motif Binding Partners Informs Design of an Antiviral from Sequence

Many RNAs cause disease; however, RNA is rarely exploited as a small-molecule drug target. Our programmatic focus is to define privileged RNA motif small-molecule interactions to enable the rational design of compounds that modulate RNA biology starting from only sequence. We completed a massive, library-versus-library screen that probed over 50 million binding events between RNA motifs and small molecules. The resulting data provide a rich encyclopedia of small-molecule RNA recognition patterns, defining chemotypes and RNA motifs that confer selective, avid binding. The resulting interaction maps were mined against the entire viral genome of hepatitis C virus (HCV). A small molecule was identified that avidly bound RNA motifs present in the HCV 30 UTR and inhibited viral replication while having no effect on host cells. Collectively, this study represents the first whole-genome pattern recognition between small molecules and RNA folds.

Discovering Highly Potent Molecules from an Initial Set of Inactives Using Iterative Screening

The versatility of similarity searching and quantitative structure-activity relationships to model the activity of compound sets within given bioactivity ranges (i.e., interpolation) is well established. However, their relative performance in the common scenario in early stage drug discovery where lots of inactive data but no active data points are available (i.e., extrapolation from the low-activity to the high-activity range) has not been thoroughly examined yet. To this aim, we have designed an iterative virtual screening strategy which was evaluated on 25 diverse bioactivity data sets from ChEMBL. We benchmark the efficiency of random forest (RF), multiple linear regression, ridge regression, similarity searching, and random selection of compounds to identify a highly active molecule in the test set among a large number of low-potency compounds. We use the number of iterations required to find this active molecule to evaluate the performance of each experimental setup. We show that linear and ridge regression often outperform RF and similarity searching, reducing the number of iterations to find an active compound by a factor of 2 or more. Even simple regression methods seem better able to extrapolate to high-bioactivity ranges than RF, which only provides output values in the range covered by the training set. In addition, examination of the scaffold diversity in the data sets used shows that in some cases similarity searching and RF require two times as many iterations as random selection depending on the chemical space covered in the initial training data. Lastly, we show using bioactivity data for COX-1 and COX-2 that our framework can be extended to multitarget drug discovery, where compounds are selected by concomitantly considering their activity against multiple targets. Overall, this study provides an approach for iterative screening where only inactive data are present in early stages of drug discovery in order to discover highly potent compounds and the best experimental set up in which to do so.

Maximal Unbiased Benchmarking Data Sets for Human Chemokine Receptors and Comparative Analysis

Chemokine receptors (CRs) have long been druggable targets for the treatment of inflammatory diseases and HIV-1 infection. As a powerful technique, virtual screening (VS) has been widely applied to identifying small molecule leads for modern drug targets including CRs. For rational selection of a wide variety of VS approaches, ligand enrichment assessment based on a benchmarking data set has become an indispensable practice. However, the lack of versatile benchmarking sets for the whole CRs family that are able to unbiasedly evaluate every single approach including both structure- and ligand-based VS somewhat hinders modern drug discovery efforts. To address this issue, we constructed Maximal Unbiased Benchmarking Data sets for human Chemokine Receptors (MUBD-hCRs) using our recently developed tools of MUBD-DecoyMaker. The MUBD-hCRs encompasses 13 subtypes out of 20 chemokine receptors, composed of 404 ligands and 15756 decoys so far and is readily expandable in the future. It had been thoroughly validated that MUBD-hCRs ligands are chemically diverse while its decoys are maximal unbiased in terms of "artificial enrichment", "analogue bias". In addition, we studied the performance of MUBD-hCRs, in particular CXCR4 and CCR5 data sets, in ligand enrichment assessments of both structure- and ligand-based VS approaches in comparison with other benchmarking data sets available in the public domain and demonstrated that MUBD-hCRs is very capable of designating the optimal VS approach. MUBD-hCRs is a unique and maximal unbiased benchmarking set that covers major CRs subtypes so far.

Adversarial Threshold Neural Computer for Molecular de Novo Design

In this article, we propose the deep neural network Adversarial Threshold Neural Computer (ATNC). The ATNC model is intended for the de novo design of novel small-molecule organic structures. The model is based on generative adversarial network architecture and reinforcement learning. ATNC uses a Differentiable Neural Computer as a generator and has a new specific block, called adversarial threshold (AT). AT acts as a filter between the agent (generator) and the environment (discriminator + objective reward functions). Furthermore, to generate more diverse molecules we introduce a new objective reward function named Internal Diversity Clustering (IDC). In this work, ATNC is tested and compared with the ORGANIC model. Both models were trained on the SMILES string representation of the molecules, using four objective functions (internal similarity, Muegge druglikeness filter, presence or absence of sp³-rich fragments, and IDC). The SMILES representations of 15K druglike molecules from the ChemDiv collection were used as a training data set. For the different functions, ATNC outperforms ORGANIC. Combined with the IDC, ATNC generates 72% of valid and 77% of unique SMILES strings, while ORGANIC generates only 7% of valid and 86% of unique SMILES strings. For each set of molecules generated by ATNC and ORGANIC, we analyzed distributions of four molecular descriptors (number of atoms, molecular weight, logP, and tpsa) and calculated five chemical statistical features (internal diversity, number of unique heterocycles, number of clusters, number of singletons, and number of compounds that have not been passed through medicinal chemistry filters). Analysis of key molecular descriptors and chemical statistical features demonstrated that the molecules generated by ATNC elicited better druglikeness properties. We also performed in vitro validation of the molecules generated by ATNC; results indicated that ATNC is an effective method for producing hit compounds.

PyBioMed: a python library for various molecular representations of chemicals, proteins and DNAs and their interactions

Background

With the increasing development of biotechnology and informatics technology, publicly available data in chemistry and biology are undergoing explosive growth. Such wealthy information in these data needs to be extracted and transformed to useful knowledge by various data mining methods. Considering the amazing rate at which data are accumulated in chemistry and biology fields, new tools that process and interpret large and complex interaction data are increasingly important. So far, there are no suitable toolkits that can effectively link the chemical and biological space in view of molecular representation. To further explore these complex data, an integrated toolkit for various molecular representation is urgently needed which could be easily integrated with data mining algorithms to start a full data analysis pipeline.

Results

Herein, the python library PyBioMed is presented, which comprises functionalities for online download for various molecular objects by providing different IDs, the pretreatment of molecular structures, the computation of various molecular descriptors for chemicals, proteins, DNAs and their interactions. PyBioMed is a feature-rich and highly customized python library used for the characterization of various complex chemical and biological molecules and interaction samples. The current version of PyBioMed could calculate 775 chemical descriptors and 19 kinds of chemical fingerprints, 9920 protein descriptors based on protein sequences, more than 6000 DNA descriptors from nucleotide sequences, and interaction descriptors from pairwise samples using three different combining strategies. Several examples and five real-life applications were provided to clearly guide the users how to use PyBioMed as an integral part of data analysis projects. By using PyBioMed, users are able to start a full pipelining from getting molecular data, pretreating molecules, molecular representation to constructing machine learning models conveniently.

Conclusion

PyBioMed provides various user-friendly and highly customized APIs to calculate various features of biological molecules and complex interaction samples conveniently, which aims at building integrated analysis pipelines from data acquisition, data checking, and descriptor calculation to modeling. PyBioMed is freely available at http://projects.scbdd.com/pybiomed.html.

Spectrophores as one-dimensional descriptors calculated from three-dimensional atomic properties: applications ranging from scaffold hopping to multi-target virtual screening

Spectrophores are novel descriptors that are calculated from the three-dimensional atomic properties of molecules. In our current implementation, the atomic properties that were used to calculate spectrophores include atomic partial charges, atomic lipophilicity indices, atomic shape deviations and atomic softness properties. This approach can easily be widened to also include additional atomic properties. Our novel methodology finds its roots in the experimental affinity fingerprinting technology developed in the 1990’s by Terrapin Technologies. Here we have translated it into a purely virtual approach using artificial affinity cages and a simplified metric to calculate the interaction between these cages and the atomic properties. A typical spectrophore consists of a vector of 48 real numbers. This makes it highly suitable for the calculation of a wide range of similarity measures for use in virtual screening and for the investigation of quantitative structure–activity relationships in combination with advanced statistical approaches such as self-organizing maps, support vector machines and neural networks. In our present report we demonstrate the applicability of our novel methodology for scaffold hopping as well as virtual screening.

Computer-Assisted Retrosynthesis Based on Molecular Similarity

We demonstrate molecular similarity to be a surprisingly effective metric for proposing and ranking one-step retrosynthetic disconnections based on analogy to precedent reactions. The developed approach mimics the retrosynthetic strategy defined implicitly by a corpus of known reactions without the need to encode any chemical knowledge. Using 40 000 reactions from the patent literature as a knowledge base, the recorded reactants are among the top 10 proposed precursors in 74.1% of 5000 test reactions, providing strong quantitative support for our methodology. Extension of the one-step strategy to multistep pathway planning is demonstrated and discussed for two exemplary drug products.

Analysing and Navigating Natural Products Space for Generating Small, Diverse, But Representative Chemical Libraries

Armed with the digital availability of two natural products libraries, amounting to some 195 885 molecular entities, we ask the question of how we can best sample from them to maximize their "representativeness" in smaller and more usable libraries of 96, 384, 1152, and 1920 molecules. The term "representativeness" is intended to include diversity, but for numerical reasons (and the likelihood of being able to perform a QSAR) it is necessary to focus on areas of chemical space that are more highly populated. Encoding chemical structures as fingerprints using the RDKit "patterned" algorithm, we first assess the granularity of the natural products space using a simple clustering algorithm, showing that there are major regions of "denseness" but also a great many very sparsely populated areas. We then apply a "hybrid" hierarchical K-means clustering algorithm to the data to produce more statistically robust clusters from which representative and appropriate numbers of samples may be chosen. There is necessarily again a trade-off between cluster size and cluster number, but within these constraints, libraries containing 384 or 1152 molecules can be found that come from clusters that represent some 18 and 30% of the whole chemical space, with cluster sizes of, respectively, 50 and 27 or above, just about sufficient to perform a QSAR. By using the online availability of molecules via the Molport system (www.molport.com), we are also able to construct (and, for the first time, provide the contents of) a small virtual library of available molecules that provided effective coverage of the chemical space described. Consistent with this, the average molecular similarities of the contents of the libraries developed is considerably smaller than is that of the original libraries. The suggested libraries may have use in molecular or phenotypic screening, including for determining possible transporter substrates.

Chemometric-assisted method development in hydrophilic interaction liquid chromatography: A review

With an enormous growth in the application of hydrophilic interaction liquid chromatography (HILIC), there has also been significant progress in HILIC method development. HILIC is a chromatographic method that utilises hydro-organic mobile phases with a high organic content, and a hydrophilic stationary phase. It has been applied predominantly in the determination of small polar compounds. Theoretical studies in computer-aided modelling tools, most importantly the predictive, quantitative structure retention relationship (QSRR) modelling methods, have attracted the attention of researchers and these approaches greatly assist the method development process. This review focuses on the application of computer-aided modelling tools in understanding the retention mechanism, the classification of HILIC stationary phases, prediction of retention times in HILIC systems, optimisation of chromatographic conditions, and description of the interaction effects of the chromatographic factors in HILIC separations. Additionally, what has been achieved in the potential application of QSRR methodology in combination with experimental design philosophy in the optimisation of chromatographic separation conditions in the HILIC method development process is communicated. Developing robust predictive QSRR models will undoubtedly facilitate more application of this chromatographic mode in a broader variety of research areas, significantly minimising cost and time of the experimental work.

The role of drug profiles as similarity metrics: applications to repurposing, adverse effects detection and drug-drug interactions

Explosion of the availability of big data sources along with the development in computational methods provides a useful framework to study drugs' actions, such as interactions with pharmacological targets and off-targets. Databases related to protein interactions, adverse effects and genomic profiles are available to be used for the construction of computational models. In this article, we focus on the description of biological profiles for drugs that can be used as a system to compare similarity and create methods to predict and analyze drugs' actions. We highlight profiles constructed with different biological data, such as target-protein interactions, gene expression measurements, adverse effects and disease profiles. We focus on the discovery of new targets or pathways for drugs already in the pharmaceutical market, also called drug repurposing, in the interaction with off-targets responsible for adverse reactions and in drug-drug interaction analysis. The current and future applications, strengths and challenges facing all these methods are also discussed. Biological profiles or signatures are an important source of data generation to deeply analyze biological actions with important implications in drug-related studies.

Analysis of drug-endogenous human metabolite similarities in terms of their maximum common substructures

In previous work, we have assessed the structural similarities between marketed drugs (‘drugs’) and endogenous natural human metabolites (‘metabolites’ or ‘endogenites’), using ‘fingerprint’ methods in common use, and the Tanimoto and Tversky similarity metrics, finding that the fingerprint encoding used had a dramatic effect on the apparent similarities observed. By contrast, the maximal common substructure (MCS), when the means of determining it is fixed, is a means of determining similarities that is largely independent of the fingerprints, and also has a clear chemical meaning. We here explored the utility of the MCS and metrics derived therefrom. In many cases, a shared scaffold helps cluster drugs and endogenites, and gives insight into enzymes (in particular transporters) that they both share. Tanimoto and Tversky similarities based on the MCS tend to be smaller than those based on the MACCS fingerprint-type encoding, though the converse is also true for a significant fraction of the comparisons. While no single molecular descriptor can account for these differences, a machine learning-based analysis of the nature of the differences (MACCS_Tanimoto vs MCS_Tversky) shows that they are indeed deterministic, although the features that are used in the model to account for this vary greatly with each individual drug. The extent of its utility and interpretability vary with the drug of interest, implying that while MCS is neither ‘better’ nor ‘worse’ for every drug–endogenite comparison, it is sufficiently different to be of value. The overall conclusion is thus that the use of the MCS provides an additional and valuable strategy for understanding the structural basis for similarities between synthetic, marketed drugs and natural intermediary metabolites.

The impact of in silico screening in the discovery of novel and safer drug candidates

Drug discovery is a multidisciplinary and multivariate optimization endeavor. As such, in silico screening tools have gained considerable importance to archive, analyze and exploit the vast and ever-increasing amount of experimental data generated throughout the process. The current review will focus on the computer-aided prediction of the numerous properties that need to be controlled during the discovery of a preliminary hit and its promotion to a viable clinical candidate. It does not pretend to the almost impossible task of an exhaustive report but will highlight a few key points that need to be collectively addressed both by chemists and biologists to fuel the drug discovery pipeline with innovative and safe drug candidates.

MayaChemTools: An Open Source Package for Computational Drug Discovery

MayaChemTools is a growing collection of Perl scripts, modules, and classes to support a variety of computational drug discovery needs, such as manipulation and analysis of data, generation of two-dimensional (2D) fingerprints, similarity searching, and calculation of physicochemical properties. MayaChemTools provides command line scripts for the following tasks: manipulation and analysis of data in SD, CSV/TSV, sequence/alignments, and PDB files; calculation of a key set of physicochemical properties, such as molecular weight, hydrogen bond donors and acceptors, logP, and topological polar surface area; generation of 2D fingerprints corresponding to atom neighborhoods, atom types, E-state indices, extended connectivity, MACCS keys, path lengths, topological atom pairs, topological atom triplets, topological atom torsions, topological pharmacophore atom pairs, and topological pharmacophore atom triplets; similarity searching and calculation of similarity matrices using available 2D fingerprints; listing properties of elements in the periodic table, amino acids, and nucleic acids; and exporting data from relational databases. An extensive set of modules and classes are also available for custom development. MayaChemTools is freely available online at www.MayaChemTools.org , under the terms of the GNU LGPL, as published by the Free Software Foundation.

Leveraging 3D chemical similarity, target and phenotypic data in the identification of drug-protein and drug-adverse effect associations

BACKGROUND

Drug-target identification is crucial to discover novel applications for existing drugs and provide more insights about mechanisms of biological actions, such as adverse drug effects (ADEs). Computational methods along with the integration of current big data sources provide a useful framework for drug-target and drug-adverse effect discovery.

RESULTS

In this article, we propose a method based on the integration of 3D chemical similarity, target and adverse effect data to generate a drug-target-adverse effect predictor along with a simple leveraging system to improve identification of drug-targets and drug-adverse effects. In the first step, we generated a system for multiple drug-target identification based on the application of 3D drug similarity into a large target dataset extracted from the ChEMBL. Next, we developed a target-adverse effect predictor combining targets from ChEMBL with phenotypic information provided by SIDER data source. Both modules were linked to generate a final predictor that establishes hypothesis about new drug-target-adverse effect candidates. Additionally, we showed that leveraging drug-target candidates with phenotypic data is very useful to improve the identification of drug-targets. The integration of phenotypic data into drug-target candidates yielded up to twofold precision improvement. In the opposite direction, leveraging drug-phenotype candidates with target data also yielded a significant enhancement in the performance.

CONCLUSIONS

The modeling described in the current study is simple and efficient and has applications at large scale in drug repurposing and drug safety through the identification of mechanism of action of biological effects.

Uniqueness: skews bit occurrence frequencies in randomly generated fingerprint libraries

Requiring that randomly generated chemical fingerprint libraries have unique fingerprints such that no two fingerprints are identical causes a systematic skew in bit occurrence frequencies, the proportion at which specified bits are set. Observed frequencies (O) at which each bit is set within the resulting libraries systematically differ from frequencies at which bits are set at fingerprint generation (E). Observed frequencies systematically skew toward 0.5, with the effect being more pronounced as library size approaches the compound space, which is the total number of unique possible fingerprints given the number of bit positions each fingerprint contains. The effect is quantified for varying library sizes as a fraction of the overall compound space, and for changes in the specified frequency E. The cause and implications for this systematic skew are subsequently discussed. When generating random libraries of chemical fingerprints, the imposition of a uniqueness requirement should either be avoided or taken into account.

Integrating computational and chemical biology tools in the discovery of antiangiogenic small molecule ligands of FGF2 derived from endogenous inhibitors

The FGFs/FGFRs system is a recognized actionable target for therapeutic approaches aimed at inhibiting tumor growth, angiogenesis, metastasis, and resistance to therapy. We previously identified a non-peptidic compound (SM27) that retains the structural and functional properties of the FGF2-binding sequence of thrombospondin-1 (TSP-1), a major endogenous inhibitor of angiogenesis. Here we identified new small molecule inhibitors of FGF2 based on the initial lead. A similarity-based screening of small molecule libraries, followed by docking calculations and experimental studies, allowed selecting 7 bi-naphthalenic compounds that bound FGF2 inhibiting its binding to both heparan sulfate proteoglycans and FGFR-1. The compounds inhibit FGF2 activity in in vitro and ex vivo models of angiogenesis, with improved potency over SM27. Comparative analysis of the selected hits, complemented by NMR and biochemical analysis of 4 newly synthesized functionalized phenylamino-substituted naphthalenes, allowed identifying the minimal stereochemical requirements to improve the design of naphthalene sulfonates as FGF2 inhibitors.