We report a Molecular Dynamics study of homogenous bubble nucleation in a Lennard-Jones fluid. The rate of bubble nucleation is estimated using forward-flux sampling (FFS). We find that cavitation starts with compact bubbles rather than with ramified structures as had been suggested by Shen and Debenedetti ( J. Chem. Phys. 111, 3581 (1999)). Our estimate of the bubblenucleation rate is higher than predicted on the basis of Classical Nucleation Theory (CNT). Our simulations show that local temperature fluctuations correlate strongly with subsequent bubble formation -this mechanism is not taken into account in CNT.2

We present an implicit solvent coarse-grained (CG) model for quantitative simulations of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. The absence of explicit solvent enables membrane simulations on large length and time scales at moderate computational expense. Despite improved computational efficiency, the model preserves chemical specificity and quantitative accuracy. The bonded and nonbonded interactions together with the effective cohesion mimicking the hydrophobic effect were systematically tuned by matching structural and mechanical properties from experiments and all-atom bilayer simulations, such as saturated area per lipid, radial distribution functions, density and pressure profiles across the bilayer, P2 order, etc. The CG lipid model is shown to self-assemble into a bilayer starting from a random dispersion. Its line tension and elastic properties, such as bending and stretching modulus, are semiquantitatively consistent with experiments. The effects of (i) reduced molecular friction and (ii) more efficient integration combine to an overall speed-up of 3−4 orders of magnitude compared to all-atom bilayer simulations. Our CG lipid model is especially useful for studies of large-scale phenomena in membranes that nevertheless require a fair description of chemical specificity, e.g., membrane patches interacting with movable and transformable membrane proteins and peptides.

Interfacial systems are at the core of fascinating phenomena in many disciplines, such as biochemistry, soft-matter physics, and food science. However, the parametrization of accurate, reliable, and consistent coarse-grained (CG) models for systems at interfaces remains a challenging endeavor. In the present work, we explore to what extent two independently developed solvent-free CG models of peptides and lipids-of different mapping schemes, parametrization methods, target functions, and validation criteria-can be combined by only tuning the cross-interactions. Our results show that the cross-parametrization can reproduce a number of structural properties of membrane peptides (for example, tilt and hydrophobic mismatch), in agreement with existing peptide-lipid CG force fields. We find encouraging results for two challenging biophysical problems: (i) membrane pore formation mediated by the cooperative action of several antimicrobial peptides, and (ii) the insertion and folding of the helix-forming peptide WALP23 in the membrane. © 2014 AIP Publishing LLC.

We present a simple, implicit-solvent model for fluid bilayer membranes. The model was designed to reproduce the elastic properties of real bilayer membranes. For this model, we observed the solid-fluid transition and studied the in-plane diffusivity of the fluid phase. As a test, we compute the elastic-bending and area-compressing moduli of fluid bilayer membranes. We find that the computed elastic properties are consistent with the available experimental data.