Evaluation of SARs for the prediction of skin irritation/corrosion potential: structural inclusion rules in the BfR decision support system

The German Federal Institute for Risk Assessment (BfR) has developed a Decision Support System (DSS) to assess certain hazardous properties of pure chemicals, including skin and eye irritation/corrosion. The BfR-DSS is a rule-based system that could be used for the regulatory classification of chemicals in the European Union. The system is based on the combined use of two predictive approaches: exclusion rules based on physicochemical cut-off values to identify chemicals that do not exhibit a certain hazard (e.g., skin irritation/corrosion), and inclusion rules based on structural alerts to identify chemicals that do show a particular toxic potential. The aim of the present study was to evaluate the structural inclusion rules implemented in the BfR-DSS for the prediction of skin irritation and corrosion. The following assessments were performed: (a) a confirmation of the structural rules by rederiving them from the original training set (1358 substances), and (b) an external validation by using a test set of 200 chemicals not used in the derivation of the rules. It was found as a result that the test data set did not match the training set relative to the inclusion of structural alerts associated with skin irritation/corrosion, albeit some skin irritants were in the test set.     

A reactivity pattern for discrimination of ER agonism and antagonism based on 3-D molecular attributes

Various models have been developed to predict the relative binding affinity (RBA) of chemicals to estrogen receptors (ER). These models can be used to prioritize chemicals for further tiered biological testing to assess the potential for endocrine disruption. One shortcoming of models predicting RBA has been the inability to distinguish potential receptor antagonism from agonism, and hence in vivo response. It has been suggested that steroid receptor antagonists are less compact than agonists; thus, ER binding of antagonists may prohibit proper alignment of receptor protein helices preventing subsequent transactivation. The current study tests the theory of chemical bulk as a defining parameter of antagonism by employing a 3-D structural approach for development of reactivity patterns for ER antagonists and agonists. Using a dataset of 23 potent ER ligands (16 agonists, 7 antagonists), molecular parameters previously found to be associated with ER binding affinity, namely global ( E HOMO ) and local (donor delocalizabilities and charges) electron donating ability of electronegative sites and steric distances between those sites, were found insufficient to discriminate ER antagonists from agonists. However, parameters related to molecular bulk, including solvent accessible surface and negatively charged Van der Waal's surface, provided reactivity patterns that were 100% successful in discriminating antagonists from agonists in the limited data set tested. The model also shows potential to discriminate pure antagonists from partial agonist/antagonist structures. Using this exploratory model it is possible to predict additional chemicals for their ability to bind but inactivate the ER, providing a further tool for hypothesis testing to elucidate chemical structural characteristics associated with estrogenicity and anti-estrogenicity.     

A reactivity pattern for discrimination of ER agonism and antagonism based on 3-D molecular attributes

Various models have been developed to predict the relative binding affinity (RBA) of chemicals to estrogen receptors (ER). These models can be used to prioritize chemicals for further tiered biological testing to assess the potential for endocrine disruption. One shortcoming of models predicting RBA has been the inability to distinguish potential receptor antagonism from agonism, and hence in vivo response. It has been suggested that steroid receptor antagonists are less compact than agonists; thus, ER binding of antagonists may prohibit proper alignment of receptor protein helices preventing subsequent transactivation. The current study tests the theory of chemical bulk as a defining parameter of antagonism by employing a 3-D structural approach for development of reactivity patterns for ER antagonists and agonists. Using a dataset of 23 potent ER ligands (16 agonists, 7 antagonists), molecular parameters previously found to be associated with ER binding affinity, namely global (E(HOMO)) and local (donor delocalizabilities and charges) electron donating ability of electronegative sites and steric distances between those sites, were found insufficient to discriminate ER antagonists from agonists. However, parameters related to molecular bulk, including solvent accessible surface and negatively charged Van der Waal's surface, provided reactivity patterns that were 100% successful in discriminating antagonists from agonists in the limited data set tested. The model also shows potential to discriminate pure antagonists from partial agonist/antagonist structures. Using this exploratory model it is possible to predict additional chemicals for their ability to bind but inactivate the ER, providing a further tool for hypothesis testing to elucidate chemical structural characteristics associated with estrogenicity and anti-estrogenicity.     

Holographic QSAR of environmental estrogens

Environmental endocrine disrupting chemicals(EDCs) or environmental estrogens pose severe healthhazard to wildlife and humans. They are believed to bethe main cause of the weakening activity and secretionof sex hormone, sperm declination, abnormal repro-d…     

Comparative PBT screening using (Q)SAR tools within REACH legislation

Small to medium sized enterprises (SMEs) in the EU are facing challenges due to the introduction of new legislation designed to protect consumers and the environment, REACH (Registration, Evaluation, Authorisation and Restriction of CHemicals). There can be high costs associated with implementing REACH because data on mammalian toxicity, environmental toxicity and environmental fate properties is required and if this data is obtained experimentally the cost is significant. These costs can be reduced if reliable quantitative structure–activity relationships ((Q)SAR) models are instead used to obtain the required information. In this paper we investigate how easily freely available (Q)SAR models can be applied for persistent, bioaccumulative and toxic (PBT) screening of 17 chemicals of interest to SMEs. In this study the PBT predictions obtained from the more user-friendly PBT Profiler and the Danish(Q)SAR database for the chemicals were compared with the results taken directly from the EPI Suite software. It was found that these widely used (Q)SAR databases might have some errors and examples are provided. It was concluded that extra care must be taken when considering the use of these databases for PBT screening. In addition, to increase the likelihood of a correct prediction, data estimates from various (Q)SAR models relevant to the PBT endpoints must be compared.     

Evaluation of SARs for the prediction of skin irritation/corrosion potential–structural inclusion rules in the BfR decision support system†

The German Federal Institute for Risk Assessment (BfR) has developed a Decision Support System (DSS) to assess certain hazardous properties of pure chemicals, including skin and eye irritation/corrosion. The BfR–DSS is a rule-based system that could be used for the regulatory classification of chemicals in the European Union. The system is based on the combined use of two predictive approaches: exclusion rules based on physicochemical cut-off values to identify chemicals that do not exhibit a certain hazard (e.g., skin irritation/corrosion), and inclusion rules based on structural alerts to identify chemicals that do show a particular toxic potential. The aim of the present study was to evaluate the structural inclusion rules implemented in the BfR–DSS for the prediction of skin irritation and corrosion. The following assessments were performed: (a) a confirmation of the structural rules by rederiving them from the original training set (1358 substances), and (b) an external validation by using a test set of 200 chemicals not used in the derivation of the rules. It was found as a result that the test data set did not match the training set relative to the inclusion of structural alerts associated with skin irritation/corrosion, albeit some skin irritants were in the test set.     

Qualitative and quantitative structure–activity relationship modelling for predicting blood-brain barrier permeability of structurally diverse chemicals

In this study, structure–activity relationship (SAR) models have been established for qualitative and quantitative prediction of the blood–brain barrier (BBB) permeability of chemicals. The structural diversity of the chemicals and nonlinear structure in the data were tested. The predictive and generalization ability of the developed SAR models were tested through internal and external validation procedures. In complete data, the QSAR models rendered ternary classification accuracy of >98.15%, while the quantitative SAR models yielded correlation (r²) of >0.926 between the measured and the predicted BBB permeability values with the mean squared error (MSE) <0.045. The proposed models were also applied to an external new in vitro data and yielded classification accuracy of >82.7% and r² > 0.905 (MSE < 0.019). The sensitivity analysis revealed that topological polar surface area (TPSA) has the highest effect in qualitative and quantitative models for predicting the BBB permeability of chemicals. Moreover, these models showed predictive performance superior to those reported earlier in the literature. This demonstrates the appropriateness of the developed SAR models to reliably predict the BBB permeability of new chemicals, which can be used for initial screening of the molecules in the drug development process.     

On some novel extended topochemical atom (ETA) parameters for effective encoding of chemical information and modelling of fundamental physicochemical properties

Extended topochemical atom (ETA) indices developed by our group have been extensively applied in our previous reports for toxicity and ecotoxicity modelling in the field of quantitative structure–activity relationships (QSARs). In the present study these indices have been further explored by defining additional novel parameters to model n-octanol–water partition coefficient (two data sets; n?=?168 and 139), water solubility (n?=?193), molar refractivity (n?=?166), and aromatic substituent constants π, MR, σ _m, and σ _p (n?=?99). All the models developed in the present study have undergone rigorous internal and external validation tests and the models have high statistical significance and prediction potential. In terms of Q ² and r ² values the models developed for the datasets of whole molecules are better than those previously reported, with topochemically arrived unique (TAU) indices on the same datasets of chemicals. An attempt has also been made to develop models using non-ETA topological and information indices. Interestingly, ETA and non-ETA models have been found to have similar predictive capacity.     

QSAR models using a large diverse set of estrogens

Endocrine disruptors (EDs) have a variety of adverse effects in humans and animals. About 58,000 chemicals, most having little safety data, must be tested in a group of tiered assays. As assays will take years, it is important to develop rapid methods to help in priority setting. For application to large data sets, we have developed an integrated system that contains sequential four phases to predict the ability of chemicals to bind to the estrogen receptor (ER), a prevalent mechanism for estrogenic EDs. Here we report the results of evaluating two types of QSAR models for inclusion in phase III to quantitatively predict chemical binding to the ER. Our data set for the relative binding affinities (RBAs) to the ER consists of 130 chemicals covering a wide range of structural diversity and a 6 orders of magnitude spread of RBAs. CoMFA and HQSAR models were constructed and compared for performance. The CoMFA model had a r2 = 0.91 and a q2LOO = 0.66. HQSAR showed reduced performance compared to CoMFA with r2 = 0.76 and q2LOO = 0.59. A number of parameters were examined to improve the CoMFA model. Of these, a phenol indicator increased the q2LOO to 0.71. When up to 50% of the chemicals were left out in the leave-N-out cross-validation, the q2 remained significant. Finally, the models were tested by using two test sets; the q2pred for these were 0.71 and 0.62, a significant result which demonstrates the utility of the CoMFA model for predicting the RBAs of chemicals not included in the training set. If used in conjunction with phases I and II, which reduced the size of the data set dramatically by eliminating most inactive chemicals, the current CoMFA model (phase III) can be used to predict the RBA of chemicals with sufficient accuracy and to provide quantitative information for priority setting.     

Modelling of ready biodegradability based on combined public and industrial data sources

ABSTRACT

The European Registration, Evaluation, Authorization and Restriction of Chemical Substances Regulation, requires marketed chemicals to be evaluated for Ready Biodegradability (RB), considering in silico prediction as valid alternative to experimental testing. However, currently available models may not be relevant to predict compounds of industrial interest, due to accuracy and applicability domain restriction issues. In this work, we present a new and extended RB dataset (2830 compounds), issued by the merging of several public data sources. It was used to train classification models, which were externally validated and benchmarked against already-existing tools on a set of 316 compounds coming from the industrial context. New models showed good performances in terms of predictive power (Balance Accuracy (BA) = 0.74–0.79) and data coverage (83–91%). The Generative Topographic Mapping approach identified several chemotypes and structural motifs unique to the industrial dataset, highlighting for which chemical classes currently available models may have less reliable predictions. Finally, public and industrial data were merged into global dataset containing 3146 compounds. This is the biggest dataset reported in the literature so far, covering some chemotypes absent in the public data. Thus, predictive model developed on the Global dataset has larger applicability domain than the existing ones.