Mechanical and Transport Properties of the Poly(ethylene oxide)−Poly(acrylic acid) Double Network Hydrogel from Molecular Dynamic Simulations

Efficient characterization of double-cross-linked networks in hydrogels using data-inspired coarse-grained molecular dynamics model

The network structure within polymers significantly influences their mechanical properties, including their strength, toughness, and fatigue resistance. All-atom molecular dynamics (AAMD) simulations offer a method to investigate the energy dissipation mechanism within polymers during deformation and fracture; Such an approach is, however, computationally inefficient when used to analyze polymers with complex network structures, such as the common chemically double-networked hydrogels. Alternatively, coarse-grained molecular dynamics (CGMD) models, which reduce the computational degrees of freedom by concentrating a set of adjacent atoms into a coarse-grained bead, can be employed. In CGMD simulations, a coarse-grained force field (CGFF) is a critical factor affecting the simulation accuracy. In this paper, we proposed a data-based method for predicting the CGFF parameters to improve the simulation efficiency of complex cross-linked network in polymers. Here, we utilized a typical chemically double-networked hydrogel as an example. An artificial neural network was selected, and it was trained with the tensile stress-strain data from the CGMD simulations using different CGFF parameters. The CGMD simulations using the predicted CGFF parameters show good agreement with the AAMD simulations and are almost fifty times faster. The data-inspired CGMD model presented here broadens the applicability of molecular dynamics simulations to cross-linked polymers and has the potential to provide insights that will aid the design of polymers with desirable mechanical properties.

Polysaccharide-Based Double-Network Hydrogels: Polysaccharide Effect, Strengthening Mechanisms, and Applications

Polysaccharides are carbohydrate polymers that are major components of plants, animals, and microorganisms, with unique properties. Biological hydrogels are polymeric networks that imbibe and retain large amounts of water and are the major components of living organisms. The mechanical properties of hydrogels are critical for their functionality and applications. Since synthetic polymeric double-network (DN) hydrogels possess unique network structures with high and tunable mechanical properties, many natural functional polysaccharides have attracted increased attention due to their rich and convenient sources, unique chemical structure and chain conformation, inherently desirable cytocompatibility, biodegradability and environmental friendliness, diverse bioactivities, and rheological properties, which rationally make them prominent constituents in designing various strong and tough polysaccharide-based DN hydrogels over the past ten years. This review focuses on the latest developments of polysaccharide-based DN hydrogels to comprehend the relationship among the polysaccharide properties, inner strengthening mechanisms, and applications. The aim of this review is to provide an insightful mechanical interpretation of the design strategy of novel polysaccharide-based DN hydrogels and their applications by introducing the correlation between performance and composition. The mechanical behavior of DN hydrogels and the roles of varieties of marine, microbial, plant, and animal polysaccharides are emphatically explained.

Mechanical deformation affects the counterion condensation in highly-swollen polyelectrolyte hydrogels

Polyelectrolyte gels can generate electric potentials under mechanical deformation. While the underlying mechanism of such a response is often attributed to changes in counterion-condensation levels or alterations in the ionic conditions in the pervaded volume of the hydrogel, the exact molecular origins are largely unknown. By using all-atom molecular dynamics simulations of a polyacrylic acid hydrogel in explicit water as a model system, we simulate the uniaxial compression and uniaxial stretching of weakly to highly swollen (i.e., between 60-90% solvent content) hydrogel networks and calculate the microscopic condensation levels of counterions around the hydrogel chains. The counterion condensation under deformation is highly non-monotonic. Ionic condensation around the constituting chains of the deformed hydrogel tends to increase as the chains are stretched. This increase reaches a maximum and decreases as the chains are strongly stretched. The condensation around the collapsed chains of the hydrogel is weakly affected by the deformation. As a result, both compressing and stretching the model hydrogel lead to an overall increase in the counterion condensation. The effect vanishes for weakly swollen hydrogels, for which most ions are already condensed. The simulations with single, stretched polyelectrolyte chains show a qualitatively similar response, suggesting the effect of chain elongation on the ionic distribution throughout the hydrogel. Notably, this deformation-induced counterion condensation phenomenon does not occur in a polyelectrolyte solution at its critical concentration, indicating the role of hydrogel topology constraining the chain ends. Our results indicate that counterion condensation in a deforming polyelectrolyte hydrogel can be highly heterogeneous and exhibit a rich behaviour of electrostatic responses.

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.

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.

Lecithin as an Effective Modifier of the Transport Properties of Variously Crosslinked Hydrogels

Transport properties are one of the most crucial assets of hydrogel samples, influencing their main application potential, i.e., as drug carriers. Depending on the type of drug or the application itself, it is very important to be able to control these transport properties in an appropriate manner. This study seeks to modify these properties by adding amphiphiles, specifically lecithin. Through its self-assembly, lecithin modifies the inner structure of the hydrogel, which affects its properties, especially the transport ones. In the proposed paper, these properties are studied mainly using various probes (organic dyes) to effectively simulate drugs in simple release diffusion experiments controlled by UV-Vis spectrophotometry. Scanning electron microscopy was used to help characterize the diffusion systems. The effects of lecithin and its concentrations, as well as the effects of variously charged model drugs, were discussed. Lecithin decreases the values of the diffusion coefficient independently of the dye used and the type of crosslinking. The ability to influence transport properties is better observed in xerogel samples. The results, complementing previously published conclusions, showed that lecithin can alter a hydrogel's structure and therefore its transport properties.

Optimizing the fracture toughness of a dual cross-linked hydrogel via molecular dynamics simulation

In this work, a coarse-grained model is adopted to explore the fracture toughness of a dual cross-linked hydrogel which consists of a physically cross-linked network and a chemically cross-linked network. By calculating the fracture energy, the optimized fracture toughness of the hydrogel appears at the intermediate content of the chemical network. To understand it, the structure change of both the chemical network and the physical network is first characterized during the tensile process. For the chemical network, the fraction and rate of broken bonds gradually improve with increasing content of the chemical network while the strain range where the bond breakage occurs is reduced. For the physical network, the number of clusters and the interaction energy first increase and then decrease with increasing strain. This reflects the breakage and reformation of the physical network, which dissipates more energy and improves the fracture energy. Furthermore, by stress decomposition, the stress is mainly borne by the physical network at small strain and the chemical network at large strain, which proves their synergistic effect in enhancing the hydrogel. Then, the number of voids is calculated as a function of strain. It is found that the voids initiate in the weak region at small strain while in the position of the bond breakage at large strain. Moreover, the number of voids decreases with increasing content of the chemical network at small strain. Finally, the effect of the strength of the chemical network or the physical network on the fracture toughness is discussed. The optimized fracture toughness of hydrogel appears at the intermediate strength.

Stretchy and disordered: Toward understanding fracture in soft network materials via mesoscopic computer simulations

Soft network materials exist in numerous forms ranging from polymer networks, such as elastomers, to fiber networks, such as collagen. In addition, in colloidal gels, an underlying network structure can be identified, and several metamaterials and textiles can be considered network materials as well. Many of these materials share a highly disordered microstructure and can undergo large deformations before damage becomes visible at the macroscopic level. Despite their widespread presence, we still lack a clear picture of how the network structure controls the fracture processes of these soft materials. In this Perspective, we will focus on progress and open questions concerning fracture at the mesoscopic scale, in which the network architecture is clearly resolved, but neither the material-specific atomistic features nor the macroscopic sample geometries are considered. We will describe concepts regarding the network elastic response that have been established in recent years and turn out to be pre-requisites to understand the fracture response. We will mostly consider simulation studies, where the influence of specific network features on the material mechanics can be cleanly assessed. Rather than focusing on specific systems, we will discuss future challenges that should be addressed to gain new fundamental insights that would be relevant across several examples of soft network materials.

Sharing the Load: Stress Redistribution Governs Fracture of Polymer Double Networks

The stress response of polymer double networks depends not only on the properties of the constituent networks but also on the interactions arising between them. Here, we demonstrate, via coarse-grained simulations, that both their global stress response and their microscopic fracture mechanics are governed by load sharing through these internetwork interactions. By comparing our results with affine predictions, where stress redistribution is by definition homogeneous, we show that stress redistribution is highly inhomogeneous. In particular, the affine prediction overestimates the fraction of broken chains by almost an order of magnitude. Furthermore, homogeneous stress distribution predicts a single fracture process, while in our simulations, fracture of sacrificial chains takes place in two steps governed by load sharing within a network and between networks, respectively. Our results thus provide a detailed microscopic picture of how inhomogeneous stress redistribution after rupture of chains governs the fracture of polymer double networks.

First-Principles-Based Machine-Learning Molecular Dynamics for Crystalline Polymers with van der Waals Interactions

Machine-learning (ML) techniques have drawn an ever-increasing focus as they enable high-throughput screening and multiscale prediction of material properties. Especially, ML force fields (FFs) of quantum mechanical accuracy are expected to play a central role for the purpose. The construction of ML-FFs for polymers is, however, still in its infancy due to the formidable configurational space of its composing atoms. Here, we demonstrate the effective development of ML-FFs using kernel functions and a Gaussian process for an organic polymer, polytetrafluoroethylene (PTFE), with a data set acquired by first-principles calculations and ab initio molecular dynamics (AIMD) simulations. Even though the training data set is sampled only with short PTFE chains, structures of longer chains optimized by our ML-FF show an excellent consistency with density functional theory calculations. Furthermore, when integrated with molecular dynamics simulations, the ML-FF successfully describes various physical properties of a PTFE bundle, such as a density, melting temperature, coefficient of thermal expansion, and Young's modulus.

Molecular Characterization of Polymer Networks

Polymer networks are complex systems consisting of molecular components. Whereas the properties of the individual components are typically well understood by most chemists, translating that chemical insight into polymer networks themselves is limited by the statistical and poorly defined nature of network structures. As a result, it is challenging, if not currently impossible, to extrapolate from the molecular behavior of components to the full range of performance and properties of the entire polymer network. Polymer networks therefore present an unrealized, important, and interdisciplinary opportunity to exert molecular-level, chemical control on material macroscopic properties. A barrier to sophisticated molecular approaches to polymer networks is that the techniques for characterizing the molecular structure of networks are often unfamiliar to many scientists. Here, we present a critical overview of the current characterization techniques available to understand the relation between the molecular properties and the resulting performance and behavior of polymer networks, in the absence of added fillers. We highlight the methods available to characterize the chemistry and molecular-level properties of individual polymer strands and junctions, the gelation process by which strands form networks, the structure of the resulting network, and the dynamics and mechanics of the final material. The purpose is not to serve as a detailed manual for conducting these measurements but rather to unify the underlying principles, point out remaining challenges, and provide a concise overview by which chemists can plan characterization strategies that suit their research objectives. Because polymer networks cannot often be sufficiently characterized with a single method, strategic combinations of multiple techniques are typically required for their molecular characterization.

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.

Study on Large Deformation Behavior of Polyacrylamide Hydrogel Using Dissipative Particle Dynamics

Meso-scale models for hydrogels are crucial to bridge the conformation change of polymer chains in micro-scale to the bulk deformation of hydrogel in macro-scale. In this study, we construct coarse-grain bead-spring models for polyacrylamide (PAAm) hydrogel and investigate the large deformation and fracture behavior by using Dissipative Particle Dynamics (DPD) to simulate the crosslinking process. The crosslinking simulations show that sufficiently large diffusion length of polymer beads is necessary for the formation of effective polymer. The constructed models show the reproducible realistic structure of PAAm hydrogel network, predict the reasonable crosslinking limit of water content and prove to be sufficiently large for statistical averaging. Incompressible uniaxial tension tests are performed in three different loading rates. From the nominal stress-stretch curves, it demonstrated that both the hyperelasticity and the viscoelasticity in our PAAm hydrogel models are reflected. The scattered large deformation behaviors of three PAAm hydrogel models with the same water content indicate that the mesoscale conformation of polymer network dominates the mechanical behavior in large stretch. This is because the effective chains with different initial length ratio stretch to straight at different time. We further propose a stretch criterion to measure the fracture stretch of PAAm hydrogel using the fracture stretch of C-C bonds. Using the stretch criterion, specific upper and lower limits of the fracture stretch are given for each PAAm hydrogel model. These ranges of fracture stretch agree quite well with experimental results. The study shows that our coarse-grain PAAm hydrogel models can be applied to numerous single network hydrogel systems.

Polymer brush in articular cartilage lubrication: Nanoscale modelling and simulation

Human knee joints move smoothly under high load conditions due to articular cartilage and synovial fluid. Much attention is paid to the role of proteoglycans. It is suggested that a part of proteoglycan forms aggregate on the cartilage surface, making a polymer brush, which has an important role in lubrication. In order to examine the lubrication mechanism in detail, we constructed a full atom model of a polymer brush system, and carried out a series of molecular dynamics simulations to analyze its frictional properties under constant shear. We use chondroitin 6-sulfate molecules grafted on resilient surface as the polymer brush and water with sodium ions as the synovial liquid. In the steady state, polymers have large deformation and the flow of synovial fluid becomes deviate from the Coutette flow, leading to a drastic reduction of friction. Longer chains have larger reduction.

Protein diffusion from microwells with contrasting hydrogel domains

Understanding and controlling molecular transport in hydrogel materials is important for biomedical tools, including engineered tissues and drug delivery, as well as life sciences tools for single-cell analysis. Here, we scrutinize the ability of microwells-micromolded in hydrogel slabs-to compartmentalize lysate from single cells. We consider both (i) microwells that are "open" to a large fluid (i.e., liquid) reservoir and (ii) microwells that are "closed," having been capped with either a slab of high-density polyacrylamide gel or an impermeable glass slide. We use numerical modeling to gain insight into the sensitivity of time-dependent protein concentration distributions on hydrogel partition and protein diffusion coefficients and open and closed microwell configurations. We are primarily concerned with diffusion-driven protein loss from the microwell cavity. Even for closed microwells, confocal fluorescence microscopy reports that a fluid (i.e., liquid) film forms between the hydrogel slabs (median thickness of 1.7 μm). Proteins diffuse from the microwells and into the fluid (i.e., liquid) layer, yet concentration distributions are sensitive to the lid layer partition coefficients and the protein diffusion coefficient. The application of a glass lid or a dense hydrogel retains protein in the microwell, increasing the protein solute concentration in the microwell by ∼7-fold for the first 15 s. Using triggered release of Protein G from microparticles, we validate our simulations by characterizing protein diffusion in a microwell capped with a high-density polyacrylamide gel lid (p > 0.05, Kolmogorov-Smirnov test). Here, we establish and validate a numerical model useful for understanding protein transport in and losses from a hydrogel microwell across a range of boundary conditions.

From Microscale to Macroscale: Nine Orders of Magnitude for a Comprehensive Modeling of Hydrogels for Controlled Drug Delivery

Because of their inherent biocompatibility and tailorable network design, hydrogels meet an increasing interest as biomaterials for the fabrication of controlled drug delivery devices. In this regard, mathematical modeling can highlight release mechanisms and governing phenomena, thus gaining a key role as complementary tool for experimental activity. Starting from the seminal contribution given by Flory-Rehner equation back in 1943 for the determination of matrix structural properties, over more than 70 years, hydrogel modeling has not only taken advantage of new theories and the increasing computational power, but also of the methods offered by computational chemistry, which provide details at the fundamental molecular level. Simulation techniques such as molecular dynamics act as a "computational microscope" and allow for obtaining a new and deeper understanding of the specific interactions between the solute and the polymer, opening new exciting possibilities for an in silico network design at the molecular scale. Moreover, system modeling constitutes an essential step within the "safety by design" paradigm that is becoming one of the new regulatory standard requirements also in the field-controlled release devices. This review aims at providing a summary of the most frequently used modeling approaches (molecular dynamics, coarse-grained models, Brownian dynamics, dissipative particle dynamics, Monte Carlo simulations, and mass conservation equations), which are here classified according to the characteristic length scale. The outcomes and the opportunities of each approach are compared and discussed with selected examples from literature.

Molecular Modeling of Complex Cross-Linked Networks of PEGDA Nanogels

Poly(ethylene glycol) (PEG)-based nanogels are attractive for biomedical applications due to their biocompatibility, versatile end group chemistry, and ability to sterically shield encapsulated drug molecules. The characteristics of a hydrogel network govern the encapsulation and efficient delivery of drug molecules for a target application. A molecular-level description of network topology can complement experimental investigations to understand its effects on the structural properties of these nanogels. In this work, atomistic molecular simulations of heterogeneous, nonideal PEG-diacrylate (PEGDA) nanogels are presented. The effects of cross-linking density and topological features on the structural properties of PEGDA nanogels were studied. The average functionality was controlled to systematically study the effect of cross-linking density on the radius of gyration, shape, and mesh size of the nanogels. For a given average functionality, the impact of distinct network topologies on the structural properties was also studied. The aspect ratios, based on the gyration tensor, were calculated to characterize the shapes of these nanogels for different topologies. Nanogel structures with higher cross-linking densities showed a globular shape, while structures with lower cross-linking density showed shape anisotropy. The distribution and connectivity of the cross-linked junctions played a key role in determining the size and shape anisotropy of PEGDA nanogels; the number of unreacted chain ends and their connectivity directly affected the anisotropy. The mesh size, denoted by the limiting "free volume element" present in the nanogel samples, does not show a significant change with increasing average functionality. This work provides insight into the structural properties of heterogeneous hydrogels that aid the design of nonideal nanogel networks for a targeted drug delivery application.

Understanding hydrogelation processes through molecular dynamics

Molecular dynamics (MD) is currently one of the preferred techniques employed to understand hydrogelation processes for its ability to include large amounts of atoms in computational calculations, since substantial amounts of solvent molecules are involved in gel formation.

Dual Physically Cross-Linked κ-Carrageenan-Based Double Network Hydrogels with Superior Self-Healing Performance for Biomedical Application

Chemically linked double network (DN) hydrogels display extraordinary mechanical attributes but mostly suffer from poor self-healing property and unsatisfactory biocompatibility due to the irreversible breaks in their chemical-linked networks and the use of toxic chemical cross-linking agents. To address these limitations, we developed a novel κ-carrageenan/polyacrylamide (KC/PAM) DN hydrogel through a dual physical-cross-linking strategy, with the ductile, hydrophobically associated PAM being the first network, and the rigid potassium ion (K⁺) cross-linked KC being the second network. The dual physically cross-linked DN (DPC-DN) hydrogels with optimized KC concentration exhibit excellent fracture tensile stress (1320 ± 46 kPa) and toughness (fracture energy: 6900 ± 280 kJ/m³), comparable to those fully chemically linked DN hydrogels and physically chemically cross-linked hybrid DN hydrogels. Moreover, because of their unique dual physical-cross-linking structures, the KC/PAM hydrogels also demonstrated rapid self-recovery, remarkable notch-insensitivity, self-healing capability, as well as excellent cytocompatibility toward stem cells. Accordingly, this work presents a new strategy toward fabricating self-repairing DPC-DN hydrogels with outstanding mechanical behaviors and biocompatibility. The new type of DN hydrogels demonstrates strong potentiality in many challenging biomedical applications such as artificial diaphragm, tendon, and cartilage.

Shell-corona microgels from double interpenetrating networks

Polymer microgels with a dense outer shell offer outstanding features as universal carriers for different guest molecules. In this paper, microgels formed by an interpenetrating network comprised of collapsed and swollen subnetworks are investigated using dissipative particle dynamics (DPD) computer simulations, and it is found that such systems can form classical core-corona structures, shell-corona structures, and core-shell-corona structures, depending on the subchain length and molecular mass of the system. The core-corona structures consisting of a dense core and soft corona are formed at small microgel sizes when the subnetworks are able to effectively separate in space. The most interesting shell-corona structures consist of a soft cavity in a dense shell surrounded with a loose corona, and are found at intermediate gel sizes; the area of their existence depends on the subchain length and the corresponding mesh size. At larger molecular masses the collapsing network forms additional cores inside the soft cavity, leading to the core-shell-corona structure.

Tailoring the mechanical properties by molecular integration of flexible and stiff polymer networks

Designing a multiple-network structure at the molecular level to tailor the mechanical properties of polymeric materials is of great scientific and technological importance. Through the coarse-grained molecular dynamics simulation, we successfully construct an interpenetrating polymer network (IPN) composed of a flexible polymer network and a stiff polymer network. First, we find that there is an optimal chain stiffness for a single network (SN) to achieve the best stress-strain behavior. Then we turn to study the mechanical behaviors of IPNs. The result shows that the stress-strain behaviors of the IPNs appreciably exceed the sum of that of the corresponding single flexible and stiff network, which highlights the advantage of the IPN structure. By systematically varying the stiffness of the stiff polymer network of the IPNs, optimal stiffness also exists to achieve the best performance. We attribute this to a much larger contribution of the non-bonded interaction energy. Last, the effect of the component concentration ratio is probed. With the increase of the concentration of the flexible network, the stress-strain behavior of the IPNs is gradually enhanced, while an optimized concentration (around 60% molar ration) of the stiff network occurs, which could result from the dominant role of the enthalpy rather than the entropy. In general, our work is expected to provide some guidelines to better tailor the mechanical properties of the IPNs made of a flexible network and a stiff network, by manipulating the stiffness of the stiff polymer network and the component concentration ratio.

Effects of thermal shrinkage temperatures and comonomers on thermal shrinkage of uniaxially-stretched PET copolymer films: a molecular dynamics simulation approach

Thermal shrinkage ratios for PET copolymer models are correlated with the conformational change of polymer chains at molecular levels.

Conducting electrospun fibres with polyanionic grafts as highly selective, label-free, electrochemical biosensor with a low detection limit for non-Hodgkin lymphoma gene

A highly selective, label-free sensor for the non-Hodgkin lymphoma gene, with an aM detection limit, utilizing electrochemical impedance spectroscopy (EIS) is presented. The sensor consists of a conducting electrospun fibre mat, surface-grafted with poly(acrylic acid) (PAA) brushes and a conducting polymer sensing element with covalently attached oligonucleotide probes. The sensor was fabricated from electrospun NBR rubber, embedded with poly(3,4-ethylenedioxythiophene) (PEDOT), followed by grafting poly(acrylic acid) brushes and then electrochemically polymerizing a conducting polymer monomer with ssDNA probe sequence pre-attached. The resulting non-Hodgkin lymphoma gene sensor showed a detection limit of 1aM (1 × 10^-18mol/L), more than 400 folds lower compared to a thin-film analogue. The sensor presented extraordinary selectivity, with only 1%, 2.7% and 4.6% of the signal recorded for the fully non-complimentary, T-A and G-C base mismatch oligonucleotide sequences, respectively. We suggest that such greatly enhanced selectivity is due to the presence of negatively charged carboxylic acid moieties from PAA grafts that electrostatically repel the non-complementary and mismatch DNA sequences, overcoming the non-specific binding.

pH-Responsive Coassembly of Oligo(ethylene glycol)-Coated Gold Nanoparticles with External Anionic Polymers via Hydrogen Bonding

Stimuli-responsive assembly of gold nanoparticles (AuNPs) with precise control of the plasmonic properties, assembly size, and stimuli responsivity has shown potential benefits with regard to biosensing devices and drug-delivery systems. Here we present a new pH-responsive coassembly system of oligo(ethylene glycol) (OEG)-coated AuNPs with anionic polymers as an external mediator via hydrogen bonding in water. Hydrogen-bond-driven coassemblies of OEG-AuNPs with poly(acrylic acid) (PAA) were confirmed by the monitoring of plasmonic peaks and hydrodynamic diameters. In this system, the protonation of anionic polymers on change in pH triggered the formation of hydrogen bond between the OEG-AuNPs and polymers, providing sensitive pH responsivity. The plasmonic properties and assembly size are affected by both the ratio of PAA to AuNPs and the molecular weight of PAAs. In addition, the attachment of hydrophobic groups to the surface ligand or anionic polymer changed the responsive pH range. These results demonstrated that the coassembly with an external mediator via hydrogen bonding provides a stimuli-responsive assembly system with tunable plasmonic properties, assembly size, and stimuli responsivity.

Structure formation in nanocomposite hydrogels

We use molecular dynamics simulations to study structure formation in physically associating nanocomposite hydrogels. Nanofillers were modeled as rigid bodies of disk-like shapes and physical crosslinks were simulated by introducing a short-range attraction between the nanofillers and polymer chain ends. The structure, dynamics and mechanics of these polymer gels were studied as a function of nanofiller volume fraction. We observe the formation of a percolated network in the hydrogels, with an ordered local structure but disordered globally, as we increase the filler fraction. This locally ordered structure was a result of the anisotropy of the disk-like fillers. The dynamics of polymers showed significant caging effects in the gel state. Stress autocorrelation and elongation results were analyzed as a function of nano-filler concentrations. Comparisons with nanofillers of different shapes showed that disk-like nanofillers are more effective in strengthening the hydrogels than spherical nanofillers.

Shape and Mechanical Control of Poly(ethylene oxide) Based Polymersome with Polyoxometalates via Hydrogen Bond

Polymersomes are self-assembled vesicles of amphiphilic block copolymers and have been explored for wide applications from drug delivery to micro/nanoreactors. As polymersomes are soft and highly deformable, their shape instability due to osmolarity difference across polymer membranes and low elasticity could conversely limit their practical use. Instead of selecting particular polymer chemical reactions to enhance the mechanical properties, we have employed inorganic polyoxometalate (POM) clusters as simple physical cross-linkers to control the shape and mechanical stability of polymersomes in aqueous suspensions. Robust spherical shape with enhanced elastic and bending moduli of POM-dressed poly(ethylene oxide) (PEO) based polymersomes is achieved. We have accounted for the hydrogen bonding between POM and PEO blocks for the adsorption and stabilization of POMS on polymersomes, whose interaction strength could also be tuned by mixing solvents of hydrogen bond donors or receptors with water. The stimuli-responsive properties of POMs are introduced in POM-dressed polymersomes upon the interaction of POMs with PEO blocks in aqueous media. As POM can be used as nanomedicines, catalysts, and other functional nanomaterials, POM-dressed polymersomes with significant shape and mechanical reinforcement could broaden the applications of PEO-based polymersomes and other PEO-tethered nanocolloids.

Hydrogen-bonded polymer complexes and nanocages of weak polyacids templated by a Pluronic® block copolymer

We investigate the phase behavior, morphology, and temperature response of hydrogen-bonded assemblies formed by a triblock copolymer Pluronic® F127 (F127) and polycarboxylic acids of varied hydrophobicity and chain lengths. As confirmed by FTIR, the complexes of poly(acrylic acid) (PAA) and poly(methacrylic acid) (PMAA) with F127 at acidic pH were stabilized by multiple hydrogen bonding between carboxylic acid groups of polyacids and ether groups of F127. The colloidal stability of the polyacid/F127 complexes (their occurrence as stable dispersions, slowly coagulating dispersions or precipitates) was dependent on the composition of complexes, polyacid molecular weight and hydrophobicity, as well as temperature. For both polyacids, complexes could not be solubilized in excess of polyacids, but excess of F127 resulted in the formation of colloidally stable nanostructured clusters whose size could be controlled from tens to hundreds of nanometers by the polyacid-to-F127 ratio, temperature, and the polyacid molecular weight. Hydrophobicity of polyacids had a dramatic effect on the temperature response of Pluronic®-enriched assemblies. While PMAA suppressed the LCST behavior of F127 due to binding within the temperature-responsive PPO core of F127, more hydrophilic PAA allowed F127 micellization and supported reversible, temperature-induced re-structuring of PAA-F127 clusters. At temperatures above the LCST of Pluronic®, low-molecular-weight PAA formed nanosized dispersed complexes, in which the polyacid chains were wrapped around individual F127 micelles. Chemical crosslinking of PAA in the shells of these complexes followed by removal of the templating F127 cores resulted in easy-to-prepare monodisperse pH-responsive polymer nanocages with controllable size and swelling amplitude.

Computing thermomechanical properties of dry homopolymers used as raw materials for formulation of biomedical hydrogels

Different static properties have been calculated with COMPASS force field for polyacrylamide, poly(2-hydroxyethylacrylate) (HEA), poly(2-hydroxyethylmethacrylate) (HEMA), poly(glycidylmethacrylate) (GMA), polyethylene glycol (PEG), and poly(2,2,2-trifluoroethylmethacrylate) (TFEM). For each polymers, the calculated values were averaged on five equilibrated configurations of amorphous cell composed of one atactic chain containing 100 repeat units. The ranking obtained from the densities calculated at 300 K is TFEM > HEA ≈ xpolycrylamide > HEMA ≈ GMA > PEG. Concerning the glass transition temperature we have obtained polyacrylamide > HEMA ≈ GMA ≈ HEA > PEG, and polyacrylamide > HEMA ≈ HEA > GMA ≈ PEG > TFEM for the bulk modulus. The calculated results, when available, have been compared with experimental data coming from literature.