Voltage-dependent hydration and conduction properties of the hydrophobic pore of the mechanosensitive channel of small conductance

A detailed picture of water and ion properties in small pores is important for understanding the behavior of biological ion channels. Several recent modeling studies have shown that small, hydrophobic pores exclude water and ions even if they are physically large enough to accommodate them, a mechanism called hydrophobic gating. This mechanism has been implicated in the gating of several channels, including the mechanosensitive channel of small conductance (MscS). Although the pore in the crystal structure of MscS is wide and was initially hypothesized to be open, it is lined by hydrophobic residues and may represent a nonconducting state. Molecular dynamics simulations were performed on MscS to determine whether or not the structure can conduct ions. Unlike previous simulations of hydrophobic nanopores, electric fields were applied to this system to model the transmembrane potential, which proved to be important. Although simulations without a potential resulted in a dehydrated, occluded pore, the application of a potential increased the hydration of the pore and resulted in current flow through the channel. The calculated channel conductance was in good agreement with experiment. Therefore, it is likely that the MscS crystal structure is closer to a conducting than a nonconducting state.     

Mechanosensitive channel MscS in the open state: modeling of the transition, explicit simulations, and experimental measurements of conductance

Mechanosensitive channels of small conductance (MscS) are ubiquitous turgor pressure regulators found in many walled cells and some intracellular organelles. Escherichia coli MscS acting as a tension-activated osmolyte release valve shows a nonsaturable conductance (1.2 nS in a 39 mS/cm electrolyte) and weak preference for anions. Pursuing the transition pathways in this channel, we applied the extrapolated motion protocol (cycles of displacements, minimizations, and short simulations) to the previously generated compact resting conformation of MscS. We observed tilting and straightening of the kinked pore-forming TM3 helices during the barrel expansion. Extended all-atom simulations confirmed the stability of the open conformation in the bilayer. A 53 degrees spontaneous axial rotation of TM3s observed after equilibration increased the width and polarity of the pore allowing for stable voltage-independent hydration and presence of both cations and anions throughout the pore. The resultant open state, characterized by a pore 1.6 nm wide, satisfied the experimental conductance and in-plane expansion. Applied transmembrane electric field (+/-100 to +/-200 mV) in simulations produced a flow of both K(+) and Cl(-), with Cl(-) current dominating at higher voltages. Electroosmotic water flux strongly correlated with the chloride current (approximately 8 waters per Cl(-)). The selectivity and rectification were in agreement with the experimental measurements performed in the same range of voltages. Among the charged residues surrounding the pore, only K169 was found to contribute noticeably in the rectification. We conclude that (a) the barrel expansion involving tilting, straightening, and rotation of TM3s provides the geometry and electrostatics that accounts for the conductive properties of the open pore; (b) the observed regimen of ion passage through the pore is similar to electrodiffusion, thus macroscopic estimations closely approximate the experimental and molecular dynamics-simulated conductances; (c) increased interaction of the opposing ionic fluxes at higher voltages may result in selectivities stronger than measured near the reversal potential.     

Ion conduction through MscS as determined by electrophysiology and simulation

The mechanosensitive channel of small conductance (MscS) is a membrane protein thought to act as a safety valve in bacteria, regulating the release of ions and small solutes when gated by membrane tension under challenging osmotic conditions. The influence of voltage on channel activation and the functional state depicted by the available crystal structure of MscS remain debated. Therefore, in an effort to relate electrophysiological measurements on MscS and properties of the MscS crystal conformation, we report here MscS's response to voltage and pressure as determined by patch-clamp experiments, as well as MscS electrostatics and transport properties as determined through all-atom molecular dynamics simulations of the protein embedded in a lipid bilayer, a 224,000-atom system. The experiments reveal that MscS is a slightly anion-selective channel with a conductance of approximately 1 ns, activated by pressure and inactivated in a voltage-dependent manner. On the other hand, the simulations, covering over 200 ns and including biasing electrostatic potentials, show that MscS restrained to the crystal conformation exhibits low conductance; unrestrained it increases the channel radius upon application of a large electrostatic bias and exhibits then ion conduction that matches experimentally determined conductances. The simulated conductance stems mainly from Cl- ions.     

Simulations of skin barrier function: free energies of hydrophobic and hydrophilic transmembrane pores in ceramide bilayers

Transmembrane pore formation is central to many biological processes such as ion transport, cell fusion, and viral infection. Furthermore, pore formation in the ceramide bilayers of the stratum corneum may be an important mechanism by which penetration enhancers such as dimethylsulfoxide (DMSO) weaken the barrier function of the skin. We have used the potential of mean constraint force (PMCF) method to calculate the free energy of pore formation in ceramide bilayers in both the innate gel phase and in the DMSO-induced fluidized state. Our simulations show that the fluid phase bilayers form archetypal water-filled hydrophilic pores similar to those observed in phospholipid bilayers. In contrast, the rigid gel-phase bilayers develop hydrophobic pores. At the relatively small pore diameters studied here, the hydrophobic pores are empty rather than filled with bulk water, suggesting that they do not compromise the barrier function of ceramide membranes. A phenomenological analysis suggests that these vapor pores are stable, below a critical radius, because the penalty of creating water-vapor and tail-vapor interfaces is lower than that of directly exposing the strongly hydrophobic tails to water. The PMCF free energy profile of the vapor pore supports this analysis. The simulations indicate that high DMSO concentrations drastically impair the barrier function of the skin by strongly reducing the free energy required for pore opening.     

Hydrophobic Gating of Ion Permeation in Magnesium Channel CorA

Ion channels catalyze ionic permeation across membranes via water-filled pores. To understand how changes in intracellular magnesium concentration regulate the influx of Mg²⁺ into cells, we examine early events in the relaxation of Mg²⁺ channel CorA toward its open state using massively-repeated molecular dynamics simulations conducted either with or without regulatory ions. The pore of CorA contains a 2-nm-long hydrophobic bottleneck which remained dehydrated in most simulations. However, rapid hydration or “wetting” events concurrent with small-amplitude fluctuations in pore diameter occurred spontaneously and reversibly. In the absence of regulatory ions, wetting transitions are more likely and include a wet state that is significantly more stable and more hydrated. The free energy profile for Mg²⁺ permeation presents a barrier whose magnitude is anticorrelated to pore diameter and the extent of hydrophobic hydration. These findings support an allosteric mechanism whereby wetting of a hydrophobic gate couples changes in intracellular magnesium concentration to the onset of ionic conduction.     

Genetic screen for potassium leaky small mechanosensitive channels (MscS) in Escherichia coli: recognition of cytoplasmic β domain as a new gating element

Mechanosensitive membrane channels in bacteria respond to the mechanical stretching of the membrane. They will open when bacteria are subjected to rapid osmotic down shock. MscS is a bacterial mechanosensitive channel of small conductance. It is a heptameric membrane protein whose transmembrane part, including the gate and its kinetics, has been well characterized. MscS has a large cytoplasmic domain of a cage-like shape that changes its conformation upon gating, but its involvement in gating is not understood. We screened MscS for mutations that cause potassium leak in Escherichia coli strains deficient in potassium transport systems. We did a phenotypic analysis of single and multiple mutants and recorded the single channel activities of some of them. After these analyses, we attributed the effects of a number of mutations to particular functional states of the channel. Our screen revealed that MscS leaks potassium in a desensitized and in an inactivated state. It also appeared that the lower part of TM3 (transmembrane, pore-forming helix) and the cytoplasmic β domain are tightly packed in the inactivated state but are dissociated in the open state. We attribute the TM3-β interaction to stabilization of the inactivated state in MscS and to the control of tight closure of its membrane pore.     

Defining the Role of the Tension Sensor in the Mechanosensitive Channel of Small Conductance

Mutations that alter the phenotypic behavior of the Escherichia coli mechanosensitive channel of small conductance (MscS) have been identified; however, most of these residues play critical roles in the transition between the closed and open states of the channel and are not directly involved in lipid interactions that transduce the tension response. In this study, we use molecular dynamic simulations to predict critical lipid interacting residues in the closed state of MscS. The physiological role of these residues was then investigated by performing osmotic downshock assays on MscS mutants where the lipid interacting residues were mutated to alanine. These experiments identified seven residues in the first and second transmembrane helices as lipid-sensing residues. The majority of these residues are hydrophobic amino acids located near the extracellular interface of the membrane. All of these residues interact strongly with the lipid bilayer in the closed state of MscS, but do not face the bilayer directly in structures associated with the open and desensitized states of the channel. Thus, the position of these residues relative to the lipid membrane appears related to the ability of the channel to sense tension in its different physiological states.     

Ion transport through membrane-spanning nanopores studied by molecular dynamics simulations and continuum electrostatics calculations

Narrow hydrophobic regions are a common feature of biological channels, with possible roles in ion-channel gating. We study the principles that govern ion transport through narrow hydrophobic membrane pores by molecular dynamics simulation of model membranes formed of hexagonally packed carbon nanotubes. We focus on the factors that determine the energetics of ion translocation through such nonpolar nanopores and compare the resulting free-energy barriers for pores with different diameters corresponding to the gating regions in closed and open forms of potassium channels. Our model system also allows us to compare the results from molecular dynamics simulations directly to continuum electrostatics calculations. Both simulations and continuum calculations show that subnanometer wide pores pose a huge free-energy barrier for ions, but a small increase in the pore diameter to approximately 1 nm nearly eliminates that barrier. We also find that in those wider channels the ion mobility is comparable to that in the bulk phase. By calculating local electrostatic potentials, we show that the long range Coulomb interactions of ions are strongly screened in the wide water-filled channels. Whereas continuum calculations capture the overall energetics reasonably well, the local water structure, which is not accounted for in this model, leads to interesting effects such as the preference of hydrated ions to move along the pore wall rather than through the center of the pore.     

Three-dimensional architecture of membrane-embedded MscS in the closed conformation

The mechanosensitive channel of small conductance (MscS) is part of a coordinated response to osmotic challenges in Escherichia coli. MscS opens as a result of membrane tension changes, thereby releasing small solutes and effectively acting as an osmotic safety valve. Both the functional state depicted by its crystal structure and its gating mechanism remain unclear. Here, we combine site-directed spin labeling, electron paramagnetic resonance spectroscopy, and molecular dynamics simulations with novel energy restraints based on experimental electron paramagnetic resonance data to investigate the native transmembrane (TM) and periplasmic molecular architecture of closed MscS in a lipid bilayer. In the closed conformation, MscS shows a more compact TM domain than in the crystal structure, characterized by a realignment of the TM segments towards the normal of the membrane. The previously unresolved NH₂-terminus forms a short helical hairpin capping the extracellular ends of TM1 and TM2 and is in close interaction with the bilayer interface. The present three-dimensional model of membrane-embedded MscS in the closed state represents a key step in determining the molecular mechanism of MscS gating.     

Gating transitions in bacterial ion channels measured at 3 microns resolution

Ion channels of high conductance (>200 pS) are widespread among prokaryotes and eukaryotes. Two examples, the Escherichia coli mechanosensitive ion channels Ec-MscS and Ec-MscL, pass currents of 125-300 pA. To resolve temporal details of conductance transitions, a patch-clamp setup was optimized for low-noise recordings at a time resolution of 3 microns (10-20 times faster than usual). Analyses of the high-resolution recordings confirm that Ec-MscL visits many subconductance states and show that most of the intersubstate transitions occur more slowly than the effective resolution of 3 micros. There is a clear trend toward longer transition times for the larger transitions. In Ec-MscS recordings, the majority of the observed full conductance transitions are also composite. We detected a short-lived (approximately 20 microns) Ec-MscS substate at 2/3 of full conductance; transitions between 2/3 and full conductance did not show fine structure and had a time course limited by the achieved resolution. Opening and closing transitions in MscS are symmetrical and are not preceded or followed by smaller, rapid currents ("anticipations" or "regrets"). Compared with other, lower-conductance channels, these measurements may detect unusually early states in the transitions from fully closed to fully open. Increased temporal resolution at the single-molecule level reveals that some elementary steps of structural transitions are composite and follow several alternative pathways, while others still escape resolution. High-bandwidth, low-noise single-channel measurements may provide details about state transitions in other high-conductance channels; and similar procedures may also be applied to channel- and nanopore-based single-molecule DNA measurements.     

The influence of geometry, surface character, and flexibility on the permeation of ions and water through biological pores

A hydrophobic constriction site can act as an efficient barrier to ion and water permeation if its diameter is less than the diameter of an ion's first hydration shell. This hydrophobic gating mechanism is thought to operate in a number of ion channels, e.g. the nicotinic receptor, bacterial mechanosensitive channels (MscL and MscS) and perhaps in some potassium channels (e.g. KcsA, MthK and KvAP). Simplified pore models allow one to investigate the primary characteristics of a conduction pathway, namely its geometry (shape, pore length, and radius), the chemical character of the pore wall surface, and its local flexibility and surface roughness. Our extended (about 0.1 micros) molecular dynamic simulations show that a short hydrophobic pore is closed to water for radii smaller than 0.45 nm. By increasing the polarity of the pore wall (and thus reducing its hydrophobicity) the transition radius can be decreased until for hydrophilic pores liquid water is stable down to a radius comparable to a water molecule's radius. Ions behave similarly but the transition from conducting to non-conducting pores is even steeper and occurs at a radius of 0.65 nm for hydrophobic pores. The presence of water vapour in a constriction zone indicates a barrier for ion permeation. A thermodynamic model can explain the behaviour of water in nanopores in terms of the surface tensions, which leads to a simple measure of 'hydrophobicity' in this context. Furthermore, increased local flexibility decreases the permeability of polar species. An increase in temperature has the same effect, and we hypothesize that both effects can be explained by a decrease in the effective solvent-surface attraction which in turn leads to an increase in the solvent-wall surface free energy.     

Structure, Dynamics, and Ion Conductance of the Phospholamban Pentamer

A 52-residue membrane protein, phospholamban (PLN) is an inhibitor of an adenosine-5′-triphosphate-driven calcium pump, the Ca²⁺-ATPase. Although the inhibition of Ca²⁺-ATPase involves PLN monomers, in a lipid bilayer membrane, PLN monomers form stable pentamers of unknown biological function. The recent NMR structure of a PLN pentamer depicts cytoplasmic helices extending normal to the bilayer in what is known as the bellflower conformation. The structure shows transmembrane helices forming a hydrophobic pore 4 Å in diameter, which is reminiscent of earlier reports of possible ion conductance through PLN pentamers. However, recent FRET measurements suggested an alternative structure for the PLN pentamer, known as the pinwheel model, which features a narrower transmembrane pore and cytoplasmic helices that lie against the bilayer. Here, we report on structural dynamics and conductance properties of the PLN pentamers from all-atom (AA) and coarse-grained (CG) molecular dynamics simulations. Our AA simulations of the bellflower model demonstrate that in a lipid bilayer membrane or a detergent micelle, the cytoplasmic helices undergo large structural fluctuations, whereas the transmembrane pore shrinks and becomes asymmetric. Similar asymmetry of the transmembrane region was observed in the AA simulations of the pinwheel model; the cytoplasmic helices remained in contact with the bilayer. Using the CG approach, structural dynamics of both models were investigated on a microsecond timescale. The cytoplasmic helices of the CG bellflower model were observed to fall against the bilayer, whereas in the CG pinwheel model the conformation of the cytoplasmic helices remained stable. Using steered molecular dynamics simulations, we investigated the feasibility of ion conductance through the pore of the bellflower model. The resulting approximate potentials of mean force indicate that the PLN pentamer is unlikely to function as an ion channel.     

Molecular dynamics simulations of water within models of ion channels.

The transbilayer pores formed by ion channel proteins contain extended columns of water molecules. The dynamic properties of such waters have been suggested to differ from those of water in its bulk state. Molecular dynamics simulations of ion channel models solvated within and at the mouths of their pores are used to investigate the dynamics and structure of intra-pore water. Three classes of channel model are investigated: a) parallel bundles of hydrophobic (Ala20) alpha-helices; b) eight-stranded hydrophobic (Ala10) antiparallel beta-barrels; and c) parallel bundles of amphipathic alpha-helices (namely, delta-toxin, alamethicin, and nicotinic acetylcholine receptor M2 helix). The self-diffusion coefficients of water molecules within the pores are reduced significantly relative to bulk water in all of the models. Water rotational reorientation rates are also reduced within the pores, particularly in those pores formed by alpha-helix bundles. In the narrowest pore (that of the Ala20 pentameric helix bundle) self-diffusion coefficients and reorientation rates of intra-pore waters are reduced by approximately an order of magnitude relative to bulk solvent. In Ala20 helix bundles the water dipoles orient antiparallel to the helix dipoles. Such dipole/dipole interaction between water and pore may explain how water-filled ion channels may be formed by hydrophobic helices. In the bundles of amphipathic helices the orientation of water dipoles is modulated by the presence of charged side chains. No preferential orientation of water dipoles relative to the pore axis is observed in the hydrophobic beta-barrel models.     

Brownian dynamics investigation into the conductance state of the MscS channel crystal structure

We suggest that the crystal structure of the mechanosensitive channel of small conductance is in a minimally conductive state rather than being fully activated. Performing Brownian dynamics simulations on the crystal structure show that no ions pass through it. When simulations are conducted on just the transmembrane domain (excluding the cytoplasmic residues 128 to 280) ions are seen to pass through the channel, but the conductance of ∼ 30 pS is well below experimentally measured values. The mutation L109S that replaces a pore lining hydrophobic residue with a polar one is found to have little effect on the conductance of the channel. Widening the hydrophobic region of the pore by 2.5 Å however, increases the channel conductance to over 200 pS suggesting that only a minimal conformational change is required to gate the pore.     

Ion channel gating: insights via molecular simulations

Ion channels are gated, i.e. they can switch conformation between a closed and an open state. Molecular dynamics simulations may be used to study the conformational dynamics of ion channels and of simple channel models. Simulations on model nanopores reveal that a narrow (<4 A) hydrophobic region can form a functionally closed gate in the channel and can be opened by either a small (approximately 1 A) increase in pore radius or an increase in polarity. Modelling and simulation studies confirm the importance of hydrophobic gating in K channels, and support a model in which hinge-bending of the pore-lining M2 (or S6 in Kv channels) helices underlies channel gating. Simulations of a simple outer membrane protein, OmpA, indicate that a gate may also be formed by interactions of charged side chains within a pore, as is also the case in ClC channels.     

The cytoplasmic cage domain of the mechanosensitive channel MscS is a sensor of macromolecular crowding

Cells actively regulate the macromolecular excluded volume of the cytoplasm to maintain the reciprocal fraction of free aqueous solution that is optimal for intracellular processes. However, the mechanisms whereby cells sense this critical parameter remain unclear. The mechanosensitive channel of small conductance (MscS channel), which is the major regulator of turgor in bacteria, mediates efflux of small osmolytes in response to increased membrane tension. At moderate sustained tensions produced by a decrease in external osmolarity, MscS undergoes slow adaptive inactivation; however, it inactivates abruptly in the presence of cytoplasmic crowding agents. To understand the mechanism underlying this rapid inactivation, we combined extrapolated and equilibrium molecular dynamics simulations with electrophysiological analyses of MscS mutants to explore possible transitions of MscS and generated models of the resting and inactivated states. Our models suggest that the coupling of the gate formed by TM3 helices to the peripheral TM1–TM2 pairs depends on the axial position of the core TM3 barrel relative to the TM1–TM2 shaft and the state of the associated hollow cytoplasmic domain (“cage”). They also indicate that the tension-driven inactivation transition separates the gate from the peripheral helices and promotes kinks in TM3s at G113 and that this conformation is stabilized by association of the TM3b segment with the β domain of the cage. We found that mutations destabilizing the TM3b–β interactions preclude inactivation and make the channel insensitive to crowding agents and voltage; mutations that strengthen this association result in a stable closed state and silent inactivation. Steered simulations showed that pressure exerted on the cage bottom in the inactivated state reduces the volume of the cage in the cytoplasm and at the same time increases the footprint of the transmembrane domain in the membrane, implying coupled sensitivity to both membrane tension and crowding pressure. The cage, therefore, provides feedback on the increasing crowding that disengages the gate and prevents excessive draining and condensation of the cytoplasm. We discuss the structural mechanics of cells surrounded by an elastic cell wall where this MscS-specific feedback mechanism may be necessary.