Search

The bioconcentration factor (BCF) is an important parameter that gives information regarding the ability of a contaminant to be taken up by organisms from the water. Per- and polyfluoroalkyl substances (PFAS) are widespread in the environment, causing concern regarding their impact on human health. Due to the lack of available bioaccumulation data for most compounds in the PFAS group, we developed a quantitative structure-property relationship (QSPR) model to predict the log BCF for fish (taxonomic class Teleostei), based on experimental data available for the most studied 33 representatives of this group of compounds. Furthermore, we implemented the developed model to predict log BCF for an external dataset of 2209 PFAS. Consequently, 1045 PFAS were found not to be bioaccumulative, 208 were classified as bioaccumulative, and 956 were predicted to be very bioaccumulative. Finally, we obtained the high correlation (R 2  = 0.844) between the log BCFs obtained in laboratory and field studies for 13 PFAS. In silico analyses indicate that PFAS bioconcentration depends on the size (chain length - number of CF 2 groups in alkyl tail/chain) of a molecule, as well as on the atomic distribution properties. In general, long-chain PFAS - above 8 and 6 carbon atoms for perfluorinated carboxylic acids (PFCAs)and perfluorinated sulfonic acids (PFSAs), respectively - tend to bioconcentrate more compared to the short-chain ones. In conclusion, predicting BCF on fish is possible for a wide range of fluorinated compounds, which can be further used for estimating PFAS behavior in the environment., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier Ltd.)

Per- and polyfluorinated alkyl substances (PFAS), known for their widespread environmental presence and slow degradation, pose significant concerns. Of the approximately 10,000 known PFAS, only a few have undergone comprehensive testing, resulting in limited experimental data. In this study, we employed a combination of physics-based methods and data-driven models to address gaps in PFAS bioaccumulation potential. Using the COnductor-like Screening MOdel for Realistic Solvents (COSMO-RS) method, we predicted n-octanol/water partition coefficients (logK OW ), crucial for PFAS bioaccumulation. Our developed Quantitative Structure-Property Relationship (QSPR) model exhibited high accuracy (R 2  = 0.95, RMSE C  = 0.75) and strong predictive ability (Q 2 LOO  = 0.93, RMSE CV  = 0.83). Leveraging the extensive NORMAN, we predicted logK OW for over 4,000 compounds, identifying 244 outliers out of 4519. Further categorizing the database into eight Organisation for Economic Co-operation and Development (OECD) categories, we confirmed fluorine atoms role in enhanced bioaccumulation. Utilizing predicted logK OW , water solubility logS W , and vapor pressure logVP values, we calculated additional physicochemical properties that are responsible for the transport and dispersion of PFAS in the environment. Parameters such as Henry's Law (k H ), air-water partition coefficient (K AW ), octanol-air coefficient (K OA ), and soil adsorption coefficient (K OC ) exhibited favorable correlations with literature data (R 2  > 0.66). Our study successfully filled data gaps, contributing to the understanding of ubiquitous PFAS in the environment and estimating missing physicochemical data for these compounds., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

The rapid development of engineered nanomaterials (ENMs) causes humans to become increasingly exposed to them. Therefore, a better understanding of the health impact of ENMs is highly demanded. Considering the 3Rs (Replacement, Reduction, and Refinement) principle, in vitro and computational methods are excellent alternatives for testing on animals. Among computational methods, nano-quantitative structure-activity relationship (nano-QSAR), which links the physicochemical and structural properties of EMNs with biological activities, is one of the leading method. The nature of toxicological experiments has evolved over the last decades; currently, one experiment can provide thousands of measurements of the organism's functioning at the molecular level. At the same time, the capacity of the in vitro systems to mimic the human organism is also improving significantly. Hence, the authors would like to discuss whether the nano-QSAR approach follows modern toxicological studies and takes full advantage of the opportunities offered by modern toxicological platforms. Challenges and possibilities for improving data integration are underlined narratively, including the need for a consensus built between the in vitro and the QSAR domains., (© 2023 Wiley-VCH GmbH.)

Nanoinformatics demands accurate predictive models to assess the potential hazards of nanomaterials (NMs). However, limited data availability and the diverse nature of NMs physicochemical properties and their interaction with biological media, hinder the development of robust nano-Quantitative Structure-Activity Relationship (QSAR) models. This article proposes an approach that combines artificially data generation techniques and topological projections to address the challenges of insufficient dataset sizes and their limited representativeness of the chemical space. By leveraging the rich information embedded in the topological features, this methodology enhances the representation of the chemical space, enabling a more an exploration of the structure-activity relationships. We demonstrate the efficacy of our approach through extensive experiments, employing various machine learning regression algorithms to validate the methodology. Finally, we compare two different resampling approaches based on different modeling scenarios. The results showcase a significant improved predictive performance of QSAR models demonstrating a promising strategy to overcome the limitations of small datasets in the field of nanoinformatics. The proposed approach offers noteworthy potential for advancing nanoinformatics research within the nanosafety domain by enabling the development of more accurate predictive models for assessing the potential hazards associated with NMs.

The acidity of Lewis-Brønsted superacids can be derived from the theoretical calculations as the Gibbs free energy of the deprotonation reaction (ΔG acid ), which describes the tendency of a studied compound to donate a proton. This paper presents the first Quantitative Structure - Property Relationship (QSPR) model that correlates the ΔG acid of superacid (HF/MeX 3 formula (X=F, Cl, Br)) with their structure. Developed model is well fitted, roubustness, has good predictive abilities, fulfills all OECD recommendation for good model. Obtained results provide the insight into the relation of structural features of superacids, which are responsible for their acid strength - the structures characterized by strong F-Me dative bond (with relatively large vibrational frequency), small positive partial atomic charge on Me central atom, possibly large polarity exhibit large acid strength. Such assumption can be used in the future as valuable information in the process of the designing new, stronger, more effective superacids., (© 2019 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Nanotechnology and the production of nanomaterials have been expanding rapidly in recent years. Since many types of engineered nanoparticles are suspected to be toxic to living organisms and to have a negative impact on the environment, the process of designing new nanoparticles and their applications must be accompanied by a thorough risk analysis. (Quantitative) Structure-Activity Relationship ([Q]SAR) modelling creates promising options among the available methods for the risk assessment. These in silico models can be used to predict a variety of properties, including the toxicity of newly designed nanoparticles. However, (Q)SAR models must be appropriately validated to ensure the clarity, consistency and reliability of predictions. This paper is a joint initiative from recently completed European research projects focused on developing (Q)SAR methodology for nanomaterials. The aim was to interpret and expand the guidance for the well-known "OECD Principles for the Validation, for Regulatory Purposes, of (Q)SAR Models", with reference to nano-(Q)SAR, and present our opinions on the criteria to be fulfilled for models developed for nanoparticles., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

Nanotechnology is one of the most important technological developments of the 21st century. In silico methods to predict toxicity, such as quantitative structure-activity relationships (QSARs), promote the safe-by-design approach for the development of new materials, including nanomaterials. In this study, a set of cytotoxicity experimental data corresponding to 19 data points for silica nanomaterials were investigated, to compare the widely employed CORAL and Random Forest approaches in terms of their usefulness for developing so-called 'nano-QSAR' models. 'External' leave-one-out cross-validation (LOO) analysis was performed, to validate the two different approaches. An analysis of variable importance measures and signed feature contributions for both algorithms was undertaken, in order to interpret the models developed. CORAL showed a more pronounced difference between the average coefficient of determination (R²) for training and for LOO (0.83 and 0.65 for training and LOO, respectively), compared to Random Forest (0.87 and 0.78 without bootstrap sampling, 0.90 and 0.78 with bootstrap sampling), which may be due to overfitting. With regard to the physicochemical properties of the nanomaterials, the aspect ratio and zeta potential were found to be the two most important variables for Random Forest, and the average feature contributions calculated for the corresponding descriptors were consistent with the clear trends observed in the data set: less negative zeta potential values and lower aspect ratio values were associated with higher cytotoxicity. In contrast, CORAL failed to capture these trends., (2016 FRAME.)

In this paper, we suggest that causal inference methods could be efficiently used in Quantitative Structure-Activity Relationships (QSAR) modeling as additional validation criteria within quality evaluation of the model. Verification of the relationships between descriptors and toxicity or other activity in the QSAR model has a vital role in understanding the mechanisms of action. The well-known phrase "correlation does not imply causation" reflects insight statistically correlated with the endpoint descriptor may not cause the emergence of this endpoint. Hence, paradigmatic shifts must be undertaken when moving from traditional statistical correlation analysis to causal analysis of multivariate data. Methods of causal discovery have been applied for broader physical insight into mechanisms of action and interpretation of the developed nano-QSAR models. Previously developed nano-QSAR models for toxicity of 17 nano-sized metal oxides towards E. coli bacteria have been validated by means of the causality criteria. Using the descriptors confirmed by the causal technique, we have developed new models consistent with the straightforward causal-reasoning account. It was proven that causal inference methods are able to provide a more robust mechanistic interpretation of the developed nano-QSAR models.

Computational techniques, such as Quantitative Structure-Property Relationship (QSPR) modeling, are very useful in predicting physicochemical properties of various chemicals. Building QSPR models requires calculating molecular descriptors and the proper choice of the geometry optimization method, which will be dedicated to specific structure of tested compounds. Herein, we examine the influence of the ionic liquids' (ILs) geometry optimization methods on the predictive ability of QSPR models by comparing three models. The models were developed based on the same experimental data on density collected for 66 ionic liquids, but with employing molecular descriptors calculated from molecular geometries optimized at three different levels of the theory, namely: (1) semi-empirical (PM7), (2) ab initio (HF/6-311+G*) and (3) density functional theory (B3LYP/6-311+G*). The model in which the descriptors were calculated by using ab initio HF/6-311+G* method indicated the best predictivity capabilities ([Formula: see text] = 0.87). However, PM7-based model has comparable values of quality parameters ([Formula: see text] = 0.84). Obtained results indicate that semi-empirical methods (faster and less expensive regarding CPU time) can be successfully employed to geometry optimization in QSPR studies for ionic liquids.

Creating suitable chemical categories and developing read-across methods, supported by quantum mechanical calculations, can be an effective solution to solving key problems related to current scarcity of data on the toxicity of various nanoparticles. This study has demonstrated that by applying a nano-read-across, the cytotoxicity of nano-sized metal oxides could be estimated with a similar level of accuracy as provided by quantitative structure-activity relationship for nanomaterials (nano-QSAR model(s)). The method presented is a suitable computational tool for the preliminary hazard assessment of nanomaterials. It also could be used for the identification of nanomaterials that may pose potential negative impact to human health and the environment. Such approaches are especially necessary when there is paucity of relevant and reliable data points to develop and validate nano-QSAR models.

Regulation for nanomaterials is urgently needed, and the drive to adopt an intelligent testing strategy is evident. Such a strategy will not only provide economic benefits but will also reduce moral and ethical concerns arising from animal testing. For regulatory purposes, such an approach is promoted by REACH, particularly the use of quantitative structure-activity relationships [(Q)SAR] as a tool for the categorisation of compounds according to their physicochemical and toxicological properties. In addition to compounds, (Q)SAR has also been applied to nanomaterials in the form of nano(Q)SAR. Although (Q)SAR in chemicals is well established, nano(Q)SAR is still in early stages of development and its successful uptake is far from reality. This article aims to identify some of the pitfalls and challenges associated with nano-(Q)SARs in relation to the categorisation of nanomaterials. Our findings show clear gaps in the research framework that must be addressed if we are to have reliable predictions from such models. Three major barriers were identified: the need to improve quality of experimental data in which the models are developed from, the need to have practical guidelines for the development of the nano(Q)SAR models and the need to standardise and harmonise activities for the purpose of regulation. Of these three, the first, i.e. the need to improve data quality requires immediate attention, as it underpins activities associated with the latter two. It should be noted that the usefulness of data in the context of nano-(Q)SAR modelling is not only about the quantity of data but also about the quality, consistency and accessibility of those data.

The development of quantitative structure-activity relationships for nanomaterials needs representation of molecular structure of extremely complex molecular systems. Obviously, various characteristics of nanomaterial could impact associated biochemical endpoints. Following features of TiO2 and ZnO nanoparticles (n=42) are considered here: (i) engineered size (nm); (ii) size in water suspension (nm); (iii) size in phosphate buffered saline (PBS, nm); (iv) concentration (mg/L); and (v) zeta potential (mV). The damage to cellular membranes (units/L) is selected as an endpoint. Quantitative features-activity relationships (QFARs) are calculated by the Monte Carlo technique for three distributions of data representing values associated with membrane damage into the training and validation sets. The obtained models are characterized by the following average statistics: 0.78

Nanotechnology has evolved as a frontrunner in the development of modern science. Current studies have established toxicity of some nanoparticles to human and environment. Lack of sufficient data and low adequacy of experimental protocols hinder comprehensive risk assessment of nanoparticles (NPs). In the present work, metal electronegativity (χ), the charge of the metal cation corresponding to a given oxide (χox), atomic number and valence electron number of the metal have been used as simple molecular descriptors to build up quantitative structure-toxicity relationship (QSTR) models for prediction of cytotoxicity of metal oxide NPs to bacteria Escherichia coli. These descriptors can be easily obtained from molecular formula and information acquired from periodic table in no time. It has been shown that a simple molecular descriptor χox can efficiently encode cytotoxicity of metal oxides leading to models with high statistical quality as well as interpretability. Based on this model and previously published experimental results, we have hypothesized the most probable mechanism of the cytotoxicity of metal oxide nanoparticles to E. coli. Moreover, the required information for descriptor calculation is independent of size range of NPs, nullifying a significant problem that various physical properties of NPs change for different size ranges., (Copyright © 2014 Elsevier Inc. All rights reserved.)

As experimental evaluation of the safety of nanoparticles (NPs) is expensive and time-consuming, computational approaches have been found to be an efficient alternative for predicting the potential toxicity of new NPs before mass production. In this background, we have developed here a regression-based nano quantitative structure-activity relationship (nano-QSAR) model to establish statistically significant relationships between the measured cellular uptakes of 109 magnetofluorescent NPs in pancreatic cancer cells with their physical, chemical, and structural properties encoded within easily computable, interpretable and reproducible descriptors. The developed model was rigorously validated internally as well as externally with the application of the principles of Organization for Economic Cooperation and Development (OECD). The test for domain of applicability was also carried out for checking reliability of the predictions. Important fragments contributing to higher/lower cellular uptake of NPs were identified through critical analysis and interpretation of the developed model. Considering all these identified structural attributes, one can choose or design safe, economical and suitable surface modifiers for NPs. The presented approach provides rich information in the context of virtual screening of relevant NP libraries., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

Quantitative structure-property/activity relationships (QSPRs/QSARs) are a tool to predict various endpoints for various substances. The "classic" QSPR/QSAR analysis is based on the representation of the molecular structure by the molecular graph. However, simplified molecular input-line entry system (SMILES) gradually becomes most popular representation of the molecular structure in the databases available on the Internet. Under such circumstances, the development of molecular descriptors calculated directly from SMILES becomes attractive alternative to "classic" descriptors. The CORAL software (http://www.insilico.eu/coral) is provider of SMILES-based optimal molecular descriptors which are aimed to correlate with various endpoints. We analyzed data set on nanoparticles uptake in PaCa2 pancreatic cancer cells. The data set includes 109 nanoparticles with the same core but different surface modifiers (small organic molecules). The concept of a QSAR as a random event is suggested in opposition to "classic" QSARs which are based on the only one distribution of available data into the training and the validation sets. In other words, five random splits into the "visible" training set and the "invisible" validation set were examined. The SMILES-based optimal descriptors (obtained by the Monte Carlo technique) for these splits are calculated with the CORAL software. The statistical quality of all these models is good., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

Nowadays nanotechnology is one of the most promising areas of science. The number and quantity of synthesized nanomaterials increase exponentially, therefore it is reasonable to expect that comprehensive risk assessment based only on empirical testing of all novel engineered nanoparticles (NPs) will very soon become impossible. Hence, the development of computational methods complementary to experimentation is very important. Quantitative structure-property relationship (QSPR) and quantitative structure-activity relationship (QSAR) models widely used in pharmaceutical chemistry and environmental science can also be modified and adopted for nanotechnology to predict physico-chemical properties and toxicity of empirically untested nanomaterials. All QSPR/QSAR modelling activities are based on experimentally derived data. It is important that, within a given data set, all values should be consistent, of high quality and measured according to a standardized protocol. Unfortunately, the amount of such data available for engineered nanoparticles in various data sources (i.e. databases and the literature) is very limited and seldom measured with a standardized protocol. Therefore, we have proposed a framework for collecting and evaluating the existing data, with the focus on possible applications for computational evaluation of properties and biological activities of nanomaterials.

Official OECD recommendations give the highest priority to application of purely empirical phys/chem data (partition coefficients, environmental half-live times etc.) in multimedia mass balance modeling of environmental overall persistence and long-range transport for potentially hazardous chemicals. We have demonstrated that the replacement of the empirical data with those predicted by employing Quantitative Structure-Property Relationships (QSPR) technique did not significantly decrease the performance of The Tool 2.0--the OECD multimedia mass balance model. To prove this, we compared each other the output results (overall persistence -P(OV); characteristic traveling distance -CTD and transport efficiency TE) obtained from 6 of multimedia models. The models utilized combinations of experimentally determined and QSPR-predicted values of the partition coefficients and half-live times. For predicting phys/chem data, we utilized 2 QSPRs developed in our laboratory and the EPI Suite package (US EPA). We did not observe any statistically significant (p<0.05) differences between the models. This conclusion is important, because it leads to reducing time and costs of laboratory studies required during the risk assessment procedure. Moreover, regarding the obtained results, we proposed to replace the single-threshold approaches established by majority of international regulations to screen substances for persistence, bioaccumulation and long-range transport potential with the approaches taking into account uncertainty of the results and/or probability of passing a given threshold., (Copyright © 2011 Elsevier B.V. All rights reserved.)

We have estimated degradation half-lives of both brominated and chlorinated dibenzo-p-dioxins (PBDDs and PCDDs), furans (PBDFs and PCDFs), biphenyls (PBBs and PCBs), naphthalenes (PBNs and PCNs), diphenyl ethers (PBDEs and PCDEs) as well as selected unsubstituted polycyclic aromatic hydrocarbons (PAHs) in air, surface water, surface soil, and sediments (in total of 1,431 compounds in four compartments). Next, we compared the persistence between chloro- (relatively well-studied) and bromo- (less studied) analogs. The predictions have been performed based on the quantitative structure-property relationship (QSPR) scheme with use of k-nearest neighbors (kNN) classifier and the semi-quantitative system of persistence classes. The classification models utilized principal components derived from the principal component analysis of a set of 24 constitutional and quantum mechanical descriptors as input variables. Accuracies of classification (based on an external validation) were 86, 85, 87, and 75% for air, surface water, surface soil, and sediments, respectively. The persistence of all chlorinated species increased with increasing halogenation degree. In the case of brominated organic pollutants (Br-OPs), the trend was the same for air and sediments. However, we noticed that the opposite trend for persistence in surface water and soil. The results suggest that, due to high photoreactivity of C-Br chemical bonds, photolytic processes occurring in surface water and soil are able to play significant role in transforming and removing Br-OPs from these compartments. This contribution is the first attempt of classifying together Br-OPs and Cl-OPs according to their persistence, in particular, environmental compartments.

The most significant achievements and challenges relating to an application of quantitative structure-activity relationship (QSAR) approach in the risk assessment of nanometer-sized materials are highlighted. Recent advances are discussed in the context of "classical" QSAR methodology. The possible ways for the structural characterization of compounds existing at the nanoscale (at least one dimension of 100 nm or less) are briefly reviewed. The applicability of the existing toxicological data for developing QSAR models is evaluated. Finally, the existing models are presented. The need to develop new interpretative descriptors for the nanosystems is also highlighted. It is suggested that, due to high variability in the molecular structures and different mechanisms of toxicity, individual classes of nanoparticles should be modeled separately.

Information about the mobility and the partitioning properties of brominated persistent organic pollutants, the environmental levels of which are sometimes higher than those of the chlorinated analogues, is limited. We estimated n-octanol/ water (log K(OW)), n-octanol/air (log K(OA)), and air/water (log K(AW)) partition coefficients for 1436 chloro- and bromo-analogues of dibenzo-p-dioxins, dibenzofurans, biphenyls, naphthalenes, diphenyl ethers, and benzenes by employing quantitative structure-property relationship (QSPR) techniques. The searches for similar partitioning patterns were performed by means of two-dimensional cluster analysis. Five classes of compounds were identified. Each of the class is characterized by similar partition coefficients and, in consequence, similar environmental properties. Finally the data was fitted into a simple multimedia model involving the partitioning map. In addition, we found thatthe changes in the partition coefficients upon the replacement of chlorine with bromine were constant: 0.11, 0.31, and -0.21 per bromine atom for log K(OW), log K(OA), and log K(AW), respectively. On the basis of this observation, a method for rapid estimation of changes in the partition coefficient upon chlorine-bromine substitution was proposed.

Molecular descriptors from calculations at the level of Density Functional Theory (B3LYP/6-311++G**) were effectively applied in QSPR estimation of supercooled liquid vapour pressures (P(L)) for individual chloronaphthalene congeners. The estimated values of log P(L) varied from 1.05 Pa to 5.6 x 10(-5) Pa, depending on the number of chlorine substituents present in the molecule and the substitution pattern. Comparison of the five chemometrical methods of modelling (approaches) led to the final conclusion, that the use of relatively simple PLS combined with one of the variable pre-selection algorithms (UVE or GA) seems to be the optimal choice in such computational studies for persistent organic pollutants. The best GA-PLS model was characterized by the value of root mean square error of prediction RMSEP = 0.108 logarithmic Pascal units.