Bruno Bienfait

Ligand-based models for the isoform specificity of cytochrome P450 3A4, 2D6, and 2C9 substrates

A data set of 379 drugs and drug analogs that are metabolized by human cytochrome P450 (CYP) isoforms 3A4, 2D6, and 2C9, respectively, was studied. A series of descriptor sets directly calculable from the constitution of these drugs was systematically investigated as to their power into classifying a compound into the CYP isoform that metabolizes it. In a four-step build-up process eventually 303 different descriptor components were investigated for 146 compounds of a training set by various model building methods, such as multinomal logistic regression, decision tree, or support vector machine (SVM). Automatic variable selection algorithms were used in order to decrease the number of descriptors. A comprehensive scheme of cross-validation (CV) experiments was applied to assess the robustness and reliability of the four models developed. In addition, the predictive power of the four models presented in this paper was inspected by predicting an external validation data set with 233 compounds. The best model has a leave-one-out (LOO) cross-validated predictivity of 89% and gives 83% correct predictions for the external validation data set. For our favored model we showed the strong influence on the predictivity of the way a data set is split into a training and test data set.

New Publicly Available Chemical Query Language, CSRML, To Support Chemotype Representations for Application to Data Mining and Modeling

Chemotypes are a new approach for representing molecules, chemical substructures and patterns, reaction rules, and reactions. Chemotypes are capable of integrating types of information beyond what is possible using current representation methods (e.g., SMARTS patterns) or reaction transformations (e.g., SMIRKS, reaction SMILES). Chemotypes are expressed in the XML-based Chemical Subgraphs and Reactions Markup Language (CSRML), and can be encoded not only with connectivity and topology but also with properties of atoms, bonds, electronic systems, or molecules. CSRML has been developed in parallel with a public set of chemotypes, i.e., the ToxPrint chemotypes, which are designed to provide excellent coverage of environmental, regulatory, and commercial-use chemical space, as well as to represent chemical patterns and properties especially relevant to various toxicity concerns. A software application, ChemoTyper has also been developed and made publicly available in order to enable chemotype searching and fingerprinting against a target structure set. The public ChemoTyper houses the ToxPrint chemotype CSRML dictionary, as well as reference implementation so that the query specifications may be adopted by other chemical structure knowledge systems. The full specifications of the XML-based CSRML standard used to express chemotypes are publicly available to facilitate and encourage the exchange of structural knowledge.

Checking the Projection Display of Multivariate Data with Colored Graphs

Projection methods such as principal component analysis (PCA), nonlinear mapping (NLM), and the self-organizing map (SOM) are valuable algorithms for visualizing multidimensional data in a two-dimensional plane. Unfortunately, the reduction of the dimensionality involves distortions. In an attempt to graphically localize the distortions of the projected data, we suggest superposing colored graphs onto the 2D plots. The color of the edges of these graphs encodes the original high-dimensional distances between the connected points. The method is applied to a cluster analysis of 37 biologically active compounds and 471 molecules represented by a structural 3D descriptor.

Combining chemoinformatics with bioinformatics: in silico prediction of bacterial flavor-forming pathways by a chemical systems biology approach "reverse pathway engineering"

The incompleteness of genome-scale metabolic models is a major bottleneck for systems biology approaches, which are based on large numbers of metabolites as identified and quantified by metabolomics. Many of the revealed secondary metabolites and/or their derivatives, such as flavor compounds, are non-essential in metabolism, and many of their synthesis pathways are unknown. In this study, we describe a novel approach, Reverse Pathway Engineering (RPE), which combines chemoinformatics and bioinformatics analyses, to predict the “missing links” between compounds of interest and their possible metabolic precursors by providing plausible chemical and/or enzymatic reactions. We demonstrate the added-value of the approach by using flavor-forming pathways in lactic acid bacteria (LAB) as an example. Established metabolic routes leading to the formation of flavor compounds from leucine were successfully replicated. Novel reactions involved in flavor formation, i.e. the conversion of alpha-hydroxy-isocaproate to 3-methylbutanoic acid and the synthesis of dimethyl sulfide, as well as the involved enzymes were successfully predicted. These new insights into the flavor-formation mechanisms in LAB can have a significant impact on improving the control of aroma formation in fermented food products. Since the input reaction databases and compounds are highly flexible, the RPE approach can be easily extended to a broad spectrum of applications, amongst others health/disease biomarker discovery as well as synthetic biology.

Applications of integrated data mining methods to exploring natural product space for acetylcholinesterase inhibitors.

Nature, especially the plant kingdom, is a rich source for novel bioactive compounds that can be used as lead compounds for drug development. In order to exploit this resource, the two neural network-based virtual screening techniques novelty detection with self-organizing maps (SOMs) and counterpropagation neural network were evaluated as tools for efficient lead structure discovery. As application scenario, significant descriptors for acetylcholinesterase (AChE) inhibitors were determined and used for model building, theoretical model validation, and virtual screening. Top-ranked virtual hits from both approaches were docked into the AChE binding site to approve the initial hits. Finally, in vitro testing of selected compounds led to the identification of forsythoside A and (+)-sesamolin as novel AChE inhibitors.

A reliable workflow for in silico assessment of genetic toxicity and application to pharmaceutical genotoxic impurities

o ChemTunes/ToxGPS o A comprehensive knowledgebase experimental toxicity data and predictive models o QSAR modeling based on biologically meaningful grouping using mechanistically selected chemotypes and molecular descriptors o Final outcome combines the evidences of QSAR models and chemotype rule-based predictions to provide good prediction performance o Robust risk assessment system providing rigorous method for quantitative weight-ofevidence o In two open challenges involving over 8,000 compounds, ToxGPS Ames mutagenicity model ranked highly for GTI relevant statistics. References

Quantifying intrinsic chemical reactivity of molecular structural features for protein binding and reactive toxicity, using the MOSES chemoinformatics system

Covalent binding of xenobiotic compounds to endogenous biomolecular sites, e.g. protein residues, leads to potentially irreversible toxic effects such as enhanced acute toxicity or skin sensitization [1]. This mechanistic knowledge provides the basis for the in silico prediction of these toxicities, as required by the EU REACH legislation and the EU Cosmetics Directive. A general toxicity prediction can be based on three consecutive steps [2]: (1.) Identification of a potential reactive protein binding mechanism via a set of molecular structural features. Those structural features can be encoded by the Chemical Subgraph Representation Markup Language (CSRML), that supports a flexible annotation of meta information, including physicochemical properties as annotations. (2.) Confirmation and quantification of (bio)chemical reactivity. The potential for a chemical to be reactive can be captured by mechanistically based QSAR models. This intrinsic reactivity is calculated rapidly by descriptors of the chemoinformatics platform Molecular Structure Encoding System (MOSES) [3]. It represents electronic, steric and resonance effects in chemical structures. The performances obtained by these reactivity models are close to time-consuming quantum chemical reactivity calculations, e.g., se = 0.44 versus 0.40 for glutathione adduct formation via Michael addition, comparing predicted values to an experimental reactivity dataset [1]. (3.) Establishing a relationship between calculated reactivity and toxicity. The predicted intrinsic reactivity values were linked to the computational prediction for different modes of toxic action, with good correlations between predicted and in vitro toxicity (up to r2=0.86). n nThe combined use of structural information and computed reactivity could assist in the non-animal based risk assessment of chemicals for regulatory purposes and in the application of integrated testing strategies (ITS). The research has received funding from the EU FP7 COSMOS Project (grant agreement n° 266835) and financing from COLIPA.

THERESA - a new reaction database-driven tool for stepwise retrosynthetic analysis

The synthesis of new compounds is a quite time consuming and cost expensive task. The need to search and evaluate alternative synthetic paths is a mandatory step before going to the lab. Searching in reaction databases may provide some information about how a compound can be synthesized but often fails if the target is not present in the database. n nThe poster presents THERESA (THE REtro-Synthesis Analyser), a novel, web-based, easy-to-use tool for the stepwise retrosynthetic analysis of a given target compound. THERESA scans reaction databases to suggest new synthetic routes and simultaneously searches in catalogs of available starting materials for the proposed precursors. Furthermore, it provides the corresponding published reaction data for each suggested synthetic step. Due to the rather general definition of reactivity, it is able to provide the chemist with new ideas for organic synthesis as it deals with a broad range and diverse chemistry, including, e.g., formation of heterocycles, pericyclic reactions, rearrangements and metathesis. n nThe poster presentation will provide insights into the general algorithms of THERESA and demonstrate its application to some simple but medicinally relevant synthetic targets.