Next Article in Journal
The Synthetic Cannabinoids THJ-2201 and 5F-PB22 Enhance In Vitro CB1 Receptor-Mediated Neuronal Differentiation at Biologically Relevant Concentrations
Previous Article in Journal
Pathophysiology of Type 2 Diabetes Mellitus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Indirect Photodegradation of Sulfamethoxazole and Trimethoprim by Hydroxyl Radicals in Aquatic Environment: Mechanisms, Transformation Products and Eco-Toxicity Evaluation

1
Environment Research Institute, Shandong University, Qingdao 266237, China
2
College of Biological and Environmental Engineering, Binzhou University, Binzhou 256600, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(17), 6276; https://doi.org/10.3390/ijms21176276
Submission received: 11 August 2020 / Revised: 26 August 2020 / Accepted: 28 August 2020 / Published: 30 August 2020
(This article belongs to the Section Physical Chemistry and Chemical Physics)

Abstract

:
The bacteriostatic antibiotics, sulfamethoxazole (SMX) and trimethoprim (TMP), have frequently been found in wastewater and surface water, which raises the concerns about their ecotoxicological effects. The indirect photochemical transformation has been proven to be an efficient way to degrade SMX and TMP. In this study, the reaction mechanisms of the degradation by SMX and TMF by OH radicals were investigated by theoretical calculations. Corresponding rate constants were determined and the eco-toxicity of SMX and TMP and its degradations products were predicted using theoretical models. The results indicate that the most favorable pathways for the transformation of SMX and TMP are both •OH-addition reaction of benzene ring site with lowest Gibbs free energy barriers (6.86 and 6.21 kcal mol−1). It was found that the overall reaction rate constants of •OH-initial reaction of SMX and TMP are 1.28 × 108 M−1 s−1 and 6.21 × 108 M−1 s−1 at 298 K, respectively. When comparing the eco-toxicity of transformation products with parent SMX and TMP, it can be concluded that the acute and chronic toxicities of the degraded products are reduced, but some products remain harmful for organisms, especially for daphnid (toxic or very toxic level). This study can give greater insight into the degradation of SMX and TMP by •OH through theoretical calculations in aquatic environment.

Graphical Abstract

1. Introduction

Pharmaceuticals and personal care products (PPCPs) as emerging organic micropollutants in the environment have been receiving increased attention due to their possible adverse effects on the aquatic organism and human health [1,2,3,4,5]. Sulfamethoxazole (SMX) and trimethoprim (TMP), as important antibiotics, have been widely used in the treatment of human and animal diseases and infections [6,7,8,9]. Excessive amounts of SMX and TMP are discharged into water for a wide range transmission to wastewater treatment plants [10], which cannot be effectively removed by conventional wastewater treatment technologies [11,12,13,14]. This has led to widespread transport of pharmaceutical contaminants in aquatic environment around the world [15].
SMX has regularly been detected in wastewater effluents and surface water with a concentration range of 100–2500 ng/L and 60–150 ng/L, even in drinking water at 12 ng/L [16,17,18]. Previous studies have reported that the residual concentration of TMP are 0.1–5 µg/L in wastewater effluents. Additionally, traces of TMP (2.2–10.9 ng/L and 0–19.8 ng/L) have been detected in surface and drinking water [19,20,21,22,23]. Based on the concentrations of SMX and TMP in the aquatic environment, it may induce the ecotoxicological concern about the chronic exposure of bacteria to trace level of antibiotics [24]. More importantly, the transformed products may retain the same level of biological properties as original compounds or even have a higher biological activity. Also, previous studies about the transformation of TMP exhibit some differences [25,26,27]. Thus, the chemical transformation of SMX and TMP deserve more attention due to their biological activity and adverse effects on aquatic ecosystems.
Much efforts have been spent on researching the degradation and environmental impact of SMX and TMP in aquatic environment. When compared with biodegradation, the photodegradation of SMX and TMP is identified as a predominant pathway in surface water [28,29,30]. Previous studies have focused on the degradation rate constant of SMX direct photolysis, which is 22.3 ± 1.86 h−1 at 3.2 pH [31]. In addition to the direct photolysis of pollutants, indirect photodegradation has been regarded as an important elimination way of organic contaminants, such as SMX and TMP in aquatic environment. In fact, advanced oxidation processes (AOPs) are broadly applied for removing organic materials of wastewaters [32,33]. Hydroxyl radical-based AOPs (such as UV/H2O2) treatments are strongly correlated with removal of pharmaceuticals with the rate constants on the order of 108–109 M−1 s−1 [34]. The experimental results shown that the rate constant for the •OH-initiated reaction of SMX is (5.8 ± 0.2) × 109 M−1 s−1 [18,35]. In addition, many researchers have studied the degradation rate constants of SMX and TMP with carbonate radical (k = 2.68 × 108 and 3.45 × 107 M−1 s−1, respectively) [36]. When combined with the experimental research, a deeper understanding of the performance of SXM/TMP transformation mechanisms and kinetics is warranted.
Therefore, the objective of this work is to make a detailed elucidation for the indirect photochemical transformation and environmental impacts of SMX and TMP with •OH in aquatic environment by quantum chemical calculation and computational toxicology method. In this paper, the transformation mechanism and kinetic calculation of •OH-initiated reaction are expounded thoroughly. Furthermore, the ecotoxicity risk of SMX and TMP and their products during indirect photolysis process is evaluated while using the ECOSAR program. The study will provide theoretical insight into the indirect photochemical behavior and ecotoxicity prediction of SMX and TMP in aquatic environments.

2. Results and Discussion

2.1. Transformation Mechanisms of SMX and TMP with •OH

In this study, the indirect photochemical transformations of SMX and TMP towards •OH in aqueous phase were investigated. For the initial reaction of SMX and TMP with •OH, the degradation mechanisms are described in Figure 1. •OH attack SMX and TMP through the three basic pathways: (i) hydrogen abstraction (Rabs), •OH extracts one H atom from SMX and TMP to form the dehydrogenated radical intermediates and H2O; (ii) hydroxylation reaction (Radd), •OH adds to the C atom of reactants to form •OH-complexes; and, (iii) electron transfer (Rset), the single-electron transfer from the reactants to •OH and generating the radical cation of reactants and •OH.

2.1.1. •OH-Initiated Reaction of SMX

For the Rset process of SMX, the radical cation species (SMX+•) and hydroxyl anion (OH) were formed via the electron transfer. The reaction of SMX with •OH is an endothermic process by Rset pathway, which absorbs 26.33 kcal mol−1 of reaction energy. The results indicate that the Rset pathway is not a spontaneous reaction, and it is not easy to occur by single-electron transfer mechanism of SMX from a thermodynamic standpoint. In addition, the Marcus theory was used to calculate the energy barrier of the Rset pathway ( Δ G S E T = 33.08 kcal mol−1), further implying that the reactant has a low reactivity during this process, and Rset process has low reaction activity during •OH-initiated degradation of SMX. In general, Rset pathway is not the primary mechanism during the •OH-initiated indirect photodegradation process in the aquatic environment. Thus, the two mechanisms (H-abstraction reaction and •OH-addition reaction) were taken into consideration, primarily in the following discussion.
•OH reacts easily with olefinic double bonds and anilines owing to the high reactivity. As shown in Figure 1, H atom of SMX can be divided to five categories (benzene ring, isoxazole, methyl, and two nitrogen-atoms), which provides the different reaction site for H-abstraction reaction (Rabs). Additionally, hydroxylation reaction (Radd) can occur at the benzene ring and isoxazole ring. The free energy profile for the •OH-initiated mechanism of the most favorable pathways of SMX and the optimized geometries of the major categories (namely, Radd3, Radd5, Radd8, Rabs1a, Rabs8a, Rabs10–12a) are provided in Figure 2. Figure S1 provides other pathways and Gibbs free energies. It is worth noting that the •OH-addition products are different at C3 site of benzene ring from other sites. The C-S bond breaks when •OH is added to C3 site, leading to the phenolic product. The Gibbs free energy barrier of the addition process (Radd3) is 7.83 kcal mol−1. For the Rabs pathways, the free energy barriers are within a range of 9.71~16.82 kcal mol−1, which is higher than the Radd pathways (6.86~9.96 kcal mol−1). The results implying that, compared with H-abstraction reactions, hydroxylation pathways of SMX have higher reactivity in aquatic environments. In addition, the reaction energies of all •OH-initiated reaction are negative (−3.11~−28.39 kcal mol−1), except a small positive number (1.48 kcal mol−1) at C8 atom of abstraction reaction. The results mean that the Radd and Rabs pathways are exothermic processes, and they can be more favorable from the standpoint of thermodynamics. Comparing the thermodynamic and kinetic data of •OH-initiated reaction at different sites of SMX, the addition reaction occurred at C5 site on the benzene ring is the most advantageous channel for the degradation of SMX. The conclusion is verified by the lowest Gibbs free energy barrier (6.86 kcal mol−1) and the highest rate constants (6.18 × 107 M−1 s−1) of Radd5 pathway. In summary, Radd and Rabs pathways predictably play an important role in the •OH-initiated photodegradation mechanisms of SMX in an aquatic environment.

2.1.2. •OH-Initiated Reaction of TMP

The indirect photolysis mechanisms of •OH-initiated reaction of TMP are similar with the photochemical transformations of SMX. For the Rset pathway, the single-electron transfer from the TMP to •OH. The radical cation intermediate (TMP+•) and hydroxyl anion (OH) are formed in this process. The Gibbs free energy of activation ( Δ G S E T ) of the SET reaction is 23.13 kcal mol−1 calculated by Marcus theory at 298.15 K. Simultaneously, the Rset mechanism of •OH-initiated reaction is endothermic with 21.88 kcal mol−1 of reaction energy. The results mean that the electron transfer reactions are not easy to occur from a thermodynamic point of view.
Figure 1 also shows that the hydroxylation reaction (Radd) and H-abstraction reaction (Rabs) of TMP could occur at different sites (including benzene ring, five-membered heterocycle, and substituent group). Figure 3 clarifies the profile of the free energy for the •OH-initiated mechanisms of the most favorable pathways of the three major categories. Additionally, other channels of •OH-initiated reactions are described in Figure S2. The four pathways (Radd1–4) of the benzene ring for the hydroxylation reactions are considered due to the symmetric structure of TMP. The lowest Gibbs free energy barrier of the addition reaction is 6.21 kcal mol−1 at C1 site, and the process releases up to 20.15 kcal mol−1 of reaction energy, as shown in Figure 3. For Rabs channels, the Gibbs free energy barrier of Rabs15a (substituent group C15 site) is the lowest. The reaction needs to overcome the transition state (TS15a) with the energy barrier of 7.07 kcal mol−1 and release heart of 22.69 kcal mol−1. Furthermore, among the overall addition and abstraction pathways, the reaction energies that are released by •OH-initiated reaction range from −0.90 to −36.53 kcal mol−1. Meanwhile, the Gibbs free energy barriers of Radd and Rabs pathways are 6.21~18.42 kcal mol−1. All of the intermediates of •OH-initiated reactions of TMP have one free electron, which can lead to the occurrence of subsequent reaction.
When considering all of the •OH-initiated photolysis reaction of SMX and TMP, it is clear that hydroxylation reaction and H-abstraction mechanisms mainly contribute to the photochemical transformation of SMX and TMP in an aquatic environment. In particularly, the most favorable pathway of the transformation is hydroxylation reaction of benzene ring site.

2.2. Kinetic Calculation

The kinetic calculation is carried out using the TST method in order to further confirm the dominant pathways at the suitable temperature ranging of 273–328 K. The rate constants of the •OH-initiated degradation of SMX and TMP are calculated based on the thermodynamic data in aquatic environment. The bimolecular rate constants of the •OH-initiated reaction are calculated according to the TST method, as shown before. The initial reaction of •OH can be described as:
SMX   ( TMP ) + OH   k 1   OH   adducts   or   dehydrogenated   intermediates
The total rate constant can be expressed using the formula:
ktotal = k1
The rate constants of •OH-addition and H-abstraction reaction are defined as kadd and kabs, respectively. The rate constants (kadd, kabs, and ktotal) of the •OH-initiated reaction of SXM and TMP at different site are displayed in Tables S1 and S2, in detail. The overall rate constants (ktotal) at 273–328 K and experimental rate constants of SMX and TMP are given in Table 1. The calculated overall rate constant of SMX and TMP are 1.28 × 108 M−1 s−1 and 6.21 × 108 M−1 s−1 at 298 K, respectively. Additionally, the rate constants of bimolecular reaction of SMX and TMP with •OH are obtained by the experimental method [18,37]. The experimental results of SMX and TMP ((5.8 ± 0.2) × 109 M−1 s−1 and 8.66 × 109 M−1 s−1, respectively.) are an order of magnitude higher than the theoretical values, which may be due to the influence of concentration of reactants, PH, aqueous medium, catalyst, and some intermediates in the experiments. In addition, as the temperature increases, the total rate constants of SMX and TMP degradation decreases. For example, the total rate constant for TMP of the •OH-initiated reaction is 9.72 × 108 M−1 s−1 at 273 K, which is 2.4 times higher than that at 328 K (4.03 × 108 M−1 s−1). The highest rate constants of •OH-initiated reaction are kadd5 (C5 site of SMX), kadd1 (C1 site of TMP), which are corresponding with the thermodynamic results, as shown in Tables S1 and S2. Thus, the initiated reaction at C5 site of SMX (or at C1 site of TMP) is most favorable.
Furthermore, the contributions of initiated reactions at different site to the total reactions of SMX are expressed as a rate branching ratio, which can be expressed as rx =kX/ktotal. The calculated results rx are shown in Figure 4. The three paths (Radd1, Radd3, and Radd5) of hydroxylation and two paths (Rabs10a and Rabs11a) of H-abstraction reactions that are most likely to occur are select to characterize the branching ratio of rate constants. It is obvious that the Radd1 and Radd5 pathways of •OH-addition reactions play a dominant role in the photochemical transformation. The rate branching ratio of these two pathways is determined as high as 89% at 298 K, and the ratio of Radd3 (C3 site) is 9% at 298 K. In a previous study, Solar et al. [38] found that the reaction of •OH with the ortho-position (C1 or C5) of aniline is predominant (account for 54%), and 10% with para-position (C3) by the absorbance characteristics experiment. The results confirmed the rate branching ratio that was carried out by theoretical calculation. Thus, the hydroxylation reaction is significant to evaluate the subsequent transformation reaction of SMX.

2.3. Subsequent Transformation Reaction of SMX and TMP

The intermediates generated from the •OH-initiated reaction of SXM and TMP contain a free radical, which can undergo the subsequent reactions due to the high reactivity. For the purpose of identifying the transformation products of SMX during the indirect photolysis process, the intermediate IM5 was selected to study the subsequent transformation, owing to the fastest reaction rate constant. Figure 5 shows the profile of free energy for subsequent reaction of IM5. The subsequent reactions of intermediate IM5 are considered in the presence of O2 in aquatic environment, as shown in Figure 5. The reaction of IM5 with O2 is exothermic, which releases the reaction energy of 3.25 kcal mol−1. In the process, the Gibbs free energy barrier of transition state TS-O2 is 11.61 kcal mol−1. Subsequently, the intermediate IM-O2 can undergo the hydrogen peroxide radical stripping reaction with the Gibbs free energy barrier of 15.54 kcal mol−1, leading to the formation of SMX transformation product (TP2). Simultaneously, the process can release the reaction energy of 28.36 kcal mol−1. The product was identified by the triple quadrupole TOF mass spectrometer in the experiment [35].
Furthermore, the intermediate IM1 of TMP is considered in the photodegradation process. Figure 6 depicts the transformation pathways of IM1 with blue color. IM1 can undergo demethylation reaction to form the IM1-1 and •OCH3 radical. The reaction requires the free energy barrier of 17.97 kcal mol−1, and the channel is endothermic with a little reaction energy of 5.08 kcal mol−1. Subsequently, the product TP01 is formed through hydroxylation reaction at C7 site. The formation process includes two steps, corresponding to the reactions IM1-1 + •OH IM1-2, and IM1-2 TP01. In this process, the intermediate IM1-2 is produced through the transition state TS1-2 with the Gibbs free energy barrier of 6.27 kcal mol−1. In addition, the next addition reaction of IM1-2 with •OH is barrier-less process and a lot of reaction energy of 67.04 kcal mol−1 is released at the same time. The ketone product TP02 is generated through oxidation reaction by losing of a molecular H2. In this process, the oxidizing agent captures the H from the carbon of C7 site and •OH, and it forms the C = O bond. The reaction is not easy to occur with a high Gibbs free energy barrier of 80.27 kcal mol−1 in theory. However, in natural water environment, the reactivity may be enhanced by the catalytic action of oxidant. Next, TP02 reacts with H2O via the hydrolysis reaction (Gibbs free energy barriers of 78.57 kcal mol−1), into TP03 (2,4-diaminopyrimidine-5-carboxylic acid, DAPC), and TP04 simultaneously with releasing energy of 6.69 kcal mol−1. The transformation products (TP01, TP02 and TP03) of TMP were verified by HRMS and MS2 method in experiment [39].
The mechanism studies elucidate that •OH can attack the organic contaminant (SMX and TMP) and the two intermediates (IM5 of SMX and IM1 of TMP) are easily formed through the •OH-addition pathway, according to above discussion. Furthermore, the transformed products of SMX (TP1 and TP2), also with the photolysis products of TMP (TP 01, TP02, TP03, and TP04), are formed in the indirect photochemical transformation process. The ecological risks of these products need to be further considered in the aquatic environment.

2.4. Eco-Toxocity Evaluation

In this work, QSAR analysis calculated by ECOSAR program was applied to predict the aquatic toxicity for fish, daphnid, and green algae, including acute toxicity and chronic toxicity. In addition, the values of 96-h median lethal concentration (LC50) for fish, 48-h LC50 for daphnid, and 96-h median effective concentration (EC50) for green algae were used to analysis the acute toxicity. Additionally ChV is defined as the geometric mean of the unobserved effect concentration and the lowest observed effect concentration. Based on the results of EC50 for D. magna and LC50 for O. latipes, the values of aniline derivative (unhindered) are selected for preliminary prediction of SMX, TMP, and transformation products [40,41].
Figure 7 illustrates the acute and chronic toxicities of SMX, TMP, and transformation products. Obviously, the acute and chronic toxicity of SMX, TP1, and TP2 to the three aquatic organisms express the same trend. The chronic toxicity of SMX and TPs for daphnid is very low, less than 0.1 mg L−1, which classified as very toxic substances according to Chinese hazard evaluation guidelines (Table S3). Additionally the product TP2 shows lower toxicity than SMX, which may be due to the increased hydrophilicity of the hydroxylated product of SMX.
The aquatic toxicity of most transformation products for fish and green alage are lower than TMP, as shown in Figure 7. The results indicate that the •OH-initiated degradation process of TMP reduced the toxicity of products. Although the acute toxicity decreases with the transformation of TMP, some of the products are classified as toxic and harmful to daphnid green alage. Similar with SMX, the daphnid is the most sensitive organism for aquatic toxicity. In summary, the eco-toxicity assessment of transformation products and parent compounds (SMX and TMP) become a necessary consideration during the degradation process of contaminants in aquatic environment.

3. Materials and Method

3.1. Mechanism Calculation

In this study, all of the quantum chemical calculations are implemented using the Gaussian 09 program [42]. The geometries of reactants, intermediates, transition states (TS), and products, as well as its vibrational frequencies are optimized with M06-2X functionals at a level 6-31 + G (d,p) basis set [43]. The M06-2X method, as a density functional theory (DFT) model, has been successfully applied to the mechanism studies of organic contaminants in the aquatic environment. [44,45,46,47,48,49]. The purpose of vibrational frequency calculation is to confirm the structure of transition states and local minimal point. Meanwhile, intrinsic reaction coordinate (IRC) calculations [49] are analyzed in order to verify the transition states that are connected the reactants with their corresponding products for the transformation pathway. Moreover, for a more accurate evaluation of energetic parameters, the single point energies of all of the geometries were carried out at the M06-2X/6-311++G (3df,3pd) level. Furthermore, the influence of solvent water is evaluated using the polarized continuum model (PCM) within a self-consistent reaction field (SCRF) theory [50,51,52]. The PCM has been proven to be flexible and accurate, in particularly, when the solute is accommodated in a cavity of realistic molecular shape and it has been widely used for the study of many chemical processes [53].
In addition, the Marcus theory [54] was used to calculate the mechanism of single-electron transfers (SET) pathway. The SET activation barrier ( Δ G S E T ) is defined relying on the reaction free energy ( Δ G S E T ) and the nuclear reorganization energy (λ):
Δ G S E T = ( λ + Δ G S E T ) 2 4 λ
The reorganization energy (λ) is determined by two energies, which has been calculated as:
λ = Δ E S E T Δ G S E T
where Δ E S E T is the energy difference between reactants and vertical products, which changes the spin multiplicity and charge at the same geometries.

3.2. Kinetics Computation

The transition state theory (TST) has been used to calculate the reaction kinetics, which has been successfully used to deal with radical reaction in former studies [55,56,57]. In this work, the rate constants of elementary reactions have been carried out with the effective method according to the activation energy barrier over a suitable temperature range. The rate constants are calculated by KiSThelP package that is based on TST with Wigner tunneling correction [58]. In order to obtain the rate constants of aqueous bimolecular reaction, the calculated rate constants need to be converted into corresponding aqueously values by dividing a coefficient (kc), kc is expressed by the following formular:
kc = RT/Pθ
where R, T, and Pθ are gas constant, temperature, and standard atmospheric pressure, respectively.

3.3. Ecotoxicity Evaluation

The ecological structure-activity relationships (ECOSAR) model, as a practical method for calculating toxicity, is employed to predict the ecotoxicity of SMX, TMP, and their transformation products [59,60]. The program can be used to estimate toxicity at the screening level of organic contaminants [61]. The eco-toxicity is evaluated by the acute toxicity (LC50 and EC50) and chronic toxicity (ChV) of SMX, TMP, and their transformation products for three aquatic organisms (green algae, daphnia, and fish).

4. Conclusions

In this paper, the transformation mechanisms, rate constants, and ecological risks for the •OH-initiated degradation process of SMX and TMP have been studied by quantum chemistry and computational toxicology methods. The lowest Gibbs free energy barrier of SMX (or TMP) initiated by •OH is 6.86 (or 6.21) kcal mol−1 through the Radd5 (or Radd1) pathway, which indicate that the most favorable pathways of the transformation are both •OH-addition reaction of benzene ring site. It is worth noting that the results shown that the hydroxylation of SMX at C3 site leading to the breaking the C-S bond and the opening of aromatic ring. This is an important finding that •OH radicals alone are able to cause the opening of aromatic ring. In addition, the rate constants of SMX and TMP were calculated for the assessment of the aquatic environmental fate. The overall rate constants of SMX and TMP reacted with •OH were 1.28 × 108 M−1 s−1 and 6.21 × 108 M−1 s−1 at 298 K, respectively. Additionally, the Radd1 and Radd5 pathways (ortho-position of aniline) of hydroxylation reactions of SMX play a dominant role (as high as 89% of branching ratio) in the photochemical transformation. Moreover, the hydroxylation products are formed through the subsequent reaction of •OH-complexes of SMX and TMP. It is noteworthy that although the aquatic toxicity decreases with the degradation of SMX and TMP, many products still maintain harmful levels to organisms, especially to daphnid (toxic or very toxic level).

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/21/17/6276/s1. Figure. S1: The pathways of •OH-initiated reaction of SMX in aquatic environment with the Gibbs free energy barriers ΔEb and reaction energies ΔEr (unit: kcal mol−1). Figure. S2: The pathways of •OH-initiated reaction of TMP in aquatic environment with the Gibbs free energy barriers ΔEb and reaction energies ΔEr (unit: kcal mol−1). Table S1: Calculated rate constants (M−1 s−1) between 273 and 328 K in the reaction of SMX and •OH. Table S2: Calculated rate constants (M−1 s−1) between 273 and 328 K in the reaction of TMP and •OH. Table S3: The acute and chronic toxic criterion (unit: mg L−1).

Author Contributions

Conceptualization, J.Y. and G.L.; methodology, J.Y.; software, Z.W.; validation, X.S.; formal analysis, J.Y.; investigation, G.L.; resources, X.S.; data curation, J.Y.; writing—original draft preparation, J.Y.; writing—review and editing, C.Z.; visualization, G.L.; supervision, X.S.; project administration, X.S.; funding acquisition, X.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (21976109), Natural Science Foundation of Shandong Province (ZR2018MB043), The Fundamental Research Funds of Shandong University (2018JC027), Shandong Province Key Research and Development Program (2019GSF109037).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Carballa, M.; Omil, F.; Lema, J.M.; Llompart, M.; García-Jares, C.; Rodríguez, I.; Gómez, M.; Ternes, T. Behavior of pharmaceuticals, cosmetics and hormones in a sewage treatment plant. Water Res. 2004, 38, 2918–2926. [Google Scholar] [CrossRef] [PubMed]
  2. Liu, J.-L.; Wong, M.-H. Pharmaceuticals and personal care products (PPCPs): A review on environmental contamination in China. Environ. Int. 2013, 59, 208–224. [Google Scholar] [CrossRef] [PubMed]
  3. Radjenović, J.; Petrović, M.; Ventura, F.; Barceló, D. Rejection of pharmaceuticals in nanofiltration and reverse osmosis membrane drinking water treatment. Water Res. 2008, 42, 3601–3610. [Google Scholar] [CrossRef] [PubMed]
  4. Zhang, X.-X.; Zhang, T.; Fang, H.H.P. Antibiotic resistance genes in water environment. Appl. Microbiol. Biotechnol. 2009, 82, 397–414. [Google Scholar] [CrossRef] [PubMed]
  5. Boyd, G.R.; Palmeri, J.M.; Zhang, S.; Grimm, D.A. Pharmaceuticals and personal care products (PPCPs) and endocrine disrupting chemicals (EDCs) in stormwater canals and Bayou St. John in New Orleans, Louisiana, USA. Sci. Total Environ. 2004, 333, 137–148. [Google Scholar] [CrossRef] [PubMed]
  6. Brogden, R.N.; Carmine, A.A.; Heel, R.C.; Speight, T.M.; Avery, G.S. Trimethoprim: A Review of its Antibacterial Activity, Pharmacokinetics and Therapeutic Use in Urinary Tract Infections. Drugs 1982, 23, 405–430. [Google Scholar] [CrossRef]
  7. Andrewes, D.A.; Chuter, P.J.; Dawson, M.J.; Eden, B.W.; Moore, R.M.A.; Freestone, D.S.; Morris, C.A. Trimethoprim and co-trimoxazole in the treatment of acute urinary tract infections: Patient compliance and efficacy. J. R. Coll. Gen. Pract. 1981, 31, 274. [Google Scholar]
  8. Murray, B.E.; Rensimer, E.R.; DuPont, H.L. Emergence of High-Level Trimethoprim Resistance in Fecal Escherichia coli during Oral Administration of Trimethoprim or Trimethoprim-Sulfamethoxazole. N. Engl. J. Med. 1982, 306, 130–135. [Google Scholar] [CrossRef]
  9. Raz, R.; Chazan, B.; Kennes, Y.; Colodner, R.; Rottensterich, E.; Dan, M.; Lavi, I.; Stamm, W.; Israeli Urinary Tract Infection Group. Empiric Use of Trimethoprim-Sulfamethoxazole (TMP-SMX) in the Treatment of Women with Uncomplicated Urinary Tract Infections, in a Geographical Area with a High Prevalence of TMP-SMX–Resistant Uropathogens. Clin. Infect. Dis. 2002, 34, 1165–1169. [Google Scholar]
  10. Boreen, A.L.; Arnold, W.A.; McNeill, K. Photodegradation of pharmaceuticals in the aquatic environment: A review. Aquat. Sci. 2003, 65, 320–341. [Google Scholar] [CrossRef]
  11. Abellán, M.N.; Giménez, J.; Esplugas, S. Photocatalytic degradation of antibiotics: The case of sulfamethoxazole and trimethoprim. Catal. Today 2009, 144, 131–136. [Google Scholar] [CrossRef]
  12. Martínez-Costa, J.I.; Rivera-Utrilla, J.; Leyva-Ramos, R.; Sánchez-Polo, M.; Velo-Gala, I. Individual and simultaneous degradation of antibiotics sulfamethoxazole and trimethoprim by UV and solar radiation in aqueous solution using bentonite and vermiculite as photocatalysts. Appl. Clay Sci. 2018, 160, 217–225. [Google Scholar] [CrossRef]
  13. Alharbi, S.K.; Kang, J.; Nghiem, L.D.; van de Merwe, J.P.; Leusch, F.D.L.; Price, W.E. Photolysis and UV/H2O2 of diclofenac, sulfamethoxazole, carbamazepine, and trimethoprim: Identification of their major degradation products by ESI–LC–MS and assessment of the toxicity of reaction mixtures. Process Saf. Environ. Protect. 2017, 112, 222–234. [Google Scholar] [CrossRef] [Green Version]
  14. Joss, A.; Zabczynski, S.; Göbel, A.; Hoffmann, B.; Löffler, D.; McArdell, C.S.; Ternes, T.A.; Thomsen, A.; Siegrist, H. Biological degradation of pharmaceuticals in municipal wastewater treatment: Proposing a classification scheme. Water Res. 2006, 40, 1686–1696. [Google Scholar] [CrossRef] [PubMed]
  15. Kolpin, D.W.; Furlong, E.T.; Meyer, M.T.; Thurman, E.M.; Zaugg, S.D.; Barber, L.B.; Buxton, H.T. Pharmaceuticals, Hormones, and Other Organic Wastewater Contaminants in U.S. Streams, 1999−2000:  A National Reconnaissance. Environ. Sci. Technol. 2002, 36, 1202–1211. [Google Scholar] [PubMed] [Green Version]
  16. Al Aukidy, M.; Verlicchi, P.; Jelic, A.; Petrovic, M.; Barcelò, D. Monitoring release of pharmaceutical compounds: Occurrence and environmental risk assessment of two WWTP effluents and their receiving bodies in the Po Valley, Italy. Sci. Total Environ. 2012, 438, 15–25. [Google Scholar] [CrossRef]
  17. Mahdi Ahmed, M.; Barbati, S.; Doumenq, P.; Chiron, S. Sulfate radical anion oxidation of diclofenac and sulfamethoxazole for water decontamination. Chem. Eng. J. 2012, 197, 440–447. [Google Scholar] [CrossRef]
  18. Boreen, A.L.; Arnold, W.A.; McNeill, K. Photochemical Fate of Sulfa Drugs in the Aquatic Environment:  Sulfa Drugs Containing Five-Membered Heterocyclic Groups. Environ. Sci. Technol. 2004, 38, 3933–3940. [Google Scholar] [CrossRef]
  19. Le-Minh, N.; Khan, S.J.; Drewes, J.E.; Stuetz, R.M. Fate of antibiotics during municipal water recycling treatment processes. Water Res. 2010, 44, 4295–4323. [Google Scholar] [CrossRef]
  20. Murata, A.; Takada, H.; Mutoh, K.; Hosoda, H.; Harada, A.; Nakada, N. Nationwide monitoring of selected antibiotics: Distribution and sources of sulfonamides, trimethoprim, and macrolides in Japanese rivers. Sci. Total Environ. 2011, 409, 5305–5312. [Google Scholar] [CrossRef]
  21. Watkinson, A.J.; Murby, E.J.; Kolpin, D.W.; Costanzo, S.D. The occurrence of antibiotics in an urban watershed: From wastewater to drinking water. Sci. Total Environ. 2009, 407, 2711–2723. [Google Scholar] [CrossRef] [PubMed]
  22. Li, B.; Zhang, T. Mass flows and removal of antibiotics in two municipal wastewater treatment plants. Chemosphere 2011, 83, 1284–1289. [Google Scholar] [CrossRef]
  23. Yang, X.; Flowers, R.C.; Weinberg, H.S.; Singer, P.C. Occurrence and removal of pharmaceuticals and personal care products (PPCPs) in an advanced wastewater reclamation plant. Water Res. 2011, 45, 5218–5228. [Google Scholar] [CrossRef]
  24. Andreozzi, R.; Caprio, V.; Ciniglia, C.; de Champdoré, M.; Lo Giudice, R.; Marotta, R.; Zuccato, E. Antibiotics in the Environment:  Occurrence in Italian STPs, Fate, and Preliminary Assessment on Algal Toxicity of Amoxicillin. Environ. Sci. Technol. 2004, 38, 6832–6838. [Google Scholar] [CrossRef] [PubMed]
  25. Brenner, C.G.B.; Mallmann, C.A.; Arsand, D.R.; Mayer, F.M.; Martins, A.F. Determination of Sulfamethoxazole and Trimethoprim and Their Metabolites in Hospital Effluent. Clean Soil Air Water 2011, 39, 28–34. [Google Scholar] [CrossRef]
  26. Eichhorn, P.; Ferguson, P.L.; Pérez, S.; Aga, D.S. Application of Ion Trap-MS with H/D Exchange and QqTOF-MS in the Identification of Microbial Degradates of Trimethoprim in Nitrifying Activated Sludge. Anal. Chem. 2005, 77, 4176–4184. [Google Scholar] [CrossRef]
  27. Yi, T.; Barr, W.; Harper, W.F., Jr. Electron density-based transformation of trimethoprim during biological wastewater treatment. Water Sci. Technol. 2012, 65, 689–696. [Google Scholar] [CrossRef]
  28. Göbel, A.; Thomsen, A.; McArdell, C.S.; Joss, A.; Giger, W. Occurrence and Sorption Behavior of Sulfonamides, Macrolides, and Trimethoprim in Activated Sludge Treatment. Environ. Sci. Technol. 2005, 39, 3981–3989. [Google Scholar] [CrossRef]
  29. Lindberg, R.H.; Olofsson, U.; Rendahl, P.; Johansson, M.I.; Tysklind, M.; Andersson, B.A.V. Behavior of Fluoroquinolones and Trimethoprim during Mechanical, Chemical, and Active Sludge Treatment of Sewage Water and Digestion of Sludge. Environ. Sci. Technol. 2006, 40, 1042–1048. [Google Scholar] [CrossRef]
  30. Batt, A.L.; Kim, S.; Aga, D.S. Enhanced Biodegradation of Iopromide and Trimethoprim in Nitrifying Activated Sludge. Environ. Sci. Technol. 2006, 40, 7367–7373. [Google Scholar] [CrossRef]
  31. Bonvin, F.; Omlin, J.; Rutler, R.; Schweizer, W.B.; Alaimo, P.J.; Strathmann, T.J.; McNeill, K.; Kohn, T. Direct Photolysis of Human Metabolites of the Antibiotic Sulfamethoxazole: Evidence for Abiotic Back-Transformation. Environ. Sci. Technol. 2013, 47, 6746–6755. [Google Scholar] [CrossRef] [PubMed]
  32. Buxton, G.V.; Greenstock, C.L.; Helman, W.P.; Ross, A.B. Critical Review of rate constants for reactions of hydrated electrons, hydrogen atoms and hydroxyl radicals (OH/O) in Aqueous Solution. J. Phys. Chem. Ref. Data 1988, 17, 513–886. [Google Scholar] [CrossRef] [Green Version]
  33. Baeza, C.; Knappe, D.R.U. Transformation kinetics of biochemically active compounds in low-pressure UV Photolysis and UV/H2O2 advanced oxidation processes. Water Res. 2011, 45, 4531–4543. [Google Scholar] [CrossRef] [PubMed]
  34. Wols, B.A.; Hofman-Caris, C.H.M. Review of photochemical reaction constants of organic micropollutants required for UV advanced oxidation processes in water. Water Res. 2012, 46, 2815–2827. [Google Scholar] [CrossRef] [PubMed]
  35. Yang, Y.; Lu, X.; Jiang, J.; Ma, J.; Liu, G.; Cao, Y.; Liu, W.; Li, J.; Pang, S.; Kong, X.; et al. Degradation of sulfamethoxazole by UV, UV/H2O2 and UV/persulfate (PDS): Formation of oxidation products and effect of bicarbonate. Water Res. 2017, 118, 196–207. [Google Scholar] [CrossRef]
  36. Zhang, R.; Sun, P.; Boyer, T.H.; Zhao, L.; Huang, C.-H. Degradation of Pharmaceuticals and Metabolite in Synthetic Human Urine by UV, UV/H2O2, and UV/PDS. Environ. Sci. Technol. 2015, 49, 3056–3066. [Google Scholar] [CrossRef]
  37. Luo, X.; Zheng, Z.; Greaves, J.; Cooper, W.J.; Song, W. Trimethoprim: Kinetic and mechanistic considerations in photochemical environmental fate and AOP treatment. Water Res. 2012, 46, 1327–1336. [Google Scholar] [CrossRef]
  38. Solar, S.; Solar, W.; Getoff, N. Resolved multisite OH-attack on aqueous aniline studied by pulse radiolysis. Int. J. Radiat. Appl. Instrum. Part C Radiat. Phys. Chem. 1986, 28, 229–234. [Google Scholar] [CrossRef]
  39. Jewell, K.S.; Castronovo, S.; Wick, A.; Falås, P.; Joss, A.; Ternes, T.A. New insights into the transformation of trimethoprim during biological wastewater treatment. Water Res. 2016, 88, 550–557. [Google Scholar] [CrossRef] [Green Version]
  40. Kim, Y.; Choi, K.; Jung, J.; Park, S.; Kim, P.-G.; Park, J. Aquatic toxicity of acetaminophen, carbamazepine, cimetidine, diltiazem and six major sulfonamides, and their potential ecological risks in Korea. Environ. Int. 2007, 33, 370–375. [Google Scholar] [CrossRef]
  41. Halling-Sørensen, B.; Lützhøft, H.C.H.; Andersen, H.R.; Ingerslev, F. Environmental risk assessment of antibiotics: Comparison of mecillinam, trimethoprim and ciprofloxacin. J. Antimicrob. Chemother. 2000, 46, 53–58. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Frisch, M.; Trucks, G.; Schlegel, H.; Scuseria, G.; Robb, M.; Cheeseman, J.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. Gaussian 09; Gaussian, Inc.: Wallingford, CT, USA, 2009; Volume 32, pp. 5648–5652. [Google Scholar]
  43. Walker, M.; Harvey, A.J.A.; Sen, A.; Dessent, C.E.H. Performance of M06, M06-2X, and M06-HF Density Functionals for Conformationally Flexible Anionic Clusters: M06 Functionals Perform Better than B3LYP for a Model System with Dispersion and Ionic Hydrogen-Bonding Interactions. J. Phys. Chem. A 2013, 117, 12590–12600. [Google Scholar] [CrossRef]
  44. Zhao, Y.; Truhlar, D.G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: Two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120, 215–241. [Google Scholar]
  45. Gao, Y.; Li, G.; Qin, Y.; Ji, Y.; Mai, B.; An, T. New theoretical insight into indirect photochemical transformation of fragrance nitro-musks: Mechanisms, eco-toxicity and health effects. Environ. Int. 2019, 129, 68–75. [Google Scholar] [CrossRef]
  46. Satheeshkumar, R.; Shankar, R.; Kaminsky, W.; Kalaiselvi, S.; Padma, V.V.; Prasad, K.J.R. Theoretical and experimental investigations on molecular structure of 7-Chloro-9-phenyl-2, 3-dihydroacridin-4(1H)-one with cytotoxic studies. J. Mol. Struct. 2016, 1109, 247–257. [Google Scholar] [CrossRef]
  47. Liu, W.; Lv, G.; Sun, X.; He, L.; Zhang, C.; Li, Z. Theoretical study on the reaction of anthracene with sulfate radical and hydroxyl radical in aqueous solution. Ecotox. Environ. Saf. 2019, 183, 109551. [Google Scholar] [CrossRef]
  48. Wang, L.; Tang, A. Oxidation mechanisms of dimethyl selenide and selenoxide in the atmosphere initiated by OH radical. A theoretical study. Chem. Phys. 2011, 382, 98–103. [Google Scholar] [CrossRef]
  49. Maeda, S.; Harabuchi, Y.; Ono, Y.; Taketsugu, T.; Morokuma, K. Intrinsic reaction coordinate: Calculation, bifurcation, and automated search. Int. J. Quantum Chem. 2015, 115, 258–269. [Google Scholar] [CrossRef]
  50. Zheng, S.; Geva, E.; Dunietz, B.D. Solvated Charge Transfer States of Functionalized Anthracene and Tetracyanoethylene Dimers: A Computational Study Based on a Range Separated Hybrid Functional and Charge Constrained Self-Consistent Field with Switching Gaussian Polarized Continuum Models. J. Chem. Theory Comput. 2013, 9, 1125–1131. [Google Scholar] [CrossRef]
  51. Fudickar, W.; Linker, T. Theoretical insights into the effect of solvents on the [4 + 2] cycloaddition of singlet oxygen to substituted anthracenes: A change from a stepwise process to a concerted process. J Phys. Org. Chem. 2019, 32, e3951. [Google Scholar] [CrossRef]
  52. Mennucci, B.; Tomasi, J.; Cammi, R.; Cheeseman, J.R.; Frisch, M.J.; Devlin, F.J.; Gabriel, S.; Stephens, P.J. Polarizable Continuum Model (PCM) Calculations of Solvent Effects on Optical Rotations of Chiral Molecules. J. Phys. Chem. A 2002, 106, 6102–6113. [Google Scholar] [CrossRef]
  53. Tomasi, J.; Cammi, R.; Mennucci, B. Medium effects on the properties of chemical systems: An overview of recent formulations in the polarizable continuum model (PCM). Int. J. Quantum Chem. 1999, 75, 783–803. [Google Scholar] [CrossRef]
  54. Mayer, J.M. Understanding Hydrogen Atom Transfer: From Bond Strengths to Marcus Theory. Acc. Chem. Res. 2011, 44, 36–46. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Vanden-Eijnden, E.; Tal, F.A. Transition state theory: Variational formulation, dynamical corrections, and error estimates. J. Chem. Phys. 2005, 123, 184103. [Google Scholar] [CrossRef] [PubMed]
  56. Sheppard, D.; Terrell, R.; Henkelman, G. Optimization methods for finding minimum energy paths. J. Chem. Phys. 2008, 128, 134106. [Google Scholar] [CrossRef] [Green Version]
  57. Georgievskii, Y.; Klippenstein, S.J. Transition State Theory for Multichannel Addition Reactions:  Multifaceted Dividing Surfaces. J. Phys. Chem. A 2003, 107, 9776–9781. [Google Scholar] [CrossRef]
  58. Canneaux, S.; Bohr, F.; Henon, E. KiSThelP: A program to predict thermodynamic properties and rate constants from quantum chemistry results. J. Comput. Chem. 2014, 35, 82–93. [Google Scholar] [CrossRef]
  59. Takata, M.; Lin, B.-L.; Xue, M.; Zushi, Y.; Terada, A.; Hosomi, M. Predicting the acute ecotoxicity of chemical substances by machine learning using graph theory. Chemosphere 2020, 238, 124604. [Google Scholar] [CrossRef]
  60. Claeys, L.; Iaccino, F.; Janssen, C.R.; Van Sprang, P.; Verdonck, F. Development and validation of a quantitative structure–activity relationship for chronic narcosis to fish. Environ. Toxicol. Chem. 2013, 32, 2217–2225. [Google Scholar] [CrossRef]
  61. Comber, M.H.I.; Walker, J.D.; Watts, C.; Hermens, J. Quantitative structure-activity relationships for predicting potential ecological hazard of organic chemicals for use in regulatory risk assessments. Environ. Toxicol. Chem. 2003, 22, 1822–1828. [Google Scholar] [CrossRef] [PubMed]
Figure 1. All the reaction pathways in the initial reaction of sulfamethoxazole (SMX) and trimethoprim (TMP) with •OH.
Figure 1. All the reaction pathways in the initial reaction of sulfamethoxazole (SMX) and trimethoprim (TMP) with •OH.
Ijms 21 06276 g001
Figure 2. Schematic diagram of free energy for the reactions of SMX with •OH at different sites at the M06-2X/6–31+G (d, p) level.
Figure 2. Schematic diagram of free energy for the reactions of SMX with •OH at different sites at the M06-2X/6–31+G (d, p) level.
Ijms 21 06276 g002
Figure 3. Schematic diagram of free energy for the reactions of TMP with •OH at different sites at the M06-2X/6–31+G (d, p) level.
Figure 3. Schematic diagram of free energy for the reactions of TMP with •OH at different sites at the M06-2X/6–31+G (d, p) level.
Ijms 21 06276 g003
Figure 4. Calculated branching ratio of each route in the initial reaction of SMX with •OH at 273–328 K.
Figure 4. Calculated branching ratio of each route in the initial reaction of SMX with •OH at 273–328 K.
Ijms 21 06276 g004
Figure 5. Schematic diagram of free energy for the subsequent reactions of IM5 during the transformation of SMX at the M06-2X/6–31 + G (d, p) level.
Figure 5. Schematic diagram of free energy for the subsequent reactions of IM5 during the transformation of SMX at the M06-2X/6–31 + G (d, p) level.
Ijms 21 06276 g005
Figure 6. Schematic diagram of free energy for the subsequent reactions of IM1 during the transformation of TMP at the M06-2X/6–31+G (d, p) level.
Figure 6. Schematic diagram of free energy for the subsequent reactions of IM1 during the transformation of TMP at the M06-2X/6–31+G (d, p) level.
Ijms 21 06276 g006
Figure 7. Acute toxicity (a) and chronic toxicity (b) (unit: mg L−1) of SMX and TMP and their transformation products.
Figure 7. Acute toxicity (a) and chronic toxicity (b) (unit: mg L−1) of SMX and TMP and their transformation products.
Ijms 21 06276 g007
Table 1. Calculated overall rate constants (M−1 s−1) between 273 and 328 K in the reaction of SMX and TMP with •OH.
Table 1. Calculated overall rate constants (M−1 s−1) between 273 and 328 K in the reaction of SMX and TMP with •OH.
T(K)SMXExperimentTMPExperiment
2731.58 × 108(5.8 ± 0.2) × 1099.72 × 1088.66 × 109
2881.39 × 1087.35 × 108
2981.28 × 1086.21 × 108
3131.15 × 1084.94 × 108
3281.05 × 1084.03 × 108

Share and Cite

MDPI and ACS Style

Yang, J.; Lv, G.; Zhang, C.; Wang, Z.; Sun, X. Indirect Photodegradation of Sulfamethoxazole and Trimethoprim by Hydroxyl Radicals in Aquatic Environment: Mechanisms, Transformation Products and Eco-Toxicity Evaluation. Int. J. Mol. Sci. 2020, 21, 6276. https://doi.org/10.3390/ijms21176276

AMA Style

Yang J, Lv G, Zhang C, Wang Z, Sun X. Indirect Photodegradation of Sulfamethoxazole and Trimethoprim by Hydroxyl Radicals in Aquatic Environment: Mechanisms, Transformation Products and Eco-Toxicity Evaluation. International Journal of Molecular Sciences. 2020; 21(17):6276. https://doi.org/10.3390/ijms21176276

Chicago/Turabian Style

Yang, Jiaoxue, Guochun Lv, Chenxi Zhang, Zehua Wang, and Xiaomin Sun. 2020. "Indirect Photodegradation of Sulfamethoxazole and Trimethoprim by Hydroxyl Radicals in Aquatic Environment: Mechanisms, Transformation Products and Eco-Toxicity Evaluation" International Journal of Molecular Sciences 21, no. 17: 6276. https://doi.org/10.3390/ijms21176276

APA Style

Yang, J., Lv, G., Zhang, C., Wang, Z., & Sun, X. (2020). Indirect Photodegradation of Sulfamethoxazole and Trimethoprim by Hydroxyl Radicals in Aquatic Environment: Mechanisms, Transformation Products and Eco-Toxicity Evaluation. International Journal of Molecular Sciences, 21(17), 6276. https://doi.org/10.3390/ijms21176276

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop