Alignment of flexible molecules at their receptor site using 3D descriptors and Hi-PCA

Integrated transcriptomics- and structure-based drug repositioning identifies drugs with proteasome inhibitor properties

Computational pharmacogenomics can potentially identify new indications for already approved drugs and pinpoint compounds with similar mechanism-of-action. Here, we used an integrated drug repositioning approach based on transcriptomics data and structure-based virtual screening to identify compounds with gene signatures similar to three known proteasome inhibitors (PIs; bortezomib, MG-132, and MLN-2238). In vitro validation of candidate compounds was then performed to assess proteasomal proteolytic activity, accumulation of ubiquitinated proteins, cell viability, and drug-induced expression in A375 melanoma and MCF7 breast cancer cells. Using this approach, we identified six compounds with PI properties ((-)-kinetin-riboside, manumycin-A, puromycin dihydrochloride, resistomycin, tegaserod maleate, and thapsigargin). Although the docking scores pinpointed their ability to bind to the β5 subunit, our in vitro study revealed that these compounds inhibited the β1, β2, and β5 catalytic sites to some extent. As shown with bortezomib, only manumycin-A, puromycin dihydrochloride, and tegaserod maleate resulted in excessive accumulation of ubiquitinated proteins and elevated HMOX1 expression. Taken together, our integrated drug repositioning approach and subsequent in vitro validation studies identified six compounds demonstrating properties similar to proteasome inhibitors.

Unlocking the Potential of Clustering and Classification Approaches: Navigating Supervised and Unsupervised Chemical Similarity

BACKGROUND

The field of toxicology has witnessed substantial advancements in recent years, particularly with the adoption of new approach methodologies (NAMs) to understand and predict chemical toxicity. Class-based methods such as clustering and classification are key to NAMs development and application, aiding the understanding of hazard and risk concerns associated with groups of chemicals without additional laboratory work. Advances in computational chemistry, data generation and availability, and machine learning algorithms represent important opportunities for continued improvement of these techniques to optimize their utility for specific regulatory and research purposes. However, due to their intricacy, deep understanding and careful selection are imperative to align the adequate methods with their intended applications.

OBJECTIVES

This commentary aims to deepen the understanding of class-based approaches by elucidating the pivotal role of chemical similarity (structural and biological) in clustering and classification approaches (CCAs). It addresses the dichotomy between general end point-agnostic similarity, often entailing unsupervised analysis, and end point-specific similarity necessitating supervised learning. The goal is to highlight the nuances of these approaches, their applications, and common misuses.

DISCUSSION

Understanding similarity is pivotal in toxicological research involving CCAs. The effectiveness of these approaches depends on the right definition and measure of similarity, which varies based on context and objectives of the study. This choice is influenced by how chemical structures are represented and the respective labels indicating biological activity, if applicable. The distinction between unsupervised clustering and supervised classification methods is vital, requiring the use of end point-agnostic vs. end point-specific similarity definition. Separate use or combination of these methods requires careful consideration to prevent bias and ensure relevance for the goal of the study. Unsupervised methods use end point-agnostic similarity measures to uncover general structural patterns and relationships, aiding hypothesis generation and facilitating exploration of datasets without the need for predefined labels or explicit guidance. Conversely, supervised techniques demand end point-specific similarity to group chemicals into predefined classes or to train classification models, allowing accurate predictions for new chemicals. Misuse can arise when unsupervised methods are applied to end point-specific contexts, like analog selection in read-across, leading to erroneous conclusions. This commentary provides insights into the significance of similarity and its role in supervised classification and unsupervised clustering approaches. https://doi.org/10.1289/EHP14001.

Block-wise Exploration of Molecular Descriptors with Multi-block Orthogonal Component Analysis (MOCA)

Data tables for machine learning and structure-activity relationship modelling (QSAR) are often naturally organized in blocks of data, where multiple molecular representations or sets of descriptors form the blocks. Multi-block Orthogonal Component Analysis (MOCA), a new analytical tool, can be used to explore such data structures in a single model, identifying principal components that are unique to a single block or joint over multiple blocks. We applied MOCA to two sets of 550 and 300 molecules and up to 9213 molecular descriptors organized in 11 blocks. The MOCA models reveal relationships between the blocks and overarching trends across the whole dataset. Based on the MOCA joint components, we propose a quantitative metric for the redundancy of blocks, useful for a priori block-wise feature selection or evaluation of new molecular representations. The second data set includes 7 ecotoxicological study endpoints for crop protection chemicals, for which we (re-)discovered some general trends and linked them to molecular properties. Using a single MOCA model we estimated the predictive potential of each block and the model-ability of the target block.

Exploring the conformational landscapes of HIV protease structural ensembles using principal component analysis

HIV protease, an essential enzyme for viral particle maturation, is an important drug target of HIV. Its structural conformation is a key determinant of both biological function as well as efficient binding of protease inhibitor molecules. In the present study we analyzed 471 crystal structures of HIV-1 protease to understand the conformational changes induced by mutations or binding of various ligands and substrates. We performed principal component analysis on the ensembles of the HIV-1 protease structures to explore the conformational landscapes. The study identified structural differences between drug resistant and drug sensitive protease structures. Conformational changes were identified in the A and B chains of homo-dimeric HIV protease structures having different combinations of mutations, and also rigidity in the binding conformation of HIV drugs within the active site of the protein.© 2018 Wiley Periodicals, Inc.

The landscape of the prion protein's structural response to mutation revealed by principal component analysis of multiple NMR ensembles

Prion Proteins (PrP) are among a small number of proteins for which large numbers of NMR ensembles have been resolved for sequence mutants and diverse species. Here, we perform a comprehensive principle components analysis (PCA) on the tertiary structures of PrP globular proteins to discern PrP subdomains that exhibit conformational change in response to point mutations and clade-specific evolutionary sequence mutation trends. This is to our knowledge the first such large-scale analysis of multiple NMR ensembles of protein structures, and the first study of its kind for PrPs. We conducted PCA on human (n = 11), mouse (n = 14), and wildtype (n = 21) sets of PrP globular structures, from which we identified five conformationally variable subdomains within PrP. PCA shows that different non-local patterns and rankings of variable subdomains arise for different pathogenic mutants. These subdomains may thus be key areas for initiating PrP conversion during disease. Furthermore, we have observed the conformational clustering of divergent TSE-non-susceptible species pairs; these non-phylogenetic clusterings indicate structural solutions towards TSE resistance that do not necessarily coincide with evolutionary divergence. We discuss the novelty of our approach and the importance of PrP subdomains in structural conversion during disease.

Prion Proteins (PrP) cause a variety of incurable TSE diseases, and are among a small number of proteins for which large numbers of NMR ensembles have been resolved for sequence mutants and diverse species. Here, we perform a comprehensive PCA study to assess conformational variation and discern the landscape of the PrP structural response to sequence mutation. This is to our knowledge the first large-scale analysis of multiple NMR ensembles for a specific protein, and the first study to perform a multivariate PCA on the native globular structures of PrP. We conducted exhaustive PCA on three PrP subsets: human and mouse subsets that include structures of sequence mutants, and the set of wild-type PrP (16 PrP species). PCA shows that different non-local patterns of variable subdomains arise for different pathogenic mutants. These subdomains may thus be key areas for initiating PrP conversion during disease. Furthermore, we observed that some evolutionarily divergent species that are non-susceptible to TSEs have surprising structural similarities in their PrPs. We discuss the novelty of our approach with respect to prions, and the advantage of this analysis as a fast, reliable starting point to identify interesting domains that may warrant further experimental and computational analysis.

Megavariate Analysis of Environmental QSAR Data. Part II – Investigating Very Complex Problem Formulations Using Hierarchical, Non-Linear and Batch-Wise Extensions of PCA and PLS

Three extensions of the basic PCA and PLS methodologies are described. These extensions are hierarchical, non-linear and batch-based in nature. The objectives of these methods are to assist in problem understanding and problem solving in very complex (QSAR) problem formulations. The method extensions are illustrated using two example QSAR data sets containing many X- and Y-variables.

Comparing the Performance of FLUFF-BALL to SEAL-CoMFA with a Large Diverse Estrogen Data Set: From Relevant Superpositions to Solid Predictions

In this work a template-based molecular mechanistic superposition algorithm FLUFF (Flexible Ligand Unified Force Field) and an accompanying local coordinate QSAR method BALL (Boundless Adaptive Localized Ligand) are validated against the benchmark techniques SEAL (Steric and Electrostatic Alignment) and CoMFA (Comparative Molecular Field Analysis) using a large diverse set of 245 xenoestrogens extracted from the EDKB (Endocrine Disruptor Knowledge Base) maintained by NCTR (National Centre for Toxicological Research). The results indicate that FLUFF is capable of generating relevant superpositions not only for BALL but also for CoMFA, as both techniques give predictive QSAR models. When the BALL and CoMFA methods are compared, it is clear that the BALL algorithm met or even exceeded the results of the standard 3D-QSAR method CoMFA using alignments either from the tailor-made superposition technique FLUFF or the reference method SEAL. The FLUFF-BALL method can be easily automated, and it is computationally light, providing thus a good computational "sieve" capable of fast screening of large molecule libraries.

Using chemometrics for navigating in the large data sets of genomics, proteomics, and metabonomics (gpm)

This article describes the applicability of multivariate projection techniques, such as principal-component analysis (PCA) and partial least-squares (PLS) projections to latent structures, to the large-volume high-density data structures obtained within genomics, proteomics, and metabonomics. PCA and PLS, and their extensions, derive their usefulness from their ability to analyze data with many, noisy, collinear, and even incomplete variables in both X and Y. Three examples are used as illustrations: the first example is a genomics data set and involves modeling of microarray data of cell cycle-regulated genes in the microorganism Saccharomyces cerevisiae. The second example contains NMR-metabonomics data, measured on urine samples of male rats treated with either of the drugs chloroquine or amiodarone. The third and last data set describes sequence-function classification studies in a set of G-protein-coupled receptors using hierarchical PCA.

Methods for reliability and uncertainty assessment and for applicability evaluations of classification- and regression-based QSARs

This article provides an overview of methods for reliability assessment of quantitative structure-activity relationship (QSAR) models in the context of regulatory acceptance of human health and environmental QSARs. Useful diagnostic tools and data analytical approaches are highlighted and exemplified. Particular emphasis is given to the question of how to define the applicability borders of a QSAR and how to estimate parameter and prediction uncertainty. The article ends with a discussion regarding QSAR acceptability criteria. This discussion contains a list of recommended acceptability criteria, and we give reference values for important QSAR performance statistics. Finally, we emphasize that rigorous and independent validation of QSARs is an essential step toward their regulatory acceptance and implementation.

Megavariate analysis of hierarchical QSAR data

Multivariate PCA- and PLS-models involving many variables are often difficult to interpret, because plots and lists of loadings, coefficients, VIPs, etc, rapidly become messy and hard to overview. There may then be a strong temptation to eliminate variables to obtain a smaller data set. Such a reduction of variables, however, often removes information and makes the modelling efforts less reliable. Model interpretation may be misleading and predictive power may deteriorate. A better alternative is usually to partition the variables into blocks of logically related variables and apply hierarchical data analysis. Such blocked data may be analyzed by PCA and PLS. This modelling forms the base-level of the hierarchical modelling set-up. On the base-level in-depth information is extracted for the different blocks. The score vectors formed on the base-level, here called 'super variables', may be linked together in new matrices on the top-level. On the top-level superficial relationships between the X- and the Y-data are investigated. In this paper the basic principles of hierarchical modelling by means of PCA and PLS are reviewed. One objective of the paper is to disseminate this concept to a broader QSAR audience. The hierarchical methods are used to analyze a set of 10 haloalkanes for which K = 30 chemical descriptors and M = 255 biological responses have been gathered. Due to the complexity of the biological data, they are sub-divided in four blocks. All the modelling steps on the base-level and the top-level are reported and the final QSAR model is interpreted thoroughly.

The effect of anions on azide binding to myoglobin: an unusual functional modulation

The effect of increasing concentrations of several anions on the azide (N(-)(3)) binding properties of sperm whale and horse ferric myoglobin has been studied. Surprisingly, a number of anions may act as heterotropic effectors, decreasing the affinity of myoglobins for N(-)(3), in the following order: ClO(-)(4)=I(-)>Br(-)>Cl(-) and SO(2-)(4), which mirrors the increase in their charge density. The largest effects were measured using ClO(-)(4) and I(-), which produce a 4-fold and 8-fold reduction of the N(-)(3) binding affinity in horse and sperm whale myoglobins, respectively. A dissociation equilibrium constant (K(d)) ranging from 150 to 250 mM was estimated for ClO(-)(4) and I(-) binding to myoglobins. In order to analyse the molecular mechanism producing the reduction of the N(-)(3) binding affinity to ferric myoglobin, the potential anionic binding sites within ferric myoglobin were investigated by a molecular modelling study using the program Grid. Analysis of the theoretical results suggests two particularly favourable binding sites: the first, next to the distal side of the haem, whose occupancy might alter the electrostatic potential surrounding the bound N(-)(3); the second, involving residues of helices B and G which are far from the haem iron atom, thus implying a long range effect on the bound N(-)(3). Based on the evidence that no significant conformational changes are found in the three-dimensional structures of N(-)(3)-free and N(-)(3)-bound myoglobin and on previous results on N(-)(3) binding to ferric myoglobin mutants in CD3 positions, we favour the first hypothesis, suggesting that the functional heterotropic modulation of monomeric myoglobin is mainly depending on a decrease of the positive charge density induced by the binding of anions to the haem distal side.