ReactionCode: format for reaction searching, analysis, classification, transform, and encoding/decoding.

In the past two decades a lot of different formats for molecules and reactions have been created. These formats were mostly developed for the purposes of identifiers, representation, classification, analysis and data exchange. A lot of efforts have been made on molecule formats but only few for reactions where the endeavors have been made mostly by companies leading to proprietary formats. Here, we present ReactionCode: a new open-source format that allows one to encode and decode a reaction into multi-layer machine readable code, which aggregates reactants and products into a condensed graph of reaction (CGR). This format is flexible and can be used in a context of reaction similarity searching and classification. It is also designed for database organization, machine learning applications and as a new transform reaction language.

SAVI, in silico generation of billions of easily synthesizable compounds through expert-system type rules.

We have made available a database of over 1 billion compounds predicted to be easily synthesizable, called Synthetically Accessible Virtual Inventory (SAVI). They have been created by a set of transforms based on an adaptation and extension of the CHMTRN/PATRAN programming languages describing chemical synthesis expert knowledge, which originally stem from the LHASA project. The chemoinformatics toolkit CACTVS was used to apply a total of 53 transforms to about 150,000 readily available building blocks (enamine.net). Only single-step, two-reactant syntheses were calculated for this database even though the technology can execute multi-step reactions. The possibility to incorporate scoring systems in CHMTRN allowed us to subdivide the database of 1.75 billion compounds in sets according to their predicted synthesizability, with the most-synthesizable class comprising 1.09 billion synthetic products. Properties calculated for all SAVI products show that the database should be well-suited for drug discovery. It is being made publicly available for free download from https://doi.org/10.35115/37n9-5738.

Adapting CHMTRN (CHeMistry TRaNslator) for a New Use.

We have adopted and extended the CHMTRN language and used it for the knowledge base of a computer program to generate a large database of synthetically accessible, drug-like chemical structures, the Synthetically Accessible Virtual Inventory (SAVI) Database. CHMTRN is a powerful language originally developed in the LHASA (Logic and Heuristics Applied to Synthetic Analysis) project at Harvard University and used together with the chemical pattern description language, PATRAN, to describe chemical retro-reactions. The languages have proven to be useful beyond the design of retrosynthetic routes and have the potential for much wider use in chemistry; this paper describes CHMTRN and PATRAN as now reimplemented for the forward-synthetic SAVI project but able to describe both forward and retro-reactions.

Toward a Comprehensive Treatment of Tautomerism in Chemoinformatics Including in InChI V2.

We have collected 86 different transforms of tautomeric interconversions. Out of those, 54 are for prototropic (non-ring-chain) tautomerism, 21 for ring-chain tautomerism, and 11 for valence tautomerism. The majority of these rules have been extracted from experimental literature. Twenty rules, covering the most well-known types of tautomerism such as keto-enol tautomerism, were taken from the default handling of tautomerism by the chemoinformatics toolkit CACTVS. The rules were analyzed against nine differerent databases totaling over 400 million (non-unique) structures as to their occurrence rates, mutual overlap in coverage, and recapitulation of the rules' enumerated tautomer sets by InChI V.1.05, both in InChI's Standard and a Nonstandard version with the increased tautomer-handling options 15T and KET turned on. These results and the background of this study are discussed in the context of the IUPAC InChI Project tasked with the redesign of handling of tautomerism for an InChI version 2. Applying the rules presented in this paper would approximately triple the number of compounds in typical small-molecule databases that would be affected by tautomeric interconversion by InChI V2. A web tool has been created to test these rules at https://cactus.nci.nih.gov/tautomerizer.

In silico prediction of Heterocyclic Aromatic Amines metabolism susceptible to form DNA adducts in humans.

Heterocyclic Aromatic Amines (HAAs) are environmental and food contaminants that are classified as probable or possible carcinogens by the International Agency for Research on Cancer. Thirty different HAAs have been identified. However the metabolism of only three of them have been fully characterized in human hepatocytes: AαC (2-amino-9H-pyrido[2,3-b]indole), MeIQx (2-amino-3,8-dimethylimidazo[4,5-f]quinoxaline) and PhIP (2-amino-1-methyl-6-phenyl-imidazo[4,5-b]pyridine). In this study, we use an integrative approach to accurately predict the biotransformation of 30 HAAs into DNA reactive and non DNA reactive compounds. We first build predicted metabolites networks by iterating a knowledge-based expert system of prediction of metabolic reactions based on fingerprint similarities. Next, we combine several methods for predicting Sites Of Metabolism (SOM) in order to reduce the metabolite reaction graphs and to predict the metabolites reactive with DNA. We validate the method by comparing the experimental versus predicted data for the known AαC, MeIQx and PhIP metabolism. 28 of the 30 experimentally determined metabolites are well predicted and 9 of the 10 metabolites known to form DNA adducts are predicted with a high probability to be reactive with DNA. Applying our approach to the 27 unknown HAAs, we generate maps for the metabolic biotransformation of each HAA, including new metabolites with a high-predicted DNA reactivity, which can be further explored through an user-friendly and interactive web interface.

A modeling approach to evaluate the balance between bioactivation and detoxification of MeIQx in human hepatocytes.

BACKGROUND: Heterocyclic aromatic amines (HAA) are environmental and food contaminants that are potentially carcinogenic for humans. 2-Amino-3,8-dimethylimidazo[4,5-f]quinoxaline (MeIQx) is one of the most abundant HAA formed in cooked meat. MeIQx is metabolized by cytochrome P450 1A2 in the human liver into detoxificated and bioactivated products. Once bioactivated, MeIQx metabolites can lead to DNA adduct formation responsible for further genome instability. METHODS: Using a computational approach, we developed a numerical model for MeIQx metabolism in the liver that predicts the MeIQx biotransformation into detoxification or bioactivation pathways according to the concentration of MeIQx. RESULTS: Our results demonstrate that (1) the detoxification pathway predominates, (2) the ratio between detoxification and bioactivation pathways is not linear and shows a maximum at 10 µM of MeIQx in hepatocyte cell models, and (3) CYP1A2 is a key enzyme in the system that regulates the balance between bioactivation and detoxification. Our analysis suggests that such a ratio could be considered as an indicator of MeIQx genotoxicity at a low concentration of MeIQx. CONCLUSIONS: Our model permits the investigation of the balance between bioactivation (i.e., DNA adduct formation pathway through the prediction of potential genotoxic compounds) and detoxification of MeIQx in order to predict the behaviour of this environmental contaminant in the human liver. It highlights the importance of complex regulations of enzyme competitions that should be taken into account in any further multi-organ models.