Accurate Specification of Molecular Structures: The Case for Zero-Order Bonds and Explicit Hydrogen Counting

PaCh (Packed Chemicals): Computationally Effective Binary Format for Chemical Structure Encoding

In this work, we propose a versatile molecule and reaction encoding binary data format that aims to bridge the gap between the advantages of SMILES, like local stereo- and implicit hydrogen encoding, and block-structured MDL MOL with a 2D layout and explicit bond encoding, while addressing their respective limitations. Our new format introduces a balance between size efficiency, processing speed, and comprehensive representation, making it well-suited for various applications in cheminformatics, including deep learning, data storage, and searching. By offering an explicit approach to store atom connectivity (including implicit hydrogens), electronic state, stereochemistry, and other crucial molecular attributes, our proposal seeks to enhance data storage efficiency and promote interoperability among different software tools.

A workflow for deriving chemical entities from crystallographic data and its application to the Crystallography Open Database

Knowledge about the 3-dimensional structure, orientation and interaction of chemical compounds is important in many areas of science and technology. X-ray crystallography is one of the experimental techniques capable of providing a large amount of structural information for a given compound, and it is widely used for characterisation of organic and metal-organic molecules. The method provides precise 3D coordinates of atoms inside crystals, however, it does not directly deliver information about certain chemical characteristics such as bond orders, delocalization, charges, lone electron pairs or lone electrons. These aspects of a molecular model have to be derived from crystallographic data using refined information about interatomic distances and atom types as well as employing general chemical knowledge. This publication describes a curated automatic pipeline for the derivation of chemical attributes of molecules from crystallographic models. The method is applied to build a catalogue of chemical entities in an open-access crystallographic database, the Crystallography Open Database (COD). The catalogue of such chemical entities is provided openly as a derived database. The content of this catalogue and the problems arising in the fully automated pipeline are discussed, along with the possibilities to introduce manual data curation into the process.

Graph isomorphism-based algorithm for cross-checking chemical and crystallographic descriptions

Published reports of chemical compounds often contain multiple machine-readable descriptions which may supplement each other in order to yield coherent and complete chemical representations. This publication presents a method to cross-check such descriptions using a canonical representation and isomorphism of molecular graphs. If immediate agreement between compound descriptions is not found, the algorithm derives the minimal set of simplifications required for both descriptions to arrive to a matching form (if any). The proposed algorithm is used to cross-check chemical descriptions from the Crystallography Open Database to identify coherently described entries as well as those requiring further curation.

SELFIES and the future of molecular string representations

Artificial intelligence (AI) and machine learning (ML) are expanding in popularity for broad applications to challenging tasks in chemistry and materials science. Examples include the prediction of properties, the discovery of new reaction pathways, or the design of new molecules. The machine needs to read and write fluently in a chemical language for each of these tasks. Strings are a common tool to represent molecular graphs, and the most popular molecular string representation, Smiles, has powered cheminformatics since the late 1980s. However, in the context of AI and ML in chemistry, Smiles has several shortcomings-most pertinently, most combinations of symbols lead to invalid results with no valid chemical interpretation. To overcome this issue, a new language for molecules was introduced in 2020 that guarantees 100% robustness: SELF-referencing embedded string (Selfies). Selfies has since simplified and enabled numerous new applications in chemistry. In this perspective, we look to the future and discuss molecular string representations, along with their respective opportunities and challenges. We propose 16 concrete future projects for robust molecular representations. These involve the extension toward new chemical domains, exciting questions at the interface of AI and robust languages, and interpretability for both humans and machines. We hope that these proposals will inspire several follow-up works exploiting the full potential of molecular string representations for the future of AI in chemistry and materials science.

TUCAN: A molecular identifier and descriptor applicable to the whole periodic table from hydrogen to oganesson

TUCAN is a canonical serialization format that is independent of domain-specific concepts of structure and bonding. The atomic number is the only chemical feature that is used to derive the TUCAN format. Other than that, the format is solely based on the molecular topology. Validation is reported on a manually curated test set of molecules as well as a library of non-chemical graphs. The serialization procedure generates a canonical “tuple-style” output which is bidirectional, allowing the TUCAN string to serve as both identifier and descriptor. Use of the Python NetworkX graph library facilitated a compact and easily extensible implementation.

Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities

Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.

Towards operando computational modeling in heterogeneous catalysis

An increased synergy between experimental and theoretical investigations in heterogeneous catalysis has become apparent during the last decade. Experimental work has extended from ultra-high vacuum and low temperature towards operando conditions. These developments have motivated the computational community to move from standard descriptive computational models, based on inspection of the potential energy surface at 0 K and low reactant concentrations (0 K/UHV model), to more realistic conditions. The transition from 0 K/UHV to operando models has been backed by significant developments in computer hardware and software over the past few decades. New methodological developments, designed to overcome part of the gap between 0 K/UHV and operando conditions, include (i) global optimization techniques, (ii) ab initio constrained thermodynamics, (iii) biased molecular dynamics, (iv) microkinetic models of reaction networks and (v) machine learning approaches. The importance of the transition is highlighted by discussing how the molecular level picture of catalytic sites and the associated reaction mechanisms changes when the chemical environment, pressure and temperature effects are correctly accounted for in molecular simulations. It is the purpose of this review to discuss each method on an equal footing, and to draw connections between methods, particularly where they may be applied in combination.

PubChem chemical structure standardization

BACKGROUND

PubChem is a chemical information repository, consisting of three primary databases: Substance, Compound, and BioAssay. When individual data contributors submit chemical substance descriptions to Substance, the unique chemical structures are extracted and stored into Compound through an automated process called structure standardization. The present study describes the PubChem standardization approaches and analyzes them for their success rates, reasons that cause structures to be rejected, and modifications applied to structures during the standardization process. Furthermore, the PubChem standardization is compared to the structure normalization of the IUPAC International Chemical Identifier (InChI) software, as manifested by conversion of the InChI back into a chemical structure.

RESULTS

The observed rejection rate for substances processed by PubChem standardization was 0.36%, which is predominantly attributed to structures with invalid atom valences that cannot be readily corrected without additional information from contributors. Of all structures that pass standardization, 44% are modified in the process, reducing the count of unique structures from 53,574,724 in substance to 45,808,881 in compound as identified by de-aromatized canonical isomeric SMILES. Even though the processing time is very low on average (only 0.4% of structures have individual standardization time above 0.1 s), total standardization time is completely dominated by edge cases: 90% of the time to standardize all structures in PubChem substance is spent on the 2.05% of structures with the highest individual standardization time. It is worth noting that 60% of the structures obtained from PubChem structure standardization are not identical to the chemical structure resulting from the InChI (primarily due to preferences for a different tautomeric form).

CONCLUSIONS

Standardization of chemical structures is complicated by the diversity of chemical information and their representations approaches. The PubChem standardization is an effective and efficient tool to account for molecular diversity and to eliminate invalid/incomplete structures. Further development will concentrate on improved tautomer consideration and an expanded stereocenter definition. Modifications are difficult to thoroughly validate, with slight changes often affecting many thousands of structures and various edge cases. The PubChem structure standardization service is accessible as a public resource ( https://pubchem.ncbi.nlm.nih.gov/standardize ), and via programmatic interfaces.

A refined QSAR model for prediction of chemical asthma hazard

BACKGROUND

A previously developed quantitative structure-activity relationship (QSAR) model has been extern ally validated as a good predictor of chemical asthma hazard (sensitivity: 79-86%, specificity: 93-99%).

AIMS

To develop and validate a second version of this model.

METHODS

Learning dataset asthmagenic chemicals with molecular weight (MW) <1 kDa were identified from reports published in the peer-reviewed literature before the end of 2012. Control chemicals for which no reported case(s) of occupational asthma had been identified were selected at random from UK and US occupational exposure limit tables. MW banding was used in an attempt to categorically match the control group for MW distribution of the asthmagens. About 10% of chemicals in each MW category were excluded for use as an external validation set. An independent researcher utilized a logistic regression approach to compare the molecular descriptors present in asthmagens and controls. The resulting equation generated a hazard index (HI), with a value between zero and one, as an estimate of the probability that the chemical had asthmagenic potential. The HI was determined for each compound in the external validation set.

RESULTS

The model development sets comprised 99 chemical asthmagens and 204 controls. The external validation showed that using a cut-point HI of 0.39, 9/10 asthmagenic (sensitivity: 90%) and 23/24 non-asthmagenic (specificity: 96%) compounds were correctly predicted. The new QSAR model showed a better receiver operating characteristic plot than the original.

CONCLUSIONS

QSAR refinement by iteration has resulted in an improved model for the prediction of chemical asthma hazard.

Machines first, humans second: on the importance of algorithmic interpretation of open chemistry data

The current rise in the use of open lab notebook techniques means that there are an increasing number of scientists who make chemical information freely and openly available to the entire community as a series of micropublications that are released shortly after the conclusion of each experiment. We propose that this trend be accompanied by a thorough examination of data sharing priorities. We argue that the most significant immediate benefactor of open data is in fact chemical algorithms, which are capable of absorbing vast quantities of data, and using it to present concise insights to working chemists, on a scale that could not be achieved by traditional publication methods. Making this goal practically achievable will require a paradigm shift in the way individual scientists translate their data into digital form, since most contemporary methods of data entry are designed for presentation to humans rather than consumption by machine learning algorithms. We discuss some of the complex issues involved in fixing current methods, as well as some of the immediate benefits that can be gained when open data is published correctly using unambiguous machine readable formats.

Graphical Abstract

New target prediction and visualization tools incorporating open source molecular fingerprints for TB Mobile 2.0

BACKGROUND

We recently developed a freely available mobile app (TB Mobile) for both iOS and Android platforms that displays Mycobacterium tuberculosis (Mtb) active molecule structures and their targets with links to associated data. The app was developed to make target information available to as large an audience as possible.

RESULTS

We now report a major update of the iOS version of the app. This includes enhancements that use an implementation of ECFP_6 fingerprints that we have made open source. Using these fingerprints, the user can propose compounds with possible anti-TB activity, and view the compounds within a cluster landscape. Proposed compounds can also be compared to existing target data, using a näive Bayesian scoring system to rank probable targets. We have curated an additional 60 new compounds and their targets for Mtb and added these to the original set of 745 compounds. We have also curated 20 further compounds (many without targets in TB Mobile) to evaluate this version of the app with 805 compounds and associated targets.

CONCLUSIONS

TB Mobile can now manage a small collection of compounds that can be imported from external sources, or exported by various means such as email or app-to-app inter-process communication. This means that TB Mobile can be used as a node within a growing ecosystem of mobile apps for cheminformatics. It can also cluster compounds and use internal algorithms to help identify potential targets based on molecular similarity. TB Mobile represents a valuable dataset, data-visualization aid and target prediction tool.

Rendering Molecular Sketches for Publication Quality Output

The creation of 2D molecular structure diagrams that make full use of the capabilities of modern display systems, using only input data expressed in file formats used for cheminformatics, is a complex task that requires a number of additional algorithms. Assuming that atom positions have been well chosen, the rendering engine is required to micromanage the precise positioning of atom labels, bonds and atom adjuncts, in such a way that the final output is correct, consistent with convention, and as pleasing to the eye as a diagram produced by a graphic designer. The techniques must be equally applicable when creating output for low-resolution screens and high resolution printed output, and make use of contemporary graphics file formats in such a way that the largest possible number of software platforms are able to display the output at any resolution without degradation or inconsistency. The main issues involved in meeting these criteria are discussed, and algorithms for satisfying them are presented.

Redefining Cheminformatics with Intuitive Collaborative Mobile Apps

The proliferation of mobile devices such as smartphones and tablet computers has recently been extended to include a growing ecosystem of increasingly sophisticated chemistry software packages, commonly known as apps. The capabilities that these apps can offer to the practicing chemist are approaching those of conventional desktop-based software, but apps tend to be focused on a relatively small range of tasks. To overcome this, chemistry apps must be able to seamlessly transfer data to other apps, and through the network to other devices, as well as to other platforms, such as desktops and servers, using documented file formats and protocols whenever possible. This article describes the development and state of the art with regard to chemistry-aware apps that make use of facile data interchange, and some of the scenarios in which these apps can be inserted into a chemical information workflow to increase productivity. A selection of contemporary apps is used to demonstrate their relevance to pharmaceutical research. Mobile apps represent a novel approach for delivery of cheminformatics tools to chemists and other scientists, and indications suggest that mobile devices represent a disruptive technology for drug discovery, as they have been to many other industries.