Molecular dynamics simulations and Kelvin probe force microscopy to study of cholesterol-induced electrostatic nanodomains in complex lipid mixtures.

The molecular arrangement of lipids and proteins within biomembranes and monolayers gives rise to complex film morphologies as well as regions of distinct electrical surface potential, topographical and electrostatic nanoscale domains. To probe these nanodomains in soft matter is a challenging task both experimentally and theoretically. This work addresses the effects of cholesterol, lipid composition, lipid charge, and lipid phase on the monolayer structure and the electrical surface potential distribution. Atomic force microscopy (AFM) was used to resolve topographical nanodomains and Kelvin probe force microscopy (KPFM) to resolve electrical surface potential of these nanodomains in lipid monolayers. Model monolayers composed of dipalmitoylphosphatidylcholine (DPPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1,2-dioleoyl-sn-glycero-3-[phospho-rac-(3-lysyl(1-glycerol))] (DOPG), and cholesterol were studied. It is shown that cholesterol changes nanoscale domain formation, affecting both topography and electrical surface potential. The molecular basis for differences in electrical surface potential was addressed with atomistic molecular dynamics (MD). MD simulations are compared the experimental results, with 100 s of mV difference in electrostatic potential between liquid-disordered bilayer (Ld, less cholesterol and lower chain order) and a liquid-ordered bilayer (Lo, more cholesterol and higher chain order). Importantly, the difference in electrostatic properties between Lo and Ld phases suggests a new mechanism by which membrane composition couples to membrane function.

Hydroxyapatite Growth Inhibition Effect of Pellicle Statherin Peptides.

In our recent studies, we have shown that in vivo-acquired enamel pellicle is a sophisticated biological structure containing a significant portion of naturally occurring salivary peptides. From a functional aspect, the identification of peptides in the acquired enamel pellicle is of interest because many salivary proteins exhibit functional domains that maintain the activities of the native protein. Among the in vivo-acquired enamel pellicle peptides that have been newly identified, 5 peptides are derived from statherin. Here, we assessed the ability of these statherin pellicle peptides to inhibit hydroxyapatite crystal growth. In addition, atomistic molecular dynamics (MD) simulations were performed to better understand the underlying physical mechanisms of hydroxyapatite growth inhibition. A microplate colorimetric assay was used to quantify hydroxyapatite growth. Statherin protein, 5 statherin-derived peptides, and a peptide lacking phosphate at residues 2 and 3 were analyzed. Statherin peptide phosphorylated on residues 2 and 3 indicated a significant inhibitory effect when compared with the 5 other peptides (P < 0.05). MD simulations showed a strong affinity and fast adsorption to hydroxyapatite for phosphopeptides, whereas unphosphorylated peptides interacted weakly with the hydroxyapatite. Our data suggest that the presence of a covalently linked phosphate group (at residues 2 and 3) in statherin peptides modulates the effect of hydroxyapatite growth inhibition. This study provides a mechanism to account for the composition and function of acquired enamel pellicle statherin peptides that will contribute as a base for the development of biologically stable and functional synthetic peptides for therapeutic use against dental caries and/or periodontal disease.

Modeling the behavior of confined colloidal particles under shear flow.

We investigate the behavior of colloidal suspensions with different volume fractions confined between parallel walls under a range of steady shears. We model the particles using molecular dynamics (MD) with full hydrodynamic interactions implemented through the use of a lattice-Boltzmann (LB) fluid. A quasi-2d ordering occurs in systems characterized by a coexistence of coupled layers with different densities, order, and granular temperature. We present a phase diagram in terms of shear and volume fraction for each layer, and demonstrate that particle exchange between layers is required for entering the disordered phase.

Effect of melatonin and cholesterol on the structure of DOPC and DPPC membranes.

The cell membrane plays an important role in the molecular mechanism of amyloid toxicity associated with Alzheimer's disease. The membrane's chemical composition and the incorporation of small molecules, such as melatonin and cholesterol, can alter its structure and physical properties, thereby affecting its interaction with amyloid peptides. Both melatonin and cholesterol have been recently linked to amyloid toxicity. Melatonin has been shown to have a protective role against amyloid toxicity. However, the underlying molecular mechanism of this protection is still not well understood, and cholesterol's role remains controversial. We used small-angle neutron diffraction (SAND) from oriented lipid multi-layers, small-angle neutron scattering (SANS) from unilamellar vesicles experiments and Molecular Dynamics (MD) simulations to elucidate non-specific interactions of melatonin and cholesterol with 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) model membranes. We conclude that melatonin decreases the thickness of both model membranes by disordering the lipid hydrocarbon chains, thus increasing membrane fluidity. This result is in stark contrast to the much accepted ordering effect induced by cholesterol, which causes membranes to thicken.

Simulations of micellization of sodium hexyl sulfate.

Micellization of the ionic surfactant sodium hexyl sulfate has been studied using atomistic explicit-solvent molecular dynamics simulations with and without excess NaCl or CaCl(2). Simulations were performed at surfactant loadings near the critical micellization concentration. Equilibrium micelle size distributions and estimates of the critical micellization concentration obtained from the simulations are in agreement with experimental data. In comparison to the sodium dodecyl sulfate surfactant, the shorter alkyl chain of sodium hexyl sulfate results in increased disorder of the micellar core and water exposure of the hydrocarbon tail groups. However, water and ions do not penetrate into the micellar core even for these weakly micellizing surfactants. Excess NaCl is observed to have a minor influence on the micelle structure but excess CaCl(2) induces drastic changes both in the structure and the dynamics of the micellar system. Furthermore, in the absence of excess salt, sodium hexyl sulfate forms predominantly spherical, disorganized aggregates but an increase in ionic strength drives an increase in aggregate size and leads to prolate aggregates.

Molecular dynamics study of prolyl oligopeptidase with inhibitor in binding cavity.

We used the crystal structure of prolyl oligopeptidase (POP) with bound Z-pro-prolinal (ZPP) inhibitor (Protein Data Bank (PDB) structure 1QFS) to perform an intensive molecular dynamics study of the POP-ZPP complex. We performed 100 ns of simulation with the hemiacetal bond, through which the ZPP is bound to the POP, removed in order to better investigate the binding cavity environment. From basic analysis, measuring the radius of gyration, root mean square deviation, solvent accessible surface area and definition of the secondary structure of protein, we determined that the protein structure is highly stable and maintains its structure over the entire simulation time. This demonstrates that such long time simulations can be performed without the protein structure losing stability. We found that water bridges and hydrogen bonds play a negligible role in binding the ZPP thus indicating the importance of the hemiacetal bond. The two domains of the protein are bound by a set of approximately 12 hydrogen bonds, specific to the particular POP protein.

Nonlinear driven response of a phase-field crystal in a periodic pinning potential.

We study numerically the phase diagram and the response under a driving force of the phase field crystal model for pinned lattice systems introduced recently for both one- and two-dimensional systems. The model describes the lattice system as a continuous density field in the presence of a periodic pinning potential, allowing for both elastic and plastic deformations of the lattice. We first present results for phase diagrams of the model in the absence of a driving force. The nonlinear response to a driving force on an initially pinned commensurate phase is then studied via overdamped dynamic equations of motion for different values of mismatch and pinning strengths. For large pinning strength the driven depinning transitions are continuous, and the sliding velocity varies with the force from the threshold with power-law exponents in agreement with analytical predictions. Transverse depinning transitions in the moving state are also found in two dimensions. Surprisingly, for sufficiently weak pinning potential we find a discontinuous depinning transition with hysteresis even in one dimension under overdamped dynamics. We also characterize structural changes of the system in some detail close to the depinning transition.

Molecular dynamics simulations of the enzyme catechol-O-methyltransferase: methodological issues.

Results from extensive 70 ns all-atom molecular dynamics simulations of catechol-O-methyltransferase (COMT) enzyme are reported. The simulations were performed with explicit TIP3P water and Mg2+ ions. Four different crystal structures of COMT, with and without different ligands, were used. These simulations are among the most extensive of their kind and as such served as a stability test for such simulations. On the methodological side we found that the initial energy minimization procedure may be a crucial step: particular hydrogen bonds may break, and this can initiate an irreversible loss of protein structure that becomes observable in longer time scales of the order of tens of nanoseconds. This has important implications for both molecular dynamics and quantum mechanics-molecular mechanics simulations.

Phase diagram and commensurate-incommensurate transitions in the phase field crystal model with an external pinning potential.

We study the phase diagram and the commensurate-incommensurate transitions in a phase field model of a two-dimensional crystal lattice in the presence of an external pinning potential. The model allows for both elastic and plastic deformations and provides a continuum description of lattice systems, such as for adsorbed atomic layers or two-dimensional vortex lattices. Analytically, a mode expansion analysis is used to determine the ground states and the commensurate-incommensurate transitions in the model as a function of the strength of the pinning potential and the lattice mismatch parameter. Numerical minimization of the corresponding free energy shows reasonable agreement with the analytical predictions and provides details on the topological defects in the transition region. We find that for small mismatch the transition is of first order, and it remains so for the largest values of mismatch studied here. Our results are consistent with results of simulations for atomistic models of adsorbed overlayers.

Cell aggregation: packing soft grains.

Cellular aggregates may be considered as collections of membrane enclosed units with a pressure difference between the internal and external liquid phases. Cells are kept together by membrane adhesion and/or confined space compression. Pattern formation and, in particular, intercellular spacing have important roles in controlling solvent diffusion within such aggregates. A physical approach is used to study generic aspects of cellular packings in a confined space. Average material properties are derived from the free energy. The appearance of penetrating intercellular void channels is found to be critically governed by the cell wall adhesion mechanisms during the formation of dense aggregates. A fully relaxed aggregate efficiently hinders solvent diffusion at high hydrostatic pressures, while a small fraction (approximately 0.1) of adhesion related packing frustration is sufficient for breaking such a blockage even at high a pressure.

Dielectrophoresis of nanocolloids: a molecular dynamics study.

Dielectrophoresis (DEP), the motion of polarizable particles in non-uniform electric fields, has become an important tool for the transport, separation, and characterization of microparticles in biomedical and nanoelectronics research. In this article we present, to our knowledge, the first molecular dynamics simulations of DEP of nanometer-sized colloidal particles. We introduce a simplified model for a polarizable nanoparticle, consisting of a large charged macroion and oppositely charged microions, in an explicit solvent. The model is then used to study DEP motion of the particle at different combinations of temperature and electric field strength. In accord with linear response theory, the particle drift velocities are shown to be proportional to the DEP force. Analysis of the colloid DEP mobility shows a clear time dependence, demonstrating the variation of friction under non-equilibrium. The time dependence of the mobility further results in an apparent weak variation of the DEP displacements with temperature.

Anomalously slow phase transitions in self-gravitating systems.

The kinetics of collapse and explosion transitions in microcanonical self-gravitating ensembles is analyzed. A system of point particles interacting via an attractive soft Coulomb potential and confined to a spherical container is considered. We observed that for 100-200 particles collapse takes 10(3) - 10(4) particle crossing times to complete; i.e., it is by two to three orders of magnitude slower than the velocity relaxation. In addition, it is found that the collapse time decreases rapidly with an increase of the soft-core radius. We found that such an anomalously long collapse time is caused by the slow energy exchange between a higher-temperature compact core and relatively cold diluted halo. The rate of energy exchange between the faster modes of the core particles and slower-moving particles of the halo is exponentially small in the ratio of the frequencies of these modes. As the soft-core radius increases and the typical core modes become slower, the ratio of core and halo frequencies decreases and the collapse accelerates. Implications for astrophysical systems and phase transition kinetics are discussed.

Collapses and explosions in self-gravitating systems.

Collapse and explosion (reverse to collapse) transitions in self-gravitating systems are studied by molecular dynamics simulations. A microcanonical ensemble of point particles confined to a spherical box is considered. The particles interact via an attractive soft Coulomb potential. It is observed that a collapse indeed takes place when the energy of the uniform state is set near or below the metastability-instability threshold (collapse energy) as predicted by the mean-field theory. Similarly, an explosion occurs when the energy of the core-halo state is increased above the explosion energy, where according to the mean-field predictions the core-halo state becomes unstable. For systems consisting of 125-500 particles, the collapse takes about 10(5) single-particle crossing times to complete, while a typical explosion is by an order of magnitude faster. A finite lifetime of metastable states is observed. It is also found that the mean-field description of the uniform and core-halo states is exact within the statistical uncertainty of the molecular dynamics data.

Molecular dynamics simulations of lipid bilayers: major artifacts due to truncating electrostatic interactions.

We study the influence of truncating the electrostatic interactions in a fully hydrated pure dipalmitoylphosphatidylcholine (DPPC) bilayer through 20 ns molecular dynamics simulations. The computations in which the electrostatic interactions were truncated are compared to similar simulations using the particle-mesh Ewald (PME) technique. All examined truncation distances (1.8-2.5 nm) lead to major effects on the bilayer properties, such as enhanced order of acyl chains together with decreased areas per lipid. The results obtained using PME, on the other hand, are consistent with experiments. These artifacts are interpreted in terms of radial distribution functions g(r) of molecules and molecular groups in the bilayer plane. Pronounced maxima or minima in g(r) appear exactly at the cutoff distance indicating that the truncation gives rise to artificial ordering between the polar phosphatidyl and choline groups of the DPPC molecules. In systems described using PME, such artificial ordering is not present.

Towards better integrators for dissipative particle dynamics simulations

Coarse-grained models that preserve hydrodynamics provide a natural approach to study collective properties of soft-matter systems. Here, we demonstrate that commonly used integration schemes in dissipative particle dynamics give rise to pronounced artifacts in physical quantities such as the compressibility and the diffusion coefficient. We assess the quality of these integration schemes, including variants based on a recently suggested self-consistent approach, and examine their relative performance. Implications of integrator-induced effects are discussed.