Chain-Length-Dependent Interactions of Polyethylene Nanoplastics with Human Catalase reveals Size based Nanoplastic Toxicity

Abstract

Nanoplastic contamination is increasingly recognized as a biochemical stressor, yet the molecular determinants governing how polymer fragment size dictates protein-level interactions remain poorly understood. Here, we combine molecular docking and all-atom molecular dynamics simulations to characterize the interaction of polyethylene (PE) oligomers of increasing length (n = 3, 5, 7, and 14 repeat units) with human catalase, a tetrameric heme enzyme central to cellular antioxidant defense. Single-chain simulations reveal that short PE fragments (3-, 5-, and 7-mer) migrate from an initial domain-connector pose toward a heme-proximal, previously characterized inhibitor-binding site, whereas the 14-mer remains kinetically confined to a peripheral surface site, implicating short-chain PE as a potential catalase inhibitor. In multichain corona systems, backbone dynamics show a non-monotonic, chain-length-dependent destabilization, with the 3- and 5-mer inducing the greatest structural perturbation and the 14-mer best preserving native-like dynamics despite forming an extensive, percolating surface network across all four subunits. MM-PBSA free energy decomposition identifies van der Waals interactions as the dominant thermodynamic driver of binding, while electrostatic surface mapping indicates that cationic surface patches instead govern long-range spatial targeting of PE chains. These findings establish a framework linking nanoplastic fragment size to protein-level toxicity, with direct implications for assessing the biochemical risk of environmentally weathered polyethylene.