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

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.

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).