Ion permeation through the alpha-hemolysin channel: theoretical studies based on Brownian dynamics and Poisson-Nernst-Plank electrodiffusion theory

Identification of the molecular interaction governing ion conduction through biological pores is one of the most important goals of modern electrophysiology. Grand canonical Monte Carlo Brownian dynamics (GCMC/BD) and three-dimensional Poisson-Nernst-Plank (3d-PNP) electrodiffusion algorithms offer powerful and general approaches to study of ion permeation through wide molecular pores. A detailed analysis of ion flows through the staphylococcal alpha-hemolysin channel based on series of simulations at different concentrations and transmembrane potentials is presented. The position-dependent diffusion coefficient is approximated on the basis of a hydrodynamic model. The channel conductance calculated by GCMC/BD is approximately 10% higher than (electrophysiologically measured) experimental values, whereas results from 3d-PNP are always 30-50% larger. Both methods are able to capture all important electrostatic interactions in equilibrium conditions. The asymmetric conductance upon the polarity of the transmembrane potential observed experimentally is reproduced by GCMC/BD and 3d-PNP. The separation of geometrical and energetic influence of the channel on ion conduction reveals that such asymmetries arise from the permanent charge distribution inside the pore. The major determinant of the asymmetry is unbalanced charge in the triad of polar residues D127, D128, and K131. The GCMC/BD or 3d-PNP calculations reproduce also experimental reversal potentials and permeability rations in asymmetric ionic solutions. The weak anionic selectivity of the channel results from the presence of the salt bridge between E111 and K147 in the constriction zone. The calculations also reproduce the experimentally derived dependence of the reversible potential to the direction of the salt gradient. The origin of such effect arises from the asymmetrical distribution of energetic barriers along the channel axis, which modulates the preferential ion passage in different directions.     

Molecular determinants of permeation through the cation channel TRPV4

We have studied the molecular determinants of ion permeation through the TRPV4 channel (VRL-2, TRP12, VR-OAC, and OTRPC4). TRPV4 is characterized by both inward and outward rectification, voltage-dependent block by Ruthenium Red, a moderate selectivity for divalent versus monovalent cations, and an Eisenman IV permeability sequence. We identify two aspartate residues, Asp(672) and Asp(682), as important determinants of the Ca(2+) sensitivity of the TRPV4 pore. Neutralization of either aspartate to alanine caused a moderate reduction of the relative permeability for divalent cations and of the degree of outward rectification. Neutralizing both aspartates simultaneously caused a much stronger reduction of Ca(2+) permeability and channel rectification and additionally altered the permeability order for monovalent cations toward Eisenman sequence II or I. Moreover, neutralizing Asp(682) but not Asp(672) strongly reduces the affinity of the channel for Ruthenium Red. Mutations to Met(680), which is located at the center of a putative selectivity filter, strongly reduced whole cell current amplitude and impaired Ca(2+) permeation. In contrast, neutralizing the only positively charged residue in the putative pore region, Lys(675), had no obvious effects on the properties of the TRPV4 channel pore. Our findings delineate the pore region of TRPV4 and give a first insight into the possible architecture of its permeation pathway.     

Molecular determinants of permeation through the cation channel TRPM6

TRPM6 and its closest relative TRPM7 are members of the Transient Receptor Potential Melastatin (TRPM) subfamily of cation channels and are known to be Mg2+ permeable. By aligning the sequence of the putative TRPM6 pore with the pore sequences of the other subfamily members, we located in the loop between the fifth and the sixth transmembrane domain, a stretch of amino acids residues, 1028GEIDVC1033, as the potential selectivity filter. Two negatively charged residues, E1024 (conserved in TRPM6, TRPM7, TRPM1 and TRPM3) and D1031 (conserved along the entire TRPM subfamily), were identified as important determinants of cation permeation through TRPM6, because neutralization of both residues into an alanine resulted in non-functional channels. Neutralization of E1029 (conserved in TRPM6, TRPM7, TRPM4 and TRPM5) resulted in channels with increased conductance for Ba2+ and Zn2+, decreased ruthenium red sensitivity and larger pore diameter compared to wild-type TRPM6. Changing the residue I1030 into methionine, resulted in channels with lower conductance for Ni2+, decreased sensitivity to ruthenium red block and reduced pore diameter. Thus, these data demonstrate that amino acid residues E1024, I1030 and D1031 are important for channel function and that subtle amino acid variation in the pore region accounts for TRPM6 permeation properties.     

Dynamics and energetics of water permeation through the aquaporin channel

Structural properties of water inside bovine aquaporin-1 are investigated by molecular simulation. The calculations, which are based on the recently determined X-ray structure at 2.2 A resolution (Sui et al., Nature 2001;414:872-878), are carried out on one monomeric subunit immersed in a water-n-octane-water bilayer. Molecular dynamics (MD) simulations suggest that His182, a fully conserved residue in the channel pore, is protonated in the delta position. Furthermore, they reveal a highly ordered water structure in the channel, induced by the electrostatic properties of the protein. Multiple-steering MD simulations are used to calculate the free-energy of water diffusion. To the best of our knowledge, this represents the first free-energy calculation based on the new, high-resolution structure of the pore. The calculated barrier is 2.5 kcal/mol, and it is associated to water permeation through the Asn-Pro-Ala (NPA) region of the pore, where water molecules are only hydrogen-bonded with themselves. These findings are fully consistent with those based on the previous MD studies on the human protein (de Groot and Grubmüller, Science 2001;294:2353-2357).     

Molecular dynamics simulation of proton transport through the influenza A virus M2 channel

The structural and dynamical properties of a solvated proton in the influenza A virus M2 channel are studied using a molecular dynamics (MD) simulation technique. The second-generation multi-state empirical valence bond (MS-EVB2) model was used to describe the interaction between the excess proton and the channel environment. Solvation structures of the excess proton and its mobility characteristics along the channel were determined. It was found that the excess proton is capable of crossing the channel gate formed by the ring of four histidine residues even though the gate was only partially open. Although the hydronium ion itself did not cross the channel gate by traditional diffusion, the excess proton was able to transport through the ring of histidine residues by hopping between two water molecules located at the opposite sides of the gate. Our data also indicate that the proton diffusion through the channel may be correlated with the changes in channel conformations. To validate this observation, a separate simulation of the proton in a "frozen" channel has been conducted, which showed that the proton mobility becomes inhibited.     

Side-chain dynamics are critical for water permeation through aquaporin-1

Molecular dynamics simulations of aquaporin-1 embedded in a solvated lipid bilayer were carried out to investigate the mechanism of water permeation. The 2.2 Å resolution crystal structure of the bovine protein was used for five independent trajectories. During the equilibration and preparatory steps in which the protein was held fixed, water molecules inside the water channel adopted the same positions as observed in the crystal structure but they did not pass through the channel, suggesting that the dynamic motion of the protein is critical for water permeation. When the protein atoms were allowed to move, the side chains of the two asparagines in the two conserved Asn-Pro-Ala motifs near the center of the channel formed hydrogen bonds with water and helped water molecules move through the channel by actively aligning them for transport. The main-chain oxygen atoms, which were exposed to the pore surface in the crystal structure, also contributed to water transfer. Besides the constriction region observed in the crystal structure (Arg¹⁹⁷, Phe⁵⁸, His¹⁸², and Cys¹⁹¹), we found that His⁷⁶ and Val¹⁵⁵ act as a valve by dynamically blocking water permeation and helping control flow.     

Molecular dynamics simulation of a synthetic ion channel.

A molecular dynamics simulation has been performed on a synthetic membrane-spanning ion channel, consisting of four alpha-helical peptides, each of which is composed of the amino acids leucine (L) and serine (S), with the sequence Ac-(LSLLLSL)3-CONH2. This four-helix bundle has been shown experimentally to act as a proton-conducting channel in a membrane environment. In the present simulation, the channel was initially assembled as a parallel bundle in the octane portion of a phase-separated water/octane system, which provided a membrane-mimetic environment. An explicit reversible multiple-time-step integrator was used to generate a dynamical trajectory, a few nanoseconds in duration for this composite system on a parallel computer, under ambient conditions. After more than 1 ns, the four helices were found to adopt an associated dimer state with twofold symmetry, which evolved into a coiled-coil tetrameric structure with a left-handed twist. In the coiled-coil state, the polar serine side chains interact to form a layered structure with the core of the bundle filled with H2O. The dipoles of these H2O molecules tended to align opposite the net dipole of the peptide bundle. The calculated dipole relaxation function of the pore H2O molecules exhibits two reorientation times. One is approximately 3.2 ps, and the other is approximately 100 times longer. The diffusion coefficient of the pore H2O is about one-third of the bulk H2O value. The total dipole moment and the inertia tensor of the peptide bundle have been calculated and reveal slow (300 ps) collective oscillatory motions. Our results, which are based on a simple united atom force-field model, suggest that the function of this synthetic ion channel is likely inextricably coupled to its dynamical behavior.     

Molecular dynamics simulation study of the vanillate transport channel of Opdk

Most water-soluble and small molecules are taken up by substrate-specific channels belonging to the Gram-negative bacteria. The protein named OpdK, a member of OprD family, plays an important role in transporting the vanillate as the only carbon source of Pseudomonas aeruginosa. P. aeruginosa infections can be serious due to the high intrinsic antibiotic resistance owing to the present of the OprD family. We applied standard molecular dynamics, steered molecular dynamics and umbrella sampling, to investigate the thermodynamics of vanillate passing through the pore of OpdK protein at physiological temperature. The results indicate that hydrogen bonds of vanillate-L3 (L3: Gly92-Gln111) and hydrophobic interactions of vanillate-L7 (Gly252-Asn278) are crucial to the transport of vanillate. Compared to L7, L3 can hardly change the shape of the pore, but its amino acids can effectively affect the transport process. The important role of charged residues in the barrel of the protein for the substrate transport has been proved in the experiment researches and our simulations also determinate that these residues may prevent the vanillate from entering the cell. These results provide detailed information that will facilitate the development of effective drugs.     

Frequency response of alternating currents through the Staphylococcus aureus alpha-hemolysin ion channel

Alternating currents were measured through transmembrane ion channels formed by Staphylococcus aureus alpha-hemolysin proteins in planar bilayer membranes as part of an investigation to determine the channel's frequency response and the appropriateness of an equivalent circuit commonly used to model electrical interactions at the surface of cells. The experimental approach includes a novel method for separating the alternating current through one or more channels, which is conductive in nature, from the capacitively coupled current through the membrane. Separation of the conductive and capacitive alternating currents made it possible to measure the frequency response of the alpha-hemolysin channels. The results of the study are consistent with an equivalent circuit of a membrane capacitor in parallel with one or more channel resistors over the frequency range 30-120 Hz. The possible usefulness of frequency response data for ion channels in cell membranes during investigations of biological effects of time-varying magnetic fields is briefly discussed.     

Control of cation permeation through the nicotinic receptor channel

We used molecular dynamics (MD) simulations to explore the transport of single cations through the channel of the muscle nicotinic acetylcholine receptor (nAChR). Four MD simulations of 16 ns were performed at physiological and hyperpolarized membrane potentials, with and without restraints of the structure, but all without bound agonist. With the structure unrestrained and a potential of −100 mV, one cation traversed the channel during a transient period of channel hydration; at −200 mV, the channel was continuously hydrated and two cations traversed the channel. With the structure restrained, however, cations did not traverse the channel at either membrane potential, even though the channel was continuously hydrated. The overall results show that cation selective transport through the nAChR channel is governed by electrostatic interactions to achieve charge selectivity, but ion translocation relies on channel hydration, facilitated by a trans-membrane field, coupled with dynamic fluctuations of the channel structure.     

Molecular dynamics simulation of the antiamoebin ion channel: linking structure and conductance

Molecular-dynamics simulations were carried out to ascertain which of the potential multimeric forms of the transmembrane peptaibol channel, antiamoebin, is consistent with its measured conductance. Estimates of the conductance obtained through counting ions that cross the channel and by solving the Nernst-Planck equation yield consistent results, indicating that the motion of ions inside the channel can be satisfactorily described as diffusive. The calculated conductance of octameric channels is markedly higher than the conductance measured in single channel recordings, whereas the tetramer appears to be nonconducting. The conductance of the hexamer was estimated to be 115 ± 34 pS and 74 ± 20 pS, at 150 mV and 75 mV, respectively, in satisfactory agreement with the value of 90 pS measured at 75 mV. On this basis, we propose that the antiamoebin channel consists of six monomers. Its pore is large enough to accommodate K⁺ and Cl⁻ with their first solvation shells intact. The free energy barrier encountered by K⁺ is only 2.2 kcal/mol whereas Cl⁻ encounters a substantially higher barrier of nearly 5 kcal/mol. This difference makes the channel selective for cations. Ion crossing events are shown to be uncorrelated and follow Poisson statistics.     

Molecular dynamics studies of ion permeation in VDAC

The voltage-dependent anion channel (VDAC) in the outer membrane of mitochondria serves an essential role in the transport of metabolites and electrolytes between the cell matrix and mitochondria. To examine its structure, dynamics, and the mechanisms underlying its electrophysiological properties, we performed a total of 1.77 μs molecular dynamics simulations of human VDAC isoform 1 in DOPE/DOPC mixed bilayers in 1 M KCl solution with transmembrane potentials of 0, ±25, ±50, ±75, and ±100 mV. The calculated conductance and ion selectivity are in good agreement with the experimental measurements. In addition, ion density distributions inside the channel reveal possible pathways for different ion species. Based on these observations, a mechanism underlying the anion selectivity is proposed; both ion species are transported across the channel, but the rate for K⁺ is smaller than that for Cl⁻ because of the attractive interactions between K⁺ and residues on the channel wall. This difference leads to the anion selectivity of VDAC.     

1,3‐propanediol binds deep inside the channel to inhibit water permeation through aquaporins

Aquaporins and aquaglyceroporins (AQPs) are membrane channel proteins responsible for transport of water and for transport of glycerol in addition to water across the cell membrane, respectively. They are expressed throughout the human body and also in other forms of life. Inhibitors of human AQPs have been sought for therapeutic treatment for various medical conditions including hypertension, refractory edema, neurotoxic brain edema, and so forth. Conducting all‐atom molecular dynamics simulations, we computed the binding affinity of acetazolamide to human AQP4 that agrees closely with in vitro experiments. Using this validated computational method, we found that 1,3‐propanediol (PDO) binds deep inside the AQP4 channel to inhibit that particular aquaporin efficaciously. Furthermore, we used the same method to compute the affinities of PDO binding to four other AQPs and one aquaglyceroporin whose atomic coordinates are available from the protein data bank (PDB). For bovine AQP1, human AQP2, AQP4, AQP5, and Plasmodium falciparum PfAQP whose structures were resolved with high resolution, we obtained definitive predictions on the PDO dissociation constant. For human AQP1 whose PDB coordinates are less accurate, we estimated the dissociation constant with a rather large error bar. Taking into account the fact that PDO is generally recognized as safe by the US FDA, we predict that PDO can be an effective diuretic which directly modulates water flow through the protein channels. It should be free from the serious side effects associated with other diuretics that change the hydro‐homeostasis indirectly by altering the osmotic gradients.     

Water transport by the bacterial channel alpha-hemolysin

This study is an investigation of the ability of the bacterial channel alpha-hemolysin to facilitate water permeation across biological membranes. alpha-Hemolysin channels were incorporated into rabbit erythrocyte ghosts at varying concentrations, and water permeation was induced by mixing the ghosts with hypertonic sucrose solutions. The resulting volume decrease of the ghosts was followed by time-resolved optical absorption at pH 5, 6, and 7. The average single-channel permeability coefficient of alpha-hemolysin for water ranged between 1.3x10-12 cm/s and 1.5x10-12 cm/s, depending on pH. The slightly increased single-channel permeability coefficient at lower pH-values was attributed to an increase in the effective pore size. The activation energy of water transport through the channel was low (Ea=5.4 kcal/mol), suggesting that the properties of water inside the alpha-hemolysin channel resemble those of bulk water. This conclusion was supported by calculations based on macroscopic hydrodynamic laws of laminar water flow. Using the known three-dimensional structure of the channel, the calculations accurately predicted the rate of water flow through the channel. The latter finding also indicated that water permeation data can provide a good estimate of the pore size for large channels.     

Molecular simulation of permeation behaviour of ethanol/water molecules with single-layer graphene oxide membranes

Graphene oxide (GO)-based materials have shown promise as water-permeating membranes in pervaporation separation. However, the feed permeation and surface affinity of single-layer nanoporous GO sheet for liquid mixtures remain unresolved. Here, the pressure-driven molecular transport of pure ethanol and pure water, as well ethanol-water mixtures, crossing through single-layer nanoporous GO sheet was studied by non-equilibrium molecular dynamics simulations. We show that single-layer GO sheet with controlled pore sizes can effectively reject ethanol and allow water permeation with high permeability. This means that porous GO sheets could act as an effective dehydration membrane, therefore providing the initial barrier for ethanol passage in GO-based membrane. The pore size effect was considered as the separation mechanism. Both ethanol and water molecules in the mixture show comparable affinity with GO surfaces. The hydrogen-bonding coupling interaction between mixture and surface functional groups provide addition influence on the molecular transport through GO pores.     

Organic cation permeation through the channel formed by polycystin-2

Polycystin-2 (PC2), a member of the transient receptor potential family of ion channels (TRPP2), forms a calcium-permeable cation channel. Mutations in PC2 lead to polycystic kidney disease. From the primary sequence and by analogy with other channels in this family, PC2 is modeled to have six transmembrane domains. However, most of the structural features of PC2, such as how large the channel is and how many subunits make up the pore of the channel, are unknown. In this study, we estimated the pore size of PC2 from the permeation properties of the channel. Organic cations of increasing size were used as current carriers through the PC2 channel after PC2 was incorporated into lipid bilayers. We found that dimethylamine, triethylamine, tetraethylammonium, tetrabutylammonium, tetrapropylammonium, and tetrapentylammonium were permeable through the PC2 channel. The slope conductance of the PC2 channel decreased as the ionic diameter of the organic cation increased. For each organic cation tested, the currents were inhibited by gadolinium and anti-PC2 antibody. Using the dimensions of the largest permeant cation, the minimum pore diameter of the PC2 channel was estimated to be at least 11 A. The large pore size suggests that the primary state of this channel found in vivo is closed to avoid rundown of cation gradients across the plasma membrane and excessive calcium leak from endoplasmic reticulum stores.     

Molecular dynamics simulations of migration of ions and molecules through the acetylcholine receptor channel

A dynamic model of the channel of an acetylcholine receptor in a closed state has been proposed. The channel is formed by five a-helices of subunit M2 and stabilized by the cyclic hydrocarbon (CH2)105. The migration of charged and unchanged van der Waals particles with a diameter of 7.72 A equivalent to the diameter of a hydrated sodium ion has been studied. The migration occurred by the action of external force applied to the complex along the channel axis. In the closed state, the inhibition of ions is due to two components: electrostatic interaction and steric constraints. The van der Waals channel gate is formed by residues 13'-A-Val255, B-Val261, C-Val269, D-Val255, and E-Ile264, and the negatively changed residues occurring in the upper part of the channel have a great effect on ion selectivity.     

Molecular dynamics simulation of deoxy and carboxy murine neuroglobin in water

Globins are respiratory proteins that reversibly bind dioxygen and other small ligands at the iron of a heme prosthetic group. Hemoglobin and myoglobin are the most prominent members of this protein family. Unexpectedly a few years ago a new member was discovered and called neuroglobin (Ngb), being predominantly expressed in the brain. Ngb is a single polypeptide of 151 amino acids and despite the small sequence similarity with other globins, it displays the typical globin fold. Oxygen, nitric oxide, or carbon monoxide can displace the distal histidine which, in ferrous Ngb as well as in ferric Ngb, is bound to the iron, yielding a reversible adduct. Recent crystallographic data on carboxy Ngb show that binding of an exogenous ligand is associated to structural changes involving heme sliding and a topological reorganization of the internal cavities; in particular, the huge internal tunnel that connects the bulk with the active site, peculiar to Ngb, is heavily reorganized. We report the results of extended (90 ns) molecular dynamics simulations in water of ferrous deoxy and carboxy murine neuroglobin, which are both coordinated on the distal site, in the latter case by CO and in the former one by the distal His(64)(E7). The long timescale of the simulations allowed us to characterize the equilibrated protein dynamics and to compare protein structure and dynamical behavior coupled to the binding of an exogenous ligand. We have characterized the heme sliding motion, the topological reorganization of the internal cavities, the dynamics of the distal histidine, and particularly the conformational change of the CD loop, whose flexibility depends ligand binding.     

Molecular dynamics simulation of TMEM16A channel: Linking structure with gating

TMEM16A, the calcium-activated chloride channel, is broadly expressed and plays pivotal roles in diverse physiological processes. To understand the structural and functional relationships of TMEM16A, it is necessary to fully clarify the structural basis of the gating of the TMEM16A channel. Herein, we performed the protein electrostatic analysis and molecular dynamics simulation on the TMEM16A in the presence and absence of Ca²⁺. Data showed that the separation of TM4 and TM6 causes pore expansion, and Q646 may be a key residue for the formation of π-helix in the middle segment of TM6. Moreover, E705 was found to form a group of H-bond interactions with D554/K588/K645 below the hydrophobic gate to stabilize the closed conformation of the pore in the Ca²⁺-free state. Interestingly, in the Ca²⁺ bound state, the E705 side chain swings 100^o to serve as Ca²⁺-binding coordination and released K645. K645 is closer to the hydrophobic gate in the calcium-bound state, which facilitates the provision of electrostatic forces for chloride ions as the ions pass through the hydrophobic gate. Our findings provide the structural-based insights to understanding the mechanisms of gating of TMEM16A.     

Simulations of ion permeation through a potassium channel: molecular dynamics of KcsA in a phospholipid bilayer

Potassium channels enable K(+) ions to move passively across biological membranes. Multiple nanosecond-duration molecular dynamics simulations (total simulation time 5 ns) of a bacterial potassium channel (KcsA) embedded in a phospholipid bilayer reveal motions of ions, water, and protein. Comparison of simulations with and without K(+) ions indicate that the absence of ions destabilizes the structure of the selectivity filter. Within the selectivity filter, K(+) ions interact with the backbone (carbonyl) oxygens, and with the side-chain oxygen of T75. Concerted single-file motions of water molecules and K(+) ions within the selectivity filter of the channel occur on a 100-ps time scale. In a simulation with three K(+) ions (initially two in the filter and one in the cavity), the ion within the central cavity leaves the channel via its intracellular mouth after approximately 900 ps; within the cavity this ion interacts with the Ogamma atoms of two T107 side chains, revealing a favorable site within the otherwise hydrophobically lined cavity. Exit of this ion from the channel is enabled by a transient increase in the diameter of the intracellular mouth. Such "breathing" motions may form the molecular basis of channel gating.