Molecular dynamics simulations of polyelectrolyte adsorption

We have performed molecular dynamics simulations of polyelectrolyte adsorption at oppositely charged surfaces from dilute polyelectrolyte solutions. In our simulations, polyelectrolytes were modeled by chains of charged Lennard-Jones particles with explicit counterions. We have studied the effects of the surface charge density, surface charge distribution, solvent quality for the polymer backbone, strength of the short-range interactions between polymers and substrates on the polymer surface coverage, and the thickness of the adsorbed layer. The polymer surface coverage monotonically increases with increasing surface charge density for almost all studied systems except for the system of hydrophilic polyelectrolytes adsorbing at hydrophilic surfaces. In this case the polymer surface coverage saturates at high surface charge densities. This is due to additional monomer-monomer repulsion between adsorbed polymer chains, which becomes important in dense polymeric layers. These interactions also preclude surface overcharging by hydrophilic polyelectrolytes at high surface charge densities. The thickness of the adsorbed layer shows monotonic dependence on the surface charge density for the systems of hydrophobic polyelectrolytes for both hydrophobic and hydrophilic surfaces. Thickness is a decreasing function of the surface charge density in the case of hydrophilic surfaces while it increases with the surface charge density for hydrophobic substrates. Qualitatively different behavior is observed for the thickness of the adsorbed layer of hydrophilic polyelectrolytes at hydrophilic surfaces. In this case, thickness first decreases with increasing surface charge density, then it begins to increase.     

Molecular dynamics simulations of water droplets on polymer surfaces

Molecular dynamics simulations were used to study the wetting of polymer surfaces with water. Contact angles of water droplets on crystalline and two amorphous polyethylene (PE) and poly(vinyl chloride) (PVC) surfaces were extracted from atomistic simulations. Crystalline surfaces were produced by duplicating the unit cell of an experimental crystal structure, and amorphous surfaces by pressing the bulk polymer step by step at elevated temperature between two repulsive grid surfaces to a target density. Different-sized water droplets on the crystalline PE surface revealed a slightly positive line tension on the order of 10(-12)-10(-11) N, whereas droplets on crystalline PVC did not yield a definite line tension. Microscopic contact angles produced by the simple point charge (SPC) water model were mostly a few degrees smaller than those produced by the extended SPC model, which, as the model with lowest bulk energy, presents an upper boundary for contact angles. The macroscopic contact angle for the SPC model was 94 degrees on crystalline PVC and 113 degrees on crystalline PE. Amorphicity of the surface increased the water contact angle on PE but decreased it on PVC, for both water models. If the simulated contact angles on crystalline and amorphous surfaces are combined in proportion to the crystallinity of the polymer in question, simulated values in relatively good agreement with measured values are obtained.     

Brownian dynamics simulations of polyelectrolyte adsorption onto charged patterned surfaces

The adsorption of single polyelectrolyte molecules onto surfaces decorated with periodic arrays of charged patches was studied using Brownian dynamics simulations. A free-draining, freely jointed bead-rod chain was used to model the polyelectrolyte, and electrostatic interactions were incorporated using a screened Coulombic potential with the excluded volume accounted for by a hard-sphere potential. The simulations predicted that the polyelectrolyte lies close to the adsorbing surface if the patch length, surface charge density, and screening length are sufficiently large. Chain conformations were found to be very sensitive to patch length, patch spacing, and the nature of the charge on adjacent patches. This is due both to the size of the polymer relative to patch length and spacing and to the structure of the electric field near the surface. In some cases, the component of the radius of gyration parallel to the surface can be made smaller than its free-solution value, which is contrary to what is observed for a uniformly charged surface. Isolated charged patches were also considered, and significant adsorption was observed above a critical surface charge density. The results demonstrate how polyelectrolyte conformations can be controlled by the design of the charged patches and may be useful for applications in which adsorbed polyelectrolyte films play a key role.     

Brownian dynamics simulations of polyelectrolyte adsorption onto topographically patterned surfaces

The effect of patterned surface topography on the adsorption of single polyelectrolyte molecules is explored using Brownian dynamics simulations. The polyelectrolyte is modeled as a free-draining, freely jointed bead-rod chain, and electrostatic interactions are incorporated using a screened Coulombic potential with excluded volume interactions accounted for by the repulsive part of a Lennard-Jones potential. Topography consisting of periodically spaced valleys of square cross section separated by flat hills is considered. Chain conformations are characterized for a wide range of valley widths, depths, and spacings as well as for several different types of surface charge distributions. Depending on the parameter values describing the topography, the chains are found to adopt conformations ranging from flat and extended to those associated with bridge-, brush-, or semi-bridge-like structures. The formation of these structures is rationalized on the basis of a free-energy model that takes into account the increase in free energy due to entropic confinement, excluded volume interactions, and chain stretching as well as the decrease in free energy due to bead-surface electrostatic attraction. The results of this work are expected to be useful in designing patterned surface topography to control the conformations of adsorbed polyelectrolyte molecules.     

Molecular dynamics simulation study of water adsorption on hydroxylated graphite surfaces

In this paper, we present results from molecular dynamic simulations devoted to the characterization of the interaction between water molecules and hydroxylated graphite surfaces considered as models for surfaces of soot emitted by aircraft. The hydroxylated graphite surfaces are modeled by anchoring several OH groups on an infinite graphite plane. The molecular dynamics simulations are based on a classical potential issued from quantum chemical calculations. They are performed at three temperatures (100, 200, and 250 K) to provide a view of the structure and dynamics of water clusters on the model soot surface. These simulations show that the water-OH sites interaction is quite weak compared to the water-water interaction. This leads to the clustering of the water molecules above the surface, and the corresponding water aggregate can only be trapped by the OH sites when the temperature is sufficiently low, or when the density of OH sites is sufficiently high.     

Fibronectin and bone morphogenetic protein-2-decorated poly(OEGMA-r-HEMA) brushes promote osseointegration of titanium surfaces

To be better used as medical implants in orthopedic and dental clinical applications, titanium and titanium-based alloys need to be capable of inducing osteogenesis. Here we describe a method that allows the facile decoration of titanium surfaces to impart an osteogenesis capacity. A Ti surface was first deposited on a poly(OEGMA-r-HEMA) film using surface-initiated atom-transfer radical polymerization (SI-ATRP) with the further step of carboxylation. The modified surfaces were resistant to cell adhesion. Fibronectin (FN) and recombinant human bone morphogenetic protein-2 (rhBMP-2) were further immobilized onto p(OEGMA-r-HEMA) matrices. Our results demonstrate that the FN- and rhBMP-2-conjugated polymer surfaces could induce the adhesion of MC3T3 cells on Ti surfaces. Moreover, the protein-tethered surface exhibited enhanced cell differentiation in terms of alkaline phosphatase activity compared to that of the pristine Ti surface at similar cell proliferation rates. This research establishes a simple modification method of Ti surfaces via Ti-thiolate self-assembled monolayers (SAMs) and SI-ATRP and identifies a dual-functional Ti surface that combines antifouling and osseointegration promotion.     

Molecular dynamics simulations of n-Alkane melts confined between solid surfaces

Molecular dynamics simulations of polyethylene chains (C_nH_2n+2 for n = 13, 16, 28, 60) have been carried out to investigate both equilibrium and dynamic properties of polymer melts confined between flat solid surfaces. We observe an oscillatory monomer density in the direction normal to the solid surfaces, which depends on the size of the monomers (i.e., on the volume density of the system). The packing manner of monomer segments, segment orientation, and local conformations of chains are found to be independent of chain length. In addition, preferential interfacial adsorption of chain ends is observed. The chains are flattened close to the surface and many molecules assume essentially two-dimensional train configurations even in the case of C₆₀H₁₂₂ melts. The apparent self-diffusivities of the centers of mass of the molecules depend on their distance from a surface. Molecules adjacent to a surface exhibit a reduced mobility perpendicular to the surface and an increased one parallel to it. The self-diffusion constant parallel to a surface depends strongly on the size of the monomers. An increase of the united atom diameter by 10% reduces the diffusion constant by a factor of three, in good agreement with the experimental value. © 1994 John Wiley & Sons, Inc.     

Molecular dynamics simulations of aggregation behavior of sodium dodecyl sulfate on SiO2 and CaCO3 surfaces

Molecular dynamics simulations were carried out to investigate the aggregation behavior of anionic surfactant sodium dodecyl sulfate (SDS) on the surfaces of SiO₂ and CaCO₃. The results indicate that SDS molecules formed a spherical micelle structure near the SiO₂ surface; moreover, there were more head groups near the SiO₂ surface. However, they could form a self‐assemble film on the CaCO₃ surface. The self‐assemble film of SDS on the CaCO₃ surface was more stable than that on the SiO₂ surface. Our simulation results have a certain significance to understand the aggregation behavior of SDS on different surfaces on molecular level.     

Molecular dynamics simulations on the adsorption and surface phenomena of simple fluid in porous media

Equilibrium molecular dynamics simulations have been used to investigate the fluid adsorption phenomena and calculate the surface tension in porous media at different temperatures, densities and pore widths. The facts that most of the fluid particles are adsorbed adherent to the pore walls and the surface forms near the walls have been found, and that the surface is not stable which means there exists an oscillation phenomenon in pores have been also found. The surface tension in pores is much bigger than that in macrovolume systems like normal liquid–liquid and liquid–vapor interfaces, and it will increase with the increase of density and pore width, but will decrease with the increase of temperature.     

Molecular dynamics simulations of the electrical double layer on smectite surfaces contacting concentrated mixed electrolyte (NaCl-CaCl2) solutions

We report new molecular dynamics results elucidating the structure of the electrical double layer (EDL) on smectite surfaces contacting mixed NaCl-CaCl(2) electrolyte solutions in the range of concentrations relevant to pore waters in geologic repositories for CO(2) or high-level radioactive waste (0.34-1.83 mol(c) dm(-3)). Our results confirm the existence of three distinct ion adsorption planes (0-, β-, and d-planes), often assumed in EDL models, but with two important qualifications: (1) the location of the β- and d-planes are independent of ionic strength or ion type and (2) "indifferent electrolyte" ions can occupy all three planes. Charge inversion occurred in the diffuse ion swarm because of the affinity of the clay surface for CaCl(+) ion pairs. Therefore, at concentrations ≥0.34 mol(c) dm(-3), properties arising from long-range electrostatics at interfaces (electrophoresis, electro-osmosis, co-ion exclusion, colloidal aggregation) will not be correctly predicted by most EDL models. Co-ion exclusion, typically neglected by surface speciation models, balanced a large part of the clay mineral structural charge in the more concentrated solutions. Water molecules and ions diffused relatively rapidly even in the first statistical water monolayer, contradicting reports of rigid "ice-like" structures for water on clay mineral surfaces.     

Molecular dynamics simulations of D2O ice photodesorption

Molecular dynamics (MD) calculations have been performed to study the ultraviolet (UV) photodissociation of D(2)O in an amorphous D(2)O ice surface at 10, 20, 60, and 90 K, in order to investigate the influence of isotope effects on the photodesorption processes. As for H(2)O, the main processes after UV photodissociation are trapping and desorption of either fragments or D(2)O molecules. Trapping mainly takes place in the deeper monolayers of the ice, whereas desorption occurs in the uppermost layers. There are three desorption processes: D atom, OD radical, and D(2)O molecule photodesorption. D(2)O desorption takes places either by direct desorption of a recombined D(2)O molecule, or when an energetic D atom produced by photodissociation kicks a surrounding D(2)O molecule out of the surface by transferring part of its momentum. Desorption probabilities are calculated for photoexcitation of D(2)O in the top four monolayers and are compared quantitatively with those for H(2)O obtained from previous MD simulations of UV photodissociation of amorphous water ice at different ice temperatures [Arasa et al., J. Chem. Phys. 132, 184510 (2010)]. The main conclusions are the same, but the average D atom photodesorption probability is smaller than that of the H atom (by about a factor of 0.9) because D has lower kinetic energy than H, whereas the average OD radical photodesorption probability is larger than that of OH (by about a factor of 2.5-2.9 depending on ice temperature) because OD has higher translational energy than OH for every ice temperature studied. The average D(2)O photodesorption probability is larger than that of H(2)O (by about a factor of 1.4-2.3 depending on ice temperature), and this is entirely due to a larger contribution of the D(2)O kick-out mechanism. This is an isotope effect: the kick-out mechanism is more efficient for D(2)O ice, because the D atom formed after D(2)O photodissociation has a larger momentum than photogenerated H atoms from H(2)O, and D transfers momentum more easily to D(2)O than H to H(2)O. The total (OD + D(2)O) yield has been compared with experiments and the total (OH + H(2)O) yield from previous simulations. We find better agreement when we compare experimental yields with calculated yields for D(2)O ice than when we compare with calculated yields for H(2)O ice.     

Molecular dynamics simulations on the effects of diameter and chirality on hydrogen adsorption in single walled carbon nanotubes

We present systematic molecular dynamics simulation studies of hydrogen storage in single walled carbon nanotubes of various diameters and chiralities using a recently developed curvature-dependent force field. Our main objective is to address the following fundamental issues: 1. For a given H2 loading and nanotube type, what is the H2 distribution in the nanotube bundle? 2. For a given nanotube type, what is the maximal loading (H2 coverage)? 3. What is the diameter range and chirality for which H2 adsorption is most energetically favorable? Our simulation results suggest strong dependence of H2 adsorption energies on the nanotube diameter but less dependence on the chirality. Substantial lattice expansion upon H2 adsorption was found. The average adsorption energy increases with the lowering of nanotube diameter (higher curvature) and decreases with higher H2 loading. The calculated H2 vibrational power spectra and radial distribution functions indicate a strong attractive interaction between H2 and nanotube walls. The calculated diffusion coefficients are much higher than what has been reported for H2 in microporous materials such as zeolites, indicating that diffusivity does not present a problem for hydrogen storage in carbon nanotubes.     

Molecular dynamics simulations of La2O3 thin films on SiO2

Classical molecular dynamics is used to investigate the equilibrium state of the surface region and interface of heteroepitaxial La2O3 thin films.Due to the lattice mismatch,heteroepitaxial thin films are subject to very large stress.For this reason the behavior of La2O3 thin films at SiO2interface becomes an important concern.Our result indicates that La2O3 can uniformly wet SiO2 surface.The properties of the simulated films are analyzed and the lack of any discernible crystalline phase in epitaxial La2O3 on SiO2 indicates that the lattice mismatch between SiO2 and La2O3 is sufficiently large to prevent the formation of even short-range orders in La2O3 film.     

Molecular dynamics simulations of La_2O_3 thin films on SiO_2

Classical molecular dynamics is used to investigate the equilibrium state of the surface region and interface of heteroepitaxial La_2O_3 thin films.Due to the lattice mismatch,heteroepitaxial thin films are subject to very large stress.For this reason the behavior of La_2O_3 thin films at SiO_2interface becomes an important concern.Our result indicates that La_2O_3 can uniformly wet SiO_2 surface.The properties of the simulated films are analyzed and the lack of any discernible crystalline phase in epitaxial La_2O_3 on SiO_2 indicates that the lattice mismatch between SiO_2 and La_2O_3 is sufficiently large to prevent the formation of even short-range orders in La_2O_3 film.     

Molecular dynamics and density functional theory simulations of cesium and strontium adsorption on illite/ smectite

Journal of Radioanalytical and Nuclear Chemistry - Based on molecular dynamics (MD) and density function theory (DFT) simulation, the adsorption mechanisms of Cs+ and Sr2+ on the...     

Comparison between heparin-conjugated fibrin and collagen sponge as bone morphogenetic protein-2 carriers for bone regeneration

Bone morphogenetic protein-2 (BMP-2) is used to promote bone regeneration. However, the bone regeneration ability of BMP-2 relies heavily on the delivery vehicle. Previously, we have developed heparin- conjugated fibrin (HCF), a vehicle for long-term delivery of BMP-2 and demonstrated that long-term delivery of BMP-2 enhanced its osteogenic efficacy as compared to short-term delivery at an equivalent dose. The aim of this study was to compare the bone-forming ability of the BMP-2 delivered by HCF to that delivered by clinically utilized BMP-2 delivery vehicle collagen sponge. An in vitro release profile of BMP-2 showed that HCF released 80% of the loaded BMP-2 within 20 days, whereas collagen sponge released the same amount within the first 6 days. Moreover, the BMP-2 released from the HCF showed significantly higher alkaline phosphatase activity than the BMP-2 released from collagen sponge at 2 weeks in vitro. Various doses of BMP-2 were delivered with HCF or collagen sponge to mouse calvarial defects. Eight weeks after the treatment, bone regeneration was evaluated by computed tomography, histology, and histomorphometric analysis. The dose of BMP-2 delivered by HCF to achieve 100% bone formation in the defects was less than half of the BMP-2 dose delivered by collagen sponge to achieve a similar level of bone formation. Additionally, bone regenerated by the HCF-BMP-2 had higher bone density than bone regenerated by the collagen sponge-BMP-2. These data demonstrate that HCF as a BMP-2 delivery vehicle exerts better osteogenic ability of BMP-2 than collagen sponge, a clinically utilized delivery vehicle.     

Molecular dynamics simulations of halide glasses

Summary Halide glasses have been extensively studied in recent years because of their potential application as infrared transmitting fibre optic materials. They are believed to be more ionic than glasses based on silica and should therefore be more amenable to molecular dynamics simulation using simple two-body potentials. Here the main features of structural models derived using such techniques are described and compared with available structural data. Possible future applications of this approach are outlined.     

Monte Carlo simulations of antibody adsorption and orientation on charged surfaces

Monte Carlo simulations were performed to study the adsorption and orientation of antibodies on charged surfaces based on both colloidal and all-atom models. The colloidal model antibody consists of 12 connected beads representing the 12 domains of an antibody molecule. The structure of the all-atom antibody model was taken from the protein databank. The effects of the surface charge sign and density, the solution pH and ionic strength on the adsorption and orientation of different colloidal model antibodies with different dipole moments were examined. Simulation results show that both the 12-bead and the all-atom models of the antibody, for which the dipole moment points from the Fc to (Fab)2 fragments, tend to have the desired "end-on" orientation on positively charged surfaces and undesired "head-on" orientation on negatively charged surfaces at high surface charge density and low solution ionic strength where electrostatic interactions dominate. At low surface charge density and high solution ionic strength where van der Waals interactions dominate, 12-bead model antibodies tend to have "lying-flat" orientation on surfaces. The orientation of adsorbed antibodies results from the compromise between electrostatic and van der Waals interactions. The dipole moment of an antibody is an important factor for antibody orientation on charged surfaces when electrostatic interactions dominate. This charge-driven protein orientation hypothesis was verified by our simulations results in this work. It was further confirmed by surface plasmon resonance biosensor and time-of-flight secondary ion mass spectrometry experiments reported elsewhere.     

Molecular dynamics simulations of peptide-surface interactions

Proteins, which are bioactive molecules, adsorb on implants placed in the body through complex and poorly understood mechanisms and directly influence biocompatibility. Molecular dynamics modeling using empirical force fields provides one of the most direct methods of theoretically analyzing the behavior of complex molecular systems and is well-suited for the simulation of protein adsorption behavior. To accurately simulate protein adsorption behavior, a force field must correctly represent the thermodynamic driving forces that govern peptide residue-surface interactions. However, since existing force fields were developed without specific consideration of protein-surface interactions, they may not accurately represent this type of molecular behavior. To address this concern, we developed a host-guest peptide adsorption model in the form of a G(4)-X-G(4) peptide (G is glycine, X is a variable residue) to enable determination of the contributions to adsorption free energy of different X residues when adsorbed to functionalized Au-alkanethiol self-assembled monolayers (SAMs). We have previously reported experimental results using surface plasmon resonance (SPR) spectroscopy to measure the free energy of peptide adsorption for this peptide model with X = G and K (lysine) on OH and COOH functionalized SAMs. The objectives of the present research were the development and assessment of methods to calculate adsorption free energy using molecular dynamics simulations with the GROMACS force field for these same peptide adsorption systems, with an oligoethylene oxide (OEG) functionalized SAM surface also being considered. By comparing simulation results to the experimental results, the accuracy of the selected force field to represent the behavior of these molecular systems can be evaluated. From our simulations, the G(4)-G-G(4) and G(4)-K-G(4) peptides showed minimal to no adsorption to the OH SAM surfaces and the G(4)-K-G(4) showed strong adsorption to the COOH SAM surface, which is in agreement with our SPR experiments. Contrary to our experimental results, however, the simulations predicted a relatively strong adsorption of G(4)-G-G(4) peptide to the COOH SAM surface. In addition, both peptides were unexpectedly predicted to adsorb to the OEG surface. These findings demonstrate the need for GROMACS force field parameters to be rebalanced for the simulation of peptide adsorption behavior on SAM surfaces. The developed methods provide a direct means of assessing, modifying, and validating force field performance for the simulation of peptide and protein adsorption to surfaces, without which little confidence can be placed in the simulation results that are generated with these types of systems.