Dataset Information

The Role of Lipid Interactions in Simulations of the ?-Hemolysin Ion-Channel-Forming Toxin.

ABSTRACT: Molecular dynamics simulations were performed to describe the function of the ion-channel-forming toxin ?-hemolysin (?HL) in lipid membranes that were composed of either 1,2-diphytanoyl-sn-glycero-3-phospho-choline or 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-choline. The simulations highlight the importance of lipid type in maintaining ?HL structure and function, enabling direct comparison to experiments for biosensing applications. We determined that although the two lipids studied are similar in structure, 1,2-diphytanoyl-sn-glycero-3-phospho-choline membranes better match the hydrophobic thickness of ?HL compared to 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-choline membranes. This hydrophobic match is essential to maintaining proper alignment of ?-sheet loops at the trans entrance of ?HL, which, when disrupted, creates an additional constriction to ion flow that decreases the channel current below experimental values and creates greater variability in channel conductance. Agreement with experiments was further improved with sufficient lipid membrane equilibration and allowed the discrimination of subtle ?HL conduction states with lipid type. Finally, we explore the effects of truncating the extramembrane cap of ?HL and its role in maintaining proper alignment of ?HL in the membrane and channel conductance. Our results demonstrate the essential role of lipid type and lipid-protein interactions in simulations of ?HL and will considerably improve the interpretation of experimental data.

SUBMITTER: Guros NB

PROVIDER: S-EPMC6225056 | biostudies-literature | 2018 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

The Role of Lipid Interactions in Simulations of the α-Hemolysin Ion-Channel-Forming Toxin.

Guros Nicholas B NB Balijepalli Arvind A Klauda Jeffery B JB

Biophysical journal 20180918 9

Molecular dynamics simulations were performed to describe the function of the ion-channel-forming toxin α-hemolysin (αHL) in lipid membranes that were composed of either 1,2-diphytanoyl-sn-glycero-3-phospho-choline or 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-choline. The simulations highlight the importance of lipid type in maintaining αHL structure and function, enabling direct comparison to experiments for biosensing applications. We determined that although the two lipids studied are similar ...[more]

PMID: 30287110

Similar Datasets

Project description:The association of the scorpion toxin Lq2 and a potassium ion (K(+)) channel has been studied using the Brownian dynamics (BD) simulation method. All of the 22 available structures of Lq2 in the Brookhaven Protein Data Bank (PDB) determined by NMR were considered during the simulation, which indicated that the conformation of Lq2 affects the binding between the two proteins significantly. Among the 22 structures of Lq2, only 4 structures dock in the binding site of the K(+) channel with a high probability and favorable electrostatic interactions. From the 4 candidates of the Lq2-K(+) channel binding models, we identified a good three-dimensional model of Lq2-K(+) channel complex through triplet contact analysis, electrostatic interaction energy estimation by BD simulation and structural refinement by molecular mechanics. Lq2 locates around the extracellular mouth of the K(+) channel and contacts the K(+) channel using its beta-sheet rather than its alpha-helix. Lys27, a conserved amino acid in the scorpion toxins, plugs the pore of the K(+) channel and forms three hydrogen bonds with the conserved residues Tyr78(A-C) and two hydrophobic contacts with Gly79 of the K(+) channel. In addition, eight hydrogen-bonds are formed between residues Arg25, Cys28, Lys31, Arg34 and Tyr36 of Lq2 and residues Pro55, Tyr78, Gly79, Asp80, and Tyr82 of K(+) channel. Many of them are formed by side chains of residues of Lq2 and backbone atoms of the K(+) channel. Thirteen hydrophobic contacts exist between residues Met29, Asn30, Lys31 and Tyr36 of Lq2 and residues Pro55, Ala58, Gly79, Asp80 and Tyr82 of the K(+) channel. These favorable interactions stabilize the association between the two proteins. These observations are in good agreement with the experimental results and can explain the binding phenomena between scorpion toxins and K(+) channels at the level of molecular structure. The consistency between the BD simulation and the experimental data indicates that our three-dimensional model of Lq2-K(+) channel complex is reasonable and can be used in further biological studies such as rational design of blocking agents of K(+) channels and mutagenesis in both toxins and K(+) channels.

Project description:Based on a homology model of the Kv1.3 potassium channel, the recognitions of the six scorpion toxins, viz. agitoxin2, charybdotoxin, kaliotoxin, margatoxin, noxiustoxin, and Pandinus toxin, to the human Kv1.3 potassium channel have been investigated by using an approach of the Brownian dynamics (BD) simulation integrating molecular dynamics (MD) simulation. Reasonable three-dimensional structures of the toxin-channel complexes have been obtained employing BD simulations and triplet contact analyses. All of the available structures of the six scorpion toxins in the Research Collaboratory for Structural Bioinformatics Protein Data Bank determined by NMR were considered during the simulation, which indicated that the conformations of the toxin significantly affect both the molecular recognition and binding energy between the two proteins. BD simulations predicted that all the six scorpion toxins in this study use their beta-sheets to bind to the extracellular entryway of the Kv1.3 channel, which is in line with the primary clues from the electrostatic interaction calculations and mutagenesis results. Additionally, the electrostatic interaction energies between the toxins and Kv1.3 channel correlate well with the binding affinities (-logK(d)s), R(2) = 0.603, suggesting that the electrostatic interaction is a dominant component for toxin-channel binding specificity. Most importantly, recognition residues and interaction contacts for the binding were identified. Lys-27 or Lys-28, residues Arg-24 or Arg-25 in the separate six toxins, and residues Tyr-400, Asp-402, His-404, Asp-386, and Gly-380 in each subunit of the Kv1.3 potassium channel, are the key residues for the toxin-channel recognitions. This is in agreement with the mutation results. MD simulations lasting 5 ns for the individual proteins and the toxin-channel complexes in a solvated lipid bilayer environment confirmed that the toxins are flexible and the channel is not flexible in the binding. The consistency between the results of the simulations and the experimental data indicated that our three-dimensional models of the toxin-channel complex are reasonable and can be used as a guide for future biological studies, such as the rational design of the blocking agents of the Kv1.3 channel and mutagenesis in both toxins and the Kv1.3 channel. Moreover, the simulation result demonstrates that the electrostatic interaction energies combined with the distribution frequencies from BD simulations might be used as criteria in ranking the binding configuration of a scorpion toxin to the Kv1.3 channel.

Project description: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.

Dataset Information

The Role of Lipid Interactions in Simulations of the ?-Hemolysin Ion-Channel-Forming Toxin.

Publications

The Role of Lipid Interactions in Simulations of the α-Hemolysin Ion-Channel-Forming Toxin.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets