- Poster presentation
- Open Access
Making the most of approximate maximum common substructure search
Journal of Cheminformaticsvolume 6, Article number: P29 (2014)
The maximum common substructure (MCS) problem is of great importance in multiple aspects of chemoinformatics. It has diverse applications ranging from lead prediction to automated reaction mapping and visual alignment of similar compounds. Many different algorithms have been developed , both exact and approximate. Since the MCS problem is NP-complete, the strict time constraints of most applications can only be realistically satisfied by fast and robust approximation methods.
We developed two efficient heuristic algorithms. One is based on the popular approach of reducing the MCS problem to finding the maximum clique in the modular product of the two molecule graphs. The other is based on a new algorithm by Kawabata, called the build-up method . We also incorporated other techniques, for example, the topological fingerprinting primarily used in substructure and similarity searching. We optimized our implementations for use in multiple applications developed at ChemAxon. In some applications, for example, hierarchical MCS-based clustering or similarity search in large databases, the algorithms are required to give close to optimal results in limited time. To meet these conflicting demands, our implementations were enhanced with strong heuristics. Upper bound calculation methods were also applied for screening and early termination purposes.
In other applications, for example, reaction mapping or visual alignment, the challenge is that topological features must also be taken into account. Apart from the size of the found common substructure, the determined one-to-one correspondence between the atoms of the molecules is also very important. Effective heuristics were developed to guide the algorithms to prefer those solutions in which the relative positions of the common fragments of the input molecules are as similar as possible.
Our implementations have been thoroughly tested and benchmarked. They have also been compared to publicly available solution methods, and integrated into different products at ChemAxon. This has shown that the presented MCS algorithms can adequately cover the conflicting requirements of typical applications. We present our methods and heuristics along with their effects on running time, memory usage, as well as the size and features of the result.
Raymond JW, Willett P: Maximum common subgraph isomorphism algorithms for the matching of chemical structures. Journal of Computer-Aided Molecular Design. 2002, 16: 521-533. 10.1023/A:1021271615909.
Kawabata T: Build-up algorithm for atomic correspondence between chemical structures. Journal of Chemical Information and Modeling. 2011, 51 (8): 1775-1787. 10.1021/ci2001023.