skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Estimating Octanol–Water Partition Coefficients of Novel Brominated Flame Retardants by Reversed-Phase High-Performance Liquid Chromatography and Computational Models
Abstract Legacy brominated flame retardants, including polybrominated diphenyl ethers (PBDEs), have been classified as persistent organic pollutants and replaced with novel brominated flame retardants (NBFRs). The octanol–water partition coefficients (log KOW) of NBFRs have been computationally estimated, but the log KOW values provided by these methods can differ by 1 to 3 orders of magnitude. Given the importance of this parameter in fate and toxicity models, we indirectly measured the log KOW values of eight NBFRs by their capacity factor (k′) on a reversed-phase high-performance liquid chromatography (HPLC) C18 column by isocratic elution and compared these measured values with those estimated by nine computational models. Log KOW values were obtained for the NBFRs 1,2-bis(2,4,6-tribromophenoxy) ethane, pentabromobenzene, pentabromoethylbenzene, pentabromotoluene, 2-ethylhexyl 2,3,4,5-tetrabromobenzoate, allyl 2,4,6-tribromophenylether, 2,3-dibromopropyl-2,4,6-tribromophenyl ether, and bis(2-ethylhexyl) tetrabromophthalate. A training set of phthalates, polychlorinated biphenyls, PBDEs, and halogenated benzenes were chosen to obtain the log k′–log KOW calibration for the NBFRs. The computational models KowWIN, XLogP3, EAS-E Suite, COSMOtherm, DirectML, and Abraham polyparameter linear free energy relationships all predicted the log KOW values of the calibration compounds to within 1 order of magnitude without significant bias. The median of these models predicted log KOW values for the calibration compounds that were close to those known in the literature with root mean square error (RMSE) = 0.224 and for the NBFRs that were close to those measured by HPLC (RMSE = 0.334). Environ Toxicol Chem 2024;43:2105–2114. © 2024 SETAC  more » « less
Award ID(s):
1804611
PAR ID:
10565003
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Environmental Toxicology and Chemistry
Volume:
43
Issue:
10
ISSN:
0730-7268
Format(s):
Medium: X Size: p. 2105-2114
Size(s):
p. 2105-2114
Sponsoring Org:
National Science Foundation
More Like this
  1. Chromatographic retention times and mass spectrometral fragmentation of per- and polyfluoroalkyl substances (PFASs) standards were determined using the optimized parameters obtained for liquid chromatography with tandem high-resolution mass spectrometry (LC-HRMS) analysis. Characteristic fragment ions obtained at various collision energies (MS2 fragmentation) were used for structural elucidation to predict the identities of newly discovered (emerging) PFASs detected in environmental samples. Moreover, the COnductor-like Screening MOdel for Realistic Solvents (COSMO-RS) was used to calculate the octanol-water partition coefficients (Kow) and mean isotropic polarizabilities of known PFASs, and the values were plotted against their chromatographic retention factors (k) to obtain a multivariable regression model that can be used to predict k values of unknown PFASs. Retention factor values of different structural isomers of the unknown PFASs were calculated and compared to the experimental k. For all the unknown PFASs, the predicted k value for the isomer that matches the corresponding MS2 fragmentation was found to be within 5% of the experimentally measured k value. This study demonstrates the applicability of a simple approach that combines the use of computationally-derived log Kow and polarizabilities, experimentally-determined k values, together with observed MS2 fragmentation patterns, in assigning the structures of emerging PFASs at environmentally relevant conditions when no reference standards are available. 
    more » « less
  2. Per-and polyfluoroalkyl substances (PFAS) are a class of contaminants of emerging concern frequently used in products like aqueous firefighting foams and non-stick coatings due to their stability and surfactant-like qualities. The lack of analytical standards for many emerging PFAS have severely limited our ability to comprehensively identify unknown PFAS contaminants in the environment, especially those that occur as isomers. Annotation of small molecules and identification of unknowns based only on elemental composition and mass fragmentation patterns remain major challenges in nontarget analysis employing liquid chromatography with high-resolution mass spectrometry (LC-HRMS). In this study, chromatographic retention factors (k) and mass spectral fragmentation patterns of 32 known PFAS were determined using our optimized parameters in LC-HRMS. The same method was then used to analyze previously unidentified PFAS in actual environmental samples. Using characteristic ions observed in the MS fragmentation of PFAS, the most probable isomeric structures of the detected PFAS were predicted. To increase confidence in the predicted molecular structure, Density Functional Theory and Conductor-like Screening Model for Realistic Solvents (COSMO-RS) calculations were used to predict physicochemical properties of different constitutional isomers. The DFT calculations facilitated geometric optimization, determination of polarizability, and calculation of the chemical potential the isomers. COSMO-RS uses the chemical potential to predict thermodynamic properties of molecules such as pKa, solubility, and Kow. These properties were then used to make a multi-variable linear regression to predict k values. The model was trained using 32 known PFAS. The properties used were log Kow of the neutral and anion species of the PFAS, and their polarizability. The model was specific enough to predict significantly different k values of unknown compounds with similar structures, which facilitated assignment of isomeric structures of PFAS. 
    more » « less
  3. In this study, we combine experiments, calculated properties, and machine learning (ML) to design new triphenylamine-based (TPA) molecules that have a high photoinduced radical (PIR) generation in crystals. A dataset of 34 crystal structures was extracted from the Cambridge Crystallographic Data Centre. Eighteen structures with experimentally reported PIR values from 0 to 0.85% were used to build an ML model trained using Random Forest that achieves an average leave-one-out test set error of 0.173% PIR. The ML model was used to screen the remaining 16 compounds, of which 4 were selected and sub-sequently compared with the experimentally measured PIR%. The predicted PIR% demonstrated good agreement with the measured values of TPA bis-urea macrocycles host-guest complexes and non-macrocyclic compounds of TPAs. Examining a broad set of molecular architectures/scaffolds allows for investigating the structural and electronic properties that lead to high PIR generation. We found very different trends for macrocycles, linear TPAs, and mono TPAs, where mono TPAs consist-ently have the lowest PIR generation. Macrocycles tend to have the highest PIR generation, especially for systems with ben-zene and fluorobenzene guests. Although linear analogs overall perform worse than macrocycles, they display clear trends with increasing excited-state dipole moment, oscillator strength and electron-hole covariance, while decreasing ionization potential and interatomic distance are generally correlated with higher PIRs. What is consistently observed is that higher PIRs are seen for brominated analogs. Our study, therefore, provides guidelines for future design strategies of TPAs for PIR generation. 
    more » « less
  4. Abstract Main‐chain boron‐containing π‐conjugated polymers are attractive for organic electronic, sensing, and imaging applications. Alternating terthiophene‐borane polymers were prepared and the effects of regioisomeric attachment of the conjugated linker and variations in the electronic effect of the pendent aryl groups (2,4,6‐tri‐tert‐butylphenyl, Mes*; 2,4,6‐tris(trifluoromethyl)phenyl, FMes) examined. Pd2dba3/P(t‐Bu)3‐catalyzed Stille polymerization of arylbis(2‐thienyl)borane and arylbis(3‐thienylborane) with 2,5‐bis(trimethylstannyl)thiophene at 120 °C gave polymers with appreciable molecular weight but MALDI‐TOF MS analyses showed evidence of unusually prominent homocoupling. These defects could be suppressed by using brominated rather than iodinated monomers, more hindered 2,5‐bis(tri‐n‐butylstannyl)thiophene as comonomer, and Pd2dba3/P(o‐tol)3as the catalyst at 100 °C. Under these conditions, macrocyclic species withn=3–10 repeating units formed preferentially according to MALDI‐TOF MS analyses. Photophysical studies revealed a prominent effect of the regiochemistry and the nature of the pendent aryl groups on the absorption and emission, giving rise to orange, yellow‐green, blue‐green, and blue emissive materials respectively. The electronic effects were rationalized through DFT calculations on bis(terthiophene) model systems. 
    more » « less
  5. Predicting adsorption of organic pollutants onto graphene nanomaterials is not only useful for exploring their potential adsorbent applications, but also helpful for better understanding their fate and risks in aquatic environments. Herein molecular dynamics (MD) simulations and theoretical linear solvation energy relationships (TLSERs) were employed to construct prediction models for adsorption of neutral organic pollutants onto graphene and graphene oxides. The MD simulations for adsorption of 43 aromatic compounds onto graphene and diverse models of graphene oxides with various functional groups (hydroxyl, epoxy and carbonyl) demonstrate that graphene has a stronger affinity for the aromatic compounds than graphene oxides. The hydroxyl and carbonyl groups of graphene oxides were found to form hydrogen bonds with the aromatic adsorbates, while epoxy groups did not. TLSER models were developed for predicting the adsorption equilibrium coefficients ( K ) onto graphene and graphene oxide nanosheets. In the graphene prediction model, H-donating ability ( ε α ) and dispersion/hydrophobic interactions ( V ) have significant effects on log  K values, while in the graphene oxide model, ε α is the most influential factor on log  K values. The models provide in silico approaches for predicting adsorption affinities onto graphenic nanomaterials. 
    more » « less