Control Strategies of Plastic Biodegradation through Adjusting Additives Ratios Using In Silico Approaches Associated with Proportional Factorial Experimental Design

Plastics, as a polymer material, have long been a source of environmental concern. This paper uses polystyrene plastics as the research object, and the relative contribution of each component of plastic additives to plastic degradation is screened using the molecular dynamics method. The factorial experimental design method is combined with molecular dynamics simulation to adjust the additive composition scheme, analyze the mechanism of interaction between the additive components, and select the plastic additive combination that is most readily absorbed and degraded by microorganisms. Seven different types of plastic additives, including plasticizers, antioxidants, light and heat stabilizers, flame retardants, lubricants, and fillers, are chosen as external stimuli affecting the biodegradability of plastics. Using molecular dynamics simulation technology, it is demonstrated that plastic additives can promote the biodegradability of plastics. The factorial experimental design analysis revealed that all plastic additives can promote plastic biodegradation and plasticizer is the most favorable factor affecting plastic degradation, that hydrophobicity interactions are the primary reason for enhancing plastic degradation, and that screening No. 116–45 (plasticizer A, light stabilizer C, flame retardant E) is the most advantageous combination of biodegradable plastic additives. The plastic biodegradation effect regulation scheme proposed in this study is based on optimizing the proportion of additive components. To continue research on aquatic biodegradable plastics, the optimal combination of plastic components that can be absorbed and degraded by microorganisms is recommended.


Introduction
Plastic, a type of polymeric material obtained through the addition or condensation polymerization of a resin monomer as a raw material [1], is widely used in all facets of life and production due to its affordability and durability [2]. Nowadays, annual plastic production can reach 350 million tons [3] and the average person discards 52 kg of plastic each year [4]. There are numerous types of plastics, including polyethylene (PE), polypropylene (PP), polystyrene (PS), poly (ethylene terephthalate) (PET) and polyvinyl chloride (PVC), among others [5], that can cause loss during marine transportation, and their microplastic particles enter the ocean as a result of randomly discarded plastic waste or factory discharge channels [6,7].PS is widely consumed, and it is one of the most frequently monitored types of plastic fragments in the environment [8]. Numerous researchers have examined the toxicity of PS plastics to marine organisms. Wan et al. [9] discovered that PS plastic could alter the metabolic spectrum of juvenile zebrafish and impair glucose, lipid, and energy Int. J. Environ. Res. Public Health 2022, 19, 5670 2 of 16 metabolism. Sussarellu et al. [10] discovered that oysters exposed to PS plastic experienced significant decreases in oocyte number, diameter, and sperm velocity of 38, 5 and 23%, respectively, and in offspring number and developmental velocity of 41 and 18%, respectively, confirming that exposure to PS microplastics causes reproductive disturbances in oysters and has serious consequences for offspring growth and development. Yu et al. [11] exposed Cryptomeria hidradiens to varying concentrations of PS microplastics to determine the toxicity and mechanism of exposure. They discovered that the degree of oxidative stress and intestinal damage were significantly correlated with the level of PS microplastic exposure. It is hypothesized that PS microplastics cause toxicity through oxidative stress and intestinal injury. As a result, consideration should be given to how to effectively manage PS plastics and mitigate their impact on the marine ecosystem.
Along with basic polymeric monomers, plastics contain a variety of plastic additives that enhance the activity, oxidation resistance, and flame retardancy of the material, among other properties [12,13]. Because the majority of plastic additives are not covalently bonded to the monomer polymer, they can be released into the environment and migrated [14]. Additionally, the additives used in plastic products have the potential to be harmful to organisms [15], but their environmental risks have been underestimated. The majority of plastic additives are endocrine disruptors, such as phthalates, brominated flame retardants, and bisphenol A [16][17][18][19], which are biotitic when released into the environment. Browne et al. [20] discovered that nonylphenol released by PVC significantly reduced Arenicolamarina's coelomic osteocytes' ability to remove pathogenic bacteria by more than 60%, and triclosan released by PVC significantly reduced its ability to produce deposits and resulted in worm death of more than 55%. Renzo et al. [21] carried out autopsies on 28 sea turtles from the Molise Adriatic coast. In all sea turtle fat and liver tissues, bisphenol A, terephthalic acid, and polyethylene terephthalate were detected. Nobre et al. [22] investigated the embryonic developmental toxicity of microplastic particles on Lytechinus variegatus and discovered that microplastic particles increased abnormal embryonic development levels by up to 66.5%, that microplastics can act as carriers of contaminants, and that plastic additives can also be leached into the environment and thus become toxic.
Because of the increasing usage of plastic products, disposing of this polymeric material has become a big concern for developed countries, as plastic pollution is persistent due to its durability [23][24][25]. Plastics are classified as thermosets and thermoplastics [26], with thermosets containing heteroatoms in the main chain, which are susceptible to hydrolytic degradation of chemical bonds such as ester or amide bonds, and thermoplastics containing long carbon chains in the main chain, which are resistant to hydrolytic degradation of chemical bonds [27]. As a thermoplastic, degradation of PS is even more challenging. Blanco et al. conducted a 13-month isothermal degradation experiment carried out at relatively low temperature (423 K) on PS and evidenced appreciable mass loss in the investigated period [28]. It was discovered that strain Rhodococcus ruber C208 could grow in the absence of a carbon source by using PS plastic as a carbon source [29]. After more than eight weeks of culture, strain C208 formed a dense biofilm and degraded a portion of PS [30], demonstrating that PS plastic can be disposed of via the biodegradation pathway.
Chen et al. [31] discovered that the aquatic toxicity of various plastic additive combinations varied and screened them for the best component combination of plastic components with the least aquatic toxicity. This paper employs keratinase from Moesziomyces antarcticus as a biodegradation protein for plastic PS and screens the major types of additives in plastic fractions that affect plastic degradation, including the best combination of plastic fractions that can be absorbed and degraded by microalgae to the greatest extent possible, in order to design a regulatory scheme for plastic biodegradation effect based on the optimization of additive fractionation ratio.

Biodegradability Assessment of the Plastic Component PS and Its Additives Using Molecular Docking and Molecular Dynamics Simulations
To investigate the biodegradability of plastic components (in this case, PS) and their additives by microorganisms, the PDB database (Protein Data Bank, PDB) was used to retrieve the keratinase (PDB ID: 7CC4) from a PS-degrading microorganism, Moesziomyces antarcticus, as a biodegrading protein for PS plastics [32].
Molecular docking is a technique for rapidly docking protein receptors with their corresponding ligand molecules. It is implemented in the DS software module Dock-Ligands (LibDock). It is primarily based on the Lock and Key principle of rigid binding between ligands and receptors, in which the interaction between ligands and receptors is modeled using the spatial shape and energy matching of the ligands and receptor macromolecules, and the complex's conformation and binding affinity can be predicted [33]. Firstly, the incomplete amino acid residues of the protein were supplemented, and the protein was hydrogenated before molecular docking. Secondly, the Powell method was used to minimize the energy of molecules in SYBYL-2.0 software (parameters are as follows: Min Energy Change: 0.005 kcal/mol; Loading charge: Gasteiger-Huckel, Max Iterations: 10,000 Times; other parameters are default). All other parameters were set to their default values [34]. Finally, the module Receptor-Ligand Interactions is used to locate and define the receptor protein's binding cavity, and the optimized ligand molecules are loaded into the DS software for docking with the receptor protein. The Libdock score is used to quantify the results of molecular docking, with higher scores indicating stronger interactions between ligand molecules and receptor proteins, more stable binding, and greater biodegradability [35].
Molecular Dynamics (MD) simulation is a molecular simulation technique based on Newtonian mechanics that is increasingly used in biochemistry, materials science, and other fields to obtain information about the thermodynamics and kinetics of receptor proteins and other macroscopic properties [36]. The dynamics of the molecular docking complex were simulated in this paper using GROMACS 4.6.5 software on a Dell PowerEdge R7525 server under the Gromos96 43A1 force field. At first, a cube box is constructed as the kinetic simulation region, into which the ligand-receptor complex is placed. Water is then added as the environmental medium in the blank group (the additive group must also add external stimuli) [37], and to balance the charge of the system and maintain it electrically neutral, Na + is added to the box in the same amount as water molecules. Second, energy minimization via the fastest descent method, where the energy converges to 1000 kJ/mol and the system energy is equal to the equilibrium state's energy. Then, in sequence, NVT temperature control (set at 300 K) and NPT pressure control (set to 1 bar to bind the ligandprotein complex site) are performed. Finally, at the MD equilibrium stage, the complex was unbound from its initial position and kinetically simulated at the initial set temperature and pressure using the frog-jump Newton integration method with a step length of 2 fs and 100,000 steps.
Additionally, the ligand's binding energy (∆G b ) after docking with the acceptor (7CC4) was calculated using the MM-PBSA method (the binding energy is the sum of van der Waals energy, electrostatic energy, non-polar solvation energy, and polar solvation energy). In this article, a PS polymer (n = 5) was docked onto the 7CC4 receptor protein as a ligand-receptor complex structure for molecular dynamics (MD) simulations using 7CC4 as the receptor and blank set. Seven different types of plastic additives were added to the blank group as external stimuli and set up as an additive group [31]. These included common plasticizers, antioxidants, flame retardants, light stabilizers, heat stabilizers, lubricants, and fillers. The kinetics of the aforementioned blank and additive group complexes were simulated using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method in GROMACS 4.6.5 software using the Gromos96 43A1 force field, and the simulation results were expressed in terms of ∆G b [38,39]. ∆G b is used to quantify the extent to which a ligand binds to a receptor. The smaller the ∆G b , the more stable the complex formed between the ligand and the receptor [40]. The binding energy of PS polymers (n = 5) and their additive to additives to receptor proteins is used to characterize the biodegradability of microplastics in this article, with a lower ∆G b indicating greater biodegradability.

Proportional Factorial Experimental Design Method for Screening the Proportioning Scheme of Plastic Components Affecting Biodegradability
Plastic, as a polymer material, is generally composed of polymer monomers and a variety of additives that vary according to the product and performance requirements. Common plastic products are primarily composed of polyolefin, polyvinyl chloride, and polyethylene terephthalate as raw materials, supplemented by plasticizers (A) (Diethylhexyl phthalate, Diisononyl phthalate, etc.), antioxidant (B) (Nonylphenol, Acetone diphenylamine, etc.), light stabilizers(C) (ultraviolet absorber UV-326, UV-328, etc.), Heat stabilizer (D) (Zinc stearate, Calcium stearate, etc.), Flame retardant (E) (commonly used brominated flame retardants: Decabromodiphenyl ether, Tetra bromo bisphenol A, etc.), Lubricant (F) (Oleamide, Stearic acid, etc.) and Filler (G) (Calcium carbonate, Calcium sulfate) and other mixtures of exogenous additives [41].The common polystyrene (PS)found in the environment is used as the plastic matrix in this paper, and seven typical plastic components (Table 1) are used as the primary research objects [31,42] to assess the biodegradability of microplastics in the environment. We adjusted the composition of plastic additives and screened for the presence of plastic components with a high potential for biodegradation. Factorial design, or full factor experimental design, is a statistical method for obtaining a large amount of information and accurately estimating the size of each experimental factor's main effect and interaction effect at all levels among factors. As a result, this method is capable of rapidly determining the optimal level of experimental factors [37]. With the assistance of Minitab 2021 software, and using the seven types of chemical additives in plastics (Plasticizer (A), Antioxidant (B), Light stabilizer (C), Heat stabilizer (D), Flame retardant (E), Lubricant (F) and Filler (G)) as the seven factors of factorial experimental design, each factor selects two representative substances as the level (Low level 1; High level 2).The "Statistics-DOE (Design of Experiment)-Factor-Create Factor Design" module generates a full factorial design experiment table of 128 (2 7 ) plastic component combinations. Second, based on the optimal combination in Table 1, the six most critical additives (Plasticizer (A), Antioxidant (B), Light stabilizer (C), Heat stabilizer (D), Flame retardant (E), and Lubricant (F)) were used as six factors in the second round of factorial experimental design (Table 2), and each factor was added as a level (Low level 0, i.e., no addition; High level 1, i.e., addition), resulting64 (2 6 ). Additionally, based on the optimal additive composition formula (Table 1 No. 116), the response value was determined using biodegradability data quantified by dynamic binding energy. According to the factorial analysis method in the "Statistics-DOE (Design of Experiment)-Factor-Analysis Factor Design" module, which is capable of analyzing the interaction between various components of plastic additives and validating the underlying mechanism.

Biodegradability Features of Plastics Using the Entire Factorial Experimental Design Group's Distribution Proportion Scheme
According to the method described in 2.2, 128 groups of various plastic additive components were created, and the binding free energy values of PS and receptor degradation protein 7CC4 in the blank (the blank group refers to the state in which no additives are present, and only polystyrene 5 polymer ligands and degraded protein receptors are present) and additive groups were calculated using a molecular docking-assisted molecular dynamics simulation method (Table 3), in order to characterize the size of plastic biodegradability. As shown in Table 3, except for the binding energies of the two combinations of No. 35 and No. 62, most of the binding energies decreased significantly when compared to the blank group (−100.374 kJ/mol) (114 groups in total). The binding energy values of the additive components were significantly reduced (at a rate of 20%), indicating that the addition and ratio of additives aided in the biodegradation of PS plastics, which could significantly improve the biodegradability of microplastics, and also demonstrate the necessity of adding exogenous additives during the plastic manufacturing process [44].  Additionally, the binding energy of additive group No. 116 (−192.483 kJ/mol) is the lowest of all additive combinations and decreases by the greatest amount (91.77%) (The percentage of change rate refers to the change rate of binding energy in the presence of additives compared to the blank group (that is, the change rate of degradability = (binding energy value of additive group-binding energy value of blank group)/binding energy value of blank group × 100%)) when compared to the binding energy of the blank group (−100.374 kJ/mol), indicating the additive combination (Diisononyl phthalate, Acetone diphenylamine, Bometriazole, Zinc stearate, Tetrabromobisphenol A, Oleamide and Calcium sulfate). Thus, using combination No. 116 as an example, the following analysis was conducted to determine the effect of various additive components and ratios on the biodegradability of plastics.

Biodegradation Characteristics of Plastics under an Additive Combination Scheme Using a Full Factorial Design
Using a fractional factorial experimental design, the effects of plastic additives on the biodegradability of plastics were studied using the optimal biodegradability combination scheme No. 116 of polymers screened in Section 3.1.
According to Table 4, the binding energies of antioxidant (B) (Acetone diphenylamine), heat stabilizer (D) (Zinc stearate) and filler G (Calcium sulfate) are −146.232 kJ/mol, −143.404 kJ/mol and −144.624 kJ/mol, respectively. The decreased range is significantly less than that of the other four additives, indicating that these three additives have a less favorable effect on plastic's biodegradability than the other four (Plasticizer (A): Isopropyl phthalate, Light stabilizer (C): Bumetrazole, Flame retardant (E): Tetrabromobisphenol A and Lubricant (F): Oleamide). Additionally, filler (Calcium sulfate) has a smaller part in plastic production than other additives, and it is mostly employed to minimize manufacturing costs [45]. Thus, to more effectively analyze the effect of plastic additives on plastic biodegradability, six different types of plastic additives were chosen (A, B, C, D, E, and F), and a full factorial experimental design of six factors and two levels was used to generate 64 groups of plastic additive combination schemes (Table 5), in order to identify the main effects of plastic additives and their interaction effects on plastic biodegradability.

The Effect of Plastic Additives on the Biodegradability of Polymers Was Investigated
The contribution of plastic additives to polymer biodegradability may not be completely independent. The DOE-analysis factor design module of Minitab software was used in this work to determine the major influence and the kind and degree of the interaction effect of additive components on the biodegradability of plastics. The primary effects of plastic additives and their interaction effects on plastic degradation were investigated, as well as the contribution rate of each order interaction effect. The results are summarized in Table 6 and Figure 1. Among the interactions of plastic additive components contributing to the biodegradability of plastics, the third-order interaction accounts for the highest proportion, followed by the second-order interaction effect, indicating that the interaction between the components of plastic additives is complex, and multi-factor synergistic or antagonistic effects dominate.
Among the major effects, plasticizer (A) has the greatest synergistic effect when microplastics and degradation proteins are combined, while antioxidant (B) has the greatest inhibitory effect when microplastics and degradation proteins are combined. In the secondorder interaction effect, the coexistence of plasticizer (A) and lubricant (F) is most favorable for the combination of microplastics and degradation proteins, followed by the synergistic effect of plasticizer (A) and antioxidant (B) (this paper uses 10% as the significant standard); the combination of light stabilizer (C) and flame retardant (E) has the most significant antagonistic effect on the combination of microplastics and degradation proteins. Additionally, antioxidant (B) and heat stabilizers (D) work in opposition to the combination of plastic additives and breakdown proteins. completely independent. The DOE-analysis factor design module of Minitab software was used in this work to determine the major influence and the kind and degree of the interaction effect of additive components on the biodegradability of plastics. The primary effects of plastic additives and their interaction effects on plastic degradation were investigated, as well as the contribution rate of each order interaction effect. The results are summarized in Table 6 and Figure 1.  Among the interactions of plastic additive components contributing to the biodegradability of plastics, the third-order interaction accounts for the highest proportion, followed by the second-order interaction effect, indicating that the interaction between the components of plastic additives is complex, and multi-factor synergistic or antagonistic effects dominate.
Among the major effects, plasticizer (A) has the greatest synergistic effect when microplastics and degradation proteins are combined, while antioxidant (B) has the greatest inhibitory effect when microplastics and degradation proteins are combined. In the second-order interaction effect, the coexistence of plasticizer (A) and lubricant (F) is most favorable for the combination of microplastics and degradation proteins, followed The effect is strong, and its magnitude exceeds 20%. As a result, it is not recommended to employ antioxidant (B) and heat stabilizer (D) as plastic additives concurrently to avoid complicating plastic biodegradation. In the third-order interaction effect, when antioxidant (B), light stabilizer (C), and heat stabilizer (D) are combined; the combination of plastic additives and degradation proteins is considerably increased, as is the combination of antioxidant (B), flame retardant (E) and lubricant (F) can be significantly increased.
Additionally, the combination of plasticizer (A), light stabilizer (C), and flame retardant (E) has a considerable synergistic impact when plastic additives and degradation proteins are combined. The combination of plasticizer (A), antioxidant (B), and flame retardant (E) is, in some ways, incompatible with the combination of plastic additives and breakdown proteins, but the effect is negligible.
The DOE module in Minitab was used to check the consistency of the degradation effect data for 64 groups of microplastics ((C 8 H 8 ) n, n = 5 ) additive combinations. Figure 2 depicts the results. All the schemes in this work follow the normal distribution trend, as seen in Figure 2. The residual distribution histogram trend is to increase first, then decline. The residual plots are dotted, indicating that the results are accurate. There is no aberrant change trend. The results are not related to the experimental order in the graph of residuals and experimental order, indicating that the experimental order has no effects on the results, demonstrating that the results are independent.

Molecular Dynamics Simulation-Based Mechanism Investigation of the Biodegradation Effect of Polystyrene 5 Polymers
The dominance analysis is primarily a statistical tool for determining the interdependence of multiple variables. As a result, this work employs the dominance analysis method to investigate the primary influencing elements of 64 different plastic additive combinations on the degradation of the polystyrene 5 polymer via protein degradation (7CC4). The van der Waals energy, electrostatic energy, polar solvation energy and non-polar solvation energy make up the binding free energy of ligand-receptor docking [46]. The lower the free energy, the better the biodegradability of the polystyrene 5 polymer. The correlation between the components of binding free energy of the polystyrene 5 polymer and receptor 7CC4 and the change value of binding free energy was analyzed under 64 additive combinations, with no external additive combination (No. 116-16) as the invariant, and which energies have a significant correlation with the binding free energy was analyzed.
the combination of antioxidant (B), flame retardant (E) and lubricant (F) can be significantly increased.
Additionally, the combination of plasticizer (A), light stabilizer (C), and flame retardant (E) has a considerable synergistic impact when plastic additives and degradation proteins are combined. The combination of plasticizer (A), antioxidant (B), and flame retardant (E) is, in some ways, incompatible with the combination of plastic additives and breakdown proteins, but the effect is negligible.
The DOE module in Minitab was used to check the consistency of the degradation effect data for 64 groups of microplastics ((C8H8)n, n = 5) additive combinations. Figure 2 depicts the results. All the schemes in this work follow the normal distribution trend, as seen in Figure 2. The residual distribution histogram trend is to increase first, then decline. The residual plots are dotted, indicating that the results are accurate. There is no aberrant change trend. The results are not related to the experimental order in the graph of residuals and experimental order, indicating that the experimental order has no effects on the results, demonstrating that the results are independent.   Figure 3 shows the 95% confidence interval in grey, with different colors representing distinct components: the bigger the correlation coefficient r, the stronger the association [47]. The correlation coefficients r of the linear regression equations between the change values of the combined free energy (the difference compared to the blank group No. 116-16) of the 64 combination schemes and van der Waals energy, non-polar solvation energy, polar solvation energy, and electrostatic energy, respectively, were 0.8185, 0.7483, 0.1010 and 0.0608 (r α = 0.2461, p = 0.05, n = 63, r > r α was significantly correlated). Furthermore, the change range of van der Waals energy is far wider than the jump range of non-polar solvation energy (SASA energy) under different additive combinations, both of which are negative, but the latter is smaller in absolute value. The solvation energy of the system is all positive, and its absolute value is larger than the absolute value of the SASA energy, indicating that the solvation effect is not conducive to the combination of the polystyrene 5 polymer and biodegrading enzymes, and the van der Waals energy is certain. Furthermore, there is a significant beneficial effect on the biodegradability of polystyrene 5 polymer. In addition, it has been found that the binding cavity of proteases is hydrophobic [47,48]. Moreover, combined with the factorial analysis results in Section 3.2, it can be found that the chemical structure of the additive with significant promoting effect mainly exhibits hydrophobic properties. Therefore, when the additive is used as an external stimulus, the additive will promote the further inward folding of the degraded protein, so that the binding of the degraded protein to the polystyrene pentamer is enhanced.

Total Factor and Amino Acid Residue Analysis Were Used to Investigate the Mechanism of Biodegradation of Polystyrene 5 Polymers
The findings of factorial analysis, when combined with Figures 1 and 4, revealed that the effect of interaction between different additives on the biodegradation of the polystyrene 5 polymer was distinct. The absolute magnitude of the biodegradation standardization effect of the polystyrene 5 polymer is depicted in the Pareto diagram ( Figure 4). (From the maximum to the minimum effect). It is statistically significant (α = 0.05) when the factor bar passes the reference line 2.074. As a result, the impact of additive A and its combination B*C*D, C*E, B*E*F, and A*C*E on the polystyrene 5 polymer's biodegrading performance was statistically significant (Figure 4a), and the addition of B*C*D, A, B*E*F, and A*C*E has a substantial favorable effect. That is, it promotes the breakdown of the polystyrene 5 polymer (Figure 4b).

Total Factor and Amino Acid Residue Analysis Were Used to Investigate the Mechanism of Biodegradation of Polystyrene 5 Polymers
The findings of factorial analysis, when combined with Figures 1 and 4, revealed that the effect of interaction between different additives on the biodegradation of the polystyrene 5 polymer was distinct. The absolute magnitude of the biodegradation standardization effect of the polystyrene 5 polymer is depicted in the Pareto diagram ( Figure 4). (From the maximum to the minimum effect). It is statistically significant (α = 0.05) when the factor bar passes the reference line 2.074. As a result, the impact of additive A and its combination B*C*D, C*E, B*E*F, and A*C*E on the polystyrene 5 polymer's biodegrading performance was statistically significant (Figure 4a), and the addition of B*C*D, A, B*E*F, and A*C*E has a substantial favorable effect. That is, it promotes the breakdown of the polystyrene 5 polymer (Figure 4b). The amino acid residues analysis of the interaction between the polystyrene 5 polymer and degradation enzyme 7CC4 is displayed in Figure 5 under the condition of an external addition with a clear favorable effect (where the green circle represented the amino acid residues involved in van der Waals interaction [49]).The number of amino acid residues around the polystyrene 5 polymer increased after adding external conditions conducive to polystyrene 5 degradation, indicating that van der Waals interaction was the main driving force for the binding of the polystyrene 5 polymer and degradation protein 7CC4. This conclusion follows the preceding one. Hydrogen bonds and hydrophobic interactions are critical in the binding of ligands to proteases, according to previous research [50]. Table 7 shows that in the absence of external additives , only four hydrophobic interactions were formed between the polystyrene 5 polymer and degradation protein 7CC4, whereas a combination of external additives  formed seven, five, five and eight hydrophobic interactions, respectively, with only a small amount of hydrogen bonds, indicating that the more hydrophobic interactions formed, the more conducive to the combination of the polystyrene 5 polymer and degradation protein 7CC4. As a result, the reason the degradation effect of polystyrene 5 polymers is stronger in the group with favorable external conditions than in the No. 116-16 group may be because the former generated more hydrophobic contacts. Combined with the factorial analysis results in Section 3.2, the No. 116-45 group contains the plasticizer with the highest proportion of plastic additives, and the factorial result of this additive combination shows a significant positive effect; the polystyrene 5 polymer formed more hydrophobic interactions with the degradation protein 7CC4, compared to the other three positive effect groups. As a result, the No. 116-45 group was recommended as the best combination scheme for the biodegradation of polystyrene 5 polymers in the study. The amino acid residues analysis of the interaction between the polystyrene 5 polymer and degradation enzyme 7CC4 is displayed in Figure 5 under the condition of an external addition with a clear favorable effect (where the green circle represented the amino acid residues involved in van der Waals interaction [49]).The number of amino acid residues around the polystyrene 5 polymer increased after adding external conditions conducive to polystyrene 5 degradation, indicating that van der Waals interaction was the main driving force for the binding of the polystyrene 5 polymer and degradation protein 7CC4. This conclusion follows the preceding one. Hydrogen bonds and hydrophobic interactions are critical in the binding of ligands to proteases, according to previous research [50]. Table 7 shows that in the absence of external additives , only four hydrophobic interactions were formed between the polystyrene 5 polymer and degradation protein 7CC4, whereas a combination of external additives  formed seven, five, five and eight hydrophobic interactions, respectively, with only a small amount of hydrogen bonds, indicating that the more hydrophobic interactions formed, the more conducive to the combination of the polystyrene 5 polymer and degradation protein 7CC4. As a result, the reason the degradation effect of polystyrene 5 polymers is stronger in the group with favorable external conditions than in the No. 116-16 group may be because the former generated more hydrophobic contacts. Combined with the factorial analysis results in Section 3.2, the No. 116-45 group contains the plasticizer with the highest proportion of plastic additives, and the factorial result of this additive combination shows a significant positive effect; the polystyrene 5 polymer formed more hydrophobic interactions with the degradation protein 7CC4, compared to the other three positive effect groups. As a result, the No. 116-45 group was recommended as the best combination scheme for the biodegradation of polystyrene 5 polymers in the study.  Zhang et al. [51] developed five environmentally friendly PAEs derivatives (DBP-CHO, DBP-COOH, DINP-NH2, DINP-NO2) using a normalization method, a pharmacophore model, and molecular docking technology. They discovered that DBP-CHO had a high biodegradability, was less toxic, and more easily degraded by light in the natural environment.
To further validate the rationality and universality of the design theory and factorial experiment design used in this study, the additive combinations with plasticizer A with an effective contribution percentage greater than 10% in the factorial results were screened and verified.   Zhang et al. [51] developed five environmentally friendly PAEs derivatives (DBP-CHO, DBP-COOH, DINP-NH 2 , DINP-NO 2 ) using a normalization method, a pharmacophore model, and molecular docking technology. They discovered that DBP-CHO had a high biodegradability, was less toxic, and more easily degraded by light in the natural environment.
To further validate the rationality and universality of the design theory and factorial experiment design used in this study, the additive combinations with plasticizer A with an effective contribution percentage greater than 10% in the factorial results were screened and verified. The plasticizers (A) in the groups  were replaced with the above five derivatives, and the replaced plasticizer was used as a new external condition (Table 8). As indicated in Figure 1, all four combinations had a positive effect value, and the contribution percentage was as follows: No. 116-28 > No. 116-26 > No. 116-45 ≈ No. 116-40. The analysis and calculation of the binding free energy data revealed that only after the DBP-CHO derivative molecule was replaced by A was the change trend of the binding free energy of the polystyrene 5 polymer and 7CC4 protease in the four groups consistent with the change trend in the contribution percentage of the effect value (the binding free energy values of the four groups were greater than those of No. 116-16 group). The screening results were consistent with those of Zhang et al., indicating that DBP-CHO had a favorable effect on the biodegradability of polystyrene 5, indicating the rationality, universality, and reliability of the design theory, the factorial experiment design, and the factorial analysis results in this study.

Conclusions
In this research, a biodegradation effect regulatory system based on the optimization of the number of plastic additives was created using factorial design and associated analysis methods supplemented by molecular docking and molecular dynamics simulation approaches. The investigation of relevant mechanisms revealed that van der Waals contact was highly associated with plastic's biodegradability. Simultaneously, it was discovered that hydrophobic interaction was the primary driving force behind the biodegradation of plastics, guiding the efficient biodegradation of plastics under the condition of a reasonable additive component combination.