MOLecular structure GENeration with MOLGEN, new features and future developments

Advanced Methods for Natural Products Discovery: Bioactivity Screening, Dereplication, Metabolomics Profiling, Genomic Sequencing, Databases and Informatic Tools, and Structure Elucidation

Natural Products (NP) are essential for the discovery of novel drugs and products for numerous biotechnological applications. The NP discovery process is expensive and time-consuming, having as major hurdles dereplication (early identification of known compounds) and structure elucidation, particularly the determination of the absolute configuration of metabolites with stereogenic centers. This review comprehensively focuses on recent technological and instrumental advances, highlighting the development of methods that alleviate these obstacles, paving the way for accelerating NP discovery towards biotechnological applications. Herein, we emphasize the most innovative high-throughput tools and methods for advancing bioactivity screening, NP chemical analysis, dereplication, metabolite profiling, metabolomics, genome sequencing and/or genomics approaches, databases, bioinformatics, chemoinformatics, and three-dimensional NP structure elucidation.

Sampling CASE Application for the Quality Control of Published Natural Product Structures

Structure elucidation with NMR correlation data is dicey, as there is no way to tell how ambiguous the data set is and how reliably it will define a constitution. Many different software tools for computer assisted structure elucidation (CASE) have become available over the past decades, all of which could ensure a better quality of the elucidation process, but their use is still not common. Since 2011, WebCocon has integrated the possibility to generate theoretical NMR correlation data, starting from an existing structural proposal, allowing this theoretical data then to be used for CASE. Now, WebCocon can also read the recently presented NMReDATA format, allowing for uncomplicated access to CASE with experimental data. With these capabilities, WebCocon presents itself as an easily accessible Web-Tool for the quality control of proposed new natural products. Results of this application to several molecules from literature are shown and demonstrate how CASE can contribute to improve the reliability of Structure elucidation with NMR correlation data.

Incorporation of 4J-HMBC and NOE Data into Computer-Assisted Structure Elucidation with WebCocon

Over the past decades, different software programs have been developed for the Computer-Assisted Structure Elucidation (CASE) with NMR data using with various approaches. WebCocon is one of them that has been continuously improved over the past 20 years. Here, we present the inclusion of 4J_CH correlations (4J-HMBC) in the HMBC interpretation of Cocon and NOE data in WebCocon. The 4J-HMBC data is used during the structure generation process, while the NOE data is used in post-processing of the results. The marine natural product oxocyclostylidol was selected to demonstrate WebCocon’s enhanced HMBC data processing capabilities. A systematic study of the 4J_CH correlations of oxocyclostylidol was performed. The application of NOEs in CASE is demonstrated using the NOE correlations of the diterpene pyrone asperginol A known from the literature. As a result, we obtained a conformation that corresponds very well to the existing X-ray structure.

Enumerating Substituted Benzene Isomers of Tree-Like Chemical Graphs

Enumeration of chemical structures is useful for drug design, which is one of the main targets of computational biology and bioinformatics. A chemical graph with no other cycles than benzene rings is called tree-like, and becomes a tree possibly with multiple edges if we contract each benzene ring into a single virtual atom of valence 6. All tree-like chemical graphs with a given tree representation are called the substituted benzene isomers of . When we replace each virtual atom in with a benzene ring to obtain a substituted benzene isomer, distinct isomers of are caused by the difference in arrangements of atom groups around a benzene ring. In this paper, we propose an efficient algorithm that enumerates all substituted benzene isomers of a given tree representation . Our algorithm first counts the number of all the isomers of the tree representation by a dynamic programming method. To enumerate all the isomers, for each , our algorithm then generates the th isomer by backtracking the counting phase of the dynamic programming. We also implemented our algorithm for computational experiments.

Computational exploration of the chemical structure space of possible reverse tricarboxylic acid cycle constituents

The reverse tricarboxylic acid (rTCA) cycle has been explored from various standpoints as an idealized primordial metabolic cycle. Its simplicity and apparent ubiquity in diverse organisms across the tree of life have been used to argue for its antiquity and its optimality. In 2000 it was proposed that chemoinformatics approaches support some of these views. Specifically, defined queries of the Beilstein database showed that the molecules of the rTCA are heavily represented in such compound databases. We explore here the chemical structure “space,” e.g. the set of organic compounds which possesses some minimal set of defining characteristics, of the rTCA cycle’s intermediates using an exhaustive structure generation method. The rTCA’s chemical space as defined by the original criteria and explored by our method is some six to seven times larger than originally considered. Acknowledging that each assumption in what is a defining criterion making the rTCA cycle special limits possible generative outcomes, there are many unrealized compounds which fulfill these criteria. That these compounds are unrealized could be due to evolutionary frozen accidents or optimization, though this optimization may also be for systems-level reasons, e.g., the way the pathway and its elements interface with other aspects of metabolism.

Computational Prediction of Electron Ionization Mass Spectra to Assist in GC/MS Compound Identification

We describe a tool, competitive fragmentation modeling for electron ionization (CFM-EI) that, given a chemical structure (e.g., in SMILES or InChI format), computationally predicts an electron ionization mass spectrum (EI-MS) (i.e., the type of mass spectrum commonly generated by gas chromatography mass spectrometry). The predicted spectra produced by this tool can be used for putative compound identification, complementing measured spectra in reference databases by expanding the range of compounds able to be considered when availability of measured spectra is limited. The tool extends CFM-ESI, a recently developed method for computational prediction of electrospray tandem mass spectra (ESI-MS/MS), but unlike CFM-ESI, CFM-EI can handle odd-electron ions and isotopes and incorporates an artificial neural network. Tests on EI-MS data from the NIST database demonstrate that CFM-EI is able to model fragmentation likelihoods in low-resolution EI-MS data, producing predicted spectra whose dot product scores are significantly better than full enumeration "bar-code" spectra. CFM-EI also outperformed previously reported results for MetFrag, MOLGEN-MS, and Mass Frontier on one compound identification task. It also outperformed MetFrag in a range of other compound identification tasks involving a much larger data set, containing both derivatized and nonderivatized compounds. While replicate EI-MS measurements of chemical standards are still a more accurate point of comparison, CFM-EI's predictions provide a much-needed alternative when no reference standard is available for measurement. CFM-EI is available at https://sourceforge.net/projects/cfm-id/ for download and http://cfmid.wishartlab.com as a web service.

Ring system-based chemical graph generation for de novo molecular design

Generating chemical graphs in silico by combining building blocks is important and fundamental in virtual combinatorial chemistry. A premise in this area is that generated structures should be irredundant as well as exhaustive. In this study, we develop structure generation algorithms regarding combining ring systems as well as atom fragments. The proposed algorithms consist of three parts. First, chemical structures are generated through a canonical construction path. During structure generation, ring systems can be treated as reduced graphs having fewer vertices than those in the original ones. Second, diversified structures are generated by a simple rule-based generation algorithm. Third, the number of structures to be generated can be estimated with adequate accuracy without actual exhaustive generation. The proposed algorithms were implemented in structure generator Molgilla. As a practical application, Molgilla generated chemical structures mimicking rosiglitazone in terms of a two dimensional pharmacophore pattern. The strength of the algorithms lies in simplicity and flexibility. Therefore, they may be applied to various computer programs regarding structure generation by combining building blocks.

Using fragmentation trees and mass spectral trees for identifying unknown compounds in metabolomics

Identification of unknown metabolites is the bottleneck in advancing metabolomics, leaving interpretation of metabolomics results ambiguous. The chemical diversity of metabolism is vast, making structure identification arduous and time consuming. Currently, comprehensive analysis of mass spectra in metabolomics is limited to library matching, but tandem mass spectral libraries are small compared to the large number of compounds found in the biosphere, including xenobiotics. Resolving this bottleneck requires richer data acquisition and better computational tools. Multi-stage mass spectrometry (MSn) trees show promise to aid in this regard. Fragmentation trees explore the fragmentation process, generate fragmentation rules and aid in sub-structure identification, while mass spectral trees delineate the dependencies in multi-stage MS of collision-induced dissociations. This review covers advancements over the past 10 years as a tool for metabolite identification, including algorithms, software and databases used to build and to implement fragmentation trees and mass spectral annotations.

Biotransformation of benzotriazoles: insights from transformation product identification and compound-specific isotope analysis

Benzotriazoles are widely used domestic and industrial corrosion inhibitors and have become omnipresent organic micropollutants in the aquatic environment. Here, the range of aerobic biological degradation mechanisms of benzotriazoles in activated sludge was investigated. Degradation pathways were elucidated by identifying transient and persistent transformation products in batch experiments using liquid chromatography-high-resolution tandem mass spectrometry (LC-HR-MS/MS). In addition, initial reactions were studied using compound-specific isotope analysis (CSIA). Biodegradation half-lives of 1.0 days for 1H-benzotriazole, 8.5 days for 4-methyl-1H-benzotriazole, and 0.9 days for 5-methyl-1H-benzotriazole with activated sludge confirmed their known partial persistence in conventional wastewater treatment. Major transformation products were identified as 4- and 5-hydroxy-1H-benzotriazole for the degradation of 1H-benzotriazole, and 1H-benzotriazole-5-carboxylic acid for the degradation of 5-methyl-1H-benzotriazole. These transformation products were found in wastewater effluents, showing their environmental relevance. Many other candidate transformation products, tentatively identified by interpretation of HR-MS/MS spectra, showed the broad range of possible reaction pathways including oxidation, alkylation, hydroxylation and indicate the significance of cometabolic processes for micropollutant degradation in biological wastewater treatment in general. The combination of evidence from product analysis with the significant carbon and nitrogen isotope fractionation suggests that aromatic monohydroxylation is the predominant step during the biotransformation of 1H-benzotriazole.

Small Molecule Identification with MOLGEN and Mass Spectrometry

This paper details the MOLGEN entries for the 2012 CASMI contest for small molecule identification to demonstrate structure elucidation using structure generation approaches. Different MOLGEN programs were used for different categories, including MOLGEN-MS/MS for Category 1, MOLGEN 3.5 and 5.0 for Category 2 and MOLGEN-MS for Categories 3 and 4. A greater focus is given to Categories 1 and 2, as most CASMI participants entered these categories. The settings used and the reasons behind them are described in detail, while various evaluations are used to put these results into perspective. As one author was also an organiser of CASMI, these submissions were not part of the official CASMI competition, but this paper provides an insight into how unknown identification could be performed using structure generation approaches. The approaches are semi-automated (category dependent) and benefit greatly from user experience. Thus, the results presented and discussed here may be better than those an inexperienced user could obtain with MOLGEN programs.

CASMI: And the Winner is

The Critical Assessment of Small Molecule Identification (CASMI) Contest was founded in 2012 to provide scientists with a common open dataset to evaluate their identification methods. In this review, we summarize the submissions, evaluate procedures and discuss the results. We received five submissions (three external, two internal) for LC–MS Category 1 (best molecular formula) and six submissions (three external, three internal) for LC–MS Category 2 (best molecular structure). No external submissions were received for the GC–MS Categories 3 and 4. The team of Dunn et al. from Birmingham had the most answers in the 1^st place for Category 1, while Category 2 was won by H. Oberacher. Despite the low number of participants, the external and internal submissions cover a broad range of identification strategies, including expert knowledge, database searching, automated methods and structure generation. The results of Category 1 show that complementing automated strategies with (manual) expert knowledge was the most successful approach, while no automated method could compete with the power of spectral searching for Category 2—if the challenge was present in a spectral library. Every participant topped at least one challenge, showing that different approaches are still necessary for interpretation diversity.

MolFind: a software package enabling HPLC/MS-based identification of unknown chemical structures

In this paper, we present MolFind, a highly multithreaded pipeline type software package for use as an aid in identifying chemical structures in complex biofluids and mixtures. MolFind is specifically designed for high-performance liquid chromatography/mass spectrometry (HPLC/MS) data inputs typical of metabolomics studies where structure identification is the ultimate goal. MolFind enables compound identification by matching HPLC/MS-based experimental data obtained for an unknown compound with computationally derived HPLC/MS values for candidate compounds downloaded from chemical databases such as PubChem. The downloaded "bins" consist of all compounds matching the monoisotopic molecular weight of the unknown. The computational HPLC/MS values predicted include retention index (RI), ECOM(50) (energy required to fragment 50% of a selected precursor ion), drift time, and collision induced dissociation (CID) spectrum. RI, ECOM(50), and drift-time models are used for filtering compounds downloaded from PubChem. The remaining candidates are then ranked based on CID spectra matching. Current RI and ECOM(50) models allow for the removal of about 28% of compounds from PubChem bins. Our estimates suggest that this could be improved to as much as 87% with additional chemical structures included in the computational models. Quantitative structure property relationship-based modeling of drift times showed a better correlation with experimentally determined drift times than did Mobcal cross-sectional areas. In 23 of 35 example cases, filtering PubChem bins with RI and ECOM(50) predictive models resulted in improved ranking of the unknown compounds compared to previous studies using CID spectra matching alone. In 19 of 35 examples, the correct candidate was ranked within the top 20 compounds in bins containing an average of 1635 compounds.

Application of preparative capillary gas chromatography (pcGC), automated structure generation and mutagenicity prediction to improve effect-directed analysis of genotoxicants in a contaminated groundwater

BACKGROUND, AIM AND SCOPE

The importance of groundwater for human life cannot be overemphasised. Besides fulfilling essential ecological functions, it is a major source of drinking water. However, in the industrial area of Bitterfeld, it is contaminated with a multitude of harmful chemicals, including genotoxicants. Therefore, recently developed methodologies including preparative capillary gas chromatography (pcGC), MOLGEN-MS structure generation and mutagenicity prediction were applied within effect-directed analysis (EDA) to reduce sample complexity and to identify candidate mutagens in the samples. A major focus was put on the added value of these tools compared to conventional EDA combining reversed-phase liquid chromatography (RP-LC) followed by GC/MS analysis and MS library search.

MATERIALS AND METHODS

We combined genotoxicity testing with umuC and RP-LC with pcGC fractionation to isolate genotoxic compounds from a contaminated groundwater sample. Spectral library information from the NIST05 database was combined with a computer-based structure generation tool called MOLGEN-MS for structure elucidation of unknowns. Finally, we applied a computer model for mutagenicity prediction (ChemProp) to identify candidate mutagens and genotoxicants.

RESULTS AND DISCUSSION

A total of 62 components were tentatively identified in genotoxic fractions. Ten of these components were predicted to be potentially mutagenic, whilst 2,4,6-trichlorophenol, 2,4-dichloro-6-methylphenol and 4-chlorobenzoic acid were confirmed as genotoxicants.

CONCLUSIONS AND PERSPECTIVES

The results suggest pcGC as a high-resolution fractionation tool and MOLGEN-MS to improve structure elucidation, whilst mutagenicity prediction failed in our study to predict identified genotoxicants. Genotoxicity, mutagenicity and carcinogenicity caused by chemicals are complex processes, and prediction from chemical structure still appears to be quite difficult. Progress in this field would significantly support EDA and risk assessment of environmental mixtures.

The use of MS classifiers and structure generation to assist in the identification of unknowns in effect-directed analysis

Structure generation and mass spectral classifiers have been incorporated into a new method to gain further information from low-resolution GC-MS spectra and subsequently assist in the identification of toxic compounds isolated using effect-directed fractionation. The method has been developed for the case where little analytical information other than the mass spectrum is available, common, for example, in effect-directed analysis (EDA), where further interpretation of the mass spectra is necessary to gain additional information about unknown peaks in the chromatogram. Structure generation from a molecular formula alone rapidly leads to enormous numbers of structures; hence reduction of these numbers is necessary to focus identification or confirmation efforts. The mass spectral classifiers and structure generation procedure in the program MOLGEN-MS was enhanced by including additional classifier information available from the NIST05 database and incorporation of post-generation 'filtering criteria'. The presented method can reduce the number of possible structures matching a spectrum by several orders of magnitude, creating much more manageable data sets and increasing the chance of identification. Examples are presented to show how the method can be used to provide 'lines of evidence' for the identity of an unknown compound. This method is an alternative to library search of mass spectra and is especially valuable for unknowns where no clear library match is available.

MOLGEN-CID--A canonizer for molecules and graphs accessible through the Internet

The MOLGEN Chemical Identifier MOLGEN-CID is a software module freely accessible via the Internet. For a molecule or graph entered in molfile format (2D) it produces, by a canonical renumbering procedure, a canonical molfile and a unique character string that is easily compared by computer to a similar string. The mode of operation of MOLGEN-CID is detailed and visualized with examples.

Manual Construction and Mathematics- and Computer-Aided Counting of Stereoisomers. The Example of Oligoinositols

Two methods to obtain numbers of stereoisomers and of achiral stereoisomers of a given molecular structure are detailed on the example of di- and triinositols. The first method is manual exhaustive construction free of redundance of all stereoisomers, which is rendered feasible by symmetry considerations despite the large number of isomeric triinositols (82176). The second method is counting without constructing, made possible by use of a mathematical tool, the Cauchy-Frobenius lemma, which actually is a formalized manner of considering symmetry. The results are compared to those obtained by computer-aided stereoisomer generation using the program MOLGEN 3.5. It is demonstrated that in their results all three methods agree.