Search

Abstract Natural products are a diverse class of compounds with promising biological properties, such as high potency and excellent selectivity. However, they have different structural motifs than typical drug-like compounds, e.g., a wider range of molecular weight, multiple stereocenters and higher fraction of sp 3 -hybridized carbons. This makes the encoding of natural products via molecular fingerprints difficult, thus restricting their use in cheminformatics studies. To tackle this issue, we explored over 30 years of research to systematically evaluate which molecular fingerprint provides the best performance on the natural product chemical space. We considered 20 molecular fingerprints from four different sources, which we then benchmarked on over 100,000 unique natural products from the COCONUT (COlleCtion of Open Natural prodUcTs) and CMNPD (Comprehensive Marine Natural Products Database) databases. Our analysis focused on the correlation between different fingerprints and their classification performance on 12 bioactivity prediction datasets. Our results show that different encodings can provide fundamentally different views of the natural product chemical space, leading to substantial differences in pairwise similarity and performance. While Extended Connectivity Fingerprints are the de-facto option to encoding drug-like compounds, other fingerprints resulted to match or outperform them for bioactivity prediction of natural products. These results highlight the need to evaluate multiple fingerprinting algorithms for optimal performance and suggest new areas of research. Finally, we provide an open-source Python package for computing all molecular fingerprints considered in the study, as well as data and scripts necessary to reproduce the results, at https://github.com/dahvida/NP_Fingerprints .

Abstract Photoemission electron microscopy (PEEM) is a powerful technique for surface characterization that provides detailed information on the chemical and structural properties of materials at the nanoscale. In this study, the potential is explored using a machine learning algorithm called self‐organizing map with a relational perspective map (SOM‐RPM) for visualizing and analyzing complex PEEM‐generated datasets. The application of SOM‐RPM is demonstrated using synchrotron‐based X‐ray magnetic circular dichroism (XMCD)‐PEEM data acquired from a pyrrhotite sample. Traditional visualization approaches for XMCD‐PEEM data may not fully capture the complexity of the sample, especially in the case of heterogeneous materials. By applying SOM‐RPM to the XMCD‐PEEM data, a colored topographic map is created that represents the spectral similarities and dissimilarities among the pixels. This approach allows for a more intuitive and easily interpretable representation of the data without the need of data binning or spectral smoothing. The results of the SOM‐RPM analysis are compared to the conventional visualization approach, highlighting the advantages of SOM‐RPM in revealing features that are not readily observable in the conventional method. This study suggests that the SOM‐RPM approach can be used complimentarily for other PEEM‐based measurements, such as core level and valence band X‐ray photoelectron spectroscopy.

Cardiovascular disease is a leading global cause of mortality. The potential cardiotoxic effects of chemicals from different classes, such as environmental contaminants, pesticides, and drugs can significantly contribute to effects on health. The same chemical can induce cardiotoxicity in different ways, following various Adverse Outcome Pathways (AOPs). In addition, the potential synergistic effects between chemicals further complicate the issue. In silico methods have become essential for tackling the problem from different perspectives, reducing the need for traditional in vivo testing, and saving valuable resources in terms of time and money. Artificial intelligence (AI) and machine learning (ML) are among today’s advanced approaches for evaluating chemical hazards. They can serve, for instance, as a first-tier component of Integrated Approaches to Testing and Assessment (IATA). This study employed ML and AI to assess interactions between chemicals and specific biological targets within the AOP networks for cardiotoxicity, starting with molecular initiating events (MIEs) and progressing through key events (KEs). We explored methods to encode chemical information in a suitable way for ML and AI. We started with commonly used approaches in Quantitative Structure–Activity Relationship (QSAR) methods, such as molecular descriptors and different types of fingerprint. We then increased the complexity of encoders, incorporating graph-based methods, auto-encoders, and character embeddings employed in neural language processing. We also developed a multimodal neural network architecture, capable of considering the complementary nature of different chemical representations simultaneously. The potential of this approach, compared to more conventional architectures designed to handle a single encoder, becomes apparent when the amount of data increases.

Abstract Biofilm formation is a major cause of hospital‐acquired infections. Research into biofilm‐resistant materials is therefore critical to reduce the frequency of these events. Polymer microarrays offer a high‐throughput approach to enable the efficient discovery of novel biofilm‐resistant polymers. Herein, bacterial attachment and surface chemistry are studied for a polymer microarray to improve the understanding of Pseudomonas aeruginosa biofilm formation on a diverse set of polymeric surfaces. The relationships between time‐of‐flight secondary ion mass spectrometry (ToF‐SIMS) data and biofilm formation are analyzed using linear multivariate analysis (partial least squares [PLS] regression) and a nonlinear self‐organizing map (SOM). The SOM models revealed several combinations of fragment ions that are positively or negatively associated with bacterial biofilm formation, which are not identified by PLS. With these insights, a second PLS model is calculated, in which interactions between key fragments (identified by the SOM) are explicitly considered. Inclusion of these terms improved the PLS model performance and shows that, without such terms, certain key fragment ions correlated with bacterial attachment may not be identified. The chemical insights provided by the combination of PLS regression and SOM will be useful for the design of materials that support negligible pathogen attachment.

Bitter pit (BP) is one of the most relevant post-harvest disorders for apple industry worldwide, which is often related to calcium (Ca) deficiency at the calyx end of the fruit. Its occurrence takes place along with an imbalance with other minerals, such as potassium (K). Although the K/Ca ratio is considered a valuable indicator of BP, a high variability in the levels of these elements occurs within the fruit, between fruits of the same plant, and between plants and orchards. Prediction systems based on the content of elements in fruit have a high variability because they are determined in samples composed of various fruits. With X-ray fluorescence (XRF) spectrometry, it is possible to characterize non-destructively the signal intensity for several mineral elements at a given position in individual fruit and thus, the complete signal of the mineral composition can be used to perform a predictive model to determine the incidence of bitter pit. Therefore, it was hypothesized that using a multivariate modeling approach, other elements beyond the K and Ca could be found that could improve the current clutter prediction capability. Two studies were carried out: on the first one an experiment was conducted to determine the K/Ca and the whole spectrum using XRF of a balanced sample of affected and non-affected ‘Granny Smith’ apples. On the second study apples of three cultivars (‘Granny Smith’, ‘Brookfield’ and ‘Fuji’), were harvested from two commercial orchards to evaluate the use of XRF to predict BP. With data from the first study a multivariate classification system was trained (balanced database of healthy and BP fruit, consisting in 176 from each group) and then the model was applied on the second study to fruit from two orchards with a history of BP. Results show that when dimensionality reduction was performed on the XRF spectra (1.5 - 8 KeV) of ‘Granny Smith’ apples, comparing fruit with and without BP, along with K and Ca, four other elements (i.e., Cl, Si, P, and S) were found to be deterministic. However, the PCA revealed that the classification between samples (BP vs. non-BP fruit) was not possible by univariate analysis (individual elements or the K/Ca ratio).Therefore, a multivariate classification approach was applied, and the classification measures (sensitivity, specificity, and balanced precision) of the PLS-DA models for all cultivars evaluated (‘Granny Smith’, ‘Fuji’ and ‘Brookfield’) on the full training samples and with both validation procedures (Venetian and Monte Carlo), ranged from 0.76 to 0.92. The results of this work indicate that using this technology at the individual fruit level is essential to understand the factors that determine this disorder and can improve BP prediction of intact fruit.

The purpose of this work is the creation of a chemical database named ChemTastesDB that includes both organic and inorganic tastants. The creation, curation pipeline and the main features of the database are described in detail. The database includes 2944 verified and curated compounds divided into nine classes, which comprise the five basic tastes (sweet, bitter, umami sour and salty) along with four additional categories: tasteless, non-sweet, multitaste and miscellaneous. ChemTastesDB provides the following information for each tastant: name, PubChem CID, CAS registry number, canonical SMILES, class taste and references to the scientific sources from which data were retrieved. The molecular structure in the HyperChem (.hin) format of each chemical is also made available. In addition, molecular fingerprints were used for characterizing and analyzing the chemical space of tastants by means of unsupervised machine learning. ChemTastesDB constitutes a useful tool to the scientific community to expand the information of taste molecules and to assist in silico studies for the taste prediction of unevaluated and as yet unsynthetized compounds, as well as the analysis of the relationships between molecular structure and taste. The database is freely accessible at https://doi.org/10.5281/zenodo.5747393.

Membrane introduction mass spectrometry (MIMS) is a direct mass spectrometry technique used to monitor online chemical systems or quickly quantify trace levels of different groups of compounds in complex matrices without extensive sample preparation steps and chromatographic separation. MIMS utilizes a thin, semi-permeable, and selective membrane that directly connects the sample and the mass spectrometer. The analytes in the sample are pre-concentrated by the membrane depending on their physicochemical properties and directly transferred, using different acceptor phases (gas, liquid or vacuum) to the mass spectrometer. Condensed phase (CP) MIMS use a liquid as a medium, extending the range to new applications to less-volatile compounds that are challenging or unsuitable to gas-phase MIMS. It directly allows the rapid quantification of selected compounds in complex matrices, the online monitoring of chemical reactions (in real-time), as well as in situ measurements. CP-MIMS has expanded beyond the measurement of several organic compounds because of the use of different types of liquid acceptor phases, geometries, dimensions, and mass spectrometers. This review surveys advancements of CP-MIMS and its applications to several molecules and matrices over the past 15 years.

According to the 2021 World Drug Report, around 275 million people use drugs of abuse, and 36 million people suffer from addiction, fostering a thriving market for illicit substances. In Italy, 30,083 people were reported to the Judicial Authority for offenses in violation of the Italian Law D.P.R. 309/1990. These offences are sentenced after a qualitative–quantitative analysis of seized materials. Given the large quantity of seized drugs and the need to perform accurate analytical determinations, Italian forensic laboratories struggle to complete analyses in a short time, delaying the entire reporting process needed to achieve sentencing. For this purpose, an UHPLC-MS/MS-based platform was developed at the University of Milano-Bicocca to support law-enforcement authorities. Software was designed to easily manage street seizure acquisition, documentation registration, and sampling. A sensitive UHPLC-MS/MS method was fully validated for the quantification of the traditional illicit substances (cocaine, heroin, 6-MAM, morphine, amphetamine, methamphetamine, MDMA, ketamine, GHB, GBL, LSD, trans-∆9-THC, and THCA) at the ppb level. The final report is relayed to the Prefecture in 3–4 days, even within 24 h for urgent requests. The platform allows for semi-automatic data handling to minimize erroneous results for an accurate report generation by standardized procedures.

The study concerns the photodegradation of the antidepressant escitalopram (ESC), the S-enantiomer of the citalopram raceme, both in ultrapure and surface water, considering the contribution of indirect photolysis through the presence of nitrate and bicarbonate. The effect of nitrate and bicarbonate concentrations was investigated by full factorial design, and only the nitrate concentration resulted in having a significant effect on the degradation. The kinetics of ESC photodegradation is the pseudo-first-order (half-life = 62.4 h in ultrapure water and 48.4 h in lake water). The generation of transformation products (TPs) was monitored through a developed and validated HPLC-MS/MS method. Fourteen TPs were identified in ultrapure water (one of them, at m/z 261, for the first time) and other two TPs at m/z 327 (found for the first time in this study) were identified only in presence of a nitrate. Several TPs were the same as those formed during the photodegradation of citalopram. The photodegradation pathway of ESC and its mechanism of degradation in water is proposed. The method was applied successfully to the analyses of surface water samples, in which a few dozen of ng L−1 of ESC was determined together with the presence of TP2, TP5 and TP12. Finally, a preliminary in silico evaluation of the toxicological profile and environmental behavior of TPs by computational models was carried out; two TPs (TP4 and TP10) were identified as of potential concern, as they were predicted mutagenic by Ames test model.

Mass spectrometry (MS) is widely used for the identification of chemical compounds by matching the experimentally acquired mass spectrum against a database of reference spectra. However, this approach suffers from a limited coverage of the existing databases causing a failure in the identification of a compound not present in the database. Among the computational approaches for mining metabolite structures based on MS data, one option is to predict molecular fingerprints from the mass spectra by means of chemometric strategies and then use them to screen compound libraries. This can be carried out by calibrating multi-task artificial neural networks from large datasets of mass spectra, used as inputs, and molecular fingerprints as outputs. In this study, we prepared a large LC-MS/MS dataset from an on-line open repository. These data were used to train and evaluate deep-learning-based approaches to predict molecular fingerprints and retrieve the structure of unknown compounds from their LC-MS/MS spectra. Effects of data sparseness and the impact of different strategies of data curing and dimensionality reduction on the output accuracy have been evaluated. Moreover, extensive diagnostics have been carried out to evaluate modelling advantages and drawbacks as a function of the explored chemical space.

Antineoplastic drugs are a wide and heterogeneous group of substances that, as universally known, can cause highly severe toxic effects to whoever is exposed. From an occupational safety point of view, surface contaminations inside preparation and administration units are a growing issue and therefore require the development and implementation of sensible and fast monitoring methods. The unlikelihood of a unique all-embracing chromatography, able to correctly retain and separate each analyte led to the need to create an orthogonal normal phase analysis, which might be able to fill the gaps in the more common reversed-phase ones. An existing hydrophilic interaction method has thus been expanded to 6 other drugs and applied to real samples after an evaluation of its performances. The experimental data were then used to evaluate the possibility of estimating reliable relationships between the chromatographic retention and the chemical-structural features of the drugs under analysis.

Neural networks are rapidly gaining popularity in chemical modeling and Quantitative Structure–Activity Relationship (QSAR) thanks to their ability to handle multitask problems. However, outcomes of neural networks depend on the tuning of several hyperparameters, whose small variations can often strongly affect their performance. Hence, optimization is a fundamental step in training neural networks although, in many cases, it can be very expensive from a computational point of view. In this study, we compared four of the most widely used approaches for tuning hyperparameters, namely, grid search, random search, tree-structured Parzen estimator, and genetic algorithms on three multitask QSAR datasets. We mainly focused on parsimonious optimization and thus not only on the performance of neural networks, but also the computational time that was taken into account. Furthermore, since the optimization approaches do not directly provide information about the influence of hyperparameters, we applied experimental design strategies to determine their effects on the neural network performance. We found that genetic algorithms, tree-structured Parzen estimator, and random search require on average 0.08% of the hours required by grid search; in addition, tree-structured Parzen estimator and genetic algorithms provide better results than random search.

The lignocellulosic materials produced after each step of a biorefinery plant using simultaneous saccharification and fermentation (SSF) technology on wheat straw (Triticum spp.) for bioethanol production were characterized by spectroscopic and chromatographic techniques in order to investigate the macromolecular interactions between the lignin and polysaccharides. In order to valorize the lignin cakes, a purification step was set up and the extraction conditions (acid pretreatment, temperature, time, and NaOH concentration) were optimized by a chemiometric approach in terms of yield and purity. Residual carbohydrate impurities, free and/or chemically bonded to lignin (lignin carbohydrates complexes), were individuate as the most critical factor for a satisfactory lignin extraction. Finally, the lignin samples collected according to the optimized extraction conditions were chemically characterized and low molecular weight, high phenols concentration, and low carboxylic acids content were recognized as interesting features for industrial applications.

This work describes a novel approach based on advanced molecular similarity to predict the sweetness of chemicals. The proposed Quantitative Structure-Taste Relationship (QSTR) model is an expert system developed keeping in mind the five principles defined by the Organization for Economic Co-operation and Development (OECD) for the validation of (Q)SARs. The 649 sweet and non-sweet molecules were described by both conformation-independent extended-connectivity fingerprints (ECFPs) and molecular descriptors. In particular, the molecular similarity in the ECFPs space showed a clear association with molecular taste and it was exploited for model development. Molecules laying in the subspaces where the taste assignation was more difficult were modeled trough a consensus between linear and local approaches (Partial Least Squares-Discriminant Analysis and N-nearest-neighbor classifier). The expert system, which was thoroughly validated through a Monte Carlo procedure and an external set, gave satisfactory results in comparison with the state-of-the-art models. Moreover, the QSTR model can be leveraged into a greater understanding of the relationship between molecular structure and sweetness, and into the design of novel sweeteners.

A series of 436 Munro database chemicals were studied with respect to their corresponding experimental LD50 values to investigate the possibility of establishing a global QSAR model for acute toxicity. Dragon molecular descriptors were used for the QSAR model development and genetic algorithms were used to select descriptors better correlated with toxicity data. Toxic values were discretized in a qualitative class on the basis of the Globally Harmonized Scheme: the 436 chemicals were divided into 3 classes based on their experimental LD50 values: highly toxic, intermediate toxic and low to non-toxic. The k-nearest neighbor (k-NN) classification method was calibrated on 25 molecular descriptors and gave a non-error rate (NER) equal to 0.66 and 0.57 for internal and external prediction sets, respectively. Even if the classification performances are not optimal, the subsequent analysis of the selected descriptors and their relationship with toxicity levels constitute a step towards the development of a global QSAR model for acute toxicity.

One of the OECD principles for model validation requires defining the Applicability Domain (AD) for the QSAR models. This is important since the reliable predictions are generally limited to query chemicals structurally similar to the training compounds used to build the model. Therefore, characterization of interpolation space is significant in defining the AD and in this study some existing descriptor-based approaches performing this task are discussed and compared by implementing them on existing validated datasets from the literature. Algorithms adopted by different approaches allow defining the interpolation space in several ways, while defined thresholds contribute significantly to the extrapolations. For each dataset and approach implemented for this study, the comparison analysis was carried out by considering the model statistics and relative position of test set with respect to the training space.