By combining large-scale dissipative particle dynamics and steered molecular dynamics simulations, we investigate the mechanochemical cellular internalization pathways of homogeneous and heterogeneous nanohydrogels and demonstrate that membrane internalization is determined by the crosslink density and encapsulation ability of nanohydrogels. The homogeneous nanohydrogels with a high crosslink density and low encapsulation ability behave as soft nanoparticles partially wrapped by the membrane, while those with a low crosslink density and high encapsulation ability permeate into the membrane. Regardless of the crosslink density, the homogeneous nanohydrogels undergo typical dual morphological deformations. The local lipid nanodomains are identified at the contacting region between the membrane and nanohydrogels because of different diffusion behaviors between lipid and receptor molecules during the internalization process. The yolk@shell heterogeneous nanohydrogels present a different mechanochemical cellular internalization pathway. The yolk with strong affinity is directly in contact with the membrane, resulting in partial membrane wrapping, and the contacting area is much reduced when compared to homogenous nanohydrogels, leading to a smaller lipid nanodomain and thus avoiding related cellular toxicity. Our findings provide a critical mechanism understanding of the biological pathways of nanohydrogels and may guide the molecular design of the hydrogel-based materials for controlled release drug delivery, tissue engineering, and cell culture.
Keywords: DPD simulations; membrane internalization; molecular encapsulation; nanohydrogel; steered molecular dynamics.