Phthalic acid esters (PAEs) are typical environmental endocrine disrupters, interfering with the endocrine system of organisms at very low concentrations. The plasma membrane is the first barrier for organic pollutants to enter the organism, so membrane permeability is a key factor affecting their biological toxicity. In this study, based on computational approaches, we investigated the permeation and intramembrane aggregation of typical PAEs (dimethyl phthalate, DMP; dibutyl phthalate, DBP; di-2-ethyl hexyl phthalate, DEHP), as well as their effects on membrane properties, and related molecular mechanisms were uncovered. Our results suggested that PAEs could enter the membrane spontaneously, preferring the headgroup-acyl chain interface of the bilayer, and the longer the side chain (DEHP > DBP > DMP), the deeper the insertion. Compared with the shortest DMP, DEHP apparently increased membrane thickness, order, and rigidity, which might be due to its stronger hydrophobicity. Potential of means force (PMF) analysis revealed the presence of an energy barrier located at the water-membrane interface, with a maximum value of 2.14 kcal mol⁻¹ obtained in the DEHP-system. Therefore, the difficulty of membrane insertion is also positively correlated with the side-chain length or hydrophobicity of PAE molecules. These findings will inspire our understanding of structure-activity relationship between PAEs and their effects on membrane properties, and provide a scientific basis for the formulation of environmental pollution standards and the prevention and control of small molecule pollutants.

Phthalate esters (PAEs) are a ubiquitous kind of environmental endocrine that disrupt chemicals, causing environmental and health issues. EstJ6 is an effective phthalate-degrading hydrolase, and its mutant with a combination of three non-conservative distal mutations has an improved activity against PAEs with unknown molecular mechanisms. Herein, we attempt to fill the significant gap between distal mutations and the activity of this enzyme using computational approaches. We found that mutations resulted in a redistribution of the enzyme’s preexisting conformational states and dynamic changes of key functional regions, especially the lid over the active site. The outward motion of the lid upon the mutations made it easier for substrates or products to enter or exit. Additionally, a stronger substrate binding affinity and conformational rearrangements of catalytic reaction-associated residues in the mutant, accompanied by the strengthened communication within the protein, could synergistically contribute to the elevated catalytic efficiency. Finally, an attempt was made to improve the thermostability of EstJ6 upon introducing a distal disulfide bond between residues A23 and A29, and the simulation results were as expected. Together, our work explored the allosteric effects caused by distal mutations, which could provide insights into the rational design of esterases for industrial applications in the future.