Molecular Design and Mechanism Analysis of Phthalic Acid Ester Substitutes: Improved Biodegradability in Processes of Sewage Treatment and Soil Remediation

Phthalic acid esters (PAEs) have the characteristics of environmental persistence. Therefore, improving the biodegradability of PAEs is the key to reducing the extent of ecological harm realized. Firstly, the scoring function values of PAEs docking with various degrading enzymes in sewage treatment were calculated. Based on this, a 3D-quantitative structure-activity relationship (3D-QSAR) model for PAE biodegradability was built, and 38 PAE substitutes were created. By predicting the endocrine-disrupting toxicity and functions of PAE substitutes, two types of PAE substitutes that are easily degraded by microorganisms, have low toxicity, and remain functional were successfully screened. Meanwhile, the differences in the mechanism of molecular degradation difference before and after PAE modification were analyzed based on the distribution characteristics of amino acid residues in the molecular docking complex. Finally, the photodegradability and microbial degradability of the PAE substitutes in the soil environment was evaluated. From the 3D-QSAR model design perspective, the modification mechanism of PAE substitutes suitable for sewage treatment and soil environment degradation was analyzed. We aim to improve the biodegradability of PAEs at the source and provide theoretical support for alleviating the environmental hazards of using PAEs.


Introduction
Phthalate esters (PAEs) are a kind of substitutes synthesized by esterification of phthalic anhydride with various alcohol compounds, which have high fat solubility and low water solubility [1,2]. PAEs can be detected in water, soil, and air environmental media as these are widely used and produced in bulk [1]. In addition, the toxic effect of ecological hormones, the toxicity of reproductive development, carcinogenesis, teratogenesis, and mutagenicity of PAEs have attracted wide attention worldwide [3]. According to earlier research, it takes approximately 400 years for plastics to degrade completely in water and soil. PAEs, the primary additives in plastics, are continuously released into the environment during plastic degradation. PAE degradation is primarily dependent on microbiological and photodegradation processes, but the degradation rate is slow, and the degree of mineralization is low [4]. It is difficult to degrade PAEs in the absence of ultraviolet light in water sediments [5]. The refractory nature of PAEs results in ecological hazards [6]. Therefore, improving the environmental degradability of PAEs while retaining their functional integrity is of great importance.
Although PAEs have low water solubility, a small part of them are still soluble in water. Due to atmospheric deposition, soil leaching, and domestic drainage, PAEs are widely lowing the process of molecular docking and the intercriteria correlation (CRITIC) (based on the results from interlayer correlation). A 3D-QSAR model was used to analyze the comprehensive biodegradability of PAEs during water treatment (W-3D-QSAR) and design substitute molecules. (2) The less toxic and easily biodegradable substitutes that exhibit good activity were screened by analyzing the score function value obtained when the PAE substitute molecules and biodegradable complex protein were docked. The QSAR model and the density functional theory (DFT) technique were used to determine the disruptive endocrine toxicity of PAEs (T-QSAR). The molecular docking technique was used to investigate the mechanisms of molecular degradation (before and after PAE modification). (3) The photodegradation and microbial degradation paths of the PAE substitutes were deduced. The energy barrier of the reaction was calculated, and the endocrine disrupting toxicity of the degradation products was predicted to evaluate the photo and biodegradability of the PAE substitutes in soil. (4) To investigate the molecular modification mechanism of PAE substitutes suitable for sewage treatment and soil degradation, the distribution characteristics of the three-dimensional isopotential map of the W-3D-QSAR were determined. The PAEs molecular soil micro-biological degradation 3D-QSAR model (S-3D-QSAR) and T-QSAR were identified. The aim of this study is to design PAE substitutes with high degradability in both sewage treatment and soil, to improve the biodegradability of PAEs themselves from the source, and to provide theoretical guidance for the preparation of environmentally friendly PAEs substitutes in the future.

Screening and Data Source of Molecular Biodegradable Enzymes for PAE Degradation during Sewage and Soil Treatment Processes
Sewage treatment processes and microorganisms in soil primarily rely on the efficiency of the secreted degradation enzymes to effectively degrade PAEs. During the process of sewage treatment, PAEs are primarily degraded by aerobic and anaerobic bacteria. Catalase (PDB ID: 1P80) [24], nitrate reductase (PDB ID: 1HZV), and esterase (PDB ID: 7CUZ), secreted by aerobic bacteria were selected as the primary degrading enzymes for aerobic degradation [25,26], while hydrolase (PDB ID: 2IOF) [27], oxidoreductase (PDB ID: 1NIC), and alcohol dehydrogenase (PDB ID: 3FMX), secreted by anaerobic bacteria, were selected as the primary degrading enzymes for anaerobic degradation [28,29]. Polyphenol peroxidase (PDB ID: 6TOY) [30], hydrolase (PDB ID: 1VA4) [27], hydrolase (PDB ID: 4PSD), and laccase (PDB ID: 6SOO) [30], which were primarily associated with PAE degradation, were chosen to study the degradation process in soil [20,31]. Du et al. identified a composite degradation enzyme (TVSP) consisting of four degradation enzymes using the proteinprotein docking technology to comprehensively evaluate the degradation effect of the four degradation enzymes in soil [32]. At the same time, the soil microbial degradation process of the PAE substitutes was evaluated by analyzing the degree of change in the docking scoring function value before and after PAE modification. The relevant data were obtained from the PDB protein structure database (https://www.rcsb.org, accessed on 25 June 2022) of the Brookhaven National Laboratory ( Figure S3).

Characterization of the Binding Ability of the PAE Molecules to Microbial Degrading Enzymes Used in the Sewage Treatment Process and Present in the Soil Environment-The Molecular Docking Method
The molecular docking method is widely used to study the interaction between molecules and receptor proteins. It can be used to consider the binding effect of ligands and receptors as a whole [32]. As research objects, we chose 28 PAE molecules (Table 1), and their molecular structures were obtained from the PubMed database ( Figure S2). We used the Define and Edit Binding Site module Receptor-Ligand Interactions of Discovery Studio 4.0 software. In the sewage treatment process, twenty-eight PAE molecules were used as ligand molecules to dock with six microbial degrading enzymes (1P80, 1HZV, 7CUZ, 2IOF, 1NIC, and 3FMX). Microbial complex degrading enzymes (TVSP) were selected to study the degradation process in the soil environment [33]. Before conducting the molecular docking process, all degrading enzymes (or the degrading enzyme complexes) were pretreated to remove the ligands, metal ions, and water molecules. Polar hydrogen and point charges were added. At the same time, the Define module was used to find and define the active sites that may bind to the systems . The active site coordinates for 1P80, 1HZV, 7CUZ,  2IOF, Figure S1. The PAE molecules were incorporated into the formed active cavity to rapidly dock with the degrading enzyme protein. Finally, the LibDock score was used to analyze the interactions between the ligands and receptors. The higher the docking score function value between the PAE molecules and the degrading enzymes (or composite degrading enzymes), the stronger the degradation ability of the enzymes or the composite degrading enzymes (secreted by microorganisms) [34].

Calculation of Comprehensive Effect Value of Microbial Degradability-The CRITIC Method
As an objective weighting method, the CRITIC method calculates the standard deviation and correlation coefficient between data, assigns weights to different data, and adds and calculates the total effect value of all data using the product of weight and each index [35]. The degradation activity data in this paper were chosen from the docking score function values of 28 PAEs molecules and six degrading enzymes secreted by microorganisms during the sewage treatment process. First, Formula (1), corresponding to the CRITIC method, was used to calculate the standard deviation between the PAE molecules and scoring function values corresponding to different degrading enzymes associated with sewage treatment. Secondly, Formula (2), corresponding to the CRITIC method, is used to calculate the conflict coefficient between other degrading enzymes associated with the sewage treatment process. Using the formula, the critical information on other degrading enzymes is calculated using the estimated standard deviation and conflict coefficient (3). Finally, the weights of the six degrading enzymes used to degrade PAEs were calculated using Formula (4), and the comprehensive effect values of the molecular ability of the six degrading enzymes were calculated using Formula (5). The specific calculation formula for the CRITIC method is as follows: The standard deviation (S j ) between the PAE molecules and the scoring function values of the different degrading enzymes associated with the sewage treatment process is calculated as follows: The correlation coefficient between the different degrading enzymes and the molecular scoring function value corresponding to the PAEs associated with the sewage treatment process (R j ) is calculated as follows: The vital information on the amount of different microbial degrading enzymes in the sewage treatment process (C j ) is calculated as follows: The weight of the six degrading enzymes (W j ) is calculated as follows: The comprehensive effect value corresponding to the simultaneous degradation of the PAE molecules degraded by the six degrading enzymes (Q) is calculated as follows: where x ij represents the value of the docking scoring function of the i PAE molecules and the j degradation enzyme, x j represents the average value of the docking scoring function of the 28th PAE molecule for the j degradation enzyme, and D j represents the value of the docking scoring function between the j degradation enzyme and 28 PAE molecules.

Construction of the 3D-QSAR Model for Comprehensive Effects of Biodegradation in PAEs Molecular Sewage Treatment Process-The SYBYL-X Software Method
SYBYL-X 2.0 was used to construct a 3D-QSAR model to study the comprehensive effects of biodegradation. Firstly, we drew the molecular structures of 28 PAEs using the Sketch Molecular module in SYBYL-X 2.0 software. In the Minimize module, we selected the Powell conjugate gradient method and the Tripos molecular force field (max iterations: 10,000; min energy change: 0.005 kcal/mol) and loaded the Gasteiger-Huckel charge. The lowest energy conformations of the PAE molecules were obtained [36]. The training set: test set ratio was set at 3:1 [15], and from the 28 PAE molecules, 21 PAE molecules were randomly selected as the training set, 7 PAE molecules were chosen as the test set, and diisopentyl phthalate (DIPP) molecules were used as the template molecules. The skeleton stacking of 28 PAE molecules was achieved using Align Database in the Applications module. Finally, the combined effect values of the PAE molecules and microbial degrading enzymes were input into the Training table. The W-3D-QSAR model of the PAE molecules was created using the comparative molecular field analysis (CoMFA) method. The types of molecular fields in the CoMFA field are electrostatic field (E) and solid field (S). The leave-one-out method was used to cross-validate the training set compounds, and the crossvalidation coefficient q2 and the optimal principal component number N were calculated. Regression analysis was performed through No Validation, and the non-cross validation coefficient r2, standard deviation SEE, and test value F were calculated. Finally, r2pred was calculated from the cross-validation of the test sets. W-3D-QSAR was completed when the above parameters met the modeling requirements [37,38]. The interface scoring function value of the PAE molecules was used to build the S-3D-QSAR, and the microbial complex degrading enzyme (TVSP) in the soil environment was the dependent variable. The scoring function values of 28 PAE molecules and six microbial degrading enzymes (1P80, 1HZV, 7CUZ, 2IOF, 1NIC, and 3FMX) associated with the sewage treatment process were calculated using Discovery Studio 4.0. The CRITIC weight method was used to calculate the comprehensive effect value (Score), as shown in Table 1. W-3D-QSAR was built using SYBYL-X and the combined effect values of the PAEs and microbial degrading enzymes as the response value. Table 1 also includes the scoring function values calculated from the interaction between the PAEs and the soil microbial complex degrading enzymes (TVSP). SYBYL-X was used to construct S-3D-QSAR using the TVSP scoring function value as the response value.

Evaluation of the 3D-QSAR Model to Study the Comprehensive Effect of Microbial Degradability of PAEs Associated with the Sewage Treatment Process
The evaluation parameters corresponding to the W-3D-QSAR model are presented in Table 2. The optimal principal component n for the W-3D-QSAR model is 8, and the crossvalidation coefficient q2 is 0.76. The value indicates that the model has a good predictive ability [39]. The non-cross validation coefficient R2 is 0.999, the standard deviation SEE is 0.009, and the F test value is 1159.179, indicating that the model has good fitting ability and stability. Overfitting was not observed [37]. The standard error obtained post crossvalidation (SEP) for the experimental and predicted values was 0.162. The cross-test coefficient corresponding to the external prediction set r2pred was 0.646, indicating high external prediction ability [38]. The contributions of the stereo (S) and electrostatic (E) fields to the combined effect of degradability were 76.5 and 23.5%, respectively, indicating that both spatial and electrical products had an impact on the microbial degradability of PAEs in the sewage treatment process. The spatial outcome was moderately significant [40]. Table 2 shows the parameters of the S-3D-QSAR model, and the evaluation of various parameters reveals that the constructed S-3D-QSAR model meets the model construction requirements. DEHP and DBP molecules were detected in sewage in high concentrations, and these were used as target molecules (Figure 1) [41]. The CoMFA-based three-dimensional equipotential diagram of the DEHP and DBP W-3D-QSAR models was examined ( Figure 2). The molecular stereoscopic field isopotential diagram was created, and the yellow area indicated that increasing the volume of the substituent groups in this area improved the biodegradability of the PAE molecules [42]. In the isophobic diagram of the molecular electrostatic field, the red region overlaps with the blue area. The blue area indicates that an increase in the positive electricity of the substituent group in this region promotes the biodegradability of PAEs [43]. DEHP and DBP molecules were detected in sewage in high concentrations, and these were used as target molecules (Figure 1) [41]. The CoMFA-based three-dimensional equipotential diagram of the DEHP and DBP W-3D-QSAR models was examined ( Figure 2). The molecular stereoscopic field isopotential diagram was created, and the yellow area indicated that increasing the volume of the substituent groups in this area improved the biodegradability of the PAE molecules [42]. In the isophobic diagram of the molecular electrostatic field, the red region overlaps with the blue area. The blue area indicates that an increase in the positive electricity of the substituent group in this region promotes the biodegradability of PAEs [43].  As can be seen from the model molecular stereo field isopotential diagram ( Figure  2), the distribution of yellow areas at the ends of sites 3, 5, and 8 of the DEHP molecule and sites 5, 6, 7, and 8 of the DBP molecule indicate that an increase in the volume of the  DEHP and DBP molecules were detected in sewage in high concentrations, and th were used as target molecules (Figure 1) [41]. The CoMFA-based three-dimensional eq potential diagram of the DEHP and DBP W-3D-QSAR models was examined (Figure The molecular stereoscopic field isopotential diagram was created, and the yellow a indicated that increasing the volume of the substituent groups in this area improved biodegradability of the PAE molecules [42]. In the isophobic diagram of the molecu electrostatic field, the red region overlaps with the blue area. The blue area indicates t an increase in the positive electricity of the substituent group in this region promotes biodegradability of PAEs [43].  As can be seen from the model molecular stereo field isopotential diagram (Figu 2), the distribution of yellow areas at the ends of sites 3, 5, and 8 of the DEHP molec and sites 5, 6, 7, and 8 of the DBP molecule indicate that an increase in the volume of  Table 3. These included 19 alternatives to DBP and 19 alternatives to DEHP (Table S1). Among the DBP substitutes that meet the increased biodegradability requirement, singlesubstituted molecules accounted for 21.05% (biodegradability increase: 0.34-0.85%), doublesubstituted molecules accounted for 52.63% (biodegradability increase: 0.17-41.98%), and triple-substituted molecules accounted for 26.32% (biodegradability increase: 2.22-21.84%) of the total molecules. All of the DEHP substitutes that can meet the biodegradability requirement are double-substituted (biodegradability increase: 0.18-17.46%). The percentage of mono-, di-, and tri-substituted molecules and the increase in microbial degradability during water treatment revealed that the rate of substituted molecules that could meet the increased degradability requirement during water treatment (following the modification of PAEs by mono-and di-substitution) was lower than that achieved for di-substituted molecules. The microbial degradability was not significantly enhanced. The molecular biodegradation effect of PAE substitutes subjected to single substitution processes was the worst, and the increase in biodegradability was ≤1.00%. When PAE molecules are modified under conditions of double substitution, the number of alternative molecules that can meet the requirement is relatively large. The improvement of microbial biodegradability is significant when the -CH3 unit present at sites 5 and 8 of DBP was replaced by -CH2CH2CH2CH2CH3, and under these conditions, the biodegradability of DBP was increased by 41.98%. Therefore, when designing degradable PAE substitutes, priority should be given to double substitution to ensure that the designed substitute molecules significantly improve degradability and design rate.

Aerobic/Anaerobic Biodegradability, Endocrine Disturbance, and Functional Evaluation of PAE Substitutes
To further verify that the PAE alternatives designed by us are readily biodegradable under both aerobic and anaerobic conditions, firstly, the degradation enzymes secreted by aerobic bacteria (1P8O, 1HZV, and 7CUZ) and the degradation enzymes secreted by anaerobic bacteria (2IOF, 1NIC, and 3FMX) were used for protein-protein docking through Discovery Studio 4.0. Aerobic biodegradable complex protein (AEDECP) and anaerobic biodegradable complex protein (ANDECP) were obtained, and the scoring function values were obtained by interacting 38 designed PAEs with the two biodegradable complex proteins [42]. The changes in the molecular scoring function values (before and after molecular modification) occurring post docking of the alternative molecules with AEDECP were analyzed. Thirteen DBP alternatives (3.64-30.34%) and nineteen DEHP alternatives were easily degraded (1.22-58.91%). Following the docking of PAE alternative molecules with ANDECP, 11 DBP alternatives (0.41-41.23%) and all DEHP alternatives were easily degraded (6.78-95.39%). This indicated that the PAE substitutes designed were easily degraded under aerobic and anaerobic conditions [44]. Aerobic microorganisms are better at degrading PAEs than anaerobic microorganisms under sewage treatment conditions [11]. The biodegradability of the PAE substitutes (recorded under anaerobic conditions) designed by us is significantly greater than that achieved in the presence of aerobic microorganisms. This can compensate (to a certain extent) for the poor degradation effect exerted by anaerobic microorganisms during the wastewater treatment process. Under these conditions, the biodegradation efficiency of the PAEs is improved as a whole during the wastewater treatment process.
The T-3D-QSAR model constructed by Han et al. was used to predict the molecular toxicity of the PAE substitutes, and it was expected that the molecular endocrine-disrupting toxicity of the PAE substitutes would not change significantly (within 10%) [45]. In addition, the energy of the stable molecular configurations was calculated using the DFT technique, and the molecular functions corresponding to the designed PAE substitutes were tested [46]. The positive frequency values of 10 DBP molecules and their substitutes were greater than 0, and positive frequency values of 15 DEHP molecules and their substitutes were greater than 0. When the molecule's positive frequency value is greater than 0, the molecular energy is locally the lowest in a particular dimension, and the molecule can exist stably under these conditions [47]. The rate of change of the energy gap for four DBP and five DEHP substitutes was less than 10. The energy gap value represents the degree of insulation. The larger the energy gap, the stronger the insulation, that is, the more stable the molecular structure [47]. Analysis of the above principles reveals that the functions of the PAE substitutes designed by us do not decrease on the whole, and the molecules can exist stably in the environment.

Analysis of the Mechanism of Differential Microbial Degradation before and after Molecular Modification of PAEs
The biodegradability of the designed DEHP alternatives is significantly improved compared to that of the DBP alternatives. Therefore, the distribution characteristics of the amino acid residues around the molecular docking complexes were used to analyze the differences in molecular biodegradation before and after the modification of PAEs. The results provide theoretical guidance for the design of PAE substitutes. The PAE molecules were docked with AEDECP and ANDECP before and after modification, respectively [42], and the distribution of the amino acid residues around the complex is shown in Figure 3. When DEHP-10 molecules bind to AEDECP and ANDECP, they form Π-Π bonds with multiple hydrophobic amino acids, unlike DEHP molecules. Therefore, the volume of the molecules exposed to solvents is reduced [48]. The hydroxyl group substituent at site 5 of the DEHP-10 molecule is a polar group with solvent affinity [49]. Therefore, the interaction force between DEHP-10 molecules and polar amino acids is strengthened, and the formation of a hydrophilic environment pro-motes the binding of molecules to complex proteins [50,51]. When DEHP-12 binds to AEDECP and ANDECP, hydrogen bonds and Π-Π bonds are formed between the hydroxyl groups of the substituent groups at site 5 of DEHP-12 and the surrounding amino acids. Therefore, the interaction between the polar amino acids and the molecules is robust [52]. This promotes the binding with complex proteins. These results indicate that the molecular biodegradability of the DEHP substitute increases significantly (40.97-93.72%) during aerobic and anaerobic biodegradation. This can be attributed primarily to interactions between molecules and amino acids. During the binding of DEHP substitutes to complex proteins, numerous Π-Π bonds are formed. Moreover, the substituent groups in the DEHP substitute molecules are primarily polar hydroxyl groups that exhibit solvent affinity. The number of polar amino acids around the molecules accounts is high. This promotes mutual binding between the DEHP substitute molecules and the complex proteins. A comparison of the degree of biodegradability of the DEHP and DBP substitutes revealed that the nature of the different substituent groups on the PAE substitutes dictated the differences in their microbial degradability. The ability to form hydrogen and Π-Π bonds with surrounding amino acids dictated the binding efficiency of the microbial degradation enzymes. Therefore, polar substituent groups should be selected while designing PAE substitute molecules. This is due to the fact that polar substituent groups are solventfriendly and can easily form hydrogen and Π-Π bonds with the amino acids around them. This improves the binding efficiency of the PAE substitutes, resulting in a greater degree of microbia degradability.

Soil Environmental Photodegradation and Microbial Degradation Pathways (before and after DEHP Molecular Modification)
The PAEs that cannot be removed during sewage treatment can enter the soil environment through secondary biochemical effluents [13], thus, affecting the soil farming environment and crop growth conditions. These can enter the human body through the food chain, posing a threat to human health [15]. Improving the degradability of PAEs in the soil is, therefore, critical. According to previous research, photodegradation and microbial degradation are the primary mechanisms for PAE degradation in soil [20,55]. Further- DBP-5 molecules bind to AEDECP and ANDECP, and the substituent ethyl group at the molecular site 6 of DBP-5 is a non-polar group that can form Π-Π bonds with PHE and ARG to improve the binding ability of the system and promote the degradation of complex proteins [49]. However, the strength of the electrostatic interactions between the DBP-5 molecule and the surrounding amino acids is reduced, and the strength of the van der Waals force is increased. The van der Waals force is an attraction between the molecules, and the hydrogen bond is the binding force of the H atom and other atoms. The hydrogen bond is stronger than the van der Waals force [53,54]. The binding ability of the DBP-5 molecule with the degradable complex protein improved, but not significantly. DBP-7 interacts with AEDECP and ANDECP, and the ethyl group at the molecular site 6 of DBP-7 is non-polar. The number of hydrophilic amino acids surrounding DBP-7 is greater than that surrounding DBP. Nonetheless, the number of interacting bonds between the DBP-7 molecule and the surrounding amino acids decreased, and hydrogen bond interactions were less than van der Waals interactions. Therefore, the binding ability of the DBP-7 molecule increased to a certain extent, but the change was not significant.
A comparison of the degree of biodegradability of the DEHP and DBP substitutes revealed that the nature of the different substituent groups on the PAE substitutes dictated the differences in their microbial degradability. The ability to form hydrogen and Π-Π bonds with surrounding amino acids dictated the binding efficiency of the microbial degradation enzymes. Therefore, polar substituent groups should be selected while designing PAE substitute molecules. This is due to the fact that polar substituent groups are solventfriendly and can easily form hydrogen and Π-Π bonds with the amino acids around them. This improves the binding efficiency of the PAE substitutes, resulting in a greater degree of microbia degradability. The PAEs that cannot be removed during sewage treatment can enter the soil environment through secondary biochemical effluents [13], thus, affecting the soil farming environment and crop growth conditions. These can enter the human body through the food chain, posing a threat to human health [15]. Improving the degradability of PAEs in the soil is, therefore, critical. According to previous research, photodegradation and microbial degradation are the primary mechanisms for PAE degradation in soil [20,55]. Furthermore, two kinds of DEHP substitutes (DEHP-10 and DEHP-12), characterized by high biodegradability, were designed and screened as examples to verify whether they are characterized by high photodegradation and microbial degradation potential in the soil environment. The DEHP molecules are cracked and hydrolyzed in the soil during photodegradation [56]. DEHP molecules are hydroxylated and cleaved by light to form 2-carboxybenzaldehyde (E-1) and octyl benzoate (E-2), respectively. The C-O bonds in E-1 and E-2 are broken to form phthalic acid (E-3) and benzoic acid (E-4), and the compounds are finally mineralized to CO2 and H2O, resulting in linear degradation [57]. The photodegradation pathway of the DEHP molecule can be used to infer the photodegradation pathway of the substitute molecule [55], as shown in Figure 4. The hydroxylation and cleavage reactions of DEHP substitutes also occur under the action of light. During the degradation of DEHP substitutes, -OH attacks the side chain and causes the C-O bond to break. This results in the formation of phthalic acid and benzoic acid, which are finally mineralized into CO2 and H2O following cleavage.

Evaluation and Mechanism
The primary microorganisms associated with microbial degradation in soil include Geobacillus, Pseudomonas adaceae, Trichoderma reesei, and Streptomyces [58], and the primary degradation mode involves the synergistic degradation of the compound under the action of mixed strains [10]. DEHP undergoes hydroxylation, oxidation, ortho-ring opening, and meta-ring opening cleavage reactions [59]. First, the DEHP molecule is hydroxylated to 2-carboxybenzaldehyde (DE-1). The C-O bond in DE-1 is then broken to form phthalic acid (DE-2), which is then oxidized to protocatechuic acid . Under epicyclic conditions, DE-5 is cleaved into acetyl-CoA, succinic acid (DE-8), pyruvic acid, and oxaloacetic acid (DE-9). It then enters the tricarboxylic acid cycle, where it is oxidized to CO2 and H2O [7]. The molecular microbial degradation path of DEHP was inferred from the molecular microbial degradation path of the substitute, as shown in Figure 4. Hydroxylation, oxidation, ortho ring opening, and meta-ring opening cracking reactions also occurred during the microbial degradation of the DEHP substitutes. During the degradation of the DEHP substitute molecules, -OH attacks the side chain and causes the C-O bond to break to form benzoic acid ester products. These are then oxidized and cracked to produce acetyl-CoA, succinic acid, pyruvate acid, and oxaloacetate. Finally, CO2 and H2O are produced under conditions of oxidation.

Evaluation of the Extent of Photodegradation and Microbial Degradation Responses Produced and Product Toxicity in Soils before and after Molecular Modification by DEHP
The photodegradation and microbial degradation processes in soil associated with DEHP substitutes, as well as the product's toxicity, were investigated. To calculate or predict the reaction energy barrier and endocrine-disrupting toxicity associated with the photodegradation and microbial degradation pathways of the DEHP substitutes, the DFT technique, and the H-QSAR model were used. The greater the extent of molecular photodegradation and microbial degradation of DEHP substitutes, the lower the reaction energy barrier [60]. It was calculated that the E10-1 reaction energy barrier (∆E = 1.05 kJ/mol) corresponding to the DEHP-10 molecular pathway A (associated with soil photodegradation) was 60.08% lower than the E-1 reaction energy barrier (∆E = 2.63 kJ/mol) recorded before modification. The E10-2 reaction energy barrier in path B (∆E = 1.84 kJ/mol) was 64.95% lower than the E-2 reaction energy barrier achieved before modification (∆E = 5.25 kJ/mol). The E12-1 reaction energy barrier (∆E = 1.58 kJ/mol) associated with the DEHP-12 molecular pathway A was 39.92% lower than the E-1 reaction energy barrier achieved before modification. The E12-2 reaction energy barrier in path B (∆E = 3.94 kJ/mol) was 24.95% lower than the E-2 reaction energy barrier achieved before modification. This indicates that the designed DEHP alternatives were more photodegradable in soil than the DEHP molecules. The toxicity of the photodegradation intermediates does not change because the same degradation intermediates are formed when DEHP substitutes and DEHP are degraded.  The reaction energy barrier corresponding to the soil microbial degradation pathway for the degradation of the DEHP substitute was calculated. It was observed that the DE10-1 reaction energy barrier (∆E = 1.05 kJ/mol) in path A corresponding to the DEHP-10 alternative molecule was 60.08% lower compared to that of the DE-1 molecule before modification (∆E = 2.63 kJ/mol). The energy barrier for the DE10-2 reaction in path B (∆E = 2.42 kJ/mol) was 91.62% lower than the energy barrier recorded for the DE-2 reaction before modification (∆E = 28.88 kJ/mol). The DE12-1 reaction energy barrier (∆E = 1.68 kJ/mol) for the DEHP-12 molecular pathway A was 36.12% lower than that recorded for the DE-1 reaction energy barrier before the modification. The energy barrier for the DE12-2 reactions in path B (∆E = 1.35 kJ/mol) was 95.33% lower than for the DE-2 reaction before the modification. This indicates that it is easier to degrade the designed DEHP alternative in the presence of microorganisms compared to the DEHP molecules. Two different intermediates (DE10-1 and DE12-1) were formed during the process of microbial degradation of the DEHP replacement molecule in soils. The T-QSAR model was used to predict the toxicity values of the DE10-1 and DE12-1 molecules [20]. The results showed that DE10-1 and DE12-1 were 4.85 and 17.26% more toxic than DE-1, respectively. Analysis revealed that DE10-1 and DE12-1 altered their toxicity by lengthening the carbon chain following substitution reactions and changing molecular weights [61]. In contrast, the intermediate DE12-1 in pathway A bears an aldehyde moiety that can covalently crosslink with proteins, nucleic acids, and lipids in organisms [62]. This results in protein glycosylation or carbonylation [63,64] and the upregulation of the DE12-1 endocrinedisrupting toxicity [7]. Therefore, given the significant increase in the toxicity value of DE12-1, an intermediate produced during the soil microbial degradation process of the DEHP-12 molecule, it is not recommended as a choice for the DEHP molecule.

Analysis of the DEHP Substitutes Suitable for Sewage Treatment and the Soil Microbial Degradation Mechanism
To further validate that DEHP-10 and DEHP-12 molecules follow easy degradation mechanisms, difference analysis methods were used from the perspective of the three-dimensional equipotential diagram of the W-3D-QSAR, S-3D-QSAR, and T-QSAR models ( Figure 5). As can be seen from Figure 5, analysis of the stereo field isopotential maps corresponding to the W-3D-QSAR and S-3D-QSAR models reveals that yellow regions are distributed at the ends of the 5 and 8 methyl groups of the DEHP molecule. It indicates that the groups located at these sites are replaced by larger groups, which promotes the degradation of molecules during the sewage treatment process in the soil environment [65]. Blue areas are distributed at the ends of the 5, 6, and 8 methyl groups of the DEHP in the isopotential diagram of the model electrostatic field. These findings indicated that the substitution of positively charged groups at these sites aided in the degradation of molecules during the sewage treatment process in the soil environment [43]. Analysis of the superimposed plot of the T-QSAR model reveals that the DEHP molecule binds to three hydrogen-bonded receptors at sites 5 and 8, suggesting that substitution at these sites can alter the endocrine-disrupting properties of the molecule [66]. Therefore, common substitution sites (5 and 8) are present in the three-dimensional equipotential map of the W-3D-QSAR, S-3D-QSAR, and T-QSAR models, and the same substituent selection rules are used to improve the degradability of the PAEs.
The DEHP-10 and DEHP-12 molecules screened in this paper were designed with a double substitution at sites 5 and 8, and the carbon chain length was extended. The substituent groups were changed to improve the microbial degradability during sewage treatment and ensure their toxicity did not increase [65]. When the DEHP-10 molecule was doubly substituted at sites 5 and 8, the -CH3 group at site 5 was replaced by the -CH2CH2OH group, and the -CH3 group at site 8 was replaced by the -CH2CH2CH3 group. The DEHP-12 molecule was doubly substituted at sites 5 and 8, whereas the -CH3 group at site 5 replaced by the -CH2OH group and the -CH3 group at site 8 replaced by the -CH2CH2CH2CH2CH2OH group. The modification pathways associated with the DEHP-10 and DEHP-12 molecules are consistent with the common distribution characteristics of the 3D isopotential maps of the W-QSAR, S-QSAR, and T-QSAR models. Therefore, the PAE replacements designed by us can satisfy both the biodegradability conditions of the sewage treatment process and the soil environment conditions while safeguarding against elevated toxicity. stituent groups were changed to improve the microbial degradability during sewage treatment and ensure their toxicity did not increase [65]. When the DEHP-10 molecule was doubly substituted at sites 5 and 8, the -CH3 group at site 5 was replaced by the -CH2CH2OH group, and the -CH3 group at site 8 was replaced by the -CH2CH2CH3 group. The DEHP-12 molecule was doubly substituted at sites 5 and 8, whereas the -CH3 group at site 5 replaced by the -CH2OH group and the -CH3 group at site 8 replaced by the -CH2CH2CH2CH2CH2OH group. The modification pathways associated with the DEHP-10 and DEHP-12 molecules are consistent with the common distribution characteristics of the 3D isopotential maps of the W-QSAR, S-QSAR, and T-QSAR models. Therefore, the PAE replacements designed by us can satisfy both the biodegradability conditions of the sewage treatment process and the soil environment conditions while safeguarding against elevated toxicity.

Conclusions
This paper proposed a molecular screening system for easily degradable substitutes based on the prediction of degradability models for water treatment processes and soil environment, and the degradation pathways were constructed for PAE molecules. Based on the W-3D-QSAR, S-3D-QSAR, and T-3D-QSAR models, two alternative molecules for PAE with unchanged endocrine disruptor toxicity and easy degradation were successfully designed. The molecular degradability of the molecules was investigated during the sewage treatment process and in the soil environment. It was discovered that the molecular degradability of the PAE substitutes screened in this paper increased by 51.31-97.13% during sewage treatment and 33.18-62.94% in the soil environment. The distribution characteristics of the amino acid residues around the PAE molecules and docking complexes (before and after modification) were analyzed, and the coupling conditions were designed based on the 3D-QSAR model. Finally, the results provide theoretical guidance for the molecular design of easily degradable PAE substitutes to solve the problem of molecular persistence of PAEs in nature.

Supplementary Materials:
The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/toxics10120783/s1, Figure S1: The docking site of proteins; Figure  S2: Molecular structure diagram of PAEs; Figure S3: The structure of proteins; Table S1: Molecular structure of DBP and DEHP substitutes.