Assessment of long-range transport potential of polychlorinated Naphthalenes based on three-dimensional QSAR models

High-throughput prediction of oral acute toxicity in Rat and Mouse of over 100,000 polychlorinated persistent organic pollutants (PC-POPs) by interpretable data fusion-driven machine learning global models

This study utilized available oral acute toxicity data in Rat and Mouse for polychlorinated persistent organic pollutants (PC-POPs) to construct data fusion-driven machine learning (ML) global models. Based on atom-centered fragments (ACFs), the collected high-throughput data overcame the applicability limitations, enabling accurate toxicity prediction for a wide range of PC-POPs series compounds using only single models. The data variances in the Rat training and test sets were 1.52 and 1.34, respectively, while for the Mouse, the values were 1.48 and 1.36, respectively. Genetic algorithm (GA) was used to build multiple linear regression (MLR) models and pre-screen descriptors, addressing the "black-box" problem prevalent in ML and enhancing model interpretability. The best ML models for Rat and Mouse achieved approximately 90 % prediction reliability for over 100,000 true untested compounds. Ultimately, a warning list of highly toxic compounds for eight categories of polychlorinated atom-centered fragments (PCACFs) was generated based on the prediction results. The analysis of descriptors revealed that dioxin analogs generally exhibited higher toxicity, because the heteroatoms and ring systems increased structural complexity and formed larger conjugated systems, contributing to greater oral acute toxicity. The present study provides valuable insights for guiding the subsequent in vivo tests, environmental risk assessment and the improvement of global governance system of pollutants.

Unveiling the Photodegradation Mechanism of Monochlorinated Naphthalenes under UV-C Irradiation: Affecting Factors Analysis, the Roles of Hydroxyl Radicals, and DFT Calculation

Polychlorinated naphthalenes (PCNs) are a new type of persistent organic pollutant (POP) characterized by persistence, bioaccumulation, dioxin-like toxicity, and long-range atmospheric transport. Focusing on one type of PCN, monochlorinated naphthalenes (CN-1, CN-2), this study aimed to examine their photodegradation in the environment. In this work, CN-1 and CN-2 were employed as the model pollutants to investigate their photodegradation process under UV-C irradiation. Factors like the pH, initial concentrations of CN-1, and inorganic anions were investigated. Next, the roles of hydroxyl radicals (•OH), superoxide anion radicals (O2•-), and singlet oxygen (1O2) in the photodegradation process were discussed and proposed via theory computation. The results show that the photodegradation of CN-1 and CN-2 follows pseudo-first-order kinetics. Acidic conditions promote the photodegradation of CN-1, while the effects of pH on the photodegradation of CN-2 are not remarkable. Cl-, NO3-, and SO32- accelerate the photodegradation of CN-1, whereas the effect of SO42- and CO32- is not significant. Additionally, the contributions of •OH and O2•- to the photodegradation of CN-1 are 20.47% and 38.80%, while, for CN-2, the contribution is 16.40% and 16.80%, respectively. Moreover, the contribution of 1O2 is 15.7%. Based on DFT calculations, C4 and C6 of the CN-1 benzene ring are prioritized attack sites for •OH, while C2 and C9 of CN-2 are prioritized attack sites.

Occurrence and Trophic Transfer of Polychlorinated Naphthalenes (PCNs) in the Arctic and Antarctic Benthic Marine Food Webs

Information about the occurrence and trophic transfer of polychlorinated naphthalenes (PCNs) in polar ecosystems is vital but scarce. In this study, PCNs were analyzed in benthic marine sediment and several biological species, collected around the Chinese polar scientific research stations in Svalbard in the Arctic and South Shetland Island in Antarctica. Total PCNs in biota ranged from 28 to 249 pg/g of lipid weight (lw) and from 11 to 284 pg/g lw in the Arctic and Antarctic regions, respectively. The concentrations and toxic equivalent (TEQ) of PCNs in polar marine matrices remained relatively low, and the compositions were dominated by lower chlorinated homologues (mono- to trichlorinated naphthalenes). Trophic magnification factors (TMFs) were calculated for congeners, homologues, and total PCNs in the polar benthic marine food webs. Opposite PCN transfer patterns were observed in the Arctic and Antarctic regions, i.e., trophic dilution and trophic magnification, respectively. This is the first comprehensive study of PCN trophic transfer behaviors in remote Arctic and Antarctic marine regions, providing support for further investigations of the biological trophodynamics and ecological risks of PCNs.

Triazine Herbicides Risk Management Strategies on Environmental and Human Health Aspects Using In-Silico Methods

As an effective herbicide, 1, 3, 5-Triazine herbicides (S-THs) are used widely in the pesticide market. However, due to their chemical properties, S-THs severely threaten the environment and human health (e.g., human lung cytotoxicity). In this study, molecular docking, Analytic Hierarchy Process—Technique for Order Preference by Similarity to the Ideal Solution (AHP-TOPSIS), and a three-dimensional quantitative structure-active relationship (3D-QSAR) model were used to design S-TH substitutes with high herbicidal functionality, high microbial degradability, and low human lung cytotoxicity. We discovered a substitute, Derivative-5, with excellent overall performance. Furthermore, Taguchi orthogonal experiments, full factorial design of experiments, and the molecular dynamics method were used to identify three chemicals (namely, the coexistence of aspartic acid, alanine, and glycine) that could promote the degradation of S-THs in maize cropping fields. Finally, density functional theory (DFT), Estimation Programs Interface (EPI), pharmacokinetic, and toxicokinetic methods were used to further verify the high microbial degradability, favorable aquatic environment, and human health friendliness of Derivative 5. This study provided a new direction for further optimizations of novel pesticide chemicals.

Prediction of adsorption capacity and biodegradability of polybrominated diphenyl ethers in soil

Polybrominated diphenyl ethers (PBDEs) are widely used brominated flame retardants with strong toxicity concerns. Understanding the behaviors of PBDEs in soil is essential to evaluate their environmental impact. However, the limited, incoherent, and inaccurate data has challenged predicting the adsorption capacity and biodegradability of all 209 PBDE congeners in the soil. Moreover, there are minimal studies regarding the interactions between adsorption and biodegradation behaviors of PBDEs in the soil. Herein, in this study, we adopted quantitative structure-property relationship (QSAR) modeling to predict the adsorption behavior of 209 PBDE congeners by estimating their organic carbon-water partition coefficient (K_OC) values. In addition, the biodegradability of commonly occurring PBDE congeners was evaluated by analyzing their affinity to extracellular enzymes responsible for biodegradation using molecular docking. The results highlight that the degree of bromination plays a significant role in both the absorption and biodegradation of PBDEs in the soil due to compound stability and molecular geometry. Our findings help to advance the knowledge on PBDE behaviors in the soil and facilitate PBDE remediation associated with a soil environment.

Modelling the octanol-air partition coefficient of aromatic pollutants based on the solvation free energy and the dimer effect

In this study, generalized predictive models were developed to estimate K_OA of four kinds of aromatic pollutants based on the calculated solvation free energy and taking the dimer effect into account. Uncorrected log K_OA values, which were directly estimated from the calculated solvation free energy of individual molecules, underestimated experimental values, and the deviation increased with increasing log K_OA. Dimers were found to greatly affect the apparent K_OA values of these aromatic pollutants, which were driven by π-π interactions. London dispersion and exchange-repulsion terms were identified to be dominant components of the underlying π-π interactions. It is interesting to find that the π-π interactions of polybrominated diphenyl ethers correlate with not only the molecular polarizability but also the size of opposing aromatic surfaces, which leads to a different trend of π-π interactions from other aromatic pollutants. A universal quantitative structure-activity relationship model was developed to estimate the proportion of dimers based on five molecular structural descriptors relevant to the π-π interactions. After calibration with the dimer effect, estimations of log K_OA were consistent with experimental values. Therefore, the dimer effect should be taken into consideration when investigating the partition behavior of aromatic pollutants, and the solvation free energy model could be an alternative method for the prediction of K_OA.

Polychlorinated naphthalenes (PCNs) in soils and plants from Svalbard, Arctic: Levels, distribution, and potential sources

As persistent organic pollutants (POPs) newly banned by the Stockholm Convention, polychlorinated naphthalenes (PCNs) have been widely detected in various environmental matrices. To date, however, the occurrence of PCNs in soils and plants in the Arctic environment has not been reported. In the current study, the concentrations and distribution of PCNs in Arctic soils and plants from Svalbard were analyzed. Total PCN concentrations ranged from 5.3 to 2550 pg/g dry weight (dw) in soils and 77 to 870 pg/g dw in plants. The higher levels of PCNs near the research stations and Longyearbyen town highlighted the significant influence of local anthropogenic emission sources. The composition of PCNs in Arctic soils and plants was dominated by lower chlorinated homologues, especially mono- to trichlorinated naphthalenes, which accounted for over 80 % of total PCNs in the soil and plant samples. The correlation analysis indicated the potential influences of total organic carbon (TOC) content on PCN concentrations in the soil, and octanol-air partition coefficients (K_OA) or octanol-water partition coefficients (K_OW) on PCN accumulation from soils to plants. To the best of our knowledge, this is the first study to report on the concentration and distribution of PCNs in Arctic soils and plants.

Theoretical design and process control of neonicotinoids insecticides suitable for synergistic degradation with the rubisco enzyme from rhizobia and carbon-fixing bacteria in soil

In this study, we studied and developed the modification schemes of environmentally friendly substitutes of neonicotinoid insecticides (NNIs) along with the regulatory measures that effectively enhanced the synergistic degradation of NNIs by soil rhizobia and carbon-fixing bacteria. Firstly, the binding ability of NNIs to the two key proteins was characterized by molecular docking; secondly, the mean square deviation decision method, which is a comprehensive evaluation method, was used to investigate the binding ability of NNI molecules with the two Rubisco rate-limiting enzymes. The three-dimensional quantitative structure-activity relationship (3D-QSAR) model was established for the synergistic degradation and single effect of rhizobia and carbon-fixing bacteria. Finally, after combining the 3D-QSAR model with a contour map analysis of the synergistic degradation effect of soil rhizobia and carbon-fixing bacteria, 102 NNI derivatives were designed. Flonicamid-36 and other four NNI derivatives passed the functional and environmentally friendly evaluation. Taguchi orthogonal experiment and factorial experiment-assisted molecular dynamics method were used to simulate the effects of 32 regulation schemes on the synergistic degradation of NNIS and its derivatives by rhizobia and carbon fixing bacteria. The synergistic degradation capacity of soil rhizobia and carbon-fixing bacteria was increased to 33.32% after right nitrogen supplementation. This indicated that supplementing the correct amount of nitrogen in the soil environment was beneficial to the microbial degradation of NNIs and their derivatives.

3D-QSAR-aided toxicity assessment of synthetic musks and their transformation by-products

Synthetic musks (SMs) are fragrance additives widely used in personal care products. SMs and their transformation by-products may reach the environment even after wastewater treatment, resulting in ecological and health concerns. The identification and toxicity assessment of SM by-products generated from different chemical and biological treatment processes have been rarely studied. This study established a 3D-QSAR model based on SMs' molecular structures (independent variable) and their lethal concentration (LC₅₀) of mysid (dependent variable). The developed model was further used to predict the LC₅₀ of SMs transformation by-products. Fifty-eight by-products of six common SMs (i.e., galaxolide (HHCB), tonalide (AHTN), phantolide (PHAN), traseolide (TRASE), celestolide (ADBI), and musk ketone (MK)) generated from biodegradation, photodegradation, advanced oxidation, and chlorination were identified through literature review and lab experiment as the model inputs. Predicted LC₅₀ results indicated that the toxicity of 40% chlorination by-products is higher than their precursors. Biodegradation is an effective method to treat AHTN. The advanced oxidation may be the best way to treat HHCB. This is the first study on biotoxicity of SM transformation by-products predicted by the 3D-QSAR model. The research outputs helped to provide valuable reference data and guidance to improve management of SMs and other emerging contaminants.

Strategies to Control Human Health Risks Arising from Antibiotics in the Environment: Molecular Modification of QNs for Enhanced Plant-Microbial Synergistic Degradation

In the present work, a comprehensive screening and evaluation system was established to improve the plant-microbial synergistic degradation effects of QNs. The study included the construction of a 3D-QSAR model, the molecular modification, environmental friendliness and functional evaluation of drugs, degradation pathway simulation, and human health risk assessment. Molecular dynamics was applied to quantify the binding capacity of QNs toward the plant degradation enzyme (peroxidase) and microbial degradation enzymes (manganese peroxidase, lignin peroxidase, and laccase). The fuzzy comprehensive evaluation method was used in combination with the weighted average method for normalization and assigning equal weights to the plant and microbial degradation effect values of the QNs. Considering the synergistic degradation effect value as the dependent variable and the molecular information of the QNs as the independent variable, a 3D-QSAR model was constructed for the plant-microbial synergistic degradation effect of QNs. The constructed model was then employed to conduct the molecular modification, environmental friendliness and functional evaluation, degradation pathway simulation, and human health risk assessment of transformation products using pharmacokinetics and toxicokinetics. The results revealed that the synergistic degradation effect 3D-QSAR (CoMSIA) model exhibited good internal and external prediction ability, fitting ability, stability, and no overfitting phenomenon. Norfloxacin (NOR) was used as the target molecule in the molecular modification. A total of 35 NOR derivatives with enhanced plant-microbial synergistic degradation effect (1.32-21.51%) were designed by introducing small-volume, strongly electronegative, and hydrophobic hydrogen bond receptor groups into the active group of the norfloxacin structure. The environment-friendliness and the functionality of NOR were evaluated prior to and after the modification, which revealed seven environment-friendly FQs derivatives exhibiting moderate improvement in stability and bactericidal efficacy. The simulation of the NOR plant and microbial degradation pathways prior to and after the modification and the calculation of the reaction energy barrier revealed Pathway A (D-17 to D-17-2) and Pathway B (D-17 to D-17-4) as the most prone degradation pathways in plants and Pathway A (D-17 to D-17-1) and Pathway B (D-17 to D-17-4) as the most prone degradation pathways in microorganisms. This demonstrated that the degradation of the modified NOR derivatives was significantly enhanced, with the hydroxylation and piperazine ring substitution reaction playing an important role in the degradation process. Finally, the parameters, including hepatotoxicity, mutagenicity, and rodent carcinogenicity, among others, predicted using the pharmacokinetics and toxicokinetics analyses revealed a significant reduction in the human health risk associated with the modified NOR, along with a considerable reduction in the toxicity of its transformation products, implying that the human health risk associated with the transformation products was reduced remarkably. The present study provides a theoretical basis for novel ideas and evaluation programs for improving the plant-microbial synergistic degradation of the QNs antibiotics for source control and drug design, thereby reducing the residues of these antibiotics and the associated hazard in the complex plant-soil environment, ultimately decreasing the potential risks to human health.

An adjusted 3D-QSAR model for the combined activity of fluoroquinolones photodegradation and microbial degradation assisted by dynamic simulation and its application in molecular modification

This study developed a comprehensive characterization method for the combined degradation effect of modified fluoroquinolones (FQs) photodegradation and microbial degradation. A combination of revised 3D-QSAR model, molecular docking, path simulation inference, pharmacokinetics, molecular dynamics (MD) simulation and toxicokinetics simulation was used to construct a systematic environment-friendly drug screening system. Five derivatives were screened with significantly improved combined degradation effect (over 20%) and functional characteristics and human health parameters through combined model verification, functional and human health risk assessment. The simulation path of photo- and microbial-degradation of gatifloxacin and new gatifloxacin molecules was derived, and the reaction energy barrier was also calculated. The ratio of the total rate-determining steps change rate of the decreased energy barrier (14.10%:26.30%) was consistent with the ratio of the increased degradation performance predicted by the model (22.87%:19.77%), demonstrating the reliability of revised 3D-QSAR model and it could be applied in molecular modification. MD and toxicokinetics simulation were used to predict the binding energy and aquatic toxicity between photo- and microbial-degradation products and the degradation enzymes, which further to screen the degradation pathways with low potential environmental risks. The findings will be helpful to screen environment-friendly drug and develop appropriate strategies for its risk management.

Environmentally Friendly Quinolones Design for a Two-Way Choice between Biotoxicity and Genotoxicity through Double-Activity 3D-QSAR Model Coupled with the Variation Weighting Method

Quinolone (QN) antibiotics are widely used, which lead to their accumulation in soil and toxic effects on ryegrass in pasture. In this study, we employed ryegrass as the research object and selected the total scores of 29 QN molecules docked with two resistant enzyme structures, superoxide dismutase (SOD, PDB ID: 1B06) and proline (Pro, PPEP-2, PDB ID: 6FPC), as dependent variables. The structural parameters of QNs were used as independent variables to construct a QN double-activity 3D-QSAR model for determining the biotoxicity on ryegrass by employing the variation weighting method. This model was constructed to determine modification sites and groups for designing QNs molecules. According to the 3D contour map of the model, by considering enrofloxacin (ENR) and sparfloxacin (SPA) as examples, 23 QN derivatives with low biotoxicity were designed, respectively. The functional properties and environmental friendliness of the QN derivatives were predicted through a two-way selection between biotoxicity and genotoxicity before and after modification; four environmentally friendly derivatives with low biotoxicity and high genotoxicity were screened out. Mixed toxicity index and molecular dynamics methods were used to verify the combined toxicity mechanism of QNs on ryegrass before and after modification. By simulating the combined pollution of ENR and its derivatives in different soils (farmland, garden, and woodland), the types of combined toxicity were determined as partial additive and synergistic. Binding energies were calculated using molecular dynamics. The designed QN derivatives with low biotoxicity, high genotoxicity, and environmental friendliness can highly reduce the combined toxicity on ryegrass and can be used as theoretic reserves to replace QN antibiotics.

A Double-Activity (Green Algae Toxicity and Bacterial Genotoxicity) 3D-QSAR Model Based on the Comprehensive Index Method and Its Application in Fluoroquinolones' Modification

The comparative molecular similarity index analysis (CoMSIA) model of double-activity quinolones targeting green algae toxicity and bacterial genotoxicity (8:2) was constructed in this paper on the basis of the comprehensive index method. The contour maps of the model were analyzed for molecular modifications with high toxicities. In the CoMSIA model, the optimum number of components n was 7, the cross-validated q² value was 0.58 (>0.5), the standard deviation standard error of estimate (SEE) was 0.02 (<0.95), F was 1265.33, and the non-cross-validated R² value was 1 (>0.9), indicating that the model had a good fit and predicting ability. The scrambling stability test parameters Q², cross-validated standard error of prediction (cSDEP), and dq²/dr²yy were 0.54, 0.25, and 0.8 (<1.2), respectively, indicating that the model had good stability. The external verification coefficient r²_pred was 0.73 (>0.6), and standard error of prediction (SEP) was 0.17, indicating that the model had a good external prediction ability. The contribution rates of the steric fields, electrostatic fields, hydrophobic fields, hydrogen bond donor, and acceptor fields were 10.9%, 19.8%, 32.7%, 13.8%, and 22.8%, respectively. Large volume groups were selected for modification of ciprofloxacin (CIP), and the derivatives with increased double-activity characterization values were screened; the increase ratio ranged from 12.31-19.09%. The frequency of derivatives were positive and total energy, bioaccumulation, and environmental persistence was reduced, indicating that the CIP derivatives had good environmental stability and friendliness. Predicted values and CoMSIA model constructed of single activities showed that the CoMSIA model of double activities had accuracy and reliability. In addition, the total scores of the derivatives docking with the D1 protein, ferredoxin-NADP (H) reductases (FNRs), and DNA gyrase increased, indicating that derivatives can be toxic to green algae by affecting the photosynthesis of green algae. The mechanism behind the bactericidal effect was also explained from a molecular perspective.

Environmental Conversion Path Inference of New Designed Fluoroquinolones and Their Potential Environmental Risk

Fluoroquinolone (FQ) derivatives with environmental friendliness regarding photodegradation, bioconcentration, and genotoxicity were selected from our previous works so that their transformation pathways of biological metabolism, photodegradation, microbial degradation, and chlorination disinfection could be studied. The pathways of these molecules and their derivatives were simulated to investigate the genotoxicity of their transformation products. The results showed that the genotoxicity of the biological metabolites, photodegradation products, and microbial degradation products of the maternal FQ derivatives partially increased, whereas the disinfection by-products exhibited lower genotoxicity than their precursors. Some designed FQ molecular derivatives still had potential environmental risks in biological metabolism, photodegradation, and microbial degradation. This study demonstrated that it is necessary to take into account the potential environmental risks of the transformed products of the modified FQs molecules during biometabolism, photodegradation, microbial degradation, and chlorination processes when designing novel FQ molecules. In future studies, assessing the potential environmental risks during various artificial or natural processes can be applied to screen environmentally friendly novel FQ molecules to avoid and or reduce their threat to environmental and human health.

Fuzzy risk assessment of modified polychlorinated naphthalenes for enhanced degradation

The three-dimensional quantitative structure-activity relationship (3D-QSAR) model is established for polychlorinated naphthalenes (PCNs) using the biological degradability (total score) results to modify CN-56 to design 37 new derivatives with higher degradability (increased by 14.55-38.79%). Furthermore, five new CN-56 derivatives are selected through evaluation of their persistent organic pollutant properties (toxicity, bioconcentration, long-range transport) and practicability (stability, insulativity, flame retardancy) using 3D-QSAR, density functional theory (DFT) and molecular docking methods. Environmental and health-based risk assessments are conducted using the multimedia fugacity model and fuzzy theory for complete screening of the new CN-56 derivatives. Whereas CN-56 is classed as high risk, three new derivatives can be classed as medium risk. The biodegradability mechanism analysis of the PCNs indicates that the electrostatic property is the main factor that affects the degradability, which provides a favorable theoretical reference to obtain environmentally friendly fire retardant and insulating materials.

Environmentally friendly polychlorinated naphthalenes (PCNs) derivatives designed using 3D-QSAR and screened using molecular docking, density functional theory and health-based risk assessment

A complete design and screening system for environmental-friendly polychlorinated naphthalene (PCN) derivatives was established through three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking, density functional theory (DFT) methods and health-based risk assessment based on dynamic multimedia fugacity model. Two types of 3D-QSAR models were established for PCNs using the experimental biological toxicity (logEC₅₀) of 14 PCNs to carry out a modification to lower the logEC₅₀ of CN-70. Consequently, 67 new monosubstituted and disubstituted derivatives with a lower biological toxicity than CN-70 were designed. Furthermore, 21 new CN-70 derivatives were selected through the evaluation of their persistent organic pollutant properties (biological toxicity, bio-concentration, long-range transport potential, biodegradability) and practicability (stability, insulativity, flame retardancy) using 3D-QSAR, molecular docking and DFT methods. Finally, the non-carcinogenic and carcinogenic risks of 19 new CN-70 derivatives in different exposure pathways were reduced, and 5 derivatives with a significant decrease both in biological toxicity (amplitude reduction: 12.73%-32.51%) and risk (amplitude reduction: 32.18%-59.19%) were selected as environmental-friendly PCN derivatives, which had been screened using the health-based risk assessment system associated with dynamic multimedia fugacity model. This study provides a theoretical basis for the design of environmental-friendly flame retardants and insulating materials.

A theoretical method for the high-sensitivity fluorescence detection of PAEs through double-substitution modification

A highly sensitive method for the detection of phthalate acid ester (PAE) derivatives by fluorescence spectroscopy based on a double-substitution modification is reported. The fluorescence intensities of 18 PAE derivatives and the template PAEs after docking with the protein bovine serum albumin (BSA) are compared; the fluorescence intensities of the PAE derivatives increase significantly (28.20 to 110.00 times) with high sensitivities, but the functionalities (stability and insulation) are not affected. The persistence, mobility, and toxicity of the PAE derivatives are reduced by varying degrees, and their amounts of bioaccumulation are not significantly changed, indicating that the derivatives are more environmentally friendly. Within the wavelength range 240-420 nm in the fluorescence spectrum, there is no interference between the PAE derivatives and polycyclic aromatic hydrocarbons (PAHs), and the derivatives can be detected with good precision. Based on the analysis of fluorescence characteristics using a 2D quantitative structure-activity relationship (QSAR) model, the occupied orbital energies and Mulliken charge numbers are the main factors that influence the spectral characteristics of the PAE derivatives.

Modified neonicotinoid insecticide with bi-directional selective toxicity and drug resistance

A three-dimensional quantitative structure-activity relationship (3D-QSAR) model was established based on the molecular structures and the negative logarithm of experimental lethal concentration 50 values (pLC₅₀) of neonicotinoid insecticides. Then, the mechanisms of bi-directional selective toxic effects and drug resistance were determined using homology modeling and molecular docking analyses. The results of the model showed that the 1-, 2-, 4-, and 12- positions of neonicotinoid insecticides strongly affected their toxicity, and that the introduction of bulky or electropositive groups at these positions could increase the pLC₅₀ values. Using Compound 19 as a template, we designed 37 derivatives with greater toxicity (increased by 0.04-11.45%). Among them, 20 derivatives had bioconcentrations lower than that of Compound 19 (reduced by 0.38-147.88%). Further screening of Compound 19 and the 20 derivatives mentioned above by homology modeling and acetylcholine receptors (AChRs) molecular docking analyses showed that 10 derivatives had bi-directional selective toxic effects against pests and bees. Further docking analyses of Compound 19 and these 10 derivatives identified that Derivative-33 showed decreased docking with superoxide dismutase (SOD) and glutathione S transferase (GST) in pests and enhanced docking with these enzymes in bees, indicating bi-directional selective resistance for pests and bees. Accordingly, Derivative-33 was selected as a new insecticide with high toxicity to pests and low toxicity to bees (bi-directional selective toxicity), low resistance in pest populations, and high resistance in bee populations. This study provides valuable reference data and will be useful for the development of strategies to produce new environmentally friendly pesticides.

Synthetic Musks: A Class of Commercial Fragrance Additives in Personal Care Products (PCPs) Causing Concern as Emerging Contaminants

Synthetic musks (SMs) are promising fragrance additives used in personal care products (PCPs). The widespread presence of SMs in environmental media remains a serious risk because of their harmful effects. Recently, the environmental hazards of SMs have been widely reported in various environmental samples including those from coastal and marine regions. This paper provides a systematic review of SMs, including their classification, synthetic routes, analysis and occurrence in environmental samples, fate and toxicity in the environment, as well as the associated risk assessment and pollution control. Research gaps and future opportunities were also identified with the hope of raising interest in this topic.

Modification of Hexachlorobenzene to Molecules with Lower Long-Range Transport Potentials Using 3D-QSAR Models with a Full Factor Experimental Design

In this study, the hexachlorobenzene molecule was modified by three-dimensional quantitative structure-activity relationship (3D-QSAR) models and a full factor experimental design to obtain new hexachlorobenzene molecules with low migration ability. The 3D-QSAR models (comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA)) were constructed by SYBLY-X 2.0 software, using experimental data of octanol-air partition coefficients (K_OA) for 12 chlorobenzenes (CBs) congeners as the dependent variable, and the structural parameters of CBs as independent variables, respectively. A target molecule (hexachlorobenzene; HCB: its long-distance migration capability leads to pollution of the marine environment in Antarctic and Arctic) was modified using the 3D-QSAR contour maps associated with resolution V of the 2^10-3 full-factorial experimental design method, and 11 modified HCB molecules were produced with a single chlorine atom (-Cl₂) and three chlorine atoms (-Cl₁, -Cl₃, and -Cl₅) replaced with electropositive groups (-COOH, -CN, -CF₃, -COF, -NO₂, -F, -CHF₂, -ONO₂, and -SiF₃) to increase the logK_OA. The new molecules had essentially similar biological enrichment functions and toxicities as HCB but were found to be more easily degraded. A 2D-QSAR model and molecular docking technology indicated that both dipole moments and highest occupied orbital energies of the substituents markedly affected migration and degradation of the new molecules. The abilities of the compounds to undergo long distance migration were assessed. The modified HCB molecules (i.e. 2-CN-HCB, 2-CF₃-HCB, 1-F-3-COOH-5-NO₂-HCB, 1-NO₂-3-CN-5-CHF₂-HCB and 1-CN-3-F-5-NO₂-HCB) moved from a long-range transport potential of the modified molecules to a relatively low mobility class, and the transport potentials of the remaining modified HCB molecules (i.e. 2-COOH-HCB, 2-COF-HCB, 1-COF-3-ONO₂-5-NO2-HCB, 1-F-3-CN-5-SiF₃-HCB, 1-F-3-COOH-5-SiF₃-HCB and 1-CN-3-SiF₃-5-ONO₂-HCB) also significantly decreased. These results provide a basic theoretical basis for designing environmentally benign molecules based on HCB.

Relationship between the binding free energy and PCBs' migration, persistence, toxicity and bioaccumulation using a combination of the molecular docking method and 3D-QSAR

The molecular docking method was used to calculate the binding free energies between biphenyl dioxygenase and 209 polychlorinated biphenyl (PCB) congeners. The relationships between the calculated binding free energies and migration (octanol-air partition coefficients, KOA), persistence (half-life, t1/2), toxicity (half maximal inhibitory concentration, IC50), and bioaccumulation (bioconcentration factor, BCF) values for the PCBs were used to gain insight into the degradation of PCBs in the presence of biphenyl dioxygenase. The relationships between the calculated binding free energies and the molecular weights, KOA, BCF, and t1/2 values for the PCBs were statistically significant (P < 0.01), whereas the relationship between the calculated binding free energies and the IC50 for the PCBs was not statistically significant (P > 0.05). The electrostatic field, derived from three-dimensional quantitative structure-activity relationship studies, was a primary factor governing the binding free energy, which agreed with literature findings for KOA, t1/2, and BCF. Comparative molecular field analysis and comparative molecular similarity indices analysis contour maps showed that the binding free energies, KOA, t1/2, and BCF values for the PCBs decreased simultaneously when substituents with electropositive groups at the 3-position or electronegative groups at the 3'-position were introduced. This indicated the binding free energy was correlated with the persistent organic pollutant characteristics of PCBs. Furthermore, low binding free energies improved the degradation of the PCBs and simultaneously decreased the KOA, t1/2, and BCF values, thereby reducing the persistent organic pollutant characteristics of PCBs in the environment. These results are expected to be beneficial in providing a theoretical foundation for further elucidation of the degradation and molecular modification of PCBs.