Comparison of predictions of developmental toxicity for compounds of solvent data set.

We have considered a series of 235 compounds technically classified as solvents. Chemically, they belong to different classes. Their potential developmental toxicity was evaluated using two models available on platform VEGA HUB; model CAESAR and the Developmental/ Reproductive Toxicity library (PG) model. Models provide beside the prediction of developmental toxicity additional information on similar compounds from models' training sets. In the report, first, we compare the predictions. Second, the sets of similar compounds have been used to implement the clustering scheme. The Kohonen artificial neural network method has been applied as a clustering method. The clusters obtained have been discussed for both models.

Chemometrical analysis of proteomics data obtained from three cell types treated with multi-walled carbon nanotubes and TiO2 nanobelts$.

Applications of nanomaterials in biomedical, industrial, and consumer goods areas are expanding rapidly because of their unique physicochemical properties. Hazard assessment of nanosubstances is necessary for the protection of human and ecological health. We studied the proteomics patterns of three cell lines: co-culture of Caco-2 and HT29-MTX cells, primary small airway epithelial cells, and THP-1macrophage-like cells. The cells were exposed at 10 µg and 100 µg concentrations for 3 and 24 hours to multi-walled carbon nanotubes and TiO2 nanobelts (TiO2-NB). The data were analysed with the hierarchical clustering method and principal components analysis. In all cases, time of exposure is the most important factor in separation and clustering of proteomics patterns. Furthermore, the sets of proteins, which are specific for long (24 hours) exposure, are identified.

Grouping of CoMPARA data with respect to compounds from the carcinogenic potency database.

Methods for clustering and measures of similarity of chemical structures have become an important supporting tool in chemoinformatics. They represent the basis for categorization of chemicals and read-across, where a molecular property is estimated from 'similar molecules'. This study proposes a clustering scheme within the given dataset with respect to a reference dataset. The scheme was applied on two datasets ToxCast_AR_Agonist and ToxCast_AR_Antagonists with 1654 and 1522 compounds, respectively. The compounds are tested to androgen receptor activity (AR) in 11 high throughput screening assays. The carcinogenic dataset was used as the reference set.

Integration of in silico methods and computational systems biology to explore endocrine-disrupting chemical binding with nuclear hormone receptors.

Thousands of potential endocrine-disrupting chemicals present difficult regulatory challenges. Endocrine-disrupting chemicals can interfere with several nuclear hormone receptors associated with a variety of adverse health effects. The U.S. Environmental Protection Agency (U.S. EPA) has released its reviews of Tier 1 screening assay results for a set of pesticides in the Endocrine Disruptor Screening Program (EDSP), and recently, the Collaborative Estrogen Receptor Activity Prediction Project (CERAPP) data. In this study, the predictive ability of QSAR and docking approaches is evaluated using these data sets. This study also presents a computational systems biology approach using carbaryl (1-naphthyl methylcarbamate) as a case study. For estrogen receptor and androgen receptor binding predictions, two commercial and two open source QSAR tools were used, as was the publicly available docking tool Endocrine Disruptome. For estrogen receptor binding predictions, the ADMET Predictor, VEGA, and OCHEM models (specificity: 0.88, 0.88, and 0.86, and accuracy: 0.81, 0.84, and 0.88, respectively) were each more reliable than the MetaDrug™ model (specificity 0.81 and accuracy 0.77). For androgen receptor binding predictions, the Endocrine Disruptome and ADMET Predictor models (specificity: 0.94 and 0.8, and accuracy: 0.78 and 0.71, respectively) were more reliable than the MetaDrug™ model (specificity 0.33 and accuracy 0.4). A consensus approach is proposed that reaches general agreement among the models (specificity 0.94 and accuracy 0.89). This study integrates QSAR, docking, and systems biology approaches as a virtual screening tool for use in risk assessment. As such, this systems biology pathways and network analysis approach provides a means to more critically assess the potential effects of endocrine-disrupting chemicals.

Robust modelling of acute toxicity towards fathead minnow (Pimephales promelas) using counter-propagation artificial neural networks and genetic algorithm.

Large worldwide use of chemicals has caused great concern about their possible adverse effects on human health, flora and fauna. Increased production of new chemicals has also increased demand for their risk assessment. Traditionally, results from animal tests have been used to assess toxicity of chemicals. However, such methods are ethically questionable since they involve killing and causing suffering of the test animals. Therefore, new in silico methods are being sought to replace the traditional in vivo and in vitro testing methods. In this article we report on one method that can be used to build robust models for the prediction of compounds' properties from their chemical structure. The method has been developed by combining a genetic algorithm, a counter-propagation artificial neural network and cross-validation. It has been tested using existing data on toxicity to fathead minnow (Pimephales promelas). The results show that the method may give reliable results for chemicals belonging to the applicability domain of the developed models. Therefore, it can aid the risk assessment of chemicals and consequently reduce demand for animal tests.

Prediction of mutagenicity and carcinogenicity using in silico modelling: A case study of polychlorinated biphenyls.

In silico modelling is an important alternative method for the evaluation of properties of chemical compounds. Basically, two concepts are used in its applications: QSAR modelling for endpoint predictions, and grouping (categorization) of large groups of chemicals. In the presented report we address both of these concepts. As a case study we present the results on a set of polychlorinated biphenyls (PCBs) and some of their metabolites. Their mutagenicity and carcinogenic potency were evaluated with CAESAR and T.E.S.T. models, which are freely available over the internet. We discuss the value and reliability of the predictions, the applicability domain of models and the ability to create prioritized groupings of PCBs as a category of chemicals.

QSTR modelling of the acute toxicity of pharmaceuticals to fish.

Extensive use of pharmaceuticals as human and veterinary medication raises concerns for their adverse effects on non-target organisms. The purpose of this study was to employ multiple linear regression (MLR) to predict the toxicities of a diverse set of pharmaceuticals to fish. The descriptor pool consisted of about 1500 descriptors calculated using Dragon 5.4, Spartan 06 and Codessa 2.2 software. Descriptor selection was made by the heuristic method available in Codessa 2.2. The data set was divided into training and test sets using Kohonen networks. The training set contained approximately 65% of the compounds of the full data set (99 compounds). The training set model contained eight descriptors from all dimensions, all of which were obtained from Dragon 5.4. The statistical parameters of the model for the training set are R(2 )= 0.664, F = 13.588, and R(cv)(2) (LOO) = 0.542 while it achieves R(2 )= 0.605 for the test set. The training, test and external sets have no response outliers considering the standardized residual greater than three. The external validation of the model was made with a set of pharmaceuticals obtained from several databases. The R(pred)(2) is 0.777, reflecting a relatively good predictive power for the external set.

Counter propagation artificial neural network categorical models for prediction of carcinogenicity for non-congeneric chemicals.

One of the main goals of the new chemical regulation REACH (Registration, Evaluation and Authorization of Chemicals) is to fill the gaps on the toxicological properties of chemicals that affect human health. Carcinogenicity is one of the endpoints under consideration. The information obtained from (quantitative) structure-activity relationship ((Q)SAR) models is accepted as an alternative solution to avoid expensive and time-consuming animal tests. The reported results were obtained within the framework of the European project 'Computer Assisted Evaluation of industrial chemical Substances According to Regulations (CAESAR)'. In this article, we demonstrate intermediate results for counter propagation artificial neural network (CP ANN) models for the prediction category of the carcinogenic potency using two-dimensional (2D) descriptors from different software programs. A total of 805 non-congeneric chemicals were extracted from the Carcinogenic Potency Database (CPDBAS). The resulting models had prediction accuracies for internal (training) and external (test) sets as high as 91-93% and 68-70%, respectively. The sensitivity and specificity of the test set were 69-73 and 63-72% correspondingly. High specificity is critical in models for regulatory use that are aimed at ensuring public safety. Thus, the errors that give rise to false negatives are much more relevant. We discuss how we can increase the number of correctly predicted carcinogens using the correlation between the threshold and the values of the sensitivity and specificity.

Assessing the reproductive toxicity of some (con)azole compounds using a structure-activity relationship approach.

The present research investigates the study of a set of 27 (con)azoles and their reproductive toxicity. (Con)azoles are used as fungicides and herbicides in agriculture for treatment of fruits, vegetables, cereals, and seeds, or as human antimycotic therapeutics. According to EEC Directive 91/414, active substances used in plant protection products must undergo reproductive toxicity testing. Reproductive toxicity is a complex biological endpoint, which includes many different biological processes and, therefore, it can only to a limited extent be assessed by a single quantitative structure-activity relationship (QSAR) model. The proposed SAR models are built using unsupervised methods, such as hierarchical clustering, principal component analysis and self-organizing maps, with the aim of studying the similarity relationships between structures. The molecular structures are represented with a set of topological and structural descriptors. The models showing clusters, closest neighbours or outliers may support the categorization and the classification of (con)azoles as potential reproductive toxicants.

Spectral representation of reduced protein models.

We consider a spectrum-like two-dimensional graphical representation of proteins based on a reduced protein model in which 20 amino acids are grouped into five classes. This particular grouping of amino acids was suggested by Riddle and co-workers in 1997. The graphical representation is based on depicting sequentially the amino acids on five horizontal lines at equal separations. One-letter codes, B, O, U, X and Y, to which numerical values 1 to 5 have been assigned, are suggested as labels for the fictional amino acids that represent all the amino acids within each group. The approach is illustrated on ND6 proteins of eight species having from 168 to 175 amino acids. While visual inspection of the novel spectral graphical representations of proteins may reveal local similarities and dissimilarities of protein sequences, arithmetic manipulations of spectra offer an elegant route to graphic visualization of the degree of similarity for selected pairs of proteins.

On novel representation of proteins based on amino acid adjacency matrix.

A novel characterization of proteins is presented based on selected properties of recently introduced 20 x 20 amino acid adjacency matrix of proteins in which matrix elements count the occurrence of all 400 possible pair-wise adjacencies obtained by reading protein primary sequence from the left to the right. In particular we consider the characterization based on the sum and the difference of the rows and the corresponding columns, which characterize proteins by a pair of 20-component vectors. The approach is illustrated on a set of ND6 proteins of eight species.

Validation of counter propagation neural network models for predictive toxicology according to the OECD principles: a case study.

The OECD has proposed five principles for validation of QSAR models used for regulatory purposes. Here we present a case study investigating how these principles can be applied to models based on Kohonen and counter propagation neural networks. The study is based on a counter propagation network model that has been built using toxicity data in fish fathead minnow for 541 compounds. The study demonstrates that most, if not all, of the OECD criteria may be met when modeling using this neural network approach.

Structure-mutagenicity study of 12 trimethylimidazopyridine isomers using orbital energies and "spectrum-like representation" as descriptors.

The set of 12 trimethylimidazopyridine isomers with mutagenic potency toward two strains of Salmonella was treated in this study. Ten isomers with known mutagenic properties were taken to build the models. Fifteen molecular orbital energies, or a "spectrum-like" representation of 3D structures, were taken as descriptors. As modeling techniques the multiple linear regression and the counter propagation neural network were applied. Models were tested with the recall ability test and the leave-one-out cross-validation tests. For two isomers, which have not been synthesized yet, we report predicted values for both mutagenic potencies obtained with different models. The best models were found when unoccupied molecular orbital energies are among the descriptors.

On 3-D graphical representation of DNA primary sequences and their numerical characterization.

In this article we (1) outline the construction of a 3-D "graphical" representation of DNA primary sequences, illustrated on a portion of the human beta globin gene; (2) describe a particular scheme that transforms the above 3-D spatial representation of DNA into a numerical matrix representation; (3) illustrate construction of matrix invariants for DNA sequences; and (4) suggest a data reduction based on statistical analysis of matrix invariants generated for DNA. Each of the four contributions represents a novel development that we hope will facilitate comparative studies of DNA and open new directions for representation and characterization of DNA primary sequences.

A study of structure-carcinogenicity relationship for 86 compounds from NTP data base using topological indices as descriptors.

86 compounds from NTP carcinogenic potency data base have been used to derive neural network models. Compounds were described with topological indices. Carcinogenicity has been given as a binary quantity--a compound is carcinogenic or non carcinogenic. Several models have been tested with a recognition ability test and with the leave-one-out cross validation method. For the best model the ratio between correct and wrong classifications was 70/30. Furthermore, the model has been used to classify 17 compounds not used for setting of the models. The predicted carcinogenic classes and the neighbors in the neural network influencing the predictions have been discussed.

On the similarity of DNA primary sequences.

We consider numerical characterization of graphical representations of DNA primary sequences. In particular we consider graphical representation of DNA of beta-globins of several species, including human, on the basis of the approach of A. Nandy in which nucleic bases are associated with a walk over integral points of a Cartesian x, y-coordinate system. With a so-generated graphical representation of DNA, we associate a distance/distance matrix, the elements of which are given by the quotient of the Euclidean and the graph theoretical distances, that is, through the space and through the bond distances for pairs of bases of graphical representation of DNA. We use eigenvalues of so-constructed matrices to characterize individual DNA sequences. The eigenvalues are used to construct numerical sequences, which are subsequently used for similarity/dissimilarity analysis. The results of such analysis have been compared and combined with similarity tables based on the frequency of occurrence of pairs of bases.

A study of the excited states in cytosine and guanine stacks in the Hartree-Fock and exciton approximations.

We report calculated exciton energies for the cytosine and guanine stacks obtained in the ab initio Hartree-Fock crystal orbital and exciton approximation, which includes the excited electron-hole interaction. This interaction plays an important role in the description of excited electron spectra in the low-energy region. The stacks were chosen as examples of polymers with helical symmetry.