Environmental impact of PFAS: Filling data gaps using theoretical quantum chemistry and QSPR modeling.

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 (logKOW), crucial for PFAS bioaccumulation. Our developed Quantitative Structure-Property Relationship (QSPR) model exhibited high accuracy (R2 = 0.95, RMSEC = 0.75) and strong predictive ability (Q2LOO = 0.93, RMSECV = 0.83). Leveraging the extensive NORMAN, we predicted logKOW 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 logKOW, water solubility logSW, 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 (kH), air-water partition coefficient (KAW), octanol-air coefficient (KOA), and soil adsorption coefficient (KOC) exhibited favorable correlations with literature data (R2 > 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.

Expanding the applicability domain of QSPRs for predicting water solubility and vapor pressure of PFAS.

This work aimed to verify whether it is possible to extend the applicability domain (AD) of existing QSPR (Quantitative Structure-Property Relationship) models by employing a strategy involving additional quantum-chemical calculations. We selected two published QSPR models: for water solubility, logSW, and vapor pressure, logVP of PFAS as case studies. We aimed to enlarge set of compounds used to build the model by applying factorial planning to plan the augmentation of the set of these compounds based on their structural features (descriptors). Next, we used the COSMO-RS model to calculate the logSW and logVP for selected chemicals. This allowed filling gaps in the experimental data for further training QSPR models. We improved the published models by significantly extending number of compounds for which theoretical predictions are reliable (i.e., extending the AD). Additionally, we performed external validation that had not been carried out in original models. To test effectiveness of the AD extension, we screened 4519 PFAS from NORMAN Database. The number of compounds outside the domain was reduced comparing the original model for both properties. Our work shows that combining physics-based methods with data-driven models can significantly improve the performance of predictions of phys-chem properties relevant for the chemical risk assessment.

How the Structure of Per- and Polyfluoroalkyl Substances (PFAS) Influences Their Binding Potency to the Peroxisome Proliferator-Activated and Thyroid Hormone Receptors-An In Silico Screening Study.

In this study, we investigated PFAS (per- and polyfluoroalkyl substances) binding potencies to nuclear hormone receptors (NHRs): peroxisome proliferator-activated receptors (PPARs) α, ß, and γ and thyroid hormone receptors (TRs) α and ß. We have simulated the docking scores of 43 perfluoroalkyl compounds and based on these data developed QSAR (Quantitative Structure-Activity Relationship) models for predicting the binding probability to five receptors. In the next step, we implemented the developed QSAR models for the screening approach of a large group of compounds (4464) from the NORMAN Database. The in silico analyses indicated that the probability of PFAS binding to the receptors depends on the chain length, the number of fluorine atoms, and the number of branches in the molecule. According to the findings, the considered PFAS group bind to the PPARα, ß, and γ only with low or moderate probability, while in the case of TR α and ß it is similar except that those chemicals with longer chains show a moderately high probability of binding.

Multivariate assessment of anticancer oleanane triterpenoids lipophilicity.

Naturally occurring molecules are excellent sources of lead compounds. A series of oleanolic acid (OA) derivatives previously synthesized in our laboratory, which show promising antitumor activity, have been analyzed in terms of lipophilicity evaluation applying chromatographic and computational approaches. Retention data obtained on three reversed-phase liquid chromatography stationary phases (RP-HPLC) and immobilized artificial membrane chromatography (IAM-HPLC) were compared with computational methods using chemometric tools such as cluster analysis, principal component analysis and sum of ranking differences. To investigate the molecular mechanism of retention quantitive structure retention relationship analysis was performed, based on the genetic algorithm coupled with multiple linear regression (GA-MLR). The obtained results suggested that the ionization potential of studied molecules significantly affects their retention in classical RP-HPLC. In IAM-HPLC additionally, polarizability-related descriptors also play an essential role in that process. The lipophilicity indices comparison shows significant differences between the computational lipophilicity and chromatographically determined ones.