Molecular Dynamics Simulation of Hydrogels Based on Phosphorylcholine-Containing Copolymers for Soft Contact Lens Applications

The structure and dynamics of copolymers of 2-hydroxyethyl methacrylate (HEMA) with 2-methacryloyloxyethyl phosphorylcholine (MPC) were studied by molecular dynamics simulations. In total, 20 systems were analyzed. They differed in numerical fractions of the MPC in the copolymer chain, equal to 0.26 and 0.74, in the sequence of mers, block and random, and the water content, from 0 to 60% by mass. HEMA side chains proved relatively rigid and stable in all considered configurations. MPC side chains, in contrast, were mobile and flexible. Water substantially influenced their dynamics. The copolymer swelling caused by water resulted in diffusion channels, pronounced in highly hydrated systems. Water in the hydrates existed in two states: those that bond to the polymer chain and the free one; the latter was similar to bulk water but with a lower self-diffusion coefficient. The results proved that molecular dynamics simulations could facilitate the preliminary selection of the polymer materials for specific purposes before their synthesis.