Prediction of contaminant persistence in aqueous phase: a quantum chemical approach

Prediction of Base-Catalyzed Hydrolysis Kinetics of Polychlorinated Dibenzo-p-Dioxins by Density Functional Theory Calculations

Polychlorinated dibenzo-p-dioxins (PCDDs), comprising 75 congeners, have gained considerable attention from the general public and the scientific community owing to their high toxic potential. The base-catalyzed hydrolysis of PCDDs is crucial for the assessment of their environmental persistence. Nonetheless, owing to the substantial number of congeners and low hydrolysis rates of PCDDs, conducting hydrolysis experiments proves to be exceedingly time-consuming and financially burdensome. Herein, density functional theory and transition state theory were employed to predict the base-catalyzed hydrolysis of PCDDs in aquatic environments. Findings reveal that PCDDs undergo base-catalyzed hydrolysis in aquatic environments with two competing pathways: prevailing dioxin ring-opening and reduced reactivity in the hydrolytic dechlorination pathway. The resultant minor products include hydroxylated PCDDs, which exhibit thermodynamic stability surpassing that of the principal product, chlorinated hydroxydiphenyl ethers. The half-lives (ranging from 17.10 to 1.33 × 1010 h at pH = 8) associated with the base-catalyzed hydrolysis of PCDDs dissolved in water were shorter compared to those within the water-sediment environmental system. This observation implies that hydroxide ions can protect aquatic environments from PCDD contamination. Notably, this study represents the first attempt to predict the base-catalyzed hydrolysis of PCDDs by using quantum chemical methods.

Prediction model on hydrolysis kinetics of phthalate monoester: A density functional theory study

As primary degradation products of phthalate esters, phthalate monoesters (MPEs) have been widely detected in various aquatic environments and drawn growing toxicological concerns. Hydrolysis kinetics that is of importance for assessing environmental persistence of chemicals remain elusive for MPEs. Herein, kinetics of base-catalyzed and neutral hydrolysis for 18 MPEs with different leaving groups was investigated by density functional theory calculation. Results indicate that MPEs with leaving groups having pKa of <10 prefer dissociative transition states. MPEs are more persistent than their parents, and their hydrolysis half-lives were calculated to vary from 3.4 min to 79.2 years (pH = 7-9). A quantitative structure-activity relationship model was developed for predicting the hydrolysis kinetics parameters. It was found that pKa of the leaving groups and electronegativity of the MPEs are key factors determining the hydrolysis kinetics. This work may lay a theoretical foundation for better understanding the chemical process that governs MPE persistence in aquatic environments.

Base-catalyzed hydrolysis of spectinomycin in aqueous solutions: Kinetics and mechanisms

Hydrolysis plays an imperative role in the abiotic transformation process of antibiotics in aqueous solutions. However, little information is available on the hydrolysis process of spectinomycin (an aminocyclitol antibiotic). This study systematically investigated the spectinomycin hydrolysis kinetics and mechanisms under different pH via experiments and density functional theory (DFT) computation. Hydrolysis was first conducted in a pure water system under pH of 4.0-9.0 and temperature of 25 °C, 50 °C and 70 °C, respectively. Results showed that hydrolysis was highly dependent on pH and temperature. When pH > 6.0, spectinomycin hydrolysis was accelerated by the catalysis of OH^-. Meanwhile, the hydrolysis rate increased with the elevation of temperature. Then, for the reference of the practical environment, the general base-catalyzed hydrolysis and mechanisms were studied under environmental pH 6.0-8.0 and 25 °C. DFT calculation demonstrated that base-catalyzed hydrolysis of spectinomycin could be more thermodynamically and kinetically favorable based on the lower Gibbs free energies of reaction and Gibbs free energies of activation. Further, instead of specific base catalysis (OH^-), the general base catalysis (e.g., phosphate buffer) was also found to promote hydrolysis efficiency. The antibacterial activity and ecotoxicities of the hydrolysis product were analyzed to be lower than the precursor, thereby decreasing the environmental impact of spectinomycin.

Quantum chemical prediction of effects of temperature on hydrolysis rate of penicillin under weakly acidic condition

Temperature and pH are important factors affecting the hydrolysis of β-lactam antibiotics in water environments. However, the determination of hydrolysis kinetics and pathways is experimentally challenging, particularly in low temperature aqueous solutions because of time and cost constraints. In this study, an equation was employed to correct the Gibbs energy calculated in aqueous solutions by density functional theory methods to predict the effect of temperature on the hydrolysis kinetics and pathways of penicillin G. The results indicate that the most likely hydrolysis mechanism involves the opening of the β-lactam ring of anionic penicillin G protonated at the β-lactam oxygen atom with the participation of the carboxyl group and a water molecule. The results also suggest that the carboxyl group of β-lactam antibiotics was crucial for the hydrogen transfer. The predicted rate constants were of the same order of magnitude as the experimental values obtained under comparable pH and temperature conditions. Therefore, the quantum chemical methodology described herein can be potentially employed to determine pH- and temperature-based two-dimensional hydrolysis rate models, which can enable the prediction of the β-lactam antibiotics persistence in frigid waters.

Temporal and Spatial Groundwater Contamination Assessment Using Geophysical and Hydrochemical Methods: The Industrial Chemical Complex of Estarreja (Portugal) Case Study

With more than a half-century in operation, the industrial chemical complex of Estarreja (ICCE) in northern Portugal has left serious environmental liabilities in the region. Although protective measures were implemented, soils, surface, and groundwater contamination caused by persistent pollutants are still prevalent. This study presents data from several geophysical and hydrochemical campaigns carried out to monitor groundwater contamination in the Estarreja region over a period of 30 years. Both geophysical and hydrochemical data showed a good agreement and revealed an important anomaly caused by groundwater contamination (high levels of Na, Cl, SO4, and Fe, among others) in 2006–2007, likely caused by the remobilization of waste pollutants (roasted pyrites, soils, and sludge) during their deposition in a sealed landfill (operating between 2003 and 2005). More recently, in 2016, this impact persists, but was more attenuated and showed a general migration pattern from E to SW according to one of the main groundwater flow paths. Groundwater flow in this region has a local radial behaviour. Drainage effluent systems, such as ditches and buried pipes formerly used by ICCE, are also likely to contribute to some contamination “hotspots”. Finally, the results obtained by the combined use of these two approaches allowed for the delineation of the contamination plume for future monitoring.

Prediction Models on pKa and Base-Catalyzed Hydrolysis Kinetics of Parabens: Experimental and Quantum Chemical Studies

Parabens for which the molecules contain hydrolytic and ionizable groups, are emerging pollutants due to their ubiquity in the environment. However, lack of pK_a and second-order base-catalyzed hydrolysis kinetics (k_B) values limits their environmental persistence assessment. Herein, six parabens were selected as reference compounds for which the pK_a and k_B values were measured experimentally. A semiempirical quantum chemical (QC) method was selected to calculate pK_a of the parabens, and density functional theory (DFT) methods were selected to calculate k_B for neutral and anionic forms of the parabens, by comparing the QC-calculated and determined values. Combining the QC-calculated and experimental pK_a and k_B values, quantitative structure-activity relationships with determination coefficients (R²) being 0.947 and 0.842 for the pK_a and k_B models, respectively, were developed, which were validated and could be employed to efficiently fill the k_B and pK_a data gaps of parabens within applicability domains. The base-catalyzed hydrolysis half-lives were estimated to range from 6 h to 1.52 × 10⁶ years (pH 7-9, 25 °C), further necessitating the in silico models due to the tedious and onerous experimental determination, and the huge number of hydrolyzable and ionizable chemicals that may be released into the environment.

Development of Prediction Models on Base-Catalyzed Hydrolysis Kinetics of Phthalate Esters with Density Functional Theory Calculation

Many phthalate esters (PAEs) are chemicals of high production volume and of toxicological concern. The second-order rate constant for base-catalyzed hydrolysis ( k_B) is a key parameter for assessing environmental persistence of PAEs. However, the k_B values for most PAEs are lacking, and the experimental determination of k_B encounters various difficulties. Herein, density functional theory (DFT) methods were selected by comparing empirical k_B values of five PAEs and five carboxylic acid esters with the DFT-calculated ones. Results indicate that PAEs with cyclic side chains are more vulnerable to base-catalyzed hydrolysis than PAEs with linear alkyl side chains, followed by PAEs with branched alkyl side chains. By combining experimental and DFT-calculated second-order rate constants for base-catalyzed hydrolysis of one side chain in PAEs ( k_{B_side chain}), quantitative structure-activity relationship models were developed. The models can differentiate PAEs with the departure of the leaving group (or the nucleophilic attack of OH^-) as the rate-determining step in the hydrolysis and estimate k_B values, which provides a promising way to predict hydrolysis kinetics of PAEs. The half-lives of the investigated PAEs were calculated and vary from 0.001 h to 558 years (pH = 7∼9), further illustrating the necessity of prediction models for hydrolysis kinetics in assessing the environmental persistence of chemicals.

The Mechanism of C-H Bond Oxidation by Aqueous Permanganate

The permanganate ion (MnO₄^-) has been widely used as a reagent for water treatment for over a century. It is a strong enough oxidant to activate carbon-hydrogen bonds, one of the most important reactions in biological and chemical systems. Our current textbook understanding of the oxidation mechanism in aqueous solution involves an initial, rate-limiting hydride abstraction by permanganate followed by reaction of the carbocation with bulk water to form an alcohol. This mechanism fits well into the classic oxidation sequence of alkane → alcohol → aldehyde → carboxylate, the central paradigm for both abiotic and biotic alkane oxidation in aqueous environments. In this study, we provide three lines of evidence through (1) a broken-symmetry density functional theory approach, (2) isotope labeling experiments, and (3) kinetic network modeling to demonstrate that aqueous permanganate can circumvent prior alcohol formation and produce aldehydes directly via a reaction path that bifurcates after the initial transition state. In contrast to classic transition state theory, the rate-limiting step is found to not determine product distribution, bearing critical implications for pathway and rate predictions. This complex reaction network provides new insights into the oxidation mechanisms of organic compounds involving transition metal complexes as well as enzyme or metal oxide surface active sites.

QSPR modeling of the logKow and logKoc of polymethoxylated, polyhydroxylated diphenyl ethers and methoxylated-, hydroxylated-polychlorinated diphenyl ethers

In the present study, the structural parameters of 209 types of polymethoxylated diphenyl ethers (PMeODEs), 209 types of polyhydroxylated diphenyl ethers (PHODEs), seven types of methoxylated-polychlorinated diphenyl ethers (MeO-PCDEs) and seven types of hydroxylated-polychlorinated diphenyl ethers (HO-PCDEs) were calculated using the Gaussian 09 program at the B3LYP/6-311G** level. Using structural and positional parameters as descriptors, quantitative structure-property relationships (QSPR) models for the prediction of n-octanol/water partition coefficient (logK_ow) and soil sorption coefficient normalized to organic carbon (logK_oc) were established and verified. The position parameters N₂₍₆₎, N₃₍₅₎ and N₄ were the main positional factors influencing logK_ow and logK_oc of PMeODEs and PHODEs. The molecular polarizability α was entered into the QSPR models of the logK_ow and logK_oc of PMeODEs, PHODEs and MeO/HO-PCDEs, indicating that the molecular volume could influence the two environment-related properties of DEs significantly. All of the established QSPR models showed good goodness-of-fit, robustness, and predictive ability. The two models for all of the tested DEs are slightly inferior compared with the models for only a class of compounds. In addition, application domain analysis indicated that the models reliably predicted the logK_ow and logK_oc of the mon- to hexa-DEs.

Effect analysis of quantum chemical descriptors and substituent characteristics on Henry's law constants of polybrominated diphenyl ethers at different temperatures

Twelve substituent descriptors, 17 quantum chemical descriptors and 1/T were selected to establish a quantitative structure-property relationship (QSPR) model of Henry's law constants for 7 polybrominated diphenyl ethers (PBDEs) at five different temperatures. Then, the lgH of 202 congeners at different temperatures were predicted. The variation rule and regulating mechanism of lgH was studied from the perspectives of both quantum chemical descriptors and substituent characteristics. The R² for modeling and testing sets of the final QSPR model are 0.977 and 0.979, respectively, thus indicating good fitness and predictive ability for Henry' law constants of PBDEs at different temperatures. The favorable hydrogen binding sites are the 5,5',6,6'-positions for high substituent congeners and the O atom of the ether bond for low substituent congeners, which affects the interaction between PBDEs and water molecules. lgH is negatively and linearly correlated with 1/T, and the variation trends of lgH with temperature are primarily regulated by individual substituent characteristics, wherein: the more substituents involved, the smaller the lgH. The significant sequence for the main effect of substituent positions is para>meta>ortho, where the ortho-positions are mainly involved in second-order interaction effect (64.01%). Having two substituents in the same ring also provides a significant effect, with 81.36% of second-order interaction effects, particularly where there is an adjacent distribution (55.02%).

In silico environmental chemical science: properties and processes from statistical and computational modelling

Quantitative structure-activity relationships (QSARs) have long been used in the environmental sciences. More recently, molecular modeling and chemoinformatic methods have become widespread. These methods have the potential to expand and accelerate advances in environmental chemistry because they complement observational and experimental data with "in silico" results and analysis. The opportunities and challenges that arise at the intersection between statistical and theoretical in silico methods are most apparent in the context of properties that determine the environmental fate and effects of chemical contaminants (degradation rate constants, partition coefficients, toxicities, etc.). The main example of this is the calibration of QSARs using descriptor variable data calculated from molecular modeling, which can make QSARs more useful for predicting property data that are unavailable, but also can make them more powerful tools for diagnosis of fate determining pathways and mechanisms. Emerging opportunities for "in silico environmental chemical science" are to move beyond the calculation of specific chemical properties using statistical models and toward more fully in silico models, prediction of transformation pathways and products, incorporation of environmental factors into model predictions, integration of databases and predictive models into more comprehensive and efficient tools for exposure assessment, and extending the applicability of all the above from chemicals to biologicals and materials.

A quantum chemical study of HOCl-induced transformations of carbamazepine

The computational chemistry approach in predicting products and recalcitrans in hypochlorous acid promoted carbamazepine degradation in the environment.

Prediction of hydrolysis pathways and kinetics for antibiotics under environmental pH conditions: a quantum chemical study on cephradine

Understanding hydrolysis pathways and kinetics of many antibiotics that have multiple hydrolyzable functional groups is important for their fate assessment. However, experimental determination of hydrolysis encounters difficulties due to time and cost restraint. We employed the density functional theory and transition state theory to predict the hydrolysis pathways and kinetics of cephradine, a model of cephalosporin with two hydrolyzable groups, two ionization states, two isomers and two nucleophilic attack directions. Results showed that the hydrolysis of cephradine at pH = 8.0 proceeds via opening of the β-lactam ring followed by intramolecular amidation. The predicted rate constants at different pH conditions are of the same order of magnitude as the experimental values, and the predicted products are confirmed by experiment. This study identified a catalytic role of the carboxyl group in the hydrolysis, and implies that the carboxyl group also plays a catalytic role in the hydrolysis of other cephalosporin and penicillin antibiotics. This is a first attempt to quantum chemically predict hydrolysis of an antibiotic with complex pathways, and indicates that to predict hydrolysis products under the environmental pH conditions, the variation of the rate constants for different pathways with pH should be evaluated.

Synthesis and QSPR study on environment-related properties of polychlorinated diphenyl sulfides (PCDPSs)

Twenty-nine congeners of polychlorinated diphenyl sulfide (PCDPS) theoretically possible were synthesized and purified, and their octanol-water partition coefficients (logK(ow)) and high performance liquid chromatography capacity factors (HPLC-logk') were determined. Then the constitutional and structural descriptors were obtained respectively to establish 2D models and comparative molecular similarity indices analysis (CoMSIA) method was used to establish 3D models. All the models were robust and predictive, indicating that the electrostatic effect may be the primary factor influencing the two properties of PCDPS. The logK(ow) and logk' values of all PCDPS congeners were predicted based on these models. The establishment of linear equation on logK(ow) between PCDPS and chlorinated or brominated diphenyl ether (PCDE, PBDE) was attempted, and the hydrophobic differences were considered to relate with the size of the molecules.