As a common substance in production and life, phthalic acid esters (PAEs), the main component of plastics, have brought more and more serious problems to the environment. This study normalized the insulation, toxicity, and bioconcentration data of 13 PAEs to eliminate the dimensional coefficients of each index, and then used the comprehensive index method to calculate the comprehensive effect value of PAEs with three properties. The comprehensive effect value was used as the data source to construct the 3D-QSAR model of PAE molecular comprehensive effect. The DAP was selected as the target molecule, the distribution of each force field in the three-dimensional equipotential map was analyzed, and 30 molecular modification schemes were created. The constructed single-effect models of insulation, toxicity, and bioconcentration of PAEs and the scoring function module of DS software were used to evaluate the stability and environmental friendliness of PAE derivative molecules. Four PAE derivatives were screened for increased comprehensive effects, enhanced insulation, and reduced toxicity and bioconcentration. By calculating the binding energy of the target molecule and the derivative molecule with the degrading enzyme under different applied electric fields, it was found that the binding energy of DAP-1-NO2-2-CH2C6H5 decreases more than DAP does when there is an applied electric field, indicating that the degradation ability of degrading enzymes on PAE derivative molecules is reduced, which indirectly proves that the insulation is enhanced. The innovation of this paper lies in the insulation, toxicity, and bioenrichment data of PAEs being processed by mathematical method for the first time, and PAEs with high insulation, low toxicity, and low bioconcentration were designed by building a comprehensive model.