The utilisation of structural descriptors to predict metabolic constants of xenobiotics in mammals

Read-across-based intelligent learning: development of a global q-RASAR model for the efficient quantitative predictions of skin sensitization potential of diverse organic chemicals

Environmental chemicals and contaminants cause a wide array of harmful implications to terrestrial and aquatic life which ranges from skin sensitization to acute oral toxicity. The current study aims to assess the quantitative skin sensitization potential of a large set of industrial and environmental chemicals acting through different mechanisms using the novel quantitative Read-Across Structure-Activity Relationship (q-RASAR) approach. Based on the identified important set of structural and physicochemical features, Read-Across-based hyperparameters were optimized using the training set compounds followed by the calculation of similarity and error-based RASAR descriptors. Data fusion, further feature selection, and removal of prediction confidence outliers were performed to generate a partial least squares (PLS) q-RASAR model, followed by the application of various Machine Learning (ML) tools to check the quality of predictions. The PLS model was found to be the best among different models. A simple user-friendly Java-based software tool was developed based on the PLS model, which efficiently predicts the toxicity value(s) of query compound(s) along with their status of Applicability Domain (AD) in terms of leverage values. This model has been developed using structurally diverse compounds and is expected to predict efficiently and quantitatively the skin sensitization potential of environmental chemicals to estimate their occupational and health hazards.

A comprehensive evaluation of liposome/water partition coefficient prediction models based on the Technique for Order Preference by Similarity to an Ideal Solution (TOPSIS) method: Challenges from different descriptor dimension reduction methods and machine learning algorithms

The liposome/water partition coefficient (K_lip/w) is a key parameter to evaluate the bioaccumulation potential of pollutants. Considering that it is difficult to determine the K_lip/w values of all pollutants through experiments, researchers gradually developed models to predict it. However, there is currently no research on how to comprehensively evaluate prediction models and recommend a compelling optimal modeling method. To remedy the defect of single parameters in a traditional model comparison, the TOPSIS evaluation method, based on entropy weight, was first proposed. We use this method to comprehensively evaluate models from multiple angles in this study. Thirty QSPR models, including 3 descriptor dimension reduction methods and 10 algorithms (belonging to 4 tribes), were used to predict K_lip/w and verify the effectiveness of the comprehensive assessment method. The results showed that RF (descriptor dimension reduction method), symbolism (tribes) and RF (algorithm) exhibited significant advantages in establishing the K_lip/w value prediction model. At present, the application of TOPSIS in environmental model evaluations is almost absent. We hope that the proposed TOPSIS evaluation method can be applied to more chemical datasets and provide a more systematic and comprehensive basis for the application of the QSPR model in environmental studies and other fields.

Machine learning and artificial intelligence in physiologically based pharmacokinetic modeling

Physiologically based pharmacokinetic (PBPK) models are useful tools in drug development and risk assessment of environmental chemicals. PBPK model development requires the collection of species-specific physiological, and chemical-specific absorption, distribution, metabolism, and excretion (ADME) parameters, which can be a time-consuming and expensive process. This raises a need to create computational models capable of predicting input parameter values for PBPK models, especially for new compounds. In this review, we summarize an emerging paradigm for integrating PBPK modeling with machine learning (ML) or artificial intelligence (AI)-based computational methods. This paradigm includes 3 steps (1) obtain time-concentration PK data and/or ADME parameters from publicly available databases, (2) develop ML/AI-based approaches to predict ADME parameters, and (3) incorporate the ML/AI models into PBPK models to predict PK summary statistics (eg, area under the curve and maximum plasma concentration). We also discuss a neural network architecture "neural ordinary differential equation (Neural-ODE)" that is capable of providing better predictive capabilities than other ML methods when used to directly predict time-series PK profiles. In order to support applications of ML/AI methods for PBPK model development, several challenges should be addressed (1) as more data become available, it is important to expand the training set by including the structural diversity of compounds to improve the prediction accuracy of ML/AI models; (2) due to the black box nature of many ML models, lack of sufficient interpretability is a limitation; (3) Neural-ODE has great potential to be used to generate time-series PK profiles for new compounds with limited ADME information, but its application remains to be explored. Despite existing challenges, ML/AI approaches will continue to facilitate the efficient development of robust PBPK models for a large number of chemicals.

IVIVE: Facilitating the Use of In Vitro Toxicity Data in Risk Assessment and Decision Making

During the past few decades, the science of toxicology has been undergoing a transformation from observational to predictive science. New approach methodologies (NAMs), including in vitro assays, in silico models, read-across, and in vitro to in vivo extrapolation (IVIVE), are being developed to reduce, refine, or replace whole animal testing, encouraging the judicious use of time and resources. Some of these methods have advanced past the exploratory research stage and are beginning to gain acceptance for the risk assessment of chemicals. A review of the recent literature reveals a burst of IVIVE publications over the past decade. In this review, we propose operational definitions for IVIVE, present literature examples for several common toxicity endpoints, and highlight their implications in decision-making processes across various federal agencies, as well as international organizations, including those in the European Union (EU). The current challenges and future needs are also summarized for IVIVE. In addition to refining and reducing the number of animals in traditional toxicity testing protocols and being used for prioritizing chemical testing, the goal to use IVIVE to facilitate the replacement of animal models can be achieved through their continued evolution and development, including a strategic plan to qualify IVIVE methods for regulatory acceptance.

Versatile modelling of polyoxymethylene-water partition coefficients for hydrophobic organic contaminants using linear and nonlinear approaches

Environmental fate or transport of hydrophobic organic contaminants (HOCs) depends on the partitioning properties of compounds within various environmental phases. Due to the wide application of polyoxymethylene (POM) in the passive sampling technique, several in silico models were developed to predict POM-water partition coefficients (K_POM-w) in accordance with the guidelines of the Organization for Economic Cooperation and Development (OECD). It is an attempt to combine conventional linear method (multiple linear regression, MLR) and popular nonlinear algorithm (artificial neural network, ANN) for estimating partition coefficients of HOCs. All models were performed on a dataset of 210 chemicals from 13 different classes. The polyparameter linear free energy relationship (pp-LFER) model included 5 molecular descriptors, namely, E, S, A, B and V, and predicted log K_POM-w with R²_adj of 0.825. The values of statistical parameters including R²_adj, Q²_ext, RMSE_tra and RMSE_ext for quantitative structure-property relationship (QSPR)-MLR and QSPR-ANN models with four descriptors (ALOGP, MeanDD, E1m and Mor24s) were: (0.928, 0.877, 0.498 and 0.649) and (0.943, 0.905, 0.443 and 0.571), with high similarity for both models, which confirmed the robustness, significance, and remarkable prediction accuracy of the QSPR models. Moreover, the mechanism interpretation revealed that the molecular volume and hydrophobicity had a major impact on distribution procedure of HOCs. The models developed herein, with the broad applicability domain (AD), provide suitable tools to fill the experimental data gap for untested chemicals and help researchers better understand the mechanistic basis of adsorption behavior of POM.

A regression-based model to predict chemical migration from packaging to food

Packaging materials can be a source of chemical contaminants in food. Process-based migration models (PMM) predict the chemical fraction transferred from packaging materials to food (F_C) for application in prioritisation tools for human exposure. These models, however, have a relatively limited applicability domain and their predictive performance is typically low. To overcome these limitations, we developed a linear mixed-effects model (LMM) to statistically relate measured F_C to properties of chemicals, food, packaging, and experimental conditions. We found a negative relationship between the molecular weight (MW) and F_C, and a positive relationship with the fat content of the food depending on the octanol-water partitioning coefficient of the migrant. We also showed that large chemicals (MW > 400 g/mol) have a higher migration potential in packaging with low crystallinity compared with high crystallinity. The predictive performance of the LMM for chemicals not included in the database in contact with untested food items but known packaging material was higher (Coefficient of Efficiency (CoE) = 0.21) compared with a recently developed PMM (CoE = -5.24). We conclude that our empirical model is useful to predict chemical migration from packaging to food and prioritise chemicals in the absence of measurements.

Synthesis, crystal structure, insecticidal activities, molecular docking and QSAR studies of some new phospho guanidines and phospho pyrazines as cholinesterase inhibitors

Novel phospho guanidine and phospho pyrazine derivatives were synthesized and characterized by ³¹P, ¹³C, ¹HNMR and IR spectroscopy to obtain novel and human-safe insecticides. Compound 35 [(C₄H₄N₂NH)₂P(O)(C₆H₆)] was investigated by X-ray crystallography. The inhibitory effects of synthesized compounds were evaluated on human and insect acetylcholinesterase (AChE) using in vitro Ellman method. A few of these compounds, which had low human toxicity, were selected for assessing the killing effects (in vivo) on the elm leaf beetle (X.luteola). The in vitro and in vivo results indicated that compounds bearing both phosphoryl groups and aromatic systems were found to possess a good selectivity for the inhibition of insect AChE over human AChE; up to 550-fold selectivity was achieved for compound 19. Docking studies were performed to explain reasons for the selective behavior of AChE inhibitors. Additionally, the quantitative structure-activity relationship (QSAR) and density functional theory (DFT) results of AChEs demonstrated that the size, shape, dipole moment, and ability to form hydrogen bond played the main role in both models. In addition, the aromatic π - π interactions and charge of the amide nitrogen had a major effect on insecticidal activity of the compounds. The present research can be helpful to gain a better understanding of the interactions between the insect AChE and its inhibitors and introduces compounds which are capable of becoming human-safe insecticides.

Chemometric modeling of Daphnia magna toxicity of agrochemicals

Over the past few years, the ecotoxicological hazard potential of agrochemicals has received much attention in the industries and regulatory agencies. In the current work, we have developed quantitative structure-activity relationship (QSAR) models for Daphnia magna toxicities of different classes of agrochemicals (fungicides, herbicides, insecticides and microbiocides) individually as well as for the combined set with the application of Organization for Economic Co-operation and Development (OECD) recommended guidelines. The models for the individual data sets as well as for the combined set were generated employing only simple and interpretable two-dimensional descriptors, and subsequently strictly validated using test set compounds. The validated individual models were used to generate consensus models, with the objective to improve the prediction quality and reduced prediction errors. All the individual models of different classes of agrochemicals as well as the global set of agrochemicals showed encouraging statistical quality and prediction ability. The general observations from the derived models suggest that the toxicity increases with lipophilicity and decreases with polarity. The generated models of different classes of agrochemicals and also for the combined set should be applicable for data gap filling for new or untested agrochemical compounds falling within the applicability domain of the developed models.

Next generation physiologically based kinetic (NG-PBK) models in support of regulatory decision making

The fields of toxicology and chemical risk assessment seek to reduce, and eventually replace, the use of animals for the prediction of toxicity in humans. In this context, physiologically based kinetic (PBK) modelling based on in vitro and in silico kinetic data has the potential to a play significant role in reducing animal testing, by providing a methodology capable of incorporating in vitro human data to facilitate the development of in vitro to in vivo extrapolation of hazard information. In the present article, we discuss the challenges in: 1) applying PBK modelling to support regulatory decision making under the toxicology and risk-assessment paradigm shift towards animal replacement; 2) constructing PBK models without in vivo animal kinetic data, while relying solely on in vitro or in silico methods for model parameterization; and 3) assessing the validity and credibility of PBK models built largely using non-animal data. The strengths, uncertainties, and limitations of PBK models developed using in vitro or in silico data are discussed in an effort to establish a higher degree of confidence in the application of such models in a regulatory context. The article summarises the outcome of an expert workshop hosted by the European Commission Joint Research Centre (EC-JRC) - European Union Reference Laboratory for Alternatives to Animal Testing (EURL ECVAM), on "Physiologically-Based Kinetic modelling in risk assessment - reaching a whole new level in regulatory decision-making" held in Ispra, Italy, in November 2016, along with results from an international survey conducted in 2017 and recently reported activities occurring within the PBK modelling field. The discussions presented herein highlight the potential applications of next generation (NG)-PBK modelling, based on new data streams.

Development of human biotransformation QSARs and application for PBT assessment refinement

Toxicokinetics heavily influence chemical toxicity as the result of Absorption, Distribution, Metabolism (Biotransformation) and Elimination (ADME) processes. Biotransformation (metabolism) reactions can lead to detoxification or, in some cases, bioactivation of parent compounds to more toxic chemicals. Moreover, biotransformation has been recognized as a key process determining chemical half-life in an organism and is thus a key determinant for bioaccumulation assessment for many chemicals. This study addresses the development of QSAR models for the prediction of in vivo whole body human biotransformation (metabolism) half-lives measured or empirically-derived for over 1000 chemicals, mainly represented by pharmaceuticals. Models presented in this study meet regulatory standards for fitting, validation and applicability domain. These QSARs were used, in combination with literature models for the prediction of biotransformation half-lives in fish, to refine the screening of the potential PBT behaviour of over 1300 Pharmaceuticals and Personal Care Products (PPCPs). The refinement of the PBT screening allowed, among others, for the identification of PPCPs, which were predicted as PBTs on the basis of their chemical structure, but may be easily biotransformed. These compounds are of lower concern in comparison to potential PBTs characterized by large predicted biotransformation half-lives.

QSARs for estimating intrinsic hepatic clearance of organic chemicals in humans

Quantitative structure-activity relationships (QSARs) were developed to predict the in vitro clearance (CLINT) of xenobiotics metabolised in human hepatocytes (118 compounds) and microsomes (115 compounds). Clearance values were gathered from the scientific literature and multiple linear models were built and validated selecting at most 6 predictors from a pool of over 2000 potential molecular descriptors. For the hepatocytes QSAR, the explained variance (Radj(2)) was 67% and the predictive ability (Rext(2)) was 62%. For the microsomes QSAR, Radj(2) was 50% and Rext(2) 30%. For both liver assays, the most important descriptor relates to electronic properties of the compound. Functional groups of fragments were useful to identify specific compounds that have a deviating reaction rate compared to the others, such as polychlorobiphenyls (PCBs) and organic amides which were poorly metabolised by hepatocytes and microsomes, respectively. For hepatocytes, clearance was predominantly determined by electronic characteristics, while size and shape characteristics were less important and partitioning properties were absent. This may suggest that uptake across the membrane and enzyme binding are not rate-limiting steps. Particularly for hepatocytes the QSAR statistics are encouraging, allowing application of the outcomes in in vitro to in vivo extrapolation.