Digging into lipid membrane permeation for cardiac ion channel blocker d-sotalol with all-atom simulations

Kevin R. DeMarco, Slava Bekker, Colleen E Clancy, Sergei Y. Noskov, Igor Vorobyov

Research output: Contribution to journalArticle

8 Citations (Scopus)

Abstract

Interactions of drug molecules with lipid membranes play crucial role in their accessibility of cellular targets and can be an important predictor of their therapeutic and safety profiles. Very little is known about spatial localization of various drugs in the lipid bilayers, their active form (ionization state) or translocation rates and therefore potency to bind to different sites in membrane proteins. All-atom molecular simulations may help to map drug partitioning kinetics and thermodynamics, thus providing in-depth assessment of drug lipophilicity. As a proof of principle, we evaluated extensively lipid membrane partitioning of d-sotalol, well-known blocker of a cardiac potassium channel Kv11.1 encoded by the hERG gene, with reported substantial proclivity for arrhythmogenesis. We developed the positively charged (cationic) and neutral d-sotalol models, compatible with the biomolecular CHARMM force field, and subjected them to all-atom molecular dynamics (MD) simulations of drug partitioning through hydrated lipid membranes, aiming to elucidate thermodynamics and kinetics of their translocation and thus putative propensities for hydrophobic and aqueous hERG access. We found that only a neutral form of d-sotalol accumulates in the membrane interior and can move across the bilayer within millisecond time scale, and can be relevant to a lipophilic channel access. The computed water-membrane partitioning coefficient for this form is in good agreement with experiment. There is a large energetic barrier for a cationic form of the drug, dominant in water, to cross the membrane, resulting in slow membrane translocation kinetics. However, this form of the drug can be important for an aqueous access pathway through the intracellular gate of hERG. This route will likely occur after a neutral form of a drug crosses the membrane and subsequently re-protonates. Our study serves to demonstrate a first step toward a framework for multi-scale in silico safety pharmacology, and identifies some of the challenges that lie therein.

Original languageEnglish (US)
Article number26
JournalFrontiers in Pharmacology
Volume9
Issue numberFEB
DOIs
StatePublished - Feb 1 2018

Fingerprint

Sotalol
Membrane Lipids
Ion Channels
Membranes
Pharmaceutical Preparations
Thermodynamics
Safety
Water
Potassium Channels
Lipid Bilayers
Molecular Dynamics Simulation
Drug Interactions
Computer Simulation
Membrane Proteins
Pharmacokinetics
Pharmacology
Genes

Keywords

  • Cardiotoxicity
  • CHARMM force field
  • hERG
  • Lipophilicity
  • Long QT syndrome
  • Molecular dynamics
  • Umbrella sampling
  • Water-membrane partitioning

ASJC Scopus subject areas

  • Pharmacology
  • Pharmacology (medical)

Cite this

Digging into lipid membrane permeation for cardiac ion channel blocker d-sotalol with all-atom simulations. / DeMarco, Kevin R.; Bekker, Slava; Clancy, Colleen E; Noskov, Sergei Y.; Vorobyov, Igor.

In: Frontiers in Pharmacology, Vol. 9, No. FEB, 26, 01.02.2018.

Research output: Contribution to journalArticle

@article{3bb888a3691b43ee89a0fd52fe49f45b,
title = "Digging into lipid membrane permeation for cardiac ion channel blocker d-sotalol with all-atom simulations",
abstract = "Interactions of drug molecules with lipid membranes play crucial role in their accessibility of cellular targets and can be an important predictor of their therapeutic and safety profiles. Very little is known about spatial localization of various drugs in the lipid bilayers, their active form (ionization state) or translocation rates and therefore potency to bind to different sites in membrane proteins. All-atom molecular simulations may help to map drug partitioning kinetics and thermodynamics, thus providing in-depth assessment of drug lipophilicity. As a proof of principle, we evaluated extensively lipid membrane partitioning of d-sotalol, well-known blocker of a cardiac potassium channel Kv11.1 encoded by the hERG gene, with reported substantial proclivity for arrhythmogenesis. We developed the positively charged (cationic) and neutral d-sotalol models, compatible with the biomolecular CHARMM force field, and subjected them to all-atom molecular dynamics (MD) simulations of drug partitioning through hydrated lipid membranes, aiming to elucidate thermodynamics and kinetics of their translocation and thus putative propensities for hydrophobic and aqueous hERG access. We found that only a neutral form of d-sotalol accumulates in the membrane interior and can move across the bilayer within millisecond time scale, and can be relevant to a lipophilic channel access. The computed water-membrane partitioning coefficient for this form is in good agreement with experiment. There is a large energetic barrier for a cationic form of the drug, dominant in water, to cross the membrane, resulting in slow membrane translocation kinetics. However, this form of the drug can be important for an aqueous access pathway through the intracellular gate of hERG. This route will likely occur after a neutral form of a drug crosses the membrane and subsequently re-protonates. Our study serves to demonstrate a first step toward a framework for multi-scale in silico safety pharmacology, and identifies some of the challenges that lie therein.",
keywords = "Cardiotoxicity, CHARMM force field, hERG, Lipophilicity, Long QT syndrome, Molecular dynamics, Umbrella sampling, Water-membrane partitioning",
author = "DeMarco, {Kevin R.} and Slava Bekker and Clancy, {Colleen E} and Noskov, {Sergei Y.} and Igor Vorobyov",
year = "2018",
month = "2",
day = "1",
doi = "10.3389/fphar.2018.00026",
language = "English (US)",
volume = "9",
journal = "Frontiers in Pharmacology",
issn = "1663-9812",
publisher = "Frontiers Media S. A.",
number = "FEB",

}

TY - JOUR

T1 - Digging into lipid membrane permeation for cardiac ion channel blocker d-sotalol with all-atom simulations

AU - DeMarco, Kevin R.

AU - Bekker, Slava

AU - Clancy, Colleen E

AU - Noskov, Sergei Y.

AU - Vorobyov, Igor

PY - 2018/2/1

Y1 - 2018/2/1

N2 - Interactions of drug molecules with lipid membranes play crucial role in their accessibility of cellular targets and can be an important predictor of their therapeutic and safety profiles. Very little is known about spatial localization of various drugs in the lipid bilayers, their active form (ionization state) or translocation rates and therefore potency to bind to different sites in membrane proteins. All-atom molecular simulations may help to map drug partitioning kinetics and thermodynamics, thus providing in-depth assessment of drug lipophilicity. As a proof of principle, we evaluated extensively lipid membrane partitioning of d-sotalol, well-known blocker of a cardiac potassium channel Kv11.1 encoded by the hERG gene, with reported substantial proclivity for arrhythmogenesis. We developed the positively charged (cationic) and neutral d-sotalol models, compatible with the biomolecular CHARMM force field, and subjected them to all-atom molecular dynamics (MD) simulations of drug partitioning through hydrated lipid membranes, aiming to elucidate thermodynamics and kinetics of their translocation and thus putative propensities for hydrophobic and aqueous hERG access. We found that only a neutral form of d-sotalol accumulates in the membrane interior and can move across the bilayer within millisecond time scale, and can be relevant to a lipophilic channel access. The computed water-membrane partitioning coefficient for this form is in good agreement with experiment. There is a large energetic barrier for a cationic form of the drug, dominant in water, to cross the membrane, resulting in slow membrane translocation kinetics. However, this form of the drug can be important for an aqueous access pathway through the intracellular gate of hERG. This route will likely occur after a neutral form of a drug crosses the membrane and subsequently re-protonates. Our study serves to demonstrate a first step toward a framework for multi-scale in silico safety pharmacology, and identifies some of the challenges that lie therein.

AB - Interactions of drug molecules with lipid membranes play crucial role in their accessibility of cellular targets and can be an important predictor of their therapeutic and safety profiles. Very little is known about spatial localization of various drugs in the lipid bilayers, their active form (ionization state) or translocation rates and therefore potency to bind to different sites in membrane proteins. All-atom molecular simulations may help to map drug partitioning kinetics and thermodynamics, thus providing in-depth assessment of drug lipophilicity. As a proof of principle, we evaluated extensively lipid membrane partitioning of d-sotalol, well-known blocker of a cardiac potassium channel Kv11.1 encoded by the hERG gene, with reported substantial proclivity for arrhythmogenesis. We developed the positively charged (cationic) and neutral d-sotalol models, compatible with the biomolecular CHARMM force field, and subjected them to all-atom molecular dynamics (MD) simulations of drug partitioning through hydrated lipid membranes, aiming to elucidate thermodynamics and kinetics of their translocation and thus putative propensities for hydrophobic and aqueous hERG access. We found that only a neutral form of d-sotalol accumulates in the membrane interior and can move across the bilayer within millisecond time scale, and can be relevant to a lipophilic channel access. The computed water-membrane partitioning coefficient for this form is in good agreement with experiment. There is a large energetic barrier for a cationic form of the drug, dominant in water, to cross the membrane, resulting in slow membrane translocation kinetics. However, this form of the drug can be important for an aqueous access pathway through the intracellular gate of hERG. This route will likely occur after a neutral form of a drug crosses the membrane and subsequently re-protonates. Our study serves to demonstrate a first step toward a framework for multi-scale in silico safety pharmacology, and identifies some of the challenges that lie therein.

KW - Cardiotoxicity

KW - CHARMM force field

KW - hERG

KW - Lipophilicity

KW - Long QT syndrome

KW - Molecular dynamics

KW - Umbrella sampling

KW - Water-membrane partitioning

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

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

U2 - 10.3389/fphar.2018.00026

DO - 10.3389/fphar.2018.00026

M3 - Article

AN - SCOPUS:85041850092

VL - 9

JO - Frontiers in Pharmacology

JF - Frontiers in Pharmacology

SN - 1663-9812

IS - FEB

M1 - 26

ER -