Properties of lipid electropores I: Molecular dynamics simulations of stabilized pores by constant charge imbalance

Molecular dynamics (MD) simulations have become a powerful tool to study electroporation (EP) in atomic detail. In the last decade, numerous MD studies have been conducted to model the effect of pulsed electric fields on membranes, providing molecular models of the EP process of lipid bilayers. Here we extend these investigations by modeling for the first time conditions comparable to experiments using long (mu s-ms) low intensity (similar to KV/cm) pulses, by studying the characteristics of pores formed in lipid bilayers maintained at a constant surface tension and subject to constant charge imbalance. This enables the evaluation of structural (size) and electrical (conductance) properties of the pores formed, providing information hardly accessible directly by experiments. Extensive simulations of EP of simple phosphatidylcholine bilayers in 1 M NaCl show that hydrophilic pores with stable radii (1-2.5 nm) form under transmembrane voltages between 420 and 630 mV, allowing for ionic conductance in the range of 6.4-29.5 nS. We discuss in particular these findings and characterize both convergence and size effects in the MD simulations. We further extend these studies in a follow-up paper (Rems et al., Bioelectrochemistry, Submitted), by proposing an improved continuum model of pore conductance consistent with the results from the MD simulations.