PFAS-Biomolecule Interactions: Case Study Using Asclepios Nodes and Automated Workflows in KNIME for Drug Discovery and Toxicology.

The Asclepios suite of KNIME nodes represents an innovative solution for conducting cheminformatics and computational chemistry tasks, specifically tailored for applications in drug discovery and computational toxicology. This suite has been developed using open-source and publicly accessible software. In this chapter, we introduce and explore the Asclepios suite through the lens of a case study. This case study revolves around investigating the interactions between per- and polyfluorinated alkyl substances (PFAS) and biomolecules, such as nuclear receptors. The objective is to characterize the potential toxicity of PFAS and gain insights into their chemical mode of action at the molecular level. The Asclepios KNIME nodes have been designed as versatile tools capable of addressing a wide range of computational toxicology challenges. Furthermore, they can be adapted and customized to accomodate the specific needs of individual users, spanning various domains such as nanoinformatics, biomedical research, and other related applications. This chapter provides an in-depth examination of the technical underpinnings and foundations of these tools. It is accompanied by a practical case study that demonstrates the utilization of Asclepios nodes in a computational toxicology investigation. This showcases the extendable functionalities that can be applied in diverse computational chemistry contexts. By the end of this chapter, we aim for readers to have a comprehensive understanding of the effectiveness of the Asclepios node functions. These functions hold significant potential for enhancing a wide spectrum of cheminformatics applications.

Data mining of PubChem bioassay records reveals diverse OXPHOS inhibitory chemotypes as potential therapeutic agents against ovarian cancer.

Focused screening on target-prioritized compound sets can be an efficient alternative to high throughput screening (HTS). For most biomolecular targets, compound prioritization models depend on prior screening data or a target structure. For phenotypic or multi-protein pathway targets, it may not be clear which public assay records provide relevant data. The question also arises as to whether data collected from disparate assays might be usefully consolidated. Here, we report on the development and application of a data mining pipeline to examine these issues. To illustrate, we focus on identifying inhibitors of oxidative phosphorylation, a druggable metabolic process in epithelial ovarian tumors. The pipeline compiled 8415 available OXPHOS-related bioassays in the PubChem data repository involving 312,093 unique compound records. Application of PubChem assay activity annotations, PAINS (Pan Assay Interference Compounds), and Lipinski-like bioavailability filters yields 1852 putative OXPHOS-active compounds that fall into 464 clusters. These chemotypes are diverse but have relatively high hydrophobicity and molecular weight but lower complexity and drug-likeness. These chemotypes show a high abundance of bicyclic ring systems and oxygen containing functional groups including ketones, allylic oxides (alpha/beta unsaturated carbonyls), hydroxyl groups, and ethers. In contrast, amide and primary amine functional groups have a notably lower than random prevalence. UMAP representation of the chemical space shows strong divergence in the regions occupied by OXPHOS-inactive and -active compounds. Of the six compounds selected for biological testing, 4 showed statistically significant inhibition of electron transport in bioenergetics assays. Two of these four compounds, lacidipine and esbiothrin, increased in intracellular oxygen radicals (a major hallmark of most OXPHOS inhibitors) and decreased the viability of two ovarian cancer cell lines, ID8 and OVCAR5. Finally, data from the pipeline were used to train random forest and support vector classifiers that effectively prioritized OXPHOS inhibitory compounds within a held-out test set (ROCAUC 0.962 and 0.927, respectively) and on another set containing 44 documented OXPHOS inhibitors outside of the training set (ROCAUC 0.900 and 0.823). This prototype pipeline is extensible and could be adapted for focus screening on other phenotypic targets for which sufficient public data are available.Scientific contributionHere, we describe and apply an assay data mining pipeline to compile, process, filter, and mine public bioassay data. We believe the procedure may be more broadly applied to guide compound selection in early-stage hit finding on novel multi-protein mechanistic or phenotypic targets. To demonstrate the utility of our approach, we apply a data mining strategy on a large set of public assay data to find drug-like molecules that inhibit oxidative phosphorylation (OXPHOS) as candidates for ovarian cancer therapies.

Workflow for predictive risk assessments of UVCBs: cheminformatics library design, QSAR, and read-across approaches applied to complex mixtures of metal naphthenates.

Substances of unknown or variable composition, complex reaction products, and biological materials (UVCBs) are commonly found in the environment. However, assessing their human toxicological risk is challenging due to their variable composition and many constituents. Metal naphthenate salts are one such category of UVCBs that are the reaction products of naphthenic acids with metals to form complex mixtures. Metal naphthenates are often found or used in household and industrial materials with potential for human exposure, but very few of these materials have been evaluated for causing human health hazards. Herein, we evaluate metal naphthenates using predictions derived from read-across and quantitative structure-activity/property relationship (QSAR/QSPR) models. Accordingly, we first built a computational chemistry library by enumerating the structures of naphthenic acids and derived 11,850 QSAR-acceptable structures; then, we used open and commercial in silico tools on these structures to predict a set of physicochemical properties and toxicity endpoints. We then compared the QSAR/QSPR predictions with available experimental data on naphthenic acids to provide a more complete picture of the contributions of the components to the toxicity profiles of metal naphthenate mixtures. The available systematic acute oral toxicity values (LD50) and QSAR LD50 predictions of all the naphthenic acid components indicated low concern for toxicity. The point of departure predictions for chronic repeated dose toxicity for the naphthenic acid components using QSAR models developed from studies on rats ranged from 25 to 50 mg/kg/day. These values are in good agreement with findings from studies on copper and zinc naphthenates, which had no observed adverse effect levels of 30 and 118 mg/kg/day, respectively. Hence, this study demonstrates how published in silico approaches can be used to identify the potential components of metal naphthenates for further testing, inform groupings of UVCBs such as naphthenates, as well as fill the data gaps using read-across and QSAR models to inform risk assessment.

Navigating a 1E+60 Chemical Space of Peptide/Peptoid Oligomers.

Herein we report a virtual library of 1E+60 members, a common estimate for the size of the drug-like chemical space. The library consists of linear or cyclic oligomers forming molecules within the size range of peptide drugs. We demonstrate ligand-based virtual screening using a genetic algorithm.

Protein language models are performant in structure-free virtual screening.

Hitherto virtual screening (VS) has been typically performed using a structure-based drug design paradigm. Such methods typically require the use of molecular docking on high-resolution three-dimensional structures of a target protein-a computationally-intensive and time-consuming exercise. This work demonstrates that by employing protein language models and molecular graphs as inputs to a novel graph-to-transformer cross-attention mechanism, a screening power comparable to state-of-the-art structure-based models can be achieved. The implications thereof include highly expedited VS due to the greatly reduced compute required to run this model, and the ability to perform early stages of computer-aided drug design in the complete absence of 3D protein structures.

Impacts of Plu kaow (Houttuynia cordata Thunb.) Ethanolic Extract on Diabetes and Dyslipidemia in STZ Induced Diabetic Rats: Phytochemical Profiling, Cheminformatics Analyses, and Molecular Docking Studies.

The increasing prevalence of diabetes and dyslipidemia poses significant health challenges, impacting millions of people globally and leading to high rates of illness and death. This study aimed to explore the potential antidiabetic and hypolipidemic effects of Plu kaow (Houttuynia cordata Thunb.) ethanolic extract (PK) in streptozotocin (STZ) induced diabetic rats, focusing on its molecular mechanisms. Diabetes was induced in fasting Long Evans rats using streptozotocin (65 mg/kg b. w.), with glibenclamide (5 mg/kg/day) used as the standard experimental drug. The treated groups received oral supplementation of PK (500 mg/kg/day) for 28 days. The study evaluated blood glucose levels, lipid status, body weight, liver, kidney, and heart function biomarkers, antioxidant activity, and histological examination of various organs. Additionally, untargeted metabolomics, cheminformatics, and molecular docking were employed to elucidate the probable mechanisms of action of PK. Based on metabolomic profiling data, the PK was found to contain various putative antidiabetic agents such as kaempferol 7-neohesperidoside, isochlorogenic acid C, rutin, datiscin, and diosmin and they have been proposed to significantly (p < 0.001) reduce blood glucose levels and modulated hyperlipidemia. PK also improved the tested liver, kidney, and heart function biomarkers and reversed damage to normal pancreatic, liver, kidney, and heart cells in histological analysis. In conclusion, PK shows promise as a potential treatment or management option for diabetes and hyperlipidemia, as well as their associated complications in diabetic rats.

Cheminformatic Identification of Tyrosyl-DNA Phosphodiesterase 1 (Tdp1) Inhibitors: A Comparative Study of SMILES-Based Supervised Machine Learning Models.

BACKGROUND: Tyrosyl-DNA phosphodiesterase 1 (Tdp1) repairs damages in DNA induced by abortive topoisomerase 1 activity; however, maintenance of genetic integrity may sustain cellular division of neoplastic cells. It follows that Tdp1-targeting chemical inhibitors could synergize well with existing chemotherapy drugs to deny cancer growth; therefore, identification of Tdp1 inhibitors may advance precision medicine in oncology. OBJECTIVE: Current computational research efforts focus primarily on molecular docking simulations, though datasets involving three-dimensional molecular structures are often hard to curate and computationally expensive to store and process. We propose the use of simplified molecular input line entry system (SMILES) chemical representations to train supervised machine learning (ML) models, aiming to predict potential Tdp1 inhibitors. METHODS: An open-sourced consensus dataset containing the inhibitory activity of numerous chemicals against Tdp1 was obtained from Kaggle. Various ML algorithms were trained, ranging from simple algorithms to ensemble methods and deep neural networks. For algorithms requiring numerical data, SMILES were converted to chemical descriptors using RDKit, an open-sourced Python cheminformatics library. RESULTS: Out of 13 optimized ML models with rigorously tuned hyperparameters, the random forest model gave the best results, yielding a receiver operating characteristics-area under curve of 0.7421, testing accuracy of 0.6815, sensitivity of 0.6444, specificity of 0.7156, precision of 0.6753, and F1 score of 0.6595. CONCLUSIONS: Ensemble methods, especially the bootstrap aggregation mechanism adopted by random forest, outperformed other ML algorithms in classifying Tdp1 inhibitors from non-inhibitors using SMILES. The discovery of Tdp1 inhibitors could unlock more treatment regimens for cancer patients, allowing for therapies tailored to the patient's condition.

Leveraging In Silico Structure-Activity Models to Predict Acute Honey Bee (Apis mellifera) Toxicity for Agrochemicals.

In the realm of crop protection products, ensuring the safety of pollinators stands as a pivotal aspect of advancing sustainable solutions. Extensive research has been dedicated to this crucial topic as well as new approach methodologies in toxicity testing. Hence, within the agricultural and chemical industries, prioritizing pollinator safety remains a constant objective during the development of predictive tools. One of these tools includes computational models like quantitative structure-activity relationships (QSARs) that are valuable in predicting the toxicity of chemicals. This research uses bee toxicity data to develop artificial neural network classification models for predicting honey bee acute toxicity. Bee toxicity data from 1542 compounds were used to develop models; the sensitivity and specificity of the best model were 0.90 and 0.91, respectively. These in silico models can aid in the discovery of next-generation crop protection products. These tools can guide the screening and selection of next-generation crop protection molecules with high margins of safety to pollinators, and candidates with favorable sustainability profiles can be identified at the early discovery stage as precursors to in vivo data generation.

Predicting ELNES/XANES spectra by machine learning with an atomic coordinate-independent descriptor and its application to ground-state electronic structures.

ELNES/XANES spectra can be observed using TEM or synchrotron radiation and can elucidate the unoccupied state electronic structures of an excited states. The computation of their features is usually demanding substantial computational resources due to the requisite structure optimization and electronic structure calculations. Herein, we leverage a machine learning technique alongside an atomic-coordinate-independent descriptor, SMILES, to yield the ELNES/XANES spectra, directly, with heightened precision. Moreover, our approach extends to obtain ground state electronic structure, namely PDOS at both occupied and unoccupied ground states, underscoring its viability for a ground-state spectroscopy. Our study revealed that incorporation of long-SMILES molecules into the training dataset enhances prediction accuracy for such molecular structures. This study's direct derivation of spectroscopy from SMILES strings holds promise for expediting spectroscopic inquiries.

Development of Drug Discovery Platforms Using Artificial Intelligence and Cheminformatics.

Recently, remarkable progress has been achieved in artificial intelligence (AI), including machine learning. Various AI models have been proposed for drug discovery, including the design of small molecules, activity prediction, and three-dimensional (3D) structure prediction of proteins. AI consists of diverse elements, including information retrieval and machine learning, and can be used in a wide range of drug discovery scenarios. In this review, we focused on AI for small-molecule drug discovery with respect to molecular design, activity prediction, and prediction of the binding poses of compounds to target molecules. We also discussed the applications of AI in academic drug discovery.

BioBricks.ai: A Versioned Data Registry for Life Sciences Data Assets.

Researchers in biomedical research, public health and the life sciences often spend weeks or months discovering, accessing, curating, and integrating data from disparate sources, significantly delaying the onset of actual analysis and innovation. Instead of countless developers creating redundant and inconsistent data pipelines, BioBricks.ai offers a centralized data repository and a suite of developer-friendly tools to simplify access to scientific data. Currently, BioBricks.ai delivers over ninety biological and chemical datasets. It provides a package manager-like system for installing and managing dependencies on data sources. Each 'brick' is a Data Version Control git repository that supports an updateable pipeline for extraction, transformation, and loading data into the BioBricks.ai backend at https://biobricks.ai. Use cases include accelerating data science workflows and facilitating the creation of novel data assets by integrating multiple datasets into unified, harmonized resources. In conclusion, BioBricks.ai offers an opportunity to accelerate access and use of public data through a single open platform.

Cheminformatics analysis of indoleamine and tryptophan 2,3-dioxygenase inhibitors: A descriptor and fingerprint based machine learning approach to disclose selectivity measures.

Indoleamine 2,3-dioxygenase (IDO) and tryptophan 2,3-dioxygenase (TDO) are attractive drug targets for cancer immunotherapy. After disappointing results of the epacadostat as a selective IDO inhibitor in phase III clinical trials, there is much interest in the development of the TDO selective inhibitors. In the current study, several data analysis methods and machine learning approaches including logistic regression, Random Forest, XGBoost and Support Vector Machines were used to model a data set of compounds retrieved from ChEMBL. Models based on the Morgan fingerprints revealed notable fragments for the selective inhibition of the IDO, TDO or both. Multiple fragment docking was performed to find the best set of bound fragments and their orientation in the space for efficient linking. Linking the fragments and optimization of the final molecules were accomplished by means of an artificial intelligence generative framework. Finally, selectivity of the optimized molecules was assessed and the top 4 lead molecules were filtered through PAINS, Brenk and NIH filters. Results indicated that phenyloxalamide, fluoroquinoline, and 3-bromo-4-fluroaniline confer selectivity towards the IDO inhibition. Correspondingly, 1-benzyl-1H-naphtho[2,3-d][1,2,3]triazole-4,9-dione was found to be an integral fragment for the selective inhibition of the TDO by constituting a coordination bond with the Fe atom of heme. In addition, furo[2,3-c]pyridine-2,3-diamine was found as a common fragment for inhibition of the both targets and can be used in the design of the dual target inhibitors of the IDO and TDO. The new fragments introduced here can be a useful building blocks for incorporation into the selective TDO or dual IDO/TDO inhibitors.

Cheminformatics-Guided Exploration of Synthetic Marine Natural Product-Inspired Brominated Indole-3-Glyoxylamides and Their Potentials for Drug Discovery.

Marine natural products (MNPs) continue to be tested primarily in cellular toxicity assays, both mammalian and microbial, despite most being inactive at concentrations relevant to drug discovery. These MNPs become missed opportunities and represent a wasteful use of precious bioresources. The use of cheminformatics aligned with published bioactivity data can provide insights to direct the choice of bioassays for the evaluation of new MNPs. Cheminformatics analysis of MNPs found in MarinLit (n = 39,730) up to the end of 2023 highlighted indol-3-yl-glyoxylamides (IGAs, n = 24) as a group of MNPs with no reported bioactivities. However, a recent review of synthetic IGAs highlighted these scaffolds as privileged structures with several compounds under clinical evaluation. Herein, we report the synthesis of a library of 32 MNP-inspired brominated IGAs (25-56) using a simple one-pot, multistep method affording access to these diverse chemical scaffolds. Directed by a meta-analysis of the biological activities reported for marine indole alkaloids (MIAs) and synthetic IGAs, the brominated IGAs 25-56 were examined for their potential bioactivities against the Parkinson's Disease amyloid protein alpha synuclein (α-syn), antiplasmodial activities against chloroquine-resistant (3D7) and sensitive (Dd2) parasite strains of Plasmodium falciparum, and inhibition of mammalian (chymotrypsin and elastase) and viral (SARS-CoV-2 3CLpro) proteases. All of the synthetic IGAs tested exhibited binding affinity to the amyloid protein α-syn, while some showed inhibitory activities against P. falciparum, and the proteases, SARS-CoV-2 3CLpro, and chymotrypsin. The cellular safety of the IGAs was examined against cancerous and non-cancerous human cell lines, with all of the compounds tested inactive, thereby validating cheminformatics and meta-analyses results. The findings presented herein expand our knowledge of marine IGA bioactive chemical space and advocate expanding the scope of biological assays routinely used to investigate NP bioactivities, specifically those more suitable for non-toxic compounds. By integrating cheminformatics tools and functional assays into NP biological testing workflows, we can aim to enhance the potential of NPs and their scaffolds for future drug discovery and development.

Multifaceted targets of cannabidiol in epilepsy: Modulating glutamate signaling and beyond.

Cannabidiol has been reported to interact with broad-spectrum biological targets with pleiotropic pharmacology including epilepsy although a cohesive mechanism is yet to be determined. Even though some studies propose that cannabidiol may manipulate glutamatergic signals, there is insufficient evidence to support cannabidiol direct effect on glutamate signaling, which is important in intervening epilepsy. Therefore, the present study aimed to analyze the epilepsy-related targets for cannabidiol, assess the differentially expressed genes with its treatment, and identify the possible glutamatergic signaling target. In this study, the epileptic protein targets of cannabidiol were identified using the Tanimoto coefficient and similarity index-based targets fishing which were later overlapped with the altered expression, epileptic biomarkers, and genetically altered proteins in epilepsy. The common proteins were then screened for possible glutamatergic signaling targets with differentially expressed genes. Later, molecular docking and simulation were performed using AutoDock Vina and GROMACS to evaluate binding affinity, ligand-protein stability, hydrophilic interaction, protein compactness, etc. Cannabidiol identified 30 different epilepsy-related targets of multiple protein classes including G-protein coupled receptors, enzymes, ion channels, etc. Glutamate receptor 2 was identified to be genetically varied in epilepsy which was targeted by cannabidiol and its expression was increased with its treatment. More importantly, cannabidiol showed a direct binding affinity with Glutamate receptor 2 forming a stable hydrophilic interaction and comparatively lower root mean squared deviation and residual fluctuations, increasing protein compactness with broad conformational changes. Based on the cheminformatic target fishing, evaluation of differentially expressed genes, molecular docking, and simulations, it can be hypothesized that cannabidiol may possess glutamate receptor 2-mediated anti-epileptic activities.

Data-science-guided calibration curve prediction of an MLCT-based ee determination assay for chiral amines.

Circular dichroism (CD) based enantiomeric excess (ee) determination assays are optical alternatives to chromatographic ee determination in high-throughput screening (HTS) applications. However, the implementation of these assays requires calibration experiments using enantioenriched materials. We present a data-driven approach that circumvents the need for chiral resolution and calibration experiments for an octahedral Fe(II) complex (1) used for the ee determination of α-chiral primary amines. By computationally parameterizing the imine ligands formed in the assay conditions, a model of the circular dichroism (CD) response of the Fe(II) assembly was developed. Using this model, calibration curves were generated for four analytes and compared to experimentally generated curves. In a single-blind ee determination study, the ee values of unknown samples were determined within 9% mean absolute error, which rivals the error using experimentally generated calibration curves.

Synthesis and Cheminformatics-Directed Antibacterial Evaluation of Echinosulfonic Acid-Inspired Bis-Indole Alkaloids.

Synthetic efforts toward complex natural product (NP) scaffolds are useful ones, particularly those aimed at expanding their bioactive chemical space. Here, we utilised an orthogonal cheminformatics-based approach to predict the potential biological activities for a series of synthetic bis-indole alkaloids inspired by elusive sponge-derived NPs, echinosulfone A (1) and echinosulfonic acids A-D (2-5). Our work includes the first synthesis of desulfato-echinosulfonic acid C, an α-hydroxy bis(3'-indolyl) alkaloid (17), and its full NMR characterisation. This synthesis provides corroborating evidence for the structure revision of echinosulfonic acids A-C. Additionally, we demonstrate a robust synthetic strategy toward a diverse range of α-methine bis(3'-indolyl) acids and acetates (11-16) without the need for silica-based purification in either one or two steps. By integrating our synthetic library of bis-indoles with bioactivity data for 2048 marine indole alkaloids (reported up to the end of 2021), we analyzed their overlap with marine natural product chemical diversity. Notably, the C-6 dibrominated α-hydroxy bis(3'-indolyl) and α-methine bis(3'-indolyl) analogues (11, 14, and 17) were found to contain significant overlap with antibacterial C-6 dibrominated marine bis-indoles, guiding our biological evaluation. Validating the results of our cheminformatics analyses, the dibrominated α-methine bis(3'-indolyl) alkaloids (11, 12, 14, and 15) were found to exhibit antibacterial activities against methicillin-sensitive and -resistant Staphylococcus aureus. Further, while investigating other synthetic approaches toward bis-indole alkaloids, 16 incorrectly assigned synthetic α-hydroxy bis(3'-indolyl) alkaloids were identified. After careful analysis of their reported NMR data, and comparison with those obtained for the synthetic bis-indoles reported herein, all of the structures have been revised to α-methine bis(3'-indolyl) alkaloids.

Exploring Tau Fibril-Disaggregating and Antioxidating Molecules Binding to Membrane-Bound Amyloid Oligomers Using Machine Learning-Enhanced Docking and Molecular Dynamics.

Intracellular tau fibrils are sources of neurotoxicity and oxidative stress in Alzheimer's. Current drug discovery efforts have focused on molecules with tau fibril disaggregation and antioxidation functions. However, recent studies suggest that membrane-bound tau-containing oligomers (mTCOs), smaller and less ordered than tau fibrils, are neurotoxic in the early stage of Alzheimer's. Whether tau fibril-targeting molecules are effective against mTCOs is unknown. The binding of epigallocatechin-3-gallate (EGCG), CNS-11, and BHT-CNS-11 to in silico mTCOs and experimental tau fibrils was investigated using machine learning-enhanced docking and molecular dynamics simulations. EGCG and CNS-11 have tau fibril disaggregation functions, while the proposed BHT-CNS-11 has potential tau fibril disaggregation and antioxidation functions like EGCG. Our results suggest that the three molecules studied may also bind to mTCOs. The predicted binding probability of EGCG to mTCOs increases with the protein aggregate size. In contrast, the predicted probability of CNS-11 and BHT-CNS-11 binding to the dimeric mTCOs is higher than binding to the tetrameric mTCOs for the homo tau but not for the hetero tau-amylin oligomers. Our results also support the idea that anionic lipids may promote the binding of molecules to mTCOs. We conclude that tau fibril-disaggregating and antioxidating molecules may bind to mTCOs, and that mTCOs may also be useful targets for Alzheimer's drug design.

Phenolic profile, cheminformatics, and antiplatelet aggregation activity of orange and purple sweet potato (Ipomoea batatas L.) storage roots.

Sweet potatoes are rich in cardioprotective phytochemicals with potential anti-platelet aggregation activity, although this benefit may vary among cultivars/genotypes. The phenolic profile [HPLC-ESI(-)-qTOF-MS2], cheminformatics (ADMET properties, affinity toward platelet proteins) and anti-PA activity of phenolic-rich hydroalcoholic extracts obtained from orange (OSP) and purple (PSP) sweet potato storage roots, was evaluated. The phenolic richness [Hydroxycinnamic acids> flavonoids> benzoic acids] was PSP > OSP. Their main chlorogenic acids could interact with platelet proteins (integrins/adhesins, kinases/metalloenzymes) but their bioavailability could be poor. Just OSP exhibited a dose-dependent anti-platelet aggregation activity [inductor (IC50, mg.ml-1): thrombin receptor activator peptide-6 (0.55) > Adenosine-5'-diphosphate (1.02) > collagen (1.56)] and reduced P-selectin expression (0.75-1.0 mg.ml-1) but not glycoprotein IIb/IIIa secretion. The explored anti-PA activity of OSP/PSP seems to be inversely related to their phenolic richness. The poor first-pass bioavailability of its chlorogenic acids (documented in silico) may represent a further obstacle for their anti-PA in vivo.

Toward the Reconciliation of Inconsistent Molecular Structures from Biochemical Databases.

Information on the structure of molecules, retrieved via biochemical databases, plays a pivotal role in various disciplines, including metabolomics, systems biology, and drug discovery. No such database can be complete and it is often necessary to incorporate data from several sources. However, the molecular structure for a given compound is not necessarily consistent between databases. This article presents StructRecon, a novel tool for resolving unique molecular structures from database identifiers. Currently, identifiers from BiGG, ChEBI, Escherichia coli Metabolome Database (ECMDB), MetaNetX, and PubChem are supported. StructRecon traverses the cross-links between entries in different databases to construct what we call identifier graphs. The goal of these graphs is to offer a more complete view of the total information available on a given compound across all the supported databases. To reconcile discrepancies met during the traversal of the databases, we develop an extensible model for molecular structure supporting multiple independent levels of detail, which allows standardization of the structure to be applied iteratively. In some cases, our standardization approach results in multiple candidate structures for a given compound, in which case a random walk-based algorithm is used to select the most likely structure among incompatible alternatives. As a case study, we applied StructRecon to the EColiCore2 model. We found at least one structure for 98.66% of its compounds, which is more than twice as many as possible when using the databases in more standard ways not considering the complex network of cross-database references captured by our identifier graphs. StructRecon is open-source and modular, which enables support for more databases in the future.