RESUMEN

Henry's law describes the vapor-liquid equilibrium for dilute gases dissolved in a liquid solvent phase. Descriptions of vapor-liquid equilibrium allow the design of improved separations in the food and beverage industry. The consumer experience of taste and odor are greatly affected by the liquid and vapor phase behavior of organic compounds. This study presents a machine learning (ML) based model that allows quick, accurate predictions of Henry's law constants (kH) for many common organic compounds. Users input only a Simplified Molecular-Input Line-Entry System (SMILES) string or a common English name, and the model returns Henry's law estimates for compounds in water and ethanol. Training was performed on 5,690 compounds. Training data were gathered from an existing database and were supplemented with quantum mechanical (QM) calculations. An extra trees regression model was generated that predicts kH with a mean absolute error of 1.3 in log space and an R2 of 0.98. The model is applied to common flavor and odor compounds in bourbon whiskey as a test case for food and beverage applications.

RESUMEN

In materials science, accurately computing properties like viscosity, melting point, and glass transition temperatures solely through physics-based models is challenging. Data-driven machine learning (ML) also poses challenges in constructing ML models, especially in the material science domain where data is limited. To address this, we integrate physics-informed descriptors from molecular dynamics (MD) simulations to enhance the accuracy and interpretability of ML models. Our current study focuses on accurately predicting viscosity in liquid systems using MD descriptors. In this work, we curated a comprehensive dataset of over 4000 small organic molecules' viscosities from scientific literature, publications, and online databases. This dataset enabled us to develop quantitative structure-property relationships (QSPR) consisting of descriptor-based and graph neural network models to predict temperature-dependent viscosities for a wide range of viscosities. The QSPR models reveal that including MD descriptors improves the prediction of experimental viscosities, particularly at the small data set scale of fewer than a thousand data points. Furthermore, feature importance tools reveal that intermolecular interactions captured by MD descriptors are most important for viscosity predictions. Finally, the QSPR models can accurately capture the inverse relationship between viscosity and temperature for six battery-relevant solvents, some of which were not included in the original data set. Our research highlights the effectiveness of incorporating MD descriptors into QSPR models, which leads to improved accuracy for properties that are difficult to predict when using physics-based models alone or when limited data is available.

RESUMEN

In this work, the vapor pressure of pesticides is employed as an indicator of their volatility potential. Quantitative Structure-Property Relationship models are established to predict the classification of compounds according to their volatility, into the high and low binary classes separated by the 1-mPa limit. A large dataset of 1005 structurally diverse pesticides with known experimental vapor pressure data at 20 °C is compiled from the publicly available Pesticide Properties DataBase (PPDB) and used for model development. The freely available PaDEL-Descriptor and ISIDA/Fragmentor molecular descriptor programs provide a large number of 19,947 non-conformational molecular descriptors that are analyzed through multivariable linear regressions and the Replacement Method technique. Through the selection of appropriate molecular descriptors of the substructure fragment type and the use of different standard classification metrics of model's quality, the classification of the structure-property relationship achieves acceptable results for discerning between the high and low volatility classes. Finally, an application of the obtained QSPR model is performed to predict the classes for 504 pesticides not having experimentally measured vapor pressures.

Asunto(s)

RESUMEN

Due to the complex composition and similar structure, the extraction denitrification of aromatic rich oil is faced with the contradiction problem of denitrification efficiency and aromatic loss which cannot be efficiently solved by experiments. However, the complex interactions involved can be analyzed from the perspective of calculation, and the prediction criteria and methods are proposed. Based on rigorous density functional theory calculation data, Simple models based on electrostatic potential (ESP) and Van der Waals potential (VdWP)-based calculations were established and validated. The twofold model provided the best prediction for interactions between extractants and nitrogen compounds and between extractants and aromatics, which determines denitrification efficiency and aromatic loss, respectively, due to the most complete description of both electrostatic and VdW force. This provides a powerful tool for evaluating the non-covalent interactions and thence tuning the efficiency of the separation process. Thus, high denitrification efficiency (43.2~66.3 %) and moderate aromatic loss (1.7~4.4 %) were obtained using screened deep eutectic solvents (DESs). This ideal observation provided the potential for mild hydrodesulfurization and manufacture of high-grade carbon materials.

RESUMEN

Physical hazards of chemical mixtures, associated for example with their fire or explosion risks, are generally characterized using experimental tools. These tests can be expensive, complex, long to perform and even dangerous for operators. Therefore, for several years and especially with the implementation of the REACH regulation, predictive methods like quantitative structure-property relationships have been encouraged as alternatives tests to determine (eco)toxicological but also physical hazards of chemical substances. Initially, these approaches were intended for pure products, by considering a molecular similarity principle. However, additional to those for pure products, QSPR models for mixtures recently appeared and represent an increasing field of research. This study proposes a state of the art of existing QSPR models specifically dedicated to the prediction of the physical hazards of mixtures. Identified models have been analysed on the key elements of model development (experimental data and fields of application, descriptors used, development and validation methods). It draws up an overview of the potential and limitations of current models as well as areas of progress towards enlarged deployment as a complement to experimental characterizations, for example in the search for safer substances (according to safety-by-design concepts).

Asunto(s)

RESUMEN

The development of biopharmaceutical downstream processes relies on exhaustive experimental studies. The root cause is the poorly understood relationship between the protein structure of monoclonal antibodies (mAbs) and their macroscopic process behavior. Especially the development of preparative chromatography processes is challenged by the increasing structural complexity of novel antibody formats and accelerated development timelines. This study introduces a multiscale in silico model consisting of homology modeling, quantitative structure-property relationships (QSPR), and mechanistic chromatography modeling leading from the amino acid sequence of a mAb to the digital representation of its cation exchange chromatography (CEX) process. The model leverages the mAbs' structural characteristics and experimental data of a diverse set of 21 therapeutic antibodies to predict elution profiles of two mAbs that were removed from the training data set. QSPR modeling identified mAb-specific protein descriptors relevant for the prediction of the thermodynamic equilibrium and the stoichiometric coefficient of the adsorption reaction. The consideration of two discrete conformational states of IgG4 mAbs enabled prediction of split-peak elution profiles. Starting from the sequence, the presented multiscale model allows in silico development of chromatography processes before protein material is available for experimental studies.

Asunto(s)

RESUMEN

Enrichment of ionic poly/perfluoroalkyl substances (PFASs) in aqueous aerosol (AA) is an important pathway for them to enter atmosphere. In this study, the enrichment behaviors of 12 legacy and emerging PFASs in AA in both single solute and mixed solutions were investigated. The enrichment factors (EF) displayed a general increasing trend with the fluorinated carbon chain length. For the first time, a robust Quantitative Structure-Property Relationship (QSPR) model coupled with partial least-square method was established with fifteen quantum chemical descriptors. Four molecular descriptors, including dipole moment (µ), molecular weight (MW), the maximal value of the molecular surface potential (Vs, max) and molecular volume (V) were identified as the key structural variables affecting the PFASs enrichment. Inorganic salts and humic acid (HA) which are common in seawater, facilitated the PFASs enrichment as a result of enhanced hydrophobicity and the bridging effect caused by divalent cations. The typical cationic and anionic surfactants, cetyltrimethylammonium bromide and sodium dodecyl sulfate, both inhibited the enrichment due to the competition between PFASs and surfactants. It is interesting that 6:2 chlorinated polyfluorinated ether sulfonate (F53B) had the highest EF among the 12 PFASs, implying its strong potential of atmosphere transport.

Asunto(s)

RESUMEN

The characterization of physical hazards of substances is a key information to manage the risks associated to their use, storage and transport. With decades of work in this area, Ineris develops and implements cutting-edge experimental facilities allowing such characterizations at different scales and under various conditions to study all of the dreaded accident scenarios. This review presents the efforts engaged by Ineris more recently in the field of chemoinformatics to develop and use new predictive methods for the anticipation and management of industrials risks associated to energetic and reactive materials as a complement to experiments. An overview of the methods used for the development of Quantitative Structure-Property Relationships for physical hazards are presented and discussed regarding the specificities associated to this class of properties. A review of models developed at Ineris is also provided from the first tentative models on the explosivity of nitro compounds to the successful application to the flammability of organic mixtures. Then, a discussion is proposed on the use of QSPR models. Good practices for robust use for QSPR models are recalled with specific comments related to physical hazards, notably for regulatory purpose. Dissemination and training efforts engaged by Ineris are also presented. The potential offered by these predictive methods in terms of in silico design and for the development of new intrinsically safer technologies in safety-by-design strategies is finally discussed. At last, challenges and perspectives to extend the application of chemoinformatics in the field of safety and in particular for the physical hazards of energetic and reactive substances are proposed.

Asunto(s)

RESUMEN

Developing models for predicting molecular properties of organic compounds is imperative for drug development and environmental safety; however, development of such models that have high predictive power and are independent of the compounds used is challenging. To overcome the challenges, we used a thermodynamics-based theoretical derivation to construct models for accurately predicting molecular properties. The free energy change that determines a property equals the sum of the free energy changes (ΔGFs) caused by the factors affecting the property. By developing or selecting molecular descriptors that are directly proportional to ΔGFs, we built a general linear free energy relationship (LFER) for predicting the property with the molecular descriptors as predictive variables. The LFER can be used to construct models for predicting various specific properties from partition coefficients. Validations show that the models constructed according to the LFER have high predictive power and their performance is independent of the compounds used, including the models for the properties having little correlation with partition coefficients. The findings in this study are highly useful for applications in drug development and environmental safety.

RESUMEN

One of the main obstacles in neurological disease treatment is the presence of the blood-brain barrier. New predictive high-throughput screening tools are essential to avoid costly failures in the advanced phases of development and to contribute to the 3 Rs policy. The objective of this work was to jointly develop a new in vitro system coupled with a physiological-based pharmacokinetic (PBPK) model able to predict brain concentration levels of different drugs in rats. Data from in vitro tests with three different cells lines (MDCK, MDCK-MDR1 and hCMEC/D3) were used together with PK parameters and three scaling factors for adjusting the model predictions to the brain and plasma profiles of six model drugs. Later, preliminary quantitative structure-property relationships (QSPRs) were constructed between the scaling factors and the lipophilicity of drugs. The predictability of the model was evaluated by internal validation. It was concluded that the PBPK model, incorporating the barrier resistance to transport, the disposition within the brain and the drug-brain binding combined with MDCK data, provided the best predictions for passive diffusion and carrier-mediated transported drugs, while in the other cell lines, active transport influence can bias predictions.

RESUMEN

The capability to predict corneal permeability based on physicochemical parameters has always been a desirable objective of ophthalmic drug development. However, previous work has been limited to cases where either the diversity of compounds used was lacking or the performance of the models was poor. Our study provides extensive quantitative structure-property relationship (QSPR) models for corneal permeability predictions. The models involved in vitro corneal permeability measurements of 189 diverse compounds. Preliminary analysis of data showed that there is no significant correlation between corneal-PAMPA (Parallel Artificial Membrane Permeability Assay) permeability values and other pharmacokinetically relevant in silico drug transport parameters like Caco-2, jejunal permeability and blood-brain partition coefficient (logBB). Two different QSPR models were developed: one for corneal permeability and one for corneal membrane retention, based on experimental corneal-PAMPA permeability data. Partial least squares regression was applied for producing the models, which contained classical molecular descriptors and ECFP fingerprints in combination. A complex validation protocol (including internal and external validation) was carried out to provide robust and appropriate predictions for the permeability and membrane retention values. Both models had an overall fit of R2 > 0.90, including R2-values not lower than 0.85 for validation runs, and provide quick and accurate predictions of corneal permeability values for a diverse set of compounds.

Asunto(s)

RESUMEN

In recent years there has been a growing interest in studying the differences between the chemical and biological space represented by natural products (NPs) of terrestrial and marine origin. In order to learn more about these two chemical spaces, marine natural products (MNPs) and terrestrial natural products (TNPs), a machine learning (ML) approach was developed in the current work to predict three classes, MNPs, TNPs and a third class of NPs that appear in both the terrestrial and marine environments. In total 22,398 NPs were retrieved from the Reaxys® database, from those 10,790 molecules are recorded as MNPs, 10,857 as TNPs, and 761 NPs appear registered as both MNPs and TNPs. Several ML algorithms such as Random Forest, Support Vector Machines, and deep learning Multilayer Perceptron networks have been benchmarked. The best performance was achieved with a consensus classification model, which predicted the external test set with an overall predictive accuracy up to 81 %. As far as we know this approach has never been intended and therefore allow to be used to better understand the chemical space defined by MNPs, TNPs or both, but also in virtual screening to define the applicability domain of QSAR models of MNPs and TNPs.

Asunto(s)

RESUMEN

A priori knowledge of physicochemical properties such as melting and boiling could expedite materials discovery. However, theoretical modeling from first principles poses a challenge for efficient virtual screening of potential candidates. As an alternative, the tools of data science are becoming increasingly important for exploring chemical datasets and predicting material properties. Herein, we extend a molecular representation, or set of descriptors, first developed for quantitative structure-property relationship modeling by Yalkowsky and coworkers known as the Unified Physicochemical Property Estimation Relationships (UPPER). This molecular representation has group-constitutive and geometrical descriptors that map to enthalpy and entropy; two thermodynamic quantities that drive thermal phase transitions. We extend the UPPER representation to include additional information about sp2-bonded fragments. Additionally, instead of using the UPPER descriptors in a series of thermodynamically-inspired calculations, as per Yalkowsky, we use the descriptors to construct a vector representation for use with machine learning techniques. The concise and easy-to-compute representation, combined with a gradient-boosting decision tree model, provides an appealing framework for predicting experimental transition temperatures in a diverse chemical space. An application to energetic materials shows that the method is predictive, despite a relatively modest energetics reference dataset. We also report competitive results on diverse public datasets of melting points (i.e., OCHEM, Enamine, Bradley, and Bergström) comprised of over 47k structures. Open source software is available at https://github.com/USArmyResearchLab/ARL-UPPER.

Asunto(s)

RESUMEN

Macrocycles constitute superior ligands for targets that have flat binding sites but often require long synthetic routes, emphasizing the need for property prediction prior to synthesis. We have investigated the scope and limitations of machine learning classification models and of regression models for predicting the cell permeability of a set of denovo-designed, drug-like macrocycles. 2D-Based classification models, which are fast to calculate, discriminated between macrocycles that had low-medium and high permeability and may be used as virtual filters in early drug discovery projects. Importantly, stereo- and regioisomer were correctly classified. QSPR studies of two small sets of comparator drugs suggested that use of 3D descriptors, calculated from biologically relevant conformations, would allow development of more precise regression models for late phase drug projects. However, a 3D permeability model could only be developed for a rigid series of macrocycles. Comparison of NMR based conformational analysis with in silico conformational sampling indicated that this shortcoming originates from the inability of the molecular mechanics force field to identify the relevant conformations for flexible macrocycles. We speculate that a Kier flexibility index of ≤10 constitutes a current upper limit for reasonably accurate 3D prediction of macrocycle cell permeability.

Asunto(s)

RESUMEN

A series of nineteen novel ring-substituted N-arylcinnamanilides was synthesized and characterized. All investigated compounds were tested against Staphylococcus aureus as the reference strain, two clinical isolates of methicillin-resistant S. aureus (MRSA), and Mycobacterium tuberculosis. (2E)-N-[3-Fluoro-4-(trifluoromethyl)phenyl]-3-phenylprop-2-enamide showed even better activity (minimum inhibitory concentration (MIC) 25.9 and 12.9 µM) against MRSA isolates than the commonly used ampicillin (MIC 45.8 µM). The screening of the cell viability was performed using THP1-Blue™ NF-κB cells and, except for (2E)-N-(4-bromo-3-chlorophenyl)-3-phenylprop-2-enamide (IC50 6.5 µM), none of the discussed compounds showed any significant cytotoxic effect up to 20 µM. Moreover, all compounds were tested for their anti-inflammatory potential; several compounds attenuated the lipopolysaccharide-induced NF-κB activation and were more potent than the parental cinnamic acid. The lipophilicity values were specified experimentally as well. In addition, in silico approximation of the lipophilicity values was performed employing a set of free/commercial clogP estimators, corrected afterwards by the corresponding pKa calculated at physiological pH and subsequently cross-compared with the experimental parameters. The similarity-driven property space evaluation of structural analogs was carried out using the principal component analysis, Tanimoto metrics, and Kohonen mapping.

Asunto(s)

RESUMEN

Sorption is one of the key process that affects the fate and mobility of pharmaceuticals in the soil environment. Several models have been developed for estimating the sorption of organic chemicals, including ionisable compounds, in soil. However, the applicability of these models to pharmaceuticals has not been extensively tested. In this study, we generated a high-quality dataset on the sorption of twenty-one pharmaceuticals in different soil types and used these data to evaluate existing models and to develop new improved models. Sorption coefficients (Kd) of the pharmaceuticals ranged from 0.2 to 1249.2 L/kg. Existing models were unable to adequately estimate the measured sorption data. Using the data, new models were developed, incorporating molecular and soil descriptors, that outperformed the published models when evaluated against external data sets. While there is a need for further evaluation of these new models against broader sorption datasets obtained at environmentally relevant concentrations, in the future they could be highly useful in supporting environmental risk assessment and prioritization efforts for pharmaceutical ingredients.

Asunto(s)

RESUMEN

The accurate prediction of solubility of drugs is still problematic. It was thought for a long time that shortfalls had been due the lack of high-quality solubility data from the chemical space of drugs. This study considers the quality of solubility data, particularly of ionizable drugs. A database is described, comprising 6355 entries of intrinsic solubility for 3014 different molecules, drawing on 1325 citations. In an earlier publication, many factors affecting the quality of the measurement had been discussed, and suggestions were offered to improve ways of extracting more reliable information from legacy data. Many of the suggestions have been implemented in this study. By correcting solubility for ionization (i.e., deriving intrinsic solubility, S0) and by normalizing temperature (by transforming measurements performed in the range 10-50 °C to 25 °C), it can now be estimated that the average interlaboratory reproducibility is 0.17 log unit. Empirical methods to predict solubility at best have hovered around the root mean square error (RMSE) of 0.6 log unit. Three prediction methods are compared here: (a) Yalkowsky's general solubility equation (GSE), (b) Abraham solvation equation (ABSOLV), and (c) Random Forest regression (RFR) statistical machine learning. The latter two methods were trained using the new database. The RFR method outperforms the other two models, as anticipated. However, the ability to predict the solubility of drugs to the level of the quality of data is still out of reach. The data quality is not the limiting factor in prediction. The statistical machine learning methodologies are probably up to the task. Possibly what's missing are solubility data from a few sparsely-covered chemical space of drugs (particularly of research compounds). Also, new descriptors which can better differentiate the factors affecting solubility between molecules could be critical for narrowing the gap between the accuracy of the prediction models and that of the experimental data.

RESUMEN

Contamination of drinking water with pharmaceuticals and personal care products (PPCPs) is an issue of health concerns. To effectively control the level of PPCPs in drinking water, a pilot study employing two parallel trains of two-stage biofiltration, i.e., a sand/anthracite (SA) biofilter coupled with a biologically-active granular activated carbon (GAC) post-filter contactor, was conducted as a post-treatment after coagulation in a drinking water treatment plant. Results showed the biofiltration process could effectively remove PPCPs with an average removal of 53.4%, where the GAC contactor played the dominant role to remove 48.1% of the total PPCPs. The molecular properties determined the removability of individual PPCPs, i.e., smaller molecules with simpler structure connectivity were more likely to be removed. Based on the quantitative structure-property relationships (QSPRs) analysis, a simple regression model was proposed to predict the removability of each PPCP across the biofiltration process. The drinking water equivalent level (DWEL) quotient method was developed to assess the health risks of detected PPCPs in water samples. The biofiltration process showed efficient capacity to reduce the health risks of PPCPs with an average removal of 79%, and the PPCPs in the effluents generally would not pose adverse health effects. Pearson correlation analysis explored the possible role of nitrogenous PPCPs (N-PPCPs) as the precursors of nitrogenous disinfection byproducts (N-DBPs) in drinking waters. Aromatic nitrogen in PPCPs was found to be a significant descriptor for the formation potential of trichloroacetonitrile (TCAN). In addition, it was found that pre-filter chlorination could slightly improve the biofiltration of PPCPs.

Asunto(s)

RESUMEN

Emerging applications of a new class of materials, sequence-defined macromolecules, are explored. Such molecularly highly defined macromolecules require stringent synthesis and purification procedures, yet offer unprecedented application possibilities. The first examples of molecular data storage and related technologies are already starting to emerge today. From a more fundamental point of view, such macromolecules offer a unique opportunity to determine quantitative structure-property relationships (QSPR), which critically aids in designing materials with applications ranging from catalysis to artificial enzymes.

Asunto(s)

RESUMEN

New Quantitative Structure-Property Relationships (QSPR) are presented to predict the flash point of binary liquid mixtures, based on more than 600 experimental flash points for 60 binary mixtures. Two models are proposed based on a GA-MLR approach that uses a genetic algorithm (GA) variable selection in multilinear regressions (MLR). In these models, mixtures were characterized by a series of mixture descriptors calculated from various mixture formula combining the molecular descriptors of the single compounds constituting the mixtures and their respective molar fractions in the mixture. The best model demonstrated good predictive capabilities with a mean absolute error of only 7.3 °C estimated for an external validation set. Moreover, this model is focused on mixture descriptors applicable to more complex mixtures, i. e. constituted of more than 2 components, and already demonstrated interesting predictions for a series of ternary mixtures.

Asunto(s)