Single chain dynamics in polymer networks: A Monte Carlo study

Molecular modeling to predict peptide accessibility for peptide-functionalized hydrogels

Peptide-functionalized (PF) hydrogels are being widely investigated by the tissue engineering and regenerative medicine communities for a broad range of applications because of their unique potential to mimic the natural extracellular matrix and promote tissue regeneration. In order for these complex material systems to perform their intended bioactive function (e.g., cell signaling), the peptides that are tethered to the hydrogel matrix must be accessible at the hydrogel surface for cell-receptor binding. The factors influencing the surface accessibility of the tethered peptide mainly include the length of the tethers, the loading (i.e., concentration) of the peptide, and the association between the tethered peptide and the hydrogel matrix. In the present work, the authors developed coarse-grained molecular models based on the all-atom polymer consistent force field for a type of poly(ethylene glycol)-based PF hydrogel and conducted molecular simulations to investigate the distribution of the peptide within the hydrogel and its surface accessibility as a function of tether length and peptide concentration. The calculated results of the effects of these design parameters on the surface accessibility of the peptide agree very well with corresponding experimental measurements in which peptide accessibility was quantified by the number of cells attached to the hydrogel surface per unit area. The developed modeling methods are able to provide unique insights into the molecular behavior of PF hydrogels and the distribution of the tethered peptides, which can serve as a guide for hydrogel design optimization.

Wavelet Monte Carlo dynamics: A new algorithm for simulating the hydrodynamics of interacting Brownian particles

We develop a new algorithm for the Brownian dynamics of soft matter systems that evolves time by spatially correlated Monte Carlo moves. The algorithm uses vector wavelets as its basic moves and produces hydrodynamics in the low Reynolds number regime propagated according to the Oseen tensor. When small moves are removed, the correlations closely approximate the Rotne-Prager tensor, itself widely used to correct for deficiencies in Oseen. We also include plane wave moves to provide the longest range correlations, which we detail for both infinite and periodic systems. The computational cost of the algorithm scales competitively with the number of particles simulated, N, scaling as N In N in homogeneous systems and as N in dilute systems. In comparisons to established lattice Boltzmann and Brownian dynamics algorithms, the wavelet method was found to be only a factor of order 1 times more expensive than the cheaper lattice Boltzmann algorithm in marginally semi-dilute simulations, while it is significantly faster than both algorithms at large N in dilute simulations. We also validate the algorithm by checking that it reproduces the correct dynamics and equilibrium properties of simple single polymer systems, as well as verifying the effect of periodicity on the mobility tensor.

Multiscale approach for the construction of equilibrated all-atom models of a poly(ethylene glycol)-based hydrogel

A multiscale modeling approach is presented for the efficient construction of an equilibrated all-atom model of a cross-linked poly(ethylene glycol) (PEG)-based hydrogel using the all-atom polymer consistent force field (PCFF). The final equilibrated all-atom model was built with a systematic simulation toolset consisting of three consecutive parts: (1) building a global cross-linked PEG-chain network at experimentally determined cross-link density using an on-lattice Monte Carlo method based on the bond fluctuation model, (2) recovering the local molecular structure of the network by transitioning from the lattice model to an off-lattice coarse-grained (CG) model parameterized from PCFF, followed by equilibration using high performance molecular dynamics methods, and (3) recovering the atomistic structure of the network by reverse mapping from the equilibrated CG structure, hydrating the structure with explicitly represented water, followed by final equilibration using PCFF parameterization. The developed three-stage modeling approach has application to a wide range of other complex macromolecular hydrogel systems, including the integration of peptide, protein, and/or drug molecules as side-chains within the hydrogel network for the incorporation of bioactivity for tissue engineering, regenerative medicine, and drug delivery applications.