ReFlex3D: Refined Flexible Alignment of Molecules Using Shape and Electrostatics

PheSA: An Open-Source Tool for Pharmacophore-Enhanced Shape Alignment

PheSA is an open-source pharmacophore- and shape-based screening and molecular alignment tool that is fully open-source as part of OpenChemLib. Supporting standard ligand-based screening, flexible refinement of alignments, and receptor-guided shape docking, PheSA is a very flexible tool and can be used for different use cases in structure-based drug design. We present the algorithm and different benchmark studies that investigate the screening performance and also the quality of the generated alignments and the pose prediction performance of the receptor-guided PheSA algorithm. An important finding is the effect of the type of similarity metric used for measuring screening enrichment (symmetric Tanimoto versus asymmetric Tversky), whereby we could observe improved enrichment rates by using Tversky. PheSA exhibits enrichments on the DUD-E that are on par with commercial methods.

SENSAAS-Flex: a joint optimization approach for aligning 3D shapes and exploring the molecular conformation space

MOTIVATION

Popular shape-based alignment methods handle molecular flexibility by utilizing conformational ensembles to select the most fitted conformer. However, the initial conformer library generation step is computationally intensive and limiting to the overall alignment process. In this work, we describe a method to perform flexible alignment of two molecular shapes by optimizing the 3D conformation. SENSAAS-Flex, an add-on to the SENSAAS tool, is able to proceed from a limited set of initial conformers through an iterative process where additional conformational optimizations are made at the substructure level and constrained by the target shape.

RESULTS

In self- and cross-alignment experiments, SENSAAS-Flex is able to reproduce the crystal structure geometry of ligands of the AstraZeneca Molecule Overlay Test set and PDBbind refined dataset. Our study shows that the point-based representation of molecular surfaces is appropriate in terms of shape constraint to sample the conformational space and perform flexible molecular alignments.

AVAILABILITY AND IMPLEMENTATION

The documentation and source code are available at https://chemoinfo.ipmc.cnrs.fr/Sensaas-flex/sensaas-flex-main.tar.gz.

Derivation of Highly Predictive 3D-QSAR Models for hERG Channel Blockers Based on the Quantum Artificial Neural Network Algorithm

The hERG potassium channel serves as an annexed target for drug discovery because the associated off-target inhibitory activity may cause serious cardiotoxicity. Quantitative structure-activity relationship (QSAR) models were developed to predict inhibitory activities against the hERG potassium channel, utilizing the three-dimensional (3D) distribution of quantum mechanical electrostatic potential (ESP) as the molecular descriptor. To prepare the optimal atomic coordinates of dataset molecules, pairwise 3D structural alignments were carried out in order for the quantum mechanical cross correlation between the template and other molecules to be maximized. This alignment method stands out from the common atom-by-atom matching technique, as it can handle structurally diverse molecules as effectively as chemical derivatives that share an identical scaffold. The alignment problem prevalent in 3D-QSAR methods was ameliorated substantially by dividing the dataset molecules into seven subsets, each of which contained molecules with similar molecular weights. Using an artificial neural network algorithm to find the functional relationship between the quantum mechanical ESP descriptors and the experimental hERG inhibitory activities, highly predictive 3D-QSAR models were derived for all seven molecular subsets to the extent that the squared correlation coefficients exceeded 0.79. Given their simplicity in model development and strong predictability, the 3D-QSAR models developed in this study are expected to function as an effective virtual screening tool for assessing the potential cardiotoxicity of drug candidate molecules.

Quantum Artificial Neural Network Approach to Derive a Highly Predictive 3D-QSAR Model for Blood-Brain Barrier Passage

A successful passage of the blood–brain barrier (BBB) is an essential prerequisite for the drug molecules designed to act on the central nervous system. The logarithm of blood–brain partitioning (LogBB) has served as an effective index of molecular BBB permeability. Using the three-dimensional (3D) distribution of the molecular electrostatic potential (ESP) as the numerical descriptor, a quantitative structure-activity relationship (QSAR) model termed AlphaQ was derived to predict the molecular LogBB values. To obtain the optimal atomic coordinates of the molecules under investigation, the pairwise 3D structural alignments were conducted in such a way to maximize the quantum mechanical cross correlation between the template and a target molecule. This alignment method has the advantage over the conventional atom-by-atom matching protocol in that the structurally diverse molecules can be analyzed as rigorously as the chemical derivatives with the same scaffold. The inaccuracy problem in the 3D structural alignment was alleviated in a large part by categorizing the molecules into the eight subsets according to the molecular weight. By applying the artificial neural network algorithm to associate the fully quantum mechanical ESP descriptors with the extensive experimental LogBB data, a highly predictive 3D-QSAR model was derived for each molecular subset with a squared correlation coefficient larger than 0.8. Due to the simplicity in model building and the high predictability, AlphaQ is anticipated to serve as an effective computational screening tool for molecular BBB permeability.

A comprehensive comparative assessment of 3D molecular similarity tools in ligand-based virtual screening

Three-dimensional (3D) molecular similarity, one major ligand-based virtual screening (VS) method, has been widely used in the drug discovery process. A variety of 3D molecular similarity tools have been developed in recent decades. In this study, we assessed a panel of 15 3D molecular similarity programs against the DUD-E and LIT-PCBA datasets, including commercial ROCS and Phase, in terms of screening power and scaffold-hopping power. The results revealed that (1) SHAFTS, LS-align, Phase Shape_Pharm and LIGSIFT showed the best VS capability in terms of screening power. Some 3D similarity tools available to academia can yield relatively better VS performance than commercial ROCS and Phase software. (2) Current 3D similarity VS tools exhibit a considerable ability to capture actives with new chemotypes in terms of scaffold hopping. (3) Multiple conformers relative to single conformations will generally improve VS performance for most 3D similarity tools, with marginal improvement observed in area under the receiving operator characteristic curve values, enrichment factor in the top 1% and hit rate in the top 1% values showed larger improvement. Moreover, redundancy and complementarity analyses of hit lists from different query seeds and different 3D similarity VS tools showed that the combination of different query seeds and/or different 3D similarity tools in VS campaigns retrieved more (and more diverse) active molecules. These findings provide useful information for guiding choices of the optimal 3D molecular similarity tools for VS practices and designing possible combination strategies to discover more diverse active compounds.

sensaas: Shape-based Alignment by Registration of Colored Point-based Surfaces

sensaas is a tool developed for aligning and comparing molecular shapes and sub-shapes. Alignment is obtained by registration of 3D point-based representations of the van der Waals surface. The method uses local properties of the shape to identify the correspondence relationships between two point clouds containing up to several thousand colored (labeled) points. Our rigid-body superimposition method follows a two-stage approach. An initial alignment is obtained by matching pose-invariant local 3D descriptors, called FPFH, of the input point clouds. This stage provides a global superimposition of the molecular surfaces, without any knowledge of their initial pose in 3D space. This alignment is then refined by optimizing the matching of colored points. In our study, each point is colored according to its closest atom, which itself belongs to a user defined physico-chemical class. Finally, sensaas provides an alignment and evaluates the molecular similarity by using Tversky coefficients. To assess the efficiency of this approach, we tested its ability to reproduce the superimposition of X-ray structures of the benchmarking AstraZeneca (AZ) data set and, compared its results with those generated by the two shape-alignment approaches shaep and shafts. We also illustrated submatching properties of our method with respect to few substructures and bioisosteric fragments. The code is available upon request from the authors (demo version at https://chemoinfo.ipmc.cnrs.fr/SENSAAS).

Development and application of a comprehensive machine learning program for predicting molecular biochemical and pharmacological properties

We establish a comprehensive quantitative structure-activity relationship (QSAR) model termed AlphaQ through the machine learning algorithm to associate the fully quantum mechanical molecular descriptors with various biochemical and pharmacological properties. Preliminarily, a novel method for molecular structural alignments was developed in such a way to maximize the quantum mechanical cross correlations among the molecules. Besides the improvement of structural alignments, three-dimensional (3D) distribution of the molecular electrostatic potential was introduced as the unique numerical descriptor for individual molecules. These dual modifications lead to a substantial accuracy enhancement in multifarious 3D-QSAR prediction models of AlphaQ. Most remarkably, AlphaQ has been proven to be applicable to structurally diverse molecules to the extent that it outperforms the conventional QSAR methods in estimating the inhibitory activity against thrombin, the water-cyclohexane distribution coefficient, the permeability across the membrane of the Caco-2 cell, and the metabolic stability in human liver microsomes. Due to the simplicity in model building and the high predictive capability for varying biochemical and pharmacological properties, AlphaQ is anticipated to serve as a valuable screening tool at both early and late stages of drug discovery.