Graph-Based Molecular Alignment (GMA)

We describe a combined 2D/3D approach for the superposition of flexible chemical structures, which is based on recent progress in the efficient identification of common subgraphs and a gradient-based torsion space optimization algorithm. The simplicity of the approach is reflected in its generality and computational efficiency: the suggested approach neither requires precalculated statistics on the conformations of the molecules nor does it make simplifying assumptions on the topology of the molecules being compared. Furthermore, graph-based molecular alignment produces alignments that are consistent with the chemistry of the molecules as well as their general structure, as it depends on both the local connectivities between atoms and the overall topology of the molecules. We validate this approach on benchmark sets taken from the literature and show that it leads to good results compared to computationally and algorithmically more involved methods. The results suggest that, for most practical purposes, graph-based molecular alignment is a viable alternative to molecular field alignment with respect to structural superposition and leads to structures of comparable quality in a fraction of the time.

[1]  Peter Willett,et al.  Maximum common subgraph isomorphism algorithms for the matching of chemical structures , 2002, J. Comput. Aided Mol. Des..

[2]  Peter Willett,et al.  RASCAL: Calculation of Graph Similarity using Maximum Common Edge Subgraphs , 2002, Comput. J..

[3]  C. Bystroff An alternative derivation of the equations of motion in torsion space for a branched linear chain. , 2001, Protein engineering.

[4]  Andreas Krämer,et al.  Fast 3D molecular superposition and similarity search in databases of flexible molecules , 2003, J. Comput. Aided Mol. Des..

[5]  Ian D. Walker,et al.  Overview of damped least-squares methods for inverse kinematics of robot manipulators , 1995, J. Intell. Robotic Syst..

[6]  Simon K. Kearsley,et al.  An alternative method for the alignment of molecular structures: Maximizing electrostatic and steric overlap , 1990 .

[7]  Peter Willett,et al.  Promoting Access to White Rose Research Papers Effectiveness of Graph-based and Fingerprint-based Similarity Measures for Virtual Screening of 2d Chemical Structure Databases , 2022 .

[8]  P. Labute,et al.  Flexible alignment of small molecules. , 2001, Journal of medicinal chemistry.

[9]  R. Cramer,et al.  Comparative molecular field analysis (CoMFA). 1. Effect of shape on binding of steroids to carrier proteins. , 1988, Journal of the American Chemical Society.

[10]  M G Rossmann,et al.  Three-dimensional structures of drug-resistant mutants of human rhinovirus 14. , 1990, Journal of molecular biology.

[11]  R. Cramer,et al.  Validation of the general purpose tripos 5.2 force field , 1989 .

[12]  W. Kabsch A solution for the best rotation to relate two sets of vectors , 1976 .

[13]  Thomas Lengauer,et al.  RigFit: A new approach to superimposing ligand molecules , 1998, German Conference on Bioinformatics.

[14]  Sung Jin Cho,et al.  FLAME: A Program to Flexibly Align Molecules. , 2006 .

[15]  Matthias Rarey,et al.  Feature trees: A new molecular similarity measure based on tree matching , 1998, J. Comput. Aided Mol. Des..

[16]  Peter Willett,et al.  Scaffold Hopping Using Clique Detection Applied to Reduced Graphs , 2006, J. Chem. Inf. Model..

[17]  Harald Mauser,et al.  A robust clustering method for chemical structures. , 2005, Journal of medicinal chemistry.

[18]  C. Lemmen,et al.  FLEXS: a method for fast flexible ligand superposition. , 1998, Journal of medicinal chemistry.