Highlight

Calculation of lipid-bilayer permeabilities using an average force. Calculations of lipid bilayer permeabilities from first principles, using molecular simulations, would be valuable to rapidly assess the bioavailability of drug candidates, as well as to decipher, at the atomic level, the mechanisms that underlie the translocation of permeants. The most common theoretical approach, the solubility–diffusion model, requires determination of the free energy and the diffusivity as functions of the position of the permeant. Quantitative predictions of permeability have, however, been stymied by acute difficulties in calculating the diffusivity, inadequate sampling, and, most insidiously, systematic biases due to imperfections in the force field, simulation parameters, and the inherent limitations of the diffusive model. In the present work, we combine importance-sampling simulations employing an adaptive biasing force with a Bayesian-inference algorithm to determine the free energy and diffusivity with noteworthy precision and spatial resolution. In multimicrosecond simulations, we probe the sensitivity of the permeability estimates to different aspects of the methodology, including the truncation of short-range interactions, the thermostat, the force-field parameters of the permeant, the time scale over which the diffusivity is estimated, and the size of the simulated system. The force-field parameters and time scale dependence of the diffusivities impose the greatest uncertainties on the permeability estimates. Our simulations highlight the importance of membrane distortion due to the presence of the permeant, which may be partially suppressed if the bilayer patch is not large enough. We suggest that improvements to force fields and more robust kinetic models may be needed to reduce systematic errors below a factor of two. J. Chem. Theory Comput. 2014.

Recent publications

Dehez, F.; Delemotte, L.; Kramar, P.; Miklavcic, D.; Tarek, M.
Evidence of conducting hydrophobic nanopores across membranes in response to an electric field
J. Phys. Chem. C

2014, 118 (13), 6752-6757.
dx.doi.org

Liu, Y.; Chipot, C.; Shao, X.; Cai, W.
Threading or tumbling? Insight into the self-inclusion mechanism of an altro-α-cyclodextrin derivative
J. Phys. Chem. C

2014, 118 (33), 19380-19386.
dx.doi.org

Liu, P.; Chipot, C.; Cai, W.; Shao, X.
Unveiling the underlying mechanism for compression and decompression strokes of a molecular engine
J. Phys. Chem. C

2014, 118 (23), 12562-12567.
dx.doi.org

News

- Two positions for postdoctoral candidates are opened. For more information.
- An update of ParseFEP is available in the latest version of VMD.
- 新的分子动力学讲义 (Dissemination).

Contact

Laboratoire International Associé
CNRS-UIUC
 
Unité mixte de recherche n°7565
 
Université de Lorraine, B.P. 70239
54506 Vandoeuvre-lès-Nancy Cedex, France

 

Phone: +33.(0)3.83.68.40.97
Fax: +33.(0)3.83.68.43.87
 
How to reach us