Diffusive models of membrane permeation with explicit orientational freedom

Jeffrey Comer, Klaus J Schulten, Christophe Chipot

Research output: Contribution to journalArticle

Abstract

Accurate calculation of permeabilities from first-principles has been a long-standing challenge for computer simulations, notably in the context of drug discovery, as a route to predict the propensity of small, organic molecules to spontaneously translocate biological membranes. Of equal importance is the understanding of the permeation process and the pathway followed by the permeant from the aqueous medium to the interior of the lipid bilayer, and back out again. A convenient framework for the computation of permeabilities is provided by the solubility-diffusion model, which requires knowledge of the underlying free-energy and diffusivity landscapes. Here, we develop a formalism that includes an explicit description of the orientational motion of the solute as it diffuses across the membrane. Toward this end, we have generalized a recently proposed method that reconciles thermodynamics and kinetics in importance-sampling simulations by means of a Bayesian-inference scheme to reverse-solve the underlying Smoluchowski equation. Performance of the proposed formalism is examined in the model cases of a water and an ethanol molecule crossing a fully hydrated lipid bilayer. Our analysis reveals a conspicuous dependence of the free-energy and rotational diffusivity on the orientation of ethanol when it lies within the headgroup region of the bilayer. Specifically, orientations for which the hydroxyl group lies among the polar lipid head groups, while the ethyl group recedes toward the hydrophobic interior are associated with free-energy minima ∼2kBT deep, as well as significantly slower orientational kinetics compared to the bulk solution or the core of the bilayer. The conspicuous orientational anisotropy of ethanol at the aqueous interface is suggestive of a complete rotation of the permeant as it crosses the hydrophobic interior of the membrane.

Original languageEnglish (US)
Pages (from-to)2710-2718
Number of pages9
JournalJournal of Chemical Theory and Computation
Volume10
Issue number7
DOIs
StatePublished - Jul 8 2014

Fingerprint

Permeation
Free energy
lipids
Lipid bilayers
Ethanol
ethyl alcohol
free energy
membranes
Membranes
diffusivity
permeability
formalism
Biological membranes
Importance sampling
Molecules
Kinetics
kinetics
inference
Hydroxyl Radical
Lipids

ASJC Scopus subject areas

  • Computer Science Applications
  • Physical and Theoretical Chemistry

Cite this

Diffusive models of membrane permeation with explicit orientational freedom. / Comer, Jeffrey; Schulten, Klaus J; Chipot, Christophe.

In: Journal of Chemical Theory and Computation, Vol. 10, No. 7, 08.07.2014, p. 2710-2718.

Research output: Contribution to journalArticle

Comer, Jeffrey ; Schulten, Klaus J ; Chipot, Christophe. / Diffusive models of membrane permeation with explicit orientational freedom. In: Journal of Chemical Theory and Computation. 2014 ; Vol. 10, No. 7. pp. 2710-2718.
@article{acc4cfd9899244e39f9950f376a20a52,
title = "Diffusive models of membrane permeation with explicit orientational freedom",
abstract = "Accurate calculation of permeabilities from first-principles has been a long-standing challenge for computer simulations, notably in the context of drug discovery, as a route to predict the propensity of small, organic molecules to spontaneously translocate biological membranes. Of equal importance is the understanding of the permeation process and the pathway followed by the permeant from the aqueous medium to the interior of the lipid bilayer, and back out again. A convenient framework for the computation of permeabilities is provided by the solubility-diffusion model, which requires knowledge of the underlying free-energy and diffusivity landscapes. Here, we develop a formalism that includes an explicit description of the orientational motion of the solute as it diffuses across the membrane. Toward this end, we have generalized a recently proposed method that reconciles thermodynamics and kinetics in importance-sampling simulations by means of a Bayesian-inference scheme to reverse-solve the underlying Smoluchowski equation. Performance of the proposed formalism is examined in the model cases of a water and an ethanol molecule crossing a fully hydrated lipid bilayer. Our analysis reveals a conspicuous dependence of the free-energy and rotational diffusivity on the orientation of ethanol when it lies within the headgroup region of the bilayer. Specifically, orientations for which the hydroxyl group lies among the polar lipid head groups, while the ethyl group recedes toward the hydrophobic interior are associated with free-energy minima ∼2kBT deep, as well as significantly slower orientational kinetics compared to the bulk solution or the core of the bilayer. The conspicuous orientational anisotropy of ethanol at the aqueous interface is suggestive of a complete rotation of the permeant as it crosses the hydrophobic interior of the membrane.",
author = "Jeffrey Comer and Schulten, {Klaus J} and Christophe Chipot",
year = "2014",
month = "7",
day = "8",
doi = "10.1021/ct500209j",
language = "English (US)",
volume = "10",
pages = "2710--2718",
journal = "Journal of Chemical Theory and Computation",
issn = "1549-9618",
publisher = "American Chemical Society",
number = "7",

}

TY - JOUR

T1 - Diffusive models of membrane permeation with explicit orientational freedom

AU - Comer, Jeffrey

AU - Schulten, Klaus J

AU - Chipot, Christophe

PY - 2014/7/8

Y1 - 2014/7/8

N2 - Accurate calculation of permeabilities from first-principles has been a long-standing challenge for computer simulations, notably in the context of drug discovery, as a route to predict the propensity of small, organic molecules to spontaneously translocate biological membranes. Of equal importance is the understanding of the permeation process and the pathway followed by the permeant from the aqueous medium to the interior of the lipid bilayer, and back out again. A convenient framework for the computation of permeabilities is provided by the solubility-diffusion model, which requires knowledge of the underlying free-energy and diffusivity landscapes. Here, we develop a formalism that includes an explicit description of the orientational motion of the solute as it diffuses across the membrane. Toward this end, we have generalized a recently proposed method that reconciles thermodynamics and kinetics in importance-sampling simulations by means of a Bayesian-inference scheme to reverse-solve the underlying Smoluchowski equation. Performance of the proposed formalism is examined in the model cases of a water and an ethanol molecule crossing a fully hydrated lipid bilayer. Our analysis reveals a conspicuous dependence of the free-energy and rotational diffusivity on the orientation of ethanol when it lies within the headgroup region of the bilayer. Specifically, orientations for which the hydroxyl group lies among the polar lipid head groups, while the ethyl group recedes toward the hydrophobic interior are associated with free-energy minima ∼2kBT deep, as well as significantly slower orientational kinetics compared to the bulk solution or the core of the bilayer. The conspicuous orientational anisotropy of ethanol at the aqueous interface is suggestive of a complete rotation of the permeant as it crosses the hydrophobic interior of the membrane.

AB - Accurate calculation of permeabilities from first-principles has been a long-standing challenge for computer simulations, notably in the context of drug discovery, as a route to predict the propensity of small, organic molecules to spontaneously translocate biological membranes. Of equal importance is the understanding of the permeation process and the pathway followed by the permeant from the aqueous medium to the interior of the lipid bilayer, and back out again. A convenient framework for the computation of permeabilities is provided by the solubility-diffusion model, which requires knowledge of the underlying free-energy and diffusivity landscapes. Here, we develop a formalism that includes an explicit description of the orientational motion of the solute as it diffuses across the membrane. Toward this end, we have generalized a recently proposed method that reconciles thermodynamics and kinetics in importance-sampling simulations by means of a Bayesian-inference scheme to reverse-solve the underlying Smoluchowski equation. Performance of the proposed formalism is examined in the model cases of a water and an ethanol molecule crossing a fully hydrated lipid bilayer. Our analysis reveals a conspicuous dependence of the free-energy and rotational diffusivity on the orientation of ethanol when it lies within the headgroup region of the bilayer. Specifically, orientations for which the hydroxyl group lies among the polar lipid head groups, while the ethyl group recedes toward the hydrophobic interior are associated with free-energy minima ∼2kBT deep, as well as significantly slower orientational kinetics compared to the bulk solution or the core of the bilayer. The conspicuous orientational anisotropy of ethanol at the aqueous interface is suggestive of a complete rotation of the permeant as it crosses the hydrophobic interior of the membrane.

UR - http://www.scopus.com/inward/record.url?scp=84904108451&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84904108451&partnerID=8YFLogxK

U2 - 10.1021/ct500209j

DO - 10.1021/ct500209j

M3 - Article

AN - SCOPUS:84904108451

VL - 10

SP - 2710

EP - 2718

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

IS - 7

ER -