A visual approach for analysis and inference of molecular activity spaces

Abstract

BackgroundMolecular space visualization can help to explore the diversity of large heterogeneous chemical data, which ultimately may increase the understanding of structure-activity relationships (SAR) in drug discovery projects. Visual SAR analysis can therefore be useful for library design, chemical classification for their biological evaluation and virtual screening for the selection of compounds for synthesis or in vitro testing. As such, computational approaches for molecular space visualization have become an important issue in cheminformatics research. The proposed approach uses molecular similarity as the sole input for computing a probabilistic surface of molecular activity (PSMA). This similarity matrix is transformed in 2D using different dimension reduction algorithms (Principal Coordinates Analysis ( PCooA), Kruskal multidimensional scaling, Sammon mapping and t-SNE). From this projection, a kernel density function is applied to compute the probability of activity for each coordinate in the new projected space.ResultsThis methodology was tested over four different quantitative structure-activity relationship (QSAR) binary classification data sets and the PSMAs were computed for each. The generated maps showed internal consistency with active molecules grouped together for all data sets and all dimensionality reduction algorithms. To validate the quality of the generated maps, the 2D coordinates of test molecules were computed into the new reference space using a data transformation matrix. In total sixteen PSMAs were built, and their performance was assessed using the Area Under Curve (AUC) and the Matthews Coefficient Correlation (MCC). For the best projections for each data set, AUC testing results ranged from 0.87 to 0.98 and the MCC scores ranged from 0.33 to 0.77, suggesting this methodology can validly capture the complexities of the molecular activity space. All four mapping functions provided generally good results yet the overall performance of PCooA and t-SNE was slightly better than Sammon mapping and Kruskal multidimensional scaling.ConclusionsOur result showed that by using an appropriate combination of metric space representation and dimensionality reduction applied over metric spaces it is possible to produce a visual PSMA for which its consistency has been validated by using this map as a classification model. The produced maps can be used as prediction tools as it is simple to project any molecule into this new reference space as long as the similarities to the molecules used to compute the initial similarity matrix can be computed.

Random-forest model for drug–target interaction prediction via Kullback–Leibler divergence

“…A well-known nonparametric density estimation method, kernel density estimation (KDE), was selected to estimate the probability density function [ 32 – 35 ]. The built-in likelihood function was maximized to estimate the probability distribution function for the data obtained.…”

Graph Convolutional Network-Enhanced Model for Screening Persistent, Mobile, and Toxic and Very Persistent and Very Mobile Substances

“…Density-based application domains were used to determine the reliable prediction regions in chemical space. 35 To define the application domains of model, the KDE method was employed, which involves the following formulas…”

“…Density-based application domains were used to determine the reliable prediction regions in chemical space . To define the application domains of model, the KDE method was employed, which involves the following formulas .25ex2ex f h ( x ) = 1 n h × ∑ i = 1 n k ( x − x i h ) k ( x , x ′ ) = 1 σ 2 π e − x − x ′ 2 / 2 σ 2 h = true( 4 σ 5 3 n true) 1 / 5 where σ represents the standard deviation of x , n represents the number of x , f h ( x ) represents the density of x in chemical space.…”

Multi-Targeting Approach in Glioblastoma Using Computer-Assisted Drug Discovery Tools to Overcome the Blood–Brain Barrier and Target EGFR/PI3Kp110β Signaling

“…Such models have limitations from the limited data from which they were fitted. So, the concept of an applicability domain (AD) defines the structural boundaries to which a model can be reliably applied [ 23 ]. Moreover, there is a trend to improve the quality of in silico predictions by allying the empirical results of QSAR models and similarity searches to 3D-docking simulations, providing a direct target–ligand interaction study to unravel the energetically most favorable binding mode [ 24 , 25 ].…”