An Insight into Structural and Mechanical Properties of Ideal‐Networked Poly(Ethylene Glycol)–Peptide Hydrogels from Molecular Dynamics Simulations

Effect of the degree of polymerization and water content on the thermal transport phenomena in PEGDA hydrogel: a molecular-dynamics-based study

A hydrogel is a 3D cross-linked polymer network that can absorb copious amounts of water or biological fluid. Due to their biocompatibility and non-toxicity, hydrogels have a wide range of applications in biomedical engineering. To develop hydrogels with superior thermal dissipation properties, atomistic-level studies are required to quantify the effect of the water content and the degree of polymerization. Classical mechanics-based non-equilibrium molecular dynamics (NEMD) simulations were performed in conjunction with a mathematical formulation developed by Müller-Plathe to explore the thermal conductivity of the poly(ethylene glycol)diacrylate (PEGDA) hydrogel. This work reveals that the thermal conductivity of the PEGDA hydrogel is enhanced with the increase in water content and approaches the value of the thermal conductivity of water at 85% water content in the hydrogel. The PEGDA-9 hydrogel, with a lower level of degree of polymerization, has a superior thermal conductivity than the PEGDA-13 and PEGDA-23 hydrogels. The lower level of degree of polymerization is associated with the higher mesh density of polymer chain network junctions that help to achieve the superior thermal conductivity at higher water contents. Increasing the water content improves the structural stability and compactness of the polymer chains, which can be further associated with the enhanced phonon transfer in PEGDA hydrogels. The work will help in the development of PEGDA-based hydrogels with superior thermal dissipation properties for tissue engineering.

PEGDA hydrogel structure from semi-dilute concentrations: insights from experiments and molecular simulations

The efficacy of hydrogel materials used in biomedical applications is dependent on polymer network topology and the structure of water-laden pore space. Hydrogel microstructure can be tuned by adjusting synthesis parameters such as macromer molar mass and concentration. Moreover, hydrogels beyond dilute conditions are needed to produce mechanically robust and dense networks for tissue engineering and/or drug delivery systems. Thus, this study utilizes a combined experimental and molecular simulation approach to characterize structural features for 4.8 and 10 kDa poly (ethylene glycol) diacrylate (PEGDA) hydrogels formed from a range of semi-dilute solution concentrations. The connection between chain-chain interactions in polymer solutions, hydrogel structure, and equilibrium swelling behavior is presented. Bulk rheology analysis revealed an entanglement concentration for PEGDA pre-gel solutions around 28 wt% for both macromers studied. A similar transition in swelling behavior was revealed around the same concentration where hydrogel capacity to retain water was drastically reduced. To understand this transition, the hydrogel structure was characterized using the swollen polymer network hypothesis and compared to pore size distributions from molecular dynamics simulations. We find in both approaches a structural transition concentration at the hydrogel swelling inflection point that is comparable to the entanglement concentration. Calculated mesh sizes from theory are compared with computationally determined average maximum pore diameters; mesh sizes from theory yielded greater feature sizes across all concentrations considered. Molecular simulations are further used to assess pore dynamics, which are shown to vary in distribution shape and number of modes compared to the time-averaged hydrogel pore features. Altogether, this work provides insights into hydrogel network features and their dynamic behavior at physiological conditions (37 °C) as a basis for hydrogel design beyond dilute conditions for biomedical applications.

Not only in silico drug discovery: Molecular modeling towards in silico drug delivery formulations

The use of methods at molecular scale for the discovery of new potential active ligands, as well as previously unknown binding sites for target proteins, is now an established reality. Literature offers many successful stories of active compounds developed starting from insights obtained in silico and approved by Food and Drug Administration (FDA). One of the most famous examples is raltegravir, a HIV integrase inhibitor, which was developed after the discovery of a previously unknown transient binding area thanks to molecular dynamics simulations. Molecular simulations have the potential to also improve the design and engineering of drug delivery devices, which are still largely based on fundamental conservation equations. Although they can highlight the dominant release mechanism and quantitatively link the release rate to design parameters (size, drug loading, et cetera), their spatial resolution does not allow to fully capture how phenomena at molecular scale influence system behavior. In this scenario, the "computational microscope" offered by simulations at atomic scale can shed light on the impact of molecular interactions on crucial parameters such as release rate and the response of the drug delivery device to external stimuli, providing insights that are difficult or impossible to obtain experimentally. Moreover, the new paradigm brought by nanomedicine further underlined the importance of such computational microscope to study the interactions between nanoparticles and biological components with an unprecedented level of detail. Such knowledge is a fundamental pillar to perform device engineering and to achieve efficient and safe formulations. After a brief theoretical background, this review aims at discussing the potential of molecular simulations for the rational design of drug delivery systems.

On shockwave propagation and attenuation in poly(ethylene glycol) diacrylate hydrogels

An analytical model is developed to predict shockwave propagation and attenuation in hydrogels by combining the classical method of shock characteristics and a solution for the shock front structure. To guide the development of the model, molecular dynamics (MD) simulations are performed. Specifically, a one-dimensional shock pulse in poly(ethylene glycol) diacrylate (PEGDA) hydrogels is simulated with the nonequilibrium MD method. The role of polymer concentration on the shock response is evaluated by constructing hydrogels with 20, 35, and 50 wt% PEGDA concentrations in an idealized crosslinked network. Steady-state pressure-density and shock-particle velocity relationships are established using the Murnaghan equation of state. Shock front structure is characterized by a power-law equation that relates the shock front thickness with shock pressure. These results are used as critical input for the shock propagation and attenuation model. The model is then evaluated via comparison with the classical method of characteristics. It shows significant improvement in accuracy and successfully captures salient features of shockwave attenuation, including the shock pressure amplitude, the velocities of the shock and release waves, and the attenuation timeline. Hydrogels with higher polymer concentrations exhibit a shorter attenuation time at all particle velocities studied. This behavior is attributed to differences in bulk properties and shock front structure in hydrogels with different polymer/water concentrations.

Molecular Simulations of PEGylated Biomolecules, Liposomes, and Nanoparticles for Drug Delivery Applications

Since the first polyethylene glycol (PEG)ylated protein was approved by the FDA in 1990, PEGylation has been successfully applied to develop drug delivery systems through experiments, but these experimental results are not always easy to interpret at the atomic level because of the limited resolution of experimental techniques. To determine the optimal size, structure, and density of PEG for drug delivery, the structure and dynamics of PEGylated drug carriers need to be understood close to the atomic scale, as can be done using molecular dynamics simulations, assuming that these simulations can be validated by successful comparisons to experiments. Starting with the development of all-atom and coarse-grained PEG models in 1990s, PEGylated drug carriers have been widely simulated. In particular, recent advances in computer performance and simulation methodologies have allowed for molecular simulations of large complexes of PEGylated drug carriers interacting with other molecules such as anticancer drugs, plasma proteins, membranes, and receptors, which makes it possible to interpret experimental observations at a nearly atomistic resolution, as well as help in the rational design of drug delivery systems for applications in nanomedicine. Here, simulation studies on the following PEGylated drug topics will be reviewed: proteins and peptides, liposomes, and nanoparticles such as dendrimers and carbon nanotubes.