One-dimensional potential of mean force underestimates activation barrier for transport across flexible lipid membranes

Abstract

Transport of a fullerene-like nanoparticle across a lipid bilayer is investigated by coarse-grained molecular dynamics (MD) simulations. Potentials of mean force (PMF) acting on the nanoparticle in a flexible bilayer suspended in water and a bilayer restrained to a flat surface are computed by constrained MD simulations. The rate of the nanoparticle transport into the bilayer interior is predicted using one-dimensional Langevin models based on these PMFs. The predictions are compared with the transport rates obtained from a series of direct (unconstrained) MD simulations of the solute transport into the flexible bilayer. It is observed that the PMF acting on the solute in the flexible membrane underestimates the transport rate by more than an order of magnitude while the PMF acting on the solute in the restrained membrane yields an accurate estimate of the activation energy for transport into the flexible membrane. This paradox is explained by a coexistence of metastable membrane configurations for a range of the solute positions inside and near the flexible membrane. This leads to a significant reduction of the contribution of the transition state to the mean force acting on the solute. Restraining the membrane shape ensures that there is only one stable membrane configuration corresponding to each solute position and thus the transition state is adequately represented in the PMF. This mechanism is quite general and thus this phenomenon is expected to occur in a wide range of interfacial systems. A simple model for the free energy landscape of the coupled solute-membrane system is proposed and validated. This model explicitly accounts for effects of the membrane deformations on the solute transport and yields an accurate prediction of the activation energy for the solute transport.