openAccessArticle: FalsePage Range: 245-245doi: 10.1016/j.parco.2011.08.005Harvest Date: 2016-01-12 15:10:37issueName:cover date: 2012-04-01pubType

We report our study of a silica-water interface using reactive molecular dynamics. This first-of-its-kind simulation achieves length and time scales required to investigate the detailed chemistry of the system. Our molecular dynamics approach is based on the ReaxFF force field of van Duin et al. ͓J. Phys. Chem. A 107, 3803 ͑2003͔͒. The specific ReaxFF implementation ͑SERIALREAX͒ and force fields are first validated on structural properties of pure silica and water systems. Chemical reactions between reactive water and dangling bonds on a freshly cut silica surface are analyzed by studying changing chemical composition at the interface. In our simulations, reactions involving silanol groups reach chemical equilibrium in ϳ250 ps. It is observed that water molecules penetrate a silica film through a proton-transfer process we call "hydrogen hopping," which is similar to the Grotthuss mechanism. In this process, hydrogen atoms pass through the film by associating and dissociating with oxygen atoms within bulk silica, as opposed to diffusion of intact water molecules. The effective diffusion constant for this process, taken to be that of hydrogen atoms within silica, is calculated to be 1.68ϫ 10 −6 cm 2 / s. Polarization of water molecules in proximity of the silica surface is also observed. The subsequent alignment of dipoles leads to an electric potential difference of ϳ10.5 V between the silica slab and water.

Molecular dynamics simulations are performed on two hydrated dipalmitoylphosphatidylcholine bilayer systems: one with pure water and one with added NaCl. Due to the rugged nature of the membrane/electrolyte interface, ion binding to the membrane surface is characterized by the loss of ion hydration. Using this structural characterization, binding of Na(+) and Cl(-) ions to the membrane is observed, although the binding of Cl(-) is seen to be slightly weaker than that of Na(+). Dehydration is seen to occur to a different extent for each type of ion. In addition, the excess binding of Na(+) gives rise to a net positive surface charge density just outside the bilayer. The positive density produces a positive electrostatic potential in this region, whereas the system without salt shows an electrostatic potential of zero.

We introduce a new force field (43A1-S3) for simulation of membranes by the Gromacs simulation package. Construction of the force fields is by standard methods of electronic structure computations for bond parameters and charge distribution and specific volumes and heats of vaporization for small-molecule components of the larger lipid molecules for van der Waals parameters. Some parameters from the earlier 43A1 force field are found to be correct in the context of these calculations, while others are modified. The validity of the force fields is demonstrated by correct replication of X-ray form factors and NMR order parameters over a wide range of membrane compositions in semi-isotropic NTP 1 atm simulations. 43-A1-S3 compares favorably with other force fields used in conjunction with the Gromacs simulation package with respect to the breadth of phenomena that it accurately reproduces.

No abstract

We have carried out an atomic-level molecular dynamics simulation of a system of nanoscopic size containing a domain of 18:0 sphingomyelin and cholesterol embedded in a fully hydrated dioleylposphatidylcholine (DOPC) bilayer. To analyze the interaction between the domain and the surrounding phospholipid, we calculate order parameters and area per molecule as a function of molecule type and proximity to the domain. We propose an algorithm based on Voronoi tessellation for the calculation of the area per molecule of various constituents in this ternary mixture. The calculated areas per sphingomyelin and cholesterol are in agreement with previous simulations. The simulation reveals that the presence of the liquid-ordered domain changes the packing properties of DOPC bilayer at a distance as large as approximately 8 nm. We calculate electron density profiles and also calculate the difference in the thickness between the domain and the surrounding DOPC bilayer. The calculated difference in thickness is consistent with data obtained in atomic force microscopy experiments.

Two mixed bilayers containing dipalmitoylphosphatidylcholine and dipalmitoylphosphatidylserine at a ratio of 5:1 are simulated in NaCl electrolyte solutions of different concentration using the molecular dynamics technique. Direct NH.O and CH.O hydrogen bonding between lipids was observed to serve as the basis of interlipid complexation. It is deduced from our results and previous studies that dipalmitoylphosphatidylcholine alone is less likely to form interlipid complexes than in the presence of bound ions or other bilayer "impurities" such as dipalmitoylphosphatidylserine. The binding of counterions is observed and quantitated. Based upon the calculated ion binding constants, the Gouy-Chapman surface potential (theta) is calculated. In addition we calculated the electrostatic potential profile (Phi) by twice integrating the system charge distribution. A large discrepancy between and the value of Phi at the membrane surface is observed. However, at "larger" distance from the bilayer surface, a qualitative similarity in the z-profiles of Phi and psi(GC) is seen. The discrepancy between the two potential profiles near the bilayer surface is attributed to the discrete and nonbulk-like nature of water in the interfacial region and to the complex geometry of this region.

It is known from experimental studies that lipid bilayers composed of unsaturated phospholipids, sphingomyelin, and cholesterol contain microdomains rich in sphingomyelin and cholesterol. These domains are similar to "rafts" isolated from cell membranes, although the latter are much smaller in lateral size. Such domain formation can be a result of very specific and subtle lipid-lipid interactions. To identify and study these interactions, we have performed two molecular dynamics simulations, of 200-ns duration, of dioleylphosphatidylcholine (DOPC), sphingomyelin (SM), and cholesterol (Chol) systems, a 1:1:1 mixture of DOPC/SM/Chol, and a 1:1 mixture of DOPC/SM. The simulations show initial stages of the onset of spontaneous phase-separated domains in the systems. On the simulation timescale cholesterol favors a position at the interface between the ordered SM region and the disordered DOPC region in the ternary system and accelerates the process of domain formation. We find that the smooth alpha-face of Chol preferentially packs next to SM molecules. Based on a comparative analysis of interaction energies, we find that Chol molecules do not show a preference for SM or DOPC. We conclude that Chol molecules assist in the process of domain formation and the process is driven by entropic factors rather than differences in interaction energies.