skip to main content


Title: Computational protocol for predicting 19F NMR chemical shifts for PFAS and connection to PFAS structure
Abstract

Per‐ and polyfluoroalkyl substances (PFAS) are robust “forever” chemicals that have become global environmental contaminants due to their inability to degrade using traditional techniques. In addition to the persistent nature of PFAS, the structural and functional diversity in PFAS creates a unique challenge in identification and remediation. Their identification is further complicated by the absence of standards for many PFAS. This work is aimed at developing a protocol for computing and establishing accurate19F NMR chemical shifts for PFAS using density functional theory (DFT), which can aid in the identification of PFAS. The impact of solvation and basis sets was evaluated by comparing the computed data with the experimental measurements. Results showed the addition of dispersion corrections in the methodology improve the accuracy of calculated NMR parameters within 4 ppm of the experimental values. Adding a second diffuse function and additional polarization did not improve the accuracy, likely because of the electronegativity of fluorine which does not allow the electron density of fluorine atoms to be polarized. The inclusion of various implicit solvation (DMSO, chloroform, and water) yielded negligible differences in accuracy, and were overall less accurate than the gas phase calculations. The most accurate methodology was then applied to more environmentally relevant PFAS, and the impact of helical nature on the NMR signatures was evaluated. The implication of this work is to be able to improve the identification of structurally diverse PFAS using the19F NMR.

 
more » « less
Award ID(s):
1710079 2109210
NSF-PAR ID:
10445223
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Journal of Computational Chemistry
Volume:
43
Issue:
20
ISSN:
0192-8651
Page Range / eLocation ID:
p. 1355-1361
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The aqueous photolysis of four pharmaceuticals with varying fluorinated functional groups was assessed under neutral, alkaline, advanced oxidation, and advanced reduction conditions with varying light sources. Solar simulator quantum yields were 2.21 × 10−1 mol Ei−1for enrofloxacin, 9.36 × 10−3 mol Ei−1for voriconazole, and 1.49 × 10−2 mol Ei−1for flecainide. Florfenicol direct photolysis was slow, taking 150 h for three degradation half‐lives. Bimolecular rate constants between pharmaceuticals and hydroxyl radicals were 109to 1010 M−1 s−1. Using a combined quantitative fluorine nuclear magnetic resonance spectroscopy (19F‐NMR) and mass spectrometry approach, fluorine mass balances and photolysis product structures were elucidated. Enrofloxacin formed a variety of short‐lived fluorinated intermediates that retained the aryl F motif. Extended photolysis time led to complete aryl F mineralization to fluoride. The aliphatic F moiety on florfenicol was also mineralized to fluoride, but the resulting product was a known antibiotic (thiamphenicol). For voriconazole, the two aryl Fs contributed more to fluoride production compared with the heteroaromatic F, indicating higher stability of the heteroaromatic F motif. The two aliphatic CF3moieties in the flecainide structure remained intact under all conditions, further supporting the stability of these moieties found in per‐ and polyfluoroalkyl substances under a variety of conditions. The advanced treatment conditions generating hydroxyl radicals or hydrated electrons accelerated the degradation, but not the defluorination, of flecainide. The combination of19F‐NMR and mass spectrometry proved powerful in allowing identification of fluorinated products and verifying the functional groups present in the intermediates and products. The results found in the present study will aid in the understanding of which fluorinated functional groups should be incorporated into pharmaceuticals to ensure organofluorine byproducts are not formed in the environment and help determine the water‐treatment processes that effectively remove specific pharmaceuticals and more generally fluorinated motifs.Environ Toxicol Chem2023;00:1–12. © 2023 The Authors.Environmental Toxicology and Chemistrypublished by Wiley Periodicals LLC on behalf of SETAC.

     
    more » « less
  2. Abstract

    Accurate estimation of solvation free energy (SFE) lays the foundation for accurate prediction of binding free energy. The Poisson‐Boltzmann (PB) or generalized Born (GB) combined with surface area (SA) continuum solvation method (PBSA and GBSA) have been widely used in SFE calculations because they can achieve good balance between accuracy and efficiency. However, the accuracy of these methods can be affected by several factors such as the charge models, polar and nonpolar SFE calculation methods and the atom radii used in the calculation. In this work, the performance of the ABCG2 (AM1‐BCC‐GAFF2) charge model as well as other two charge models, that is, RESP (Restrained Electrostatic Potential) and AM1‐BCC (Austin Model 1‐bond charge corrections), on the SFE prediction of 544 small molecules in water by PBSA/GBSA was evaluated. In order to improve the performance of the PBSA prediction based on the ABCG2 charge, we further explored the influence of atom radii on the prediction accuracy and yielded a set of atom radius parameters for more accurate SFE prediction using PBSA based on the ABCG2/GAFF2 by reproducing the thermodynamic integration (TI) calculation results. The PB radius parameters of carbon, oxygen, sulfur, phosphorus, chloride, bromide and iodine, were adjusted. New atom types,on,oi,hn1,hn2,hn3, were introduced to further improve the fitting performance. Then, we tuned the parameters in the nonpolar SFE model using the experimental SFE data and the PB calculation results. By adopting the new radius parameters and new nonpolar SFE model, the root mean square error (RMSE) of the SFE calculation for the 544 molecules decreased from 2.38 to 1.05 kcal/mol. Finally, the new radius parameters were applied in the prediction of protein‐ligand binding free energies using the MM‐PBSA method. For the eight systems tested, we could observe higher correlation between the experiment data and calculation results and smaller prediction errors for the absolute binding free energies, demonstrating that our new radius parameters can improve the free energy calculation using the MM‐PBSA method.

     
    more » « less
  3. Abstract

    We report the first spectroscopic evidence for a [C−F−C]+fluoronium ion in solution. Extensive NMR studies (19F,1H,13C) characterize a symmetric cage‐like species in which fluorine exhibits substantial covalent bonding to each of the two carbon atoms involved in the three‐center interaction. Experimental NMR data comport well with calculated values to lend credence to the structural assignment. As the culminating experiment, a Saunders isotopic perturbation test confirmed the symmetric structure. Congruent with the trend in other types of onium ions, the calculated charge at fluorine moves in a more positive (less negative) direction from the neutral. It is this important trend that explains in part the extraordinary historical difficulty in making theoretical predictions of fluoronium ions come true in solution, and why it takes fluorine captured in a cage to produce, finally, a stable ion and complete the historical arc of the organic halonium ion story.

     
    more » « less
  4. Abstract

    We report the first spectroscopic evidence for a [C−F−C]+fluoronium ion in solution. Extensive NMR studies (19F,1H,13C) characterize a symmetric cage‐like species in which fluorine exhibits substantial covalent bonding to each of the two carbon atoms involved in the three‐center interaction. Experimental NMR data comport well with calculated values to lend credence to the structural assignment. As the culminating experiment, a Saunders isotopic perturbation test confirmed the symmetric structure. Congruent with the trend in other types of onium ions, the calculated charge at fluorine moves in a more positive (less negative) direction from the neutral. It is this important trend that explains in part the extraordinary historical difficulty in making theoretical predictions of fluoronium ions come true in solution, and why it takes fluorine captured in a cage to produce, finally, a stable ion and complete the historical arc of the organic halonium ion story.

     
    more » « less
  5. Abstract

    19F NMR spectroscopy is an attractive and growing area of research with broad applications in biochemistry, chemical biology, medicinal chemistry, and materials science. We have explored fast magic angle spinning (MAS)19F solid‐state NMR spectroscopy in assemblies of HIV‐1 capsid protein. Tryptophan residues with fluorine substitution at the 5‐position of the indole ring were used as the reporters. The19F chemical shifts for the five tryptophan residues are distinct, reflecting differences in their local environment. Spin‐diffusion and radio‐frequency‐driven‐recoupling experiments were performed at MAS frequencies of 35 kHz and 40–60 kHz, respectively. Fast MAS frequencies of 40–60 kHz are essential for consistently establishing19F–19F correlations, yielding interatomic distances of the order of 20 Å. Our results demonstrate the potential of fast MAS19F NMR spectroscopy for structural analysis in large biological assemblies.

     
    more » « less