Abstract
BACKGROUND: There is a need for computational approaches to compare small organic molecules based on chemical similarity or for evaluating biochemical transformations. No tool currently exists to generate global molecular alignments for small organic molecules. The study introduces a new approach to molecular alignment in the Simplified Molecular Input Line Entry System (SMILES) format. This method leverages programming and scoring alignments to minimize differences in electronegativity, here using a measure of atomic partial charges to address the challenge of understanding structural transformations in reaction pathways. This can be applied to study transitions from linear to cyclical pathways. RESULTS: The proposed method is based on the Needleman-Wunsch algorithm for sequence alignment, but it uses a modified scoring function for different input data. Validation against a benchmarked dataset from the Krebs cycle, based on the known chemical transformations in the pathway, confirmed the efficacy of the approach in aligning atoms that are known to be the same across the transformation. The algorithm also quantified each transformation of metabolites in the Pentose Phosphate Pathway and in Glycolysis. The method was used to study the difference in chemical similarity over transformations between linear and cyclical pathways. The study found a midpoint dissimilarity peak in cyclical pathways (particularly the Krebs Cycle) and a progressive decrease in molecular similarity in linear pathways, consistent with expectations. CONCLUSIONS: The study introduces an algorithm that quantifies molecular transformations in metabolic pathways. The algorithm effectively highlights structural changes and was applied to a hypothesis about the transition from linear to cyclical structures. The software, which provides valuable insights into molecular transformations, is available at: https://github.com/24atang/SMILES-Alignment.git.