Finiteness effects in computer simulation of fluids with spherical boundary conditions

Molecular dynamics simulations of a reversibly folding beta-heptapeptide in methanol: influence of the treatment of long-range electrostatic interactions

Eight 100-ns molecular dynamics simulations of a beta-heptapeptide in methanol at 340 K (within cubic periodic computational boxes of about 6-nm edge) are reported and compared. These simulations were performed with three different charge-state combinations at the peptide termini, one of them with or without a neutralizing chloride counterion, and using either the lattice-sum (LS) or reaction-field (RF) scheme to handle electrostatic interactions. The choice of the electrostatic scheme has essentially no influence on the folding-unfolding equilibrium when the peptide termini are uncharged and only a small influence when the peptide is positively charged at its N-terminus (with or without inclusion of a neutralizing chloride counterion). However, when the peptide is zwitterionic, the LS scheme leads to preferential sampling of the high-dipole folded helical state, whereas the RF scheme leads to preferential sampling of a low-dipole unfolded salt-bridged state. A continuum electrostatics analysis based on the sampled configurations (zwitterionic case) suggests that the LS scheme stabilizes the helical state through artificial periodicity, but that the magnitude of this perturbation is essentially negligible (compared to the thermal energy) for the large box size and relatively polar solvent considered. The results thus provide clear evidence (continuum electrostatics analysis) for the absence of LS artifacts and some indications (still not definitive because of the limited sampling of the folding-unfolding transition) for the presence of RF artifacts in this specific system.

Effect of geometrical confinement on the interaction between charged colloidal suspensions

The effective interaction between charged colloidal particles confined between two planar like-charged walls is investigated using computer simulations of the primitive model describing asymmetric electrolytes. In detail, we calculate the effective force acting onto a single macroion and onto a macroion pair in the presence of slitlike confinement. For moderate Coulomb coupling, we find that this force is repulsive. Under strong-coupling conditions, however, the sign of the force depends on the distance to the plates and on the interparticle distance. In particular, the particle-plate interaction becomes strongly attractive for small distances which may explain the occurrence of colloidal crystalline layers near the plates observed in recent experiments.