Reaction Decoder Tool (RDT): extracting features from chemical reactions.

UNLABELLED: Extracting chemical features like Atom-Atom Mapping (AAM), Bond Changes (BCs) and Reaction Centres from biochemical reactions helps us understand the chemical composition of enzymatic reactions. Reaction Decoder is a robust command line tool, which performs this task with high accuracy. It supports standard chemical input/output exchange formats i.e. RXN/SMILES, computes AAM, highlights BCs and creates images of the mapped reaction. This aids in the analysis of metabolic pathways and the ability to perform comparative studies of chemical reactions based on these features.
AVAILABILITY AND IMPLEMENTATION: This software is implemented in Java, supported on Windows, Linux and Mac OSX, and freely available at https://github.com/asad/ReactionDecoder CONTACT: : asad@ebi.ac.uk or s9asad@gmail.com.

1 Introduction

Large-scale chemical reaction databases such as KEGG (Kanehisa ), BRENDA (Chang ), Rhea (Alcántara ) and MetaCyc (Latendresse ) link reactions to enzymes and provide data-mining opportunities for novel pathways (Hatzimanikatis ; Rahman, 2007), and the discovery of drugs, natural products and green chemistry. One of the primary bottlenecks for automated analyses of these chemical reactions comes from the realizations of the imperfect quality of data, such as unmapped or unbalanced reactions. Accurate Atom–Atom Mapping (AAM)—the one-to-one correspondence between the substrate and product atoms (Gasteiger, 2003), will lead to correct prediction of bond changes (BCs) (Rahman ) and the ability to locate the fate of interesting atoms or substructure across metabolic networks (Faulon and Bender, 2010) etc. Linking novel pathways or optimizing pathways of biological/commercial relevance demands better understanding of metabolic routes (Rahman, 2007) and pathway annotation (May ).

We present Reaction Decoder Tool (RDT), a robust open source software for mining reaction features, i.e. BCs, reaction centres and to calculate similarity between reactions etc. The algorithm and competing tools have been published in our previous article (Rahman ). Here we present the source code with relevant changes and below mentioned features.

2 Features

BCs in chemical reactions refers to the cleavage and formation of chemical bonds, changes in bond order and stereo changes, which are due to chemical processes such as chiral inversions or cis-trans isomerization(s). In the chemical reaction diagram and tables atoms are connected by a bond i.e. single ‘–‘, double ‘=’or ring ‘%’ etc. (Fig. 1a), for instance C–C means a single carbon–carbon bond that is cleaved (‘|’) or formed in the reaction (‘‖’). Bond order changes are represented by star ‘*’, e.g. C–C * C = C means a single carbon–carbon bond turning into double carbon-carbon bond or vice versa. Stereo changes are represented as atoms that change their absolute configuration, for instance C(R/S) means a carbon atom that changes from R to S configuration. A reaction centre is the collection of atoms and bonds that are changed during the reaction (Warr, 2014), also known as the local atomic environment around the atoms involved in BCs.

Fig. 1

(a) AAM performed by RDT and the resulting BCs (i.e. formed/cleaved: C%C (Ring), Order Change: C%C*C = C, C-C*C = C) and (b) reaction centres are highlighted in pink colour in the image, where as neighbouring atoms are highlighted in green

The key features of this tool are:

Ability to perform AAM on chemical reactions catalyzed by enzymes (Fig. 1).

In a balanced reaction, the total number of atoms on the left side of the equation (Reactant), equals the total number of atoms on the right (Product). In unbalanced reactions a best ‘guess’ is made.

Generates images of the mapped reactions where matching substructures are highlighted.

Generates reaction patterns and BCs for input reactions.

The input format and resulting mapped reaction with AAM information can be SMILES or RXN file (Gasteiger, 2003).

The SMSD (Rahman ) and CDK (Steinbeck ) are used to process chemical information.

3 Usage and applications

Tools like EC-Blast (Rahman ), FunTree (Sillitoe and Furnham, 2016), MACiE (Holliday ) etc. use RDT in the background to mine and extract chemical information from thousands of enzyme reactions. The success rate of mapping is >99% when compared with manual AAM mappings (Rahman ). Originally developed to explore enzyme reactions, the tool is also useful to explore any kind of organic chemical reaction (Martínez Cuesta ).

4 Conclusion

Reaction Decoder is a robust tool to compute AAM and extract chemical features and calculate similarity between chemical reactions. This is coded in java and optimized to run as a computationally asynchronous process. It is distributed under GNU-GPL V3 license.