Prediction of octanol-air partition coefficients for polychlorinated biphenyls (PCBs) using 3D-QSAR models

The q-RASPR approach for predicting the property and fate of persistent organic pollutants

This study presents a quantitative read-across structure-property relationship (q-RASPR) approach that integrates the chemical similarity information used in read-across with traditional quantitative structure-property relationship (QSPR) models. This novel framework is applied to predict the physicochemical properties and environmental behaviors of persistent organic pollutants, specifically polychlorinated biphenyls (PCBs) and polybrominated diphenyl ethers (PBDEs). By utilizing a curated dataset and incorporating similarity-based descriptors, the q-RASPR approach improves the accuracy of predictions, particularly for compounds with limited experimental data. The models' performances were assessed using internal cross-validation and external testing, demonstrating significant enhancements in predictive reliability compared to conventional QSPR models. The findings highlight the potential of q-RASPR for use in regulatory risk assessments and optimizing remediation strategies by providing more precise insights into the environmental fate of these contaminants.

Red Emission Carbon Nanoparticles Which Can Simultaneously Responding to Hypochlorite and pH

Multi-targets detection has obtained much attention because this sensing mode can realize the detection of multi-targets simultaneously, which is helpful for biomedical analysis. Carbon nanoparticles have attracted extensive attention due to their superior optical and chemical properties, but there are few reports about red emission carbon nanoparticles for simultaneous detection of multi-targets. In this paper, a red emission fluorescent carbon nanoparticles were prepared by 1, 2, 4-triaminobenzene dihydrochloride at room temperature. The as-prepared red emission fluorescent carbon nanoparticles exhibited strong emission peak located at 635 nm with an absolute quantum yield up to 24%. They showed excellent solubility, high photostability and good biocompatibility. Furthermore, it could sensitively and selectively response to hypochlorite and pH, thus simultaneous detection of hypochlorite and pH was achieved by combining the red emission fluorescent carbon nanoparticles with computational chemistry. The formation mechanisms of red emission fluorescent carbon nanoparticles and their response to hypochlorite and pH were investigated, respectively.

Multi-process regulation of novel brominated flame retardants: Environmentally friendly substitute design, screening and environmental risk regulation

Novel brominated flame retardants (NBFRs), one of the most widely used synthetic flame-retardant materials, have been considered as a new group of pollutants that potentially affect human health. To overcome the adverse effects of NBFRs, a systematic approach for molecular design, screening, and performance evaluation was developed to generate environmentally friendly NBFR derivatives with unaltered functionality. In the present study, the features of NBFRs (long-distance migration, biotoxicity, bioenrichment, and environmental persistence) were determined and characterized by the multifactor comprehensive characterization method with equal weight addition, and the similarity index analysis (CoMSIA) model was constructed. Based on the three-dimensional equipotential diagram of the target molecule 2-ethylhexyl tetrabromobenzoic acid (TBB), 23 TBB derivatives were designed. Of these, 22 derivatives with decreased environmental impact and unaltered functional properties (i.e., flame retardancy and stability) were selected using 3D-QSAR models and density functional theory methods. The health risks of these derivatives to humans were assessed by toxicokinetic analysis; the results narrowed down the number of candidates to three (Derivative-7, Derivative-10, and Derivative-15). The environmental impact of these candidates was further evaluated and regulated in the real-world environment by using molecular dynamics simulation assisted by the Taguchi experimental design method. The relationship between the binding effects and the nonbonding interaction resultant force (TBB derivatives-receptor proteins) was also studied, and it was found that the larger the modulus of the binding force, the stronger the binding ability of the two. This finding indicated that the environmental impact of the designed NBFR derivatives was decreased. The present study aimed to provide a new idea and method for designing NBFR substitutes and to provide theoretical support for restraining the potential environmental risks of NBFRs.

Human Endocrine-Disrupting Effects of Phthalate Esters through Adverse Outcome Pathways: A Comprehensive Mechanism Analysis

Phthalate esters (PAEs) are widely exposed in the environment as plasticizers in plastics, and they have been found to cause significant environmental and health hazards, especially in terms of endocrine disruption in humans. In order to investigate the processes underlying the endocrine disruption effects of PAEs, three machine learning techniques were used in this study to build an adverse outcome pathway (AOP) for those effects on people. According to the results of the three machine learning techniques, the random forest and XGBoost models performed well in terms of prediction. Subsequently, sensitivity analysis was conducted to identify the initial events, key events, and key features influencing the endocrine disruption effects of PAEs on humans. Key features, such as Mol.Wt, Q+, QH+, ELUMO, minHCsats, MEDC-33, and EG, were found to be closely related to the molecular structure. Therefore, a 3D-QSAR model for PAEs was constructed, and, based on the three-dimensional potential energy surface information, it was discovered that the hydrophobic, steric, and electrostatic fields of PAEs significantly influence their endocrine disruption effects on humans. Lastly, an analysis of the contributions of amino acid residues and binding energy (BE) was performed, identifying and confirming that hydrogen bonding, hydrophobic interactions, and van der Waals forces are important factors affecting the AOP of PAEs' molecular endocrine disruption effects. This study defined and constructed a comprehensive AOP for the endocrine disruption effects of PAEs on humans and developed a method based on theoretical simulation to characterize the AOP, providing theoretical guidance for studying the mechanisms of toxicity caused by other pollutants.

Multi-Dimensional Elimination of β-Lactams in the Rural Wetland: Molecule Design and Screening for More Antibacterial and Degradable Substitutes

Restricted economic conditions and limited sewage treatment facilities in rural areas lead to the discharge of small-scale breeding wastewater containing higher values of residual beta-lactam antibiotics (β-lactams), which seriously threatens the aquatic environment. In this paper, molecular docking and a comprehensive method were performed to quantify and fit the source modification for the combined biodegradation of β-lactams. Using penicillin (PNC) as the target molecule, combined with contour maps for substitute modification, a three-dimensional quantitative structure-activity relationship (3D-QSAR) model was constructed for the high-performance combined biodegradation of β-lactams. The selected candidate with better environmental friendliness, functionality, and high performance was screened. By using the homology modeling algorithms, the mutant penicillin-binding proteins (PBPs) of Escherichia coli were constructed to have antibacterial resistance against β-lactams. The molecular docking was applied to obtain the target substitute by analyzing the degree of antibacterial resistance of β-lactam substitute. The combined biodegradation of β-lactams and substitute in the constructed wetland (CW) by different wetland plant root secretions was studied using molecular dynamics simulations. The result showed a 49.28% higher biodegradation of the substitutes than PNC when the combined wetland plant species of Eichhornia crassipes, Phragmites australis, and Canna indica L. were employed.

Changes and distribution of gas-phase polycyclic aromatic hydrocarbons and organochlorine compounds in a high-mountain gradient over a three-year period (Pyrenees, 2017-2020)

The atmospheric gas-phase concentrations of several polychlorinated biphenyls (PCBs), polycyclic aromatic hydrocarbons (PAHs), hexachlorobenzene (HCB), and pentachlorobenzene (PeCB) were measured in six high-mountain sites in the Pyrenees (1619-2453 m). Polyurethane foam passive air samplers were used for this purpose, providing continuous records spanning over three years (2017-2020). The mean concentrations of ∑PCBs, HCB, and PeCB, 13 ± 4 pg m^-3, 44 ± 18 pg m^-3, and 23 ± 20 pg m^-3, respectively, were of the order of those reported in other mountain sites and similar to those measured 20 years ago in the same area, evidencing the persistence of these compounds despite the international regulatory actions. The mean concentration of ∑PAHs was 631 ± 238 pg m^-3, representing between two- and three-times lower values than 20 years ago in the same area, but still in the range of other mountain regions. Statistically significant increases in gas-phase concentrations at higher temperatures were observed for most compounds. The experimental phase-change pseudo-enthalpies calculated from the slopes of the regressions between the natural logarithm of the concentrations and the reciprocal of temperature were lower than the reference values for nearly all compounds. This difference suggested a main contribution of long-range atmospheric transport of the gas-phase PAH and organochlorine concentrations in this mountain area. However, the less volatile compounds such as benz[a]anthracene, PCB138, and PCB180 showed a closer similarity between experimental and laboratory enthalpies, indicating that a significant portion of the variations in concentration of these compounds originated from temperature-dependent diffusive exchange by re-volatilization from local surfaces. The concentrations found in these sentinel ecosystems demonstrate that long-range transport of organic pollutants remains a risk in remote continental environments.

Field comparison of passive polyurethane foam and active air sampling techniques for analysis of gas-phase semi-volatile organic compounds at a remote high-mountain site

Polyurethane foam passive air samplers (PUF-PAS) are good candidates for the determination of gas-phase semi-volatile organic compound (SVOC) air concentrations in high-mountain areas over long periods because they do not require an energy supply. However, the harsh meteorological conditions present in such locations can increase the uncertainties inherently associated to PAS sampling rates due to the many variables involved in their calculation and to the assumptions made regarding PUF diffusive uptake mechanics, which can considerably bias the resulting concentrations. Therefore, we studied the performance of PUF-PASs in a remote location in the Pyrenees mountain range for the analysis of several SVOCs in air, including polychlorobiphenyls (PCBs), hexachlorobenzene, pentachlorobenzene, polycyclic aromatic hydrocarbons (PAHs), and the less studied emerging organophosphate flame retardants (OPFRs). An in-situ PUF-PAS calibration using Performance Reference Compounds (PRCs) provided compound- and sampler-specific sampling rates, showing mean experimental errors (12%) that adequately conformed to an estimate of their expanded theoretical uncertainties (15%). This showcases the suitability of this calibration strategy in an area with conditions beyond those typically considered in calibration efforts available to date. Moreover, gas-phase concentrations of the studied pollutants from PUF-PAS samples showed very good agreement (R² up to 0.91, p < 0.01) when compared to those obtained using a conventional high-volume active air sampler (PUF-AAS), with some minor deviations observed for PAHs caused by the seasonality in their atmospheric concentrations. No relevant levels of pollutants preferentially bound to the particle phase were detected in the PUF-PASs, the particle infiltration efficiency of the sampler configuration used was found to be low, and compounds typically distributed between the gas and particle phases of AAS samples revealed profiles consistent with their vapor pressures, except for some OPFRs.

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.

Reliable Prediction of the Octanol-Air Partition Ratio

The octanol-air equilibrium partition ratio (K_OA ) is frequently used to describe the volatility of organic chemicals, whereby n-octanol serves as a substitute for a variety of organic phases ranging from organic matter in atmospheric particles and soils, to biological tissues such as plant foliage, fat, blood, and milk, and to polymeric sorbents. Because measured K_OA values exist for just over 500 compounds, most of which are nonpolar halogenated aromatics, there is a need for tools that can reliably predict this parameter for a wide range of organic molecules, ideally at different temperatures. The ability of five techniques, specifically polyparameter linear free energy relationships (ppLFERs) with either experimental or predicted solute descriptors, EPISuite's KOAWIN, COSMOtherm, and OPERA, to predict the K_OA of organic substances, either at 25 °C or at any temperature, was assessed by comparison with all K_OA values measured to date. In addition, three different ppLFER equations for K_OA were evaluated, and a new modified equation is proposed. A technique's performance was quantified with the mean absolute error (MAE), the root mean square error (RMSE), and the estimated uncertainty of future predicted values, that is, the prediction interval. We also considered each model's applicability domain and accessibility. With an RMSE of 0.37 and a MAE of 0.23 for predictions of log K_OA at 25 °C and RMSE of 0.32 and MAE of 0.21 for predictions made at any temperature, the ppLFER equation using experimental solute descriptors predicted the K_OA the best. Even if solute descriptors must be predicted in the absence of experimental values, ppLFERs are the preferred method, also because they are easy to use and freely available. Environ Toxicol Chem 2021;40:3166-3180. © 2021 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Potential Environmental Risk Characteristics of PCB Transformation Products in the Environmental Medium

The complementary construction of polychlorinated biphenyl (PCB) phytotoxicity and the biotoxicity 3D-QSAR model, combined with the constructed PCB environmental risk characterization model, was carried out to evaluate the persistent organic pollutant (POP) properties (toxicity (phytotoxicity and biotoxicity), bioconcentration, migration, and persistence) of PCBs and their corresponding transformation products (phytodegradation, microbial degradation, biometabolism, and photodegradation). The transformation path with a significant increase in environmental risks was analyzed. Some environmentally friendly PCB derivatives, exhibiting a good modification effect, and their parent molecules were selected as precursor molecules. Their transformation processes were simulated and evaluated for assessing the environmental risks. Some transformation products displayed increased environmental risks. The environmental risks of plant degradation products of the PCBs in the environmental media showed the maximum risk, indicating that the potential risks of the transformation products of the PCBs and their environmentally friendly derivatives could not be neglected. It is essential to further improve the ability of plants to degrade their transformation products. The improvement of some degradation products for environmentally friendly PCB derivatives indicates that the theoretical modification of a single environmental feature cannot completely control the potential environmental risks of molecules. In addition, this method can be used to analyze and evaluate environmentally friendly PCB derivatives to avoid and reduce the potential environmental and human health risks caused by environmentally friendly PCB derivatives.

Design of a Microbial Remediation Inoculation Program for Petroleum Hydrocarbon Contaminated Sites Based on Degradation Pathways

This paper analyzed the degradation pathways of petroleum hydrocarbon degradation bacteria, screened the main degradation pathways, and found the petroleum hydrocarbon degradation enzymes corresponding to each step of the degradation pathway. Through the Copeland method, the best inoculation program of petroleum hydrocarbon degradation bacteria in a polluted site was selected as follows: single oxygenation path was dominated by Streptomyces avermitilis, hydroxylation path was dominated by Methylosinus trichosporium OB3b, secondary oxygenation path was dominated by Pseudomonas aeruginosa, secondary hydroxylation path was dominated by Methylococcus capsulatus, double oxygenation path was dominated by Acinetobacter baylyi ADP1, hydrolysis path was dominated by Rhodococcus erythropolis, and CoA path was dominated by Geobacter metallireducens GS-15 to repair petroleum hydrocarbon contaminated sites. The Copeland method score for this solution is 22, which is the highest among the 375 solutions designed in this paper, indicating that it has the best degradation effect. Meanwhile, we verified its effect by the Cdocker method, and the Cdocker energy of this solution is −285.811 kcal/mol, which has the highest absolute value. Among the inoculation programs of the top 13 petroleum hydrocarbon degradation bacteria, the effect of the best inoculation program of petroleum hydrocarbon degradation bacteria was 18% higher than that of the 13th group, verifying that this solution has the best overall degradation effect. The inoculation program of petroleum hydrocarbon degradation bacteria designed in this paper considered the main pathways of petroleum hydrocarbon pollutant degradation, especially highlighting the degradability of petroleum hydrocarbon intermediate degradation products, and enriching the theoretical program of microbial remediation of petroleum hydrocarbon contaminated sites.

Enhanced plant-microbe remediation of PCBs in soil using enzyme modification technique combined with molecular docking and molecular dynamics

The study on the enhanced mechanisms of the enzymes involved in plant absorption, plant degradation, and microbial mineralization in the remediation of soils contaminated with polychlorinated biphenyls (PCBs) is of great significance for the application of plant-microbe combined remediation technique in PCB-contaminated soils. The present study first used a combination of molecular docking and molecular dynamics methods to calculate the effects of the plant absorption enzyme, plant degradation enzyme, and microbial mineralization enzyme on the PCBs in the soil environment. A multifunctional plant degradation enzyme was constructed with three functional roles of absorption, degradation, and mineralization using amino acid sequence recombination and site-directed mutagenesis to modify the template of plant degradation enzyme. Finally, using the Taguchi experimental design-assisted molecular dynamics simulation method, the suitable external environmental conditions of plant-microbe combined remediation of the PCB-contaminated soil were determined. In total, six multifunctional plant degradation enzymes were designed, which exhibited a significantly improved efficiency of PCB degradation. In comparison to the complex of plant absorption enzyme, plant degradation enzyme, and microorganism mineralization enzyme (6QIM-3GZX-1B85), the six multifunctional plant degradation enzymes exhibited significantly higher efficiency (2.10-2.38 times) in degrading the PCBs, with a maximum of 2.69 times under suitable external environmental conditions.

Mitigating the Adverse Effects of Polychlorinated Biphenyl Derivatives on Estrogenic Activity via Molecular Modification Techniques

The aim of this paper is to explore the mechanism of the change in oestrogenic activity of PCBs molecules before and after modification by designing new PCBs derivatives in combination with molecular docking techniques through the constructed model of oestrogenic activity of PCBs molecules. We found that the weakened hydrophobic interaction between the hydrophobic amino acid residues and hydrophobic substituents at the binding site of PCB derivatives and human oestrogen receptor alpha (hERα) was the main reason for the weakened binding force and reduced anti-oestrogenic activity. It was consistent with the information that the hydrophobic field displayed by the 3D contour maps in the constructed oestrogen activity CoMSIA model was one of the main influencing force fields. The hydrophobic interaction between PCB derivatives and oestrogen-active receptors was negatively correlated with the average distance between hydrophobic substituents and hydrophobic amino acid residues at the hERα-binding site, and positively correlated with the number of hydrophobic amino acid residues. In other words, the smaller the average distance between the hydrophobic amino acid residues at the binding sites between the two and the more the number of them, and the stronger the oestrogen activity expression degree of PCBS derivative molecules. Therefore, hydrophobic interactions between PCB derivatives and the oestrogen receptor can be reduced by altering the microenvironmental conditions in humans. This reduces the ability of PCB derivatives to bind to the oestrogen receptor and can effectively modulate the risk of residual PCB derivatives to produce oestrogenic activity in humans.

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.

Combined remediation of polychlorinated naphthalene-contaminated soil under multiple scenarios: An integrated method of genetic engineering and environmental remediation technology

This study explored the types of polychlorinated naphthalene (PCN)-contaminated soil and determined the practicable scheme of combined remediation using an integrated method of genetic engineering and environmental remediation technology. A multi-scenario comprehensive evaluation system of a plant-microbial combined bioremediation of PCN-contaminated soil was established using the intelligent integration of analytic hierarchy process and formula evaluation methods based on the current situation of PCN contamination in China, which showed the bioremediation of PCN-contaminated soil by the plant-microbial system could be divided into four scenarios. QSAR models were constructed to quantify the remediation mechanism that electronic parameter ∆E was the key factor changing the efficiency of combined bioremediation. Moreover, the macro-control scheme of PCN-contaminated soil was established, which indicated that four new multifunctional proteins promoted the absorption, degradation, and mineralization of PCNs in specific soil pollution types significantly, were obtained through cross gene recombination. The molecular dynamics (MD) simulation results showed the efficiency of the plant-microbial combined bioremediation were increased by 15.45% (Scenario 1, 2, 3) and 20.02% (Scenario 4) under the optimal regulation scheme. The findings will be helpful to realize the regional control of PCN-contaminated soil.

QSAR models for the acute toxicity of 1,2,4-triazole fungicides to zebrafish (Danio rerio) embryos

In recent decades, the 1,2,4-triazole fungicides are widely used for crop diseases control, and their toxicity to wild lives and pollution to ecosystem have attracted more and more attention. However, how to quickly and efficiently evaluate the toxicity of these compounds to environmental organisms is still a challenge. In silico method, such like Quantitative Structure-Activity Relationship (QSAR), provides a good alternative to evaluate the environmental toxicity of a large number of chemicals. At the present study, the acute toxicity of 23 1,2,4-triazole fungicides to zebrafish (Danio rerio) embryos was firstly tested, and the LC₅₀ (median lethal concentration) values were used as the bio-activity endpoint to conduct QSAR modelling for these triazoles. After the comparative study of several QSAR models, the 2D-QSAR model was finally constructed using the stepwise multiple linear regression algorithm combining with two physicochemical parameters (logD and μ), an electronic parameter (Q_N1) and a topological parameter (X^vPC₄). The optimal model could be mathematically described as following: pLC₅₀ = -7.24-0.30X^vPC₄ + 0.76logD - 26.15Q_N1 - 0.08μ. The internal validation by leave-one-out (LOO) cross-validation showed that the R²_adj (adjusted noncross-validation squared correlation coefficient), Q² (cross-validation correlation coefficient) and RMSD (root-mean-square error) was 0.88, 0.84 and 0.17, respectively. The external validation indicated the model had a robust predictability with the q² (predictive squared correlation coefficient) of 0.90 when eliminated tricyclazole. The present study provided a potential tool for predicting the acute toxicity of new 1,2,4-triazole fungicides which contained an independent triazole ring group in their molecules to zebrafish embryos, and also provided a reference for the development of more environmentally-friendly 1,2,4-triazole pesticides in the future.

Environmentally Friendly Fluoroquinolone Derivatives with Lower Plasma Protein Binding Rate Designed Using 3D-QSAR, Molecular Docking and Molecular Dynamics Simulation

Comparative molecular similarity index analysis (CoMSIA) was used to establish a three-dimensional quantitative structure–activity relationship (3D-QSAR) model with structural parameters of quinolones as the independent variables and plasma protein binding rate (logf_b) as the dependent variable to predict the logf_b values of remaining quinolones in this study. In addition, the mono-substituted and bis-substituted reaction schemes that significantly influenced the plasma protein binding rate of quinolones were determined through an analysis of the 3D-QSAR contour maps. It was found that the replacement of small groups, hydrophobic groups, electronegative groups, or hydrogen bond acceptor groups at the substitution sites significantly reduce the logf_b values of quinolone derivatives. Furthermore, the mechanism of decrease in binding rate between trovafloxacin (TRO) derivatives and plasma protein was revealed qualitatively and quantitatively based on molecular docking and molecular dynamics simulation. After modification of the target molecule, 11 TRO derivatives with low plasma protein binding rates were screened (reduced by 0.50–24.18%). Compared with the target molecule, the molecular genotoxicity and photodegradability of the TRO derivatives was higher (genotoxicity increased by 4.89–21.36%, and photodegradability increased by 9.04–20.56%), and their bioconcentration was significantly lower (by 36.90–61.41%).

Enhanced Biodegradation of Phthalic Acid Esters' Derivatives by Plasticizer-Degrading Bacteria (Burkholderia cepacia, Archaeoglobus fulgidus, Pseudomonas aeruginosa) Using a Correction 3D-QSAR Model

A phthalic acid ester’s (PAEs) comprehensive biodegradability three-dimensional structure-activity relationship (3D-QSAR) model was established, to design environmentally friendly PAE derivatives, which could be simultaneously degraded by plasticizer-degrading bacteria, such as Burkholderia cepacia, Archaeoglobus fulgidus, and Pseudomonas aeruginosa. Only three derivatives of diethyl phthalate (DEP (DEP-27, DEP-28 and DEP-29)) were suited for their functionality and environmental friendliness, which had an improved stability in the environment and improved the characteristics (bio-toxicity, bioaccumulation, persistence, and long-range migration) of the persistent organic pollutants (POPs). The simulation inference of the microbial degradation path before and after DEP modification and the calculation of the reaction energy barrier exhibited the energy barrier for degradation being reduced after DEP modification and was consistent with the increased ratio of comprehensive biodegradability. This confirmed the effectiveness of the comparative molecular similarity index analysis (CoMSIA) model of the PAE’s comprehensive biodegradability. In addition, a molecular dynamics simulation revealed that the binding of the DEP-29 derivative with the three plasticizer-degradation enzymes increased significantly. DEP-29 could be used as a methyl phthalate derivative that synergistically degrades with microplastics, providing directional selection and theoretical designing for plasticizer replacement.

Integration of Fuzzy Matter-Element Method and 3D-QSAR Model for Generation of Environmentally Friendly Quinolone Derivatives

The environmental pollution of quinolone antibiotics (QAs) has caused rising public concern due to their widespread usage. In this study, Gaussian 09 software was used to obtain the infrared spectral intensity (IRI) and ultraviolet spectral intensity (UVI) of 24 QAs based on the Density Functional Theory (DFT). Rather than using two single-factor inputs, a fuzzy matter-element method was selected to calculate the combined effects of infrared and ultraviolet spectra (CI). The Comparative Molecular Field Analysis (CoMFA) was then used to construct a three-dimensional quantitative structure–activity relationship (3D-QSAR) with QAs’ molecular structure as the independent variable and CI as the dependent variable. Using marbofloxacin and levofloxacin as target molecules, the molecular design of 87 QA derivatives was carried out. The developed models were further used to determine the stability, functionality (genetic toxicity), and the environmental effects (bioaccumulation, biodegradability) of these designed QA derivatives. Results indicated that all QA derivatives are stable in the environment with their IRI, UVI, and CI enhanced. Meanwhile, the genetic toxicity of the 87 QA derivatives increased by varying degrees (0.24%–29.01%), among which the bioaccumulation and biodegradability of 43 QA derivatives were within the acceptable range. Through integration of fuzzy matter-element method and 3D-QSAR, this study advanced the QAs research with the enhanced CI and helped to generate the proposed environmentally friendly quinolone derivatives so as to aid the management of this class of antibiotics.

Potential of topological descriptors to model the retention of polychlorinated biphenyls in different gas chromatography stationary phases, including ionic liquid-based columns

The aim of this study was to develop a statistical model based on a set of intuitive topological descriptors that will help to determine the influence of the polychlorinated biphenyls (PCBs) structural features on the chromatographic behavior of these analytes in a variety of gas chromatographic stationary phases, including the highly polar ionic liquid (IL)-based SLB-IL76 and SLB-IL60 columns. The model was developed using the stepwise multiple linear regression method, and constructed through several levels of increasing complexity to make evident the relative influence of the selected descriptors. The proposed model was easy to implement and provided similar satisfactory results irrespective of the dependent variables used (i.e., retention index or retention time) or the chromatographic conditions applied (i.e., pseudo-isotherm and programmed temperature) for IL-based phases. The model also allowed the correct prediction of the elution order of selected PCBs in these and other less polar phases evaluated (i.e., SW-10, DB-17, ZB-5 and HT-8). To our knowledge, this is the first models based on topological descriptors described in the literature that provided a satisfactory fitting of the PCB behavior in IL-based phases. Our results indicated that the mechanism governing the chromatographic separation of PCBs in these highly polar columns showed significant differences compared with those observed in other less polar stationary phases.

QSPR versus fragment-based methods to predict octanol-air partition coefficients: Revisiting a recent comparison of both approaches

The octanol-air partition coefficient (K_OA) is useful to assess the fate of organic chemicals in the environment. Very recently, an interesting comparison of current methods to predict this property (Chemosphere 148 (2016) 118-125) highlighted a newly introduced Quantitative Structure-Property Relationship (QSPR), as a group-contribution (GC) method and a quantum chemical solvation model were reported to yield significantly less accurate results. Based on the observation that the so-called GC method investigated in this earlier study was inconsistent with the temperature dependence of K_OA, the previously recommended QSPR is presently compared to the geometrical fragment (GF) additivity scheme. In addition to providing some improvement in terms of accuracy, this fragment-based procedure exhibits many advantages in terms of simplicity, interpretability, applicability and availability.

Prediction of octanol-air partition coefficients for PCBs at different ambient temperatures based on the solvation free energy and the dimer ratio

Temperature-dependent octanol-air partition coefficients (K_OA) are of great importance in assessing the environmental behavior and fate of persistent organic pollutants including polychlorinated biphenyls (PCBs). Due to the tremendous amounts of time, effort and cost needed for the experimental determination of K_OA, it is desirable to develop a rapid and precise predictive method to estimate K_OA just based on molecular structure. In the present study, a predictive model for log K_OA of PCBs at ambient temperatures was developed based on the thermodynamic relationship between K_OA and the solvation free energy from air to octanol (ΔG_OA). For the calculation of ΔG_OA of PCBs, the optimal combination of theoretical method and basis-set was identified to be HF/MIDI!6D for both geometry optimization and energy calculation. Dimer formation could affect the partition behavior and promote the apparent K_OA values of PCBs. After taking the effect of dimer formation into account, the goodness-of-fit, predictive ability, and robustness of the predictive model were significantly improved. Apparent log K_OA values of PCBs at different ambient temperatures ranging from 283.15 to 303.15 K were predicted. Compared with other reported models, the model developed in the present study had not only comparable goodness-of-fit and predictive ability, but also a universal application domain and the relative independency of experimental data. Therefore, the solvation free energy method could be a promising method for the prediction of K_OA.

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.

Combined 2D-QSAR, Principal Component Analysis and Sensitivity Analysis Studies on Fluoroquinolones' Genotoxicity

In this paper, two-dimensional quantitative structure-activity relationship (2D-QSAR) and principal component analysis (PCA) methods were employed to screen the main parameters affecting the genotoxicity of fluoroquinolones (FQs), and the rules affecting the genetic toxicity of FQs were investigated by combining 2D-QSAR and PCA with the sensitivity analysis method. First, four types of parameters were calculated, namely, the geometric parameters (7), electronic parameters (5), physical and chemical parameters (8), and spectral parameters (7), but the physical and chemical parameters heat of formation (HF) and critical volume (CV) were excluded after the establishment of the 2D-QSAR model. Then, after PCA, it was found that the first principal component represented the main driving factors affecting the molecular genetic toxicity of FQs. In addition, after comprehensive analysis of the factor loading of the first, second, and third principal components, seven parameters affecting the genotoxicity of the FQs were screened out, namely, total energy (TE), critical temperature (CT), and molecular weight (Mol Wt) (increased with increasing genotoxicity of the FQs) and steric parameter (MR), quadrupole moment QXX (QXX), quadrupole moment QYY (QYY), and boiling point (BP) (decreased with increasing genotoxicity of the FQs); the above key parameters were also verified by sensitivity analysis. The obtained rules could be used to determine the substitution sites and the substitution groups associated with higher genotoxicity in the process of FQ modification, and these rules agreed well with the hologram quantitative structure-activity relationship (HQSAR) model. Finally, it was also found through SPSS analysis that the parameters screened in this paper were significantly correlated with FQ derivatives' genetic toxicity.

Fuzzy Comprehensive Evaluation Assistant 3D-QSAR of Environmentally Friendly FQs to Reduce ADRs

Most studies on adverse drug reactions (ADRs) of fluoroquinolones (FQs) have focused on the mechanisms of single ADRs, and no quantitative structure–activity relationship (QSAR) method studies have been carried out that combine several ADRs of FQs. In this study, an improved three-dimensional (3D) QSAR method was established using fuzzy comprehensive evaluation. This method could simultaneously consider three common ADRs of FQs using molecular parameters. The improved method could comprehensively predict three ADRs of FQs and provide direction for the development of new drugs with lower ADRs than the originals. According to the improved method, 48 derivatives with lower ADRs (decreased by 4.86% to 50.92%) were designed from pazufloxacin. Three derivatives with a higher genotoxicity, higher photodegradation, and lower bioconcentration than pazufloxacin were selected using the constructed QSAR methods of the FQs. Finally, three traditional 3D-QSAR methods of single ADR were constructed to validate the improved method. The improved method was reasonable, with a relative error range of 0.96% to 4.30%. This study provides valuable reference data and will be useful for the development of strategies to produce new drugs with few ADRs. In the absence of complementary biological studies of these adverse drug reactions, the results reported here may be quite divergent from those found in humans or experimental animals in vivo. One major reason for this is that many adverse drug reactions are dependent upon enzyme-catalyzed metabolic activation (toxication) or on non-enzymatic conversion to toxic products and are not due to the parent drug moiety.

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.

Biological enrichment prediction of polychlorinated biphenyls and novel molecular design based on 3D-QSAR/HQSAR associated with molecule docking

Based on the experimental data of octanol-water partition coefficients (Kow, represents bioaccumulation) for 13 polychlorinated biphenyl (PCB) congeners, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were used to establish 3D-QSAR models, combined with the hologram quantitative structure-activity relationship (HQSAR), the substitution sites (mono-substituted and bis-substituted) and substituent groups (electron-withdrawing hydrophobic groups) that significantly affect the octanol-water partition coefficients values of PCBs were identified, a total of 63 monosubstituted and bis-substituted were identified. Compared with using 3D-QSAR model alone, the coupling of 3D-QSAR and HQSAR models greatly increased the number of newly designed bis-substituted molecules, and the logKow reduction in newly designed bis-substituted molecules was larger than that of monosubstituted molecules. This was established to predict the Kow values of 196 additional PCBs and carry out a modification of target molecular PCB-207 to lower its Kow (biological enrichment) significantly, simultaneously maintaining the flame retardancy and insulativity after calculation by using Gaussian09. Simultaneously, molecular docking could further screen out three more environmental friendly low biological enrichment newly designed PCB-207 molecules (5-methyl-PCB-207, 5-amino-PCB-207, and 4-amino-5-ethyl-PCB-207).

Biomagnification and risk assessment of polychlorinated biphenyls in food web components from Zhoushan fishing ground, China

Trophodynamics and risks of polychlorinated biphenyls (PCBs) in organisms from the Zhoushan fishing ground (ZFG), China were studied. Σ₂₂PCBs varied from 1.36 to 36.75 ng/g wet weight, which were far below the maximum residue levels allowed in fishery products. However, estimated daily intake and hazard ratio calculations present possible adverse effects due to PCB pollution. Significantly positive correlations appeared between wet-weighted concentrations of target chemicals and trophic levels (TLs) of the organisms, with trophic magnification factors (TMFs) from 1.15 to 9.72. The TMF values first increased with an increase of compound's K_OW values, and then decreased, with log K_OW around 7.0 as an inflection point. TL is suggested as the key factor controlling contaminant burden among the species for only PCBs 105, 138, 153, and 171. For the remaining PCBs, lipid content of the organism or metabolite capacity of the compound may be more important influence on their bioaccumulation.

Combined QSAR/QSPR and molecular docking study on fluoroquinolones to reduce biological enrichment

With the aim of reducing the adverse effects of fluoroquinolones in the environment, a complete design and screening system for the low biological enrichment and high photodegradabilities of 29 fluoroquinolones was established through a three-dimensional quantitative structure-activity relationship (3D-QSAR) model and molecular docking. The interaction mechanisms of the fluoroquinolones with Gram-negative bacteria (DNA gyrase in Escherichia coli) and Gram-positive bacteria (Topoisomerase IV in Staphylococcus aureus) were also evaluated. Consequently, the 3D-QSAR model showed that the 3- and 18-positions of the fluoroquinolones strongly affected their biological enrichment, and that the introduction of electropositive or hydrophobic groups at these positions reduced the logarithm of the octanol-water partition coefficient. Using nadifloxacin as a template, 23 derivatives with lower biological enrichment than nadifloxacin (decreased by 30.12%-94.18%) were designed. Meanwhile, the photodegradabilities of 15 derivatives were increased compared with nadifloxacin. Finally, the further screening by molecular docking of nadifloxacin and the above 15 derivatives with DNA gyrase and Topoisomerase IV showed that 13 of the derivatives had lower biological enrichment (decreased by 0.30%-16.76%) than nadifloxacin in the bacteria.

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.

Elucidation of endocrine-disrupting polychlorinated biphenyls binding potency with steroidogenic genes: Integration of in silico methods and ensemble docking approaches

A myriad of polychlorinated biphenyls (PCBs) may have potential to reproductive axis and endocrine disruptions. PCBs mostly breach the cholesterol biotransformation in mitochondria through interfering with steroidogenic genes and lead to adverse consequences in steroidogenesis; however, studies are scanty. In this examination, the combinations of quantitative structure-activity relationship (QSAR) modeling and ensemble docking approaches was performed to envisage structural properties of PCBs that influence the developmental toxicity, estrogen receptor-mediated impacts, and to provide a better comprehension of binding levels between PCBs, steroidogenic acute regulatory protein (StAR) and cholesterol side-chain cleavage enzyme (CYP11A1). Prognostic QSAR models were illustrated with good robustness and predictive ability. Models provide extensive data on applicability domain and similarities between PCBs and training set compounds was used to implement for clustering and classification of toxic PCBs by employing Self-Organizing Maps. Docking and interaction profiles interpretations provided various insights into the structural features of PCBs that influence the cholesterol binding to StAR and CYP11A1 domains. The non-polar, atomic π - stacking and halogen bonds of PCBs with novel hotspots residues of StAR and CYP11A1 was indicated subtle conformational changes that barrier to cholesterol binding and/or locks to cholesterol ejection from Ω1-loop of StAR, and inhibits cholesterol to pregnenolone biosynthesis by CYP11A1; thus, these are probably revealed as block-cluster mechanisms. These outcomes are potential to be useful to predict developmental toxicity, endocrine disruption potencies and anti-steroidogenic activities of other environmental pollutants and provided sorted pinpoints for further evaluation of interaction mechanisms of PCBs with other sterodogenic genes.

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.

Quantitative structure-property relationship (QSPR) models for predicting the physicochemical properties of polychlorinated biphenyls (PCBs) using deep belief network

Octanol/water partition coefficient (log P), octanol/air partition coefficient (log K_OA) and bioconcentration factor (log BCF) are important physiochemical properties of organic substances. Quantitative structure-property relationship (QSPR) models are a promising alternative method of reducing and replacing experimental steps in determination of log P, log K_OA and log BCF. In the current study, we propose a new QSPR model based on a deep belief network (DBN) to predict the physicochemical properties of polychlorinated biphenyls (PCBs). The prediction accuracy of the proposed model was compared to the results of previous reported models. The predictive ability of the DBN model, validated with a test set, is clearly superior to the other models. All results showed that the proposed model is robust and satisfactory, and can effectively predict the physiochemical properties of PCBs without highly reliable experimental values.

Molecular docking, 3D-QSAR and structural optimization on imidazo-pyridine derivatives dually targeting AT1 and PPARg

Telmisartan, a bifunctional agent of blood pressure lowering and glycemia reduction, was previously reported to antagonize angiotensin II type 1 (AT1) receptor and partially activate peroxisome proliferator-activated receptor γ (PPARγ) simultaneously. Through the modification to telmisartan, researchers designed and obtained imidazo-\pyridine derivatives with the IC₅₀s of 0.49∼94.1 nM against AT1 and EC₅₀s of 20∼3640 nM towards PPARγ partial activation. For minutely inquiring the interaction modes with the relevant receptor and analyzing the structure-activity relationships, molecular docking and 3D-QSAR (Quantitative structure-activity relationships) analysis of these imidazo-\pyridines on dual targets were conducted in this work. Docking approaches of these derivatives with both receptors provided explicit interaction behaviors and excellent matching degree with the binding pockets. The best CoMFA (Comparative Molecular Field Analysis) models exhibited predictive results of q²=0.553, r²=0.954, SEE=0.127, r²_pred=0.779 for AT1 and q²=0.503, r²=1.00, SEE=0.019, r²_pred=0.604 for PPARγ, respectively. The contour maps from the optimal model showed detailed information of structural features (steric and electrostatic fields) towards the biological activity. Combining the bioisosterism with the valuable information from above studies, we designed six molecules with better predicted activities towards AT1 and PPARγ partial activation. Overall, these results could be useful for designing potential dual AT1 antagonists and partial PPARγ agonists.

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.

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

Experimentally determined octanol-air partition coefficients (K _OA) for 43 polychlorinated naphthalene (PCN) congeners and experimentally determined subcooled liquid vapor pressures (P _L) for 17 PCN congeners were used with comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) to generate three-dimensional quantitative structure-activity relationship (3D-QSAR) models. The data were used to predict K _OA values for the other 32 congeners and P _L values for the other 58 congeners. The CoMFA and CoMSIA model contour maps showed that the electrostatic fields of the PCN molecules are the most important factors affecting the K _OA and P _L values. The long-range transport potentials of several PCN homologs were assessed using the following grading system: high mobility (MoCNs), relatively high mobility (DiCNs to TeCNs), relatively low mobility (PeCNs to HeCNs) and low mobility (HeCNs and OCN). The PCN-2 molecule was modified using the contour maps of the two models, and the results showed that introducing an electronegative R₁ substituent increased the K _OA value but introducing an electropositive R₆ substituent decreased the P _L value. PCN-2 was in the high mobility class, but introducing these substituents moved the long-range transport potentials of the modified molecules to the relatively high mobility class.