The Construction and Validation of All-Atom Bulk-Phase Models of Amorphous Polymers Using the TIGER2/TIGER3 Empirical Sampling Method

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.

Perspectives on the simulation of protein-surface interactions using empirical force field methods

Protein-surface interactions are of fundamental importance for a broad range of applications in the fields of biomaterials and biotechnology. Present experimental methods are limited in their ability to provide a comprehensive depiction of these interactions at the atomistic level. In contrast, empirical force field based simulation methods inherently provide the ability to predict and visualize protein-surface interactions with full atomistic detail. These methods, however, must be carefully developed, validated, and properly applied before confidence can be placed in results from the simulations. In this perspectives paper, I provide an overview of the critical aspects that I consider being of greatest importance for the development of these methods, with a focus on the research that my combined experimental and molecular simulation groups have conducted over the past decade to address these issues. These critical issues include the tuning of interfacial force field parameters to accurately represent the thermodynamics of interfacial behavior, adequate sampling of these types of complex molecular systems to generate results that can be comparable with experimental data, and the generation of experimental data that can be used for simulation results evaluation and validation.

The Structure of Hydrated Poly (D, L - Lactic Acid) Studied With X-Ray Diffraction and Molecular Simulation Methods

The effect of hydration on the molecular structure of amorphous poly (D, L-lactic acid) (PDLLA) with 50:50 L-to-D ratio has been studied by combining experiments with molecular simulations. X-ray diffraction measurements revealed significant changes upon hydration in the structure functions of the copolymer. Large changes in the structure functions at ~ 10 days of incubation coincided with the large increase in the water uptake from ~1 to ~40% and the formation of voids in the film. Computer modeling based on the recently developed TIGER2/TIGER3 mixed sampling scheme was used to interpret these changes by efficiently equilibrating both dry and hydrated models of PDLLA. Realistic models of bulk amorphous PDLLA structure were generated as demonstrated by close agreement between the calculated and the experimental structure functions. These molecular simulations were used to identify the interactions between water and the polymer at the atomic level including the change of positional order between atoms in the polymer due to hydration. Changes in the partial O-O structure functions, about 95% of which were due to water-polymer interactions, were apparent in the radial distribution functions. These changes, and somewhat smaller changes in the C-C and C-O partial structure functions, clearly demonstrated the ability of the model to capture the hydrogen bonding interactions between water and the polymer, with the probability of water forming hydrogen bonds with the carbonyl oxygen of the ester group being about four times higher than with its ether oxygen.