Gain-of-function mutations reveal expanded intermediate states and a sequential action of two gates in MscL

The tension-driven gating transition in the large mechanosensitive channel MscL proceeds through detectable states of intermediate conductance. Gain-of-function (GOF) mutants with polar or charged substitutions in the main hydrophobic gate display altered patterns of subconducting states, providing valuable information about gating intermediates. Here we present thermodynamic analysis of several GOF mutants to clarify the nature and position of low-conducting conformations in the transition pathway. Unlike wild-type (WT) MscL, which predominantly occupies the closed and fully open states with very brief substates, the mild V23T GOF mutant frequently visits a multitude of short-lived subconducting states. Severe mutants V23D and G22N open in sequence: closed (C) --> low-conducting substate (S) --> open (O), with the first subtransition occurring at lower tensions. Analyses of equilibrium state occupancies as functions of membrane tension show that the C-->S subtransition in WT MscL is associated with only a minor conductance increment, but the largest in-plane expansion and free energy change. The GOF substitutions strongly affect the first subtransition by reducing area ((Delta)A) and energy ((Delta)E) changes between C and S states commensurably with the severity of mutation. GOF mutants also exhibited a considerably larger (Delta)E associated with the second (S-->O) subtransition, but a (Delta)A similar to WT. The area changes indicate that closed conformations of GOF mutants are physically preexpanded. The tension dependencies of rate constants for channel closure (k(off)) predict different positions of rate-limiting barriers on the energy-area profiles for WT and GOF MscL. The data support the two-gate mechanism in which the first subtransition (C-->S) can be viewed as opening of the central (M1) gate, resulting in an expanded water-filled "leaky" conformation. Strong facilitation of this step by polar GOF substitutions suggests that separation of M1 helices associated with hydration of the pore in WT MscL is the major energetic barrier for opening. Mutants with a stabilized S1 gate demonstrate impeded transitions from low-conducting substates to the fully open state, whereas extensions of S1-M1 linkers result in a much higher probability of reverse O-->S transitions. These data strongly suggest that the bulk of conductance gain in the second subtransition (S-->O) occurs through the opening of the NH2-terminal (S1) gate and the linkers are coupling elements between the M1 and S1 gates.     

Defining the physical gate of a mechanosensitive channel, MscL, by engineering metal-binding sites

The mechanosensitive channel of large conductance, MscL, of Escherichia coli is one of the best-studied mechanosensitive proteins. Although the structure of the closed or "nearly-closed" state of the Mycobacterium tuberculosis ortholog has been solved and mechanisms of gating have been proposed, the transition from the closed to the open states remains controversial. Here, we probe the relative position of specific residues predicted to line the pore of MscL in either the closed state or during the closed-to-open transition by engineering single-site histidine substitutions and assessing the ability of Ni2+, Cd2+ or Zn2+ ions to affect channel activity. All residues predicted to be within the pore led to a change in channel threshold pressure, although the direction and extent of this change were dependent upon the mutation and metal used. One of the MscL mutants, L19H, exhibited gating that was inhibited by Cd2+ but stimulated by Ni2+, suggesting that these metals bind to and influence different states of the channel. Together, the results derived from this study support the hypotheses that the crystal structure depicts a "nearly closed" rather than a "fully closed" state of MscL, and that a clockwise rotation of transmembrane domain 1 occurs early in the gating process.     

CFTR displays voltage dependence and two gating modes during stimulation

The patch-clamp technique in conjunction with current noise analysis was employed to clarify the events underlying the regulation of the CFTR (cystic fibrosis transmembrane conductance regulator) during cAMP- dependent stimulation. 3T3 fibroblast cells expressing the CFTR were stimulated in cell-attached mode with forskolin. The number (N) of activated channels per patch ranged from 1 to approximately 100. In true single-channel recordings, CFTR's gating was best described by two open states (approximately 5 and approximately 100 ms) and three closed states (< or = 5, approximately 100, and approximately 1,000 ms). Current noise analysis resulted in spectra containing two distinct Lorentzian noise components with corner frequencies of 1.3 Hz and approximately 50 Hz, respectively. Single-channel time constants were dependent on voltage. The fastest closed state increased its contribution from 48% at +100 mV to 87% at -100 mV, and the medium open state reduced its length to one half, resulting in gating dominated by fast events. Similarly, the fast Lorentzian increased its amplitude, and its corner frequency increased from 44 Hz at +100 mV to 91 Hz at - 100 mV, while the slow Lorentzian was voltage independent. In multi- channel recordings N.Po (i.e., N times open probability) increased significantly, on average by 52% between -90 and +90 mV. Stimulation with forskolin increased Po of CFTR to approximately 0.5, which resulted from a decrease of the longest closed state while the faster open and closed states were unaffected. Neither corner frequency was affected during stimulation. Recordings from multichannel patches revealed in addition, unique, very long channel openings (high Po mode, average 13 s). Channels exhibiting high Po (i.e., Po approximately 1.0) or low Po (i.e., Po approximately 0.5) gating modes were both present in multichannel recordings, and CFTRs switched modes during stimulation. In addition, the switch to the high Po mode appeared to be a cooperative event for channel pairs. High forskolin concentration (i.e., 10 microM) favored transition into the high Po mode, suggesting a cellularly mediated regulation of model switching due to a fundamental change in configuration of the CFTR. Thus, during stimulation the CFTR increased its activity through two distinct effects: the reduction of the long closed state and modal switching to the high Po mode.     

Two distinct gating mechanisms in gap junction channels: CO2-sensitive and voltage-sensitive.

The chemical gating of single-gap junction channels was studied by the dual whole-cell voltage-clamp method in HeLa cells transfected with connexin43 (HeLa43) and in fibroblasts from sciatic nerves. Junctional current (Ij), single-channel conductance, and Ij kinetics were studied in cell pairs during CO2 uncoupling and recoupling at small transjunctional voltages (Vj < 35 mV: Vj gating absent) and at high Vj (Vj > 40 mV: Vj gating strongly activated). In the absence of Vj gating, CO2 exclusively caused Ij slow transitions from open to closed channel states (mean transition time: approximately 10 ms), corresponding to a single-channel conductance of approximately 120 pS. At Vj > 40 mV, Vj gating induced fast Ij flickering between open, gamma j(main state), and residual, gamma j(residual), states (transition time: approximately 2 ms). The ratio gamma j(main state)/gamma j(residual) was approximately 4-5. No obvious correlation between Ij fast flickering and CO2 treatment was noticed. At high Vj, in addition to slow Ij transitions between open and closed states, CO2 induced slow transitions between residual and closed states. During recoupling, each channel reopened by a slow transition (mean transition time: approximately 10 ms) from closed to open state (rarely from closed to residual state). Fast Ij flickering between open and residual states followed. The data are in agreement with the hypothesis that gap junction channels possess two gating mechanisms, and indicate that CO2 induces channel gating exclusively by the slow gating mechanism.     

Gating of MscL studied by steered molecular dynamics

Steered molecular dynamics simulations of the mechanosensitive channel of large conductance, MscL, were used to investigate how forces arising from membrane tension induce gating of the channel. A homology model of the closed form of MscL from Escherichia coli was subjected to external forces of 35-70 pN applied to residues near the membrane-water interface. The magnitude and location of these forces corresponded to those determined from the lateral pressure profile computed from a lipid bilayer simulation. A fully expanded state was obtained on the 10-ns timescale that revealed the mechanism for transducing membrane forces into channel opening. The expanded state agrees well with proposed models of MscL gating, in that it entails an irislike expansion of the pore accompanied by tilting of the transmembrane helices. The channel was most easily opened when force was applied predominantly on the cytoplasmic side of MscL. Comparison of simulations in which gating progressed to varying degrees identified residues that pose steric hindrance to channel opening.     

Structural determinants of MscL gating studied by molecular dynamics simulations

The mechanosensitive channel of large conductance (MscL) in prokaryotes plays a crucial role in exocytosis as well as in the response to osmotic downshock. The channel can be gated by tension in the membrane bilayer. The determination of functionally important residues in MscL, patch-clamp studies of pressure-conductance relationships, and the recently elucidated crystal structure of MscL from Mycobacterium tuberculosis have guided the search for the mechanism of MscL gating. Here, we present a molecular dynamics study of the MscL protein embedded in a fully hydrated POPC bilayer. Simulations totaling 3 ns in length were carried out under conditions of constant temperature and pressure using periodic boundary conditions and full electrostatics. The protein remained in the closed state corresponding to the crystal structure, as evidenced by its impermeability to water. Analysis of equilibrium fluctuations showed that the protein was least mobile in the narrowest part of the channel. The gating process was investigated through simulations of the bare protein under conditions of constant surface tension. Under a range of conditions, the transmembrane helices flattened as the pore widened. Implications for the gating mechanism in light of these and experimental results are discussed.     

Hydrophilicity of a single residue within MscL correlates with increased channel mechanosensitivity.

Mechanosensitive channel large (MscL) encodes the large conductance mechanosensitive channel of the Escherichia coli inner membrane that protects bacteria from lysis upon osmotic shock. To elucidate the molecular mechanism of MscL gating, we have comprehensively substituted Gly(22) with all other common amino acids. Gly(22) was highlighted in random mutagenesis screens of E. coli MscL (, Proc. Nat. Acad. Sci. USA. 95:11471-11475). By analogy to the recently published MscL structure from Mycobacterium tuberculosis (, Science. 282:2220-2226), Gly(22) is buried within the constriction that closes the pore. Substituting Gly(22) with hydrophilic residues decreased the threshold pressure at which channels opened and uncovered an intermediate subconducting state. In contrast, hydrophobic substitutions increased the threshold pressure. Although hydrophobic substitutions had no effect on growth, similar to the effect of an MscL deletion, channel hyperactivity caused by hydrophilic substitutions correlated with decreased proliferation. These results suggest a model for gating in which Gly(22) moves from a hydrophobic, and through a hydrophilic, environment upon transition from the closed to open conformation.     

Gating of the mechanosensitive channel protein MscL: the interplay of membrane and protein

The mechanosensitive channel of large conductance (MscL) belongs to a family of transmembrane channel proteins in bacteria and functions as a safety valve that relieves the turgor pressure produced by osmotic downshock. MscL gating can be triggered solely by stretching of the membrane. This work reports an effort to understand this mechanotransduction by means of molecular dynamics (MD) simulation on the MscL of mycobacterium tuberculosis embedded in a palmitoyloleoylphosphatidylethanolamine membrane. Equilibrium MD under zero membrane tension produced a more compact protein structure, as measured by its radii of gyration, compared to the crystal structure, in agreement with previous experimental findings. Even under a large applied tension up to 1000 dyn/cm, the MscL lateral dimension largely remained unchanged after up to 20 ns of simulation. A nonequilibrium MD simulation of 3% membrane expansion showed a significant increase in membrane rigidity upon MscL inclusion, which can contribute to efficient mechanotransduction. Direct observation of channel opening was possible only when an explicit lateral bias force was applied to each of the five subunits of MscL in the radially outward direction. Using this force, open structures with a large pore of radius 10 Å could be obtained. The channel opening takes place in a stepwise manner and concurrently with the water chain formation across the channel, which occurs without direct involvement of protein hydrophilic residues. The N-terminal S1 helices stabilize the open structure, and the membrane asymmetry (different lipid density on the two leaflets of membrane) promotes channel opening.     

Effect of ATP concentration on CFTR Cl- channels: a kinetic analysis of channel regulation.

Phosphorylated cystic fibrosis transmembrane conductance regulator (CFTR) Cl- channels require nucleoside triphosphates, such as ATP, to open. As the concentration of intracellular ATP increases, the probability of the channel being open (Po) increases. To better understand how ATP regulates the channel, we studied excised inside-out membrane patches that contained single, phosphorylated CFTR Cl- channels and examined the kinetics of gating at different concentrations of ATP. As the ATP concentration increased from 0.1 to 3 mM the mean closed time decreased, but mean open time did not change. Analysis of the data using histograms of open- and closed-state durations, the maximum likelihood method, and the log-likelihood ratio test suggested that channel behavior could be described by a model containing one open and two closed states (C1<==>C2<==>O). ATP regulated phosphorylated channels at the transition between the closed states C1 and C2: as the concentration of ATP increased, the rate of transition from C1 to C2 (C1-->C2) increased. In contrast, transitions from C2 to C1 and between C2 and the open state (O) were not significantly altered by ATP. Addition of ADP in the presence of ATP decreased the transition rate from C1 to C2 without affecting other transition rates. These data suggest that ATP regulates CFTR Cl- channels through an interaction that increases the rate of transition from the closed state to a bursting state in which the channel flickers back and forth between an open and a closed state (C2). This transition may reflect ATP binding or perhaps a step subsequent to binding.     

Mechanical Activation of MscL Revealed by a Locally Distributed Tension Molecular Dynamics Approach

Membrane tension perceived by mechanosensitive (MS) proteins mediates cellular responses to mechanical stimuli and osmotic stresses, and it also guides multiple biological functions including cardiovascular control and development. In bacteria, MS channels function as tension-activated pores limiting excessive turgor pressure, with MS channel of large conductance (MscL) acting as an emergency release valve preventing cell lysis. Previous attempts to simulate gating transitions in MscL by either directly applying steering forces to the protein or by increasing the whole-system tension were not fully successful and often disrupted the integrity of the system. We present a novel, to our knowledge, locally distributed tension molecular dynamics (LDT-MD) simulation method that allows application of forces continuously distributed among lipids surrounding the channel using a specially constructed collective variable. We report reproducible and reversible transitions of MscL to the open state with measured parameters of lateral expansion and conductivity that exactly satisfy experimental values. The LDT-MD method enables exploration of the MscL-gating process with different pulling velocities and variable tension asymmetry between the inner and outer membrane leaflets. We use LDT-MD in combination with well-tempered metadynamics to reconstruct the tension-dependent free-energy landscape for the opening transition in MscL. The flexible definition of the LDT collective variable allows general application of our method to study mechanical activation of any membrane-embedded protein.     

A large iris-like expansion of a mechanosensitive channel protein induced by membrane tension

MscL, a bacterial mechanosensitive channel of large conductance, is the first structurally characterized mechanosensor protein. Molecular models of its gating mechanisms are tested here. Disulfide crosslinking shows that M1 transmembrane alpha-helices in MscL of resting Escherichia coli are arranged similarly to those in the crystal structure of MscL from Mycobacterium tuberculosis. An expanded conformation was trapped in osmotically shocked cells by the specific bridging between Cys 20 and Cys 36 of adjacent M1 helices. These bridges stabilized the open channel. Disulfide bonds engineered between the M1 and M2 helices of adjacent subunits (Cys 32-Cys 81) do not prevent channel gating. These findings support gating models in which interactions between M1 and M2 of adjacent subunits remain unaltered while their tilts simultaneously increase. The MscL barrel, therefore, undergoes a large concerted iris-like expansion and flattening when perturbed by membrane tension.     

Structural Investigation of MscL Gating Using Experimental Data and Coarse Grained MD Simulations

The mechanosensitive channel of large conductance (MscL) has become a model system in which to understand mechanosensation, a process involved in osmoregulation and many other physiological functions. While a high resolution closed state structure is available, details of the open structure and the gating mechanism remain unknown. In this study we combine coarse grained simulations with restraints from EPR and FRET experiments to study the structural changes involved in gating with much greater level of conformational sampling than has previously been possible. We generated a set of plausible open pore structures that agree well with existing open pore structures and gating models. Most interestingly, we found that membrane thinning induces a kink in the upper part of TM1 that causes an outward motion of the periplasmic loop away from the pore centre. This previously unobserved structural change might present a new mechanism of tension sensing and might be related to a functional role in osmoregulation.     

The gating mechanism of the bacterial mechanosensitive channel MscL revealed by molecular dynamics simulations

One of the ultimate goals of the study on mechanosensitive (MS) channels is to understand the biophysical mechanisms of how the MS channel protein senses forces and how the sensed force induces channel gating. The bacterial MS channel MscL is an ideal subject to reach this goal owing to its resolved 3D protein structure in the closed state on the atomic scale and large amounts of electrophysiological data on its gating kinetics. However, the structural basis of the dynamic process from the closed to open states in MscL is not fully understood. In this study, we performed molecular dynamics (MD) simulations on the initial process of MscL opening in response to a tension increase in the lipid bilayer. To identify the tension-sensing site(s) in the channel protein, we calculated interaction energy between membrane lipids and candidate amino acids (AAs) facing the lipids. We found that Phe78 has a conspicuous interaction with the lipids, suggesting that Phe78 is the primary tension sensor of MscL. Increased membrane tension by membrane stretch dragged radially the inner (TM1) and outer (TM2) helices of MscL at Phe78, and the force was transmitted to the pentagon-shaped gate that is formed by the crossing of the neighboring TM1 helices in the inner leaflet of the bilayer. The radial dragging force induced radial sliding of the crossing portions, leading to a gate expansion. Calculated energy for this expansion is comparable to an experimentally estimated energy difference between the closed and the first subconductance state, suggesting that our model simulates the initial step toward the full opening of MscL. The model also successfully mimicked the behaviors of a gain of function mutant (G22N) and a loss of function mutant (F78N), strongly supporting that our MD model did simulate some essential biophysical aspects of the mechano-gating in MscL.     

The gating mechanism of the bacterial mechanosensitive channel MscL revealed by molecular dynamics simulations: From tension sensing to channel opening

One of the ultimate goals of the study on mechanosensitive (MS) channels is to understand the biophysical mechanisms of how the MS channel protein senses forces and how the sensed force induces channel gating. The bacterial MS channel MscL is an ideal subject to reach this goal owing to its resolved 3D protein structure in the closed state on the atomic scale and large amounts of electrophysiological data on its gating kinetics. However, the structural basis of the dynamic process from the closed to open states in MscL is not fully understood. In this study, we performed molecular dynamics (MD) simulations on the initial process of MscL opening in response to a tension increase in the lipid bilayer. To identify the tension-sensing site(s) in the channel protein, we calculated interaction energy between membrane lipids and candidate amino acids (AAs) facing the lipids. We found that Phe78 has a conspicuous interaction with the lipids, suggesting that Phe78 is the primary tension sensor of MscL. Increased membrane tension by membrane stretch dragged radially the inner (TM1) and outer (TM2) helices of MscL at Phe78, and the force was transmitted to the pentagon-shaped gate that is formed by the crossing of the neighboring TM1 helices in the inner leaflet of the bilayer. The radial dragging force induced radial sliding of the crossing portions, leading to a gate expansion. Calculated energy for this expansion is comparable to an experimentally estimated energy difference between the closed and the first subconductance state, suggesting that our model simulates the initial step toward the full opening of MscL. The model also successfully mimicked the behaviors of a gain of function mutant (G22N) and a loss of function mutant (F78N), strongly supporting that our MD model did simulate some essential biophysical aspects of the mechano-gating in MscL.     

New insights of membrane environment effects on MscL channel mechanics from theoretical approaches

The prokaryotic mechanosensitive channel of large conductance (MscL) is a remarkable integral membrane protein. During hypo-osmotic shock, it responses to membrane tension through large conformational changes, that lead to an open state of the pore. The structure of the channel from Mycobacterium tuberculosis has been resolved in the closed state. Numerous experiments have attempted to trap the channel in its open state but they did not succeed in obtaining a structure. A gating mechanism has been proposed based on different experimental data but there is no experimental technique available to follow this process in atomic details. In addition, it has been shown that a decrease of the lipid bilayer thickness lowered MscL activation energy and stabilized a structurally distinct closed channel intermediate. Here, we use atomistic molecular dynamics simulations to investigate the effect of the lipid bilayer thinning on our model of the structure of the Escherichia coli. We thoroughly analyze simulations of the channel embedded in two pre-equilibrated membranes differing by their hydrophobic tail length (DMPE and POPE). The MscL structure remains stable in POPE, whereas a distinct structural state is obtained in DMPE in response to hydrophobic mismatch. This latter is obtained by tilts and kinks of the transmembrane helices, leading to a widening and a diminution of the channel height. Part of these motions is guided by a competition between solvent and lipids for the interaction with the periplasmic loops. We finally conduct a principal component analysis of the simulation and compare anharmonic motions with harmonic ones, previously obtained from a coarse-grained normal mode analysis performed on the same structural model. Significant similarities exist between low-frequency harmonic motions and those observed with essential dynamics in DMPE. In summary, change in membrane thickness permits to accelerate the conformational changes involved in the mechanics of the E. coli channel, providing a closed structural intermediate en route to the open state. These results give clues for better understanding why the channel activation energy is lowered in a thinner membrane.     

On the structure of the N-terminal domain of the MscL channel: helical bundle or membrane interface

The mechanosensitive channel of large conductance, MscL, serves as a biological emergency release valve protecting bacteria from acute osmotic downshock and is to date the best characterized mechanosensitive channel. A well-recognized and supported model for Escherichia coli MscL gating proposes that the N-terminal 11 amino acids of this protein form a bundle of amphipathic helices in the closed state that functionally serves as a cytoplasmic second gate. However, a recently reexamined crystal structure of a closed state of the Mycobacterium tuberculosis MscL shows these helices running along the cytoplasmic surface of the membrane. Thus, it is unclear if one structural model is correct or if they both reflect valid closed states. Here, we have systematically reevaluated this region utilizing cysteine-scanning, in vivo functional characterization, in vivo SCAM, electrophysiological studies, and disulfide-trapping experiments. The disulfide-trapping pattern and functional studies do not support the helical bundle and second-gate hypothesis but correlate well with the proposed structure for M. tuberculosis MscL. We propose a functional model that is consistent with the collective data.     

Gating kinetics of batrachotoxin-modified Na+ channels in the squid giant axon. Voltage and temperature effects.

The gating kinetics of batrachotoxin-modified Na+ channels were studied in outside-out patches of axolemma from the squid giant axon by means of the cut-open axon technique. Single channel kinetics were characterized at different membrane voltages and temperatures. The probability of channel opening (Po) as a function of voltage was well described by a Boltzmann distribution with an equivalent number of gating particles of 3.58. The voltage at which the channel was open 50% of the time was a function of [Na+] and temperature. A decrease in the internal [Na+] induced a shift to the right of the Po vs. V curve, suggesting the presence of an integral negative fixed charge near the activation gate. An increase in temperature decreased Po, indicating a stabilization of the closed configuration of the channel and also a decrease in entropy upon channel opening. Probability density analysis of dwell times in the closed and open states of the channel at 0 degrees C revealed the presence of three closed and three open states. The slowest open kinetic component constituted only a small fraction of the total number of transitions and became negligible at voltages greater than -65 mV. Adjacent interval analysis showed that there is no correlation in the duration of successive open and closed events. Consistent with this analysis, maximum likelihood estimation of the rate constants for nine different single-channel models produced a preferred model (model 1) having a linear sequence of closed states and two open states emerging from the last closed state. The effect of temperature on the rate constants of model 1 was studied. An increase in temperature increased all rate constants; the shift in Po would be the result of an increase in the closing rates predominant over the change in the opening rates. The temperature study also provided the basis for building an energy diagram for the transitions between channel states.     

Kinetic analysis of channel gating. Application to the cholinergic receptor channel and the chloride channel from Torpedo californica.

Identification of the minimum number of ways in which open and closed states communicate is a crucial step in defining the gating kinetics of multistate channels. We used certain correlation functions to extract information about the pathways connecting the open and closed states of the cation channel of the purified nicotinic acetylcholine receptor and of the chloride channel of Torpedo californica electroplax membranes. Single channel currents were recorded from planar lipid bilayers containing the membrane channel proteins under investigation. The correlation functions are conveniently computed from single channel current records and yield information on E, the minimum number of entry/exit states into the open or closed aggregates. E gives a lower limit on the numbers of transition pathways between open and closed states. For the acetylcholine receptor, the autocorrelation analysis shows that there are at least two entry/exit states through which the open and closed aggregates communicate. The chloride channel fluctuates between three conductance substates, here indentified as C, M, and H for closed, intermediate, and high conductance, respectively. Correlation analysis shows that E is greater than or equal to 2 for the M aggregate, indicating that there are at least two distinct entry/exit states in the M aggregate. In contrast, there is no evidence for the existence of more than one entry/exit state in the C or H aggregates. Thus, these correlation functions provide a simple and general strategy to extract information on channel gating kinetics.     

Dynamical properties of the MscL of Escherichia coli: a normal mode analysis

The mechanosensitive channel (MscL) is an integral membrane protein which gates in response to membrane tension. Physiological data have shown that the gating transition involves a very large change in the conformation, and that the open state of the channel forms a large non-specific pore with a high conductance. The Escherichia coli channel structure was first modeled by homology modeling, starting with the X-ray structure of the homologous from Mycobacterium tuberculosis. Then, the dynamical and conformational properties of the channel were explored, using normal mode analysis. Such an analysis was also performed with the different structures proposed recently by Sukharev and co-workers. Similar dynamical behaviors are observed, which are characteristic of the channel architecture, subtle differences being due to the different relative positioning of the structural elements. The ability of particular regions of the channel to deform is discussed with respect to the functional and structural properties, implied in the gating process. Our results show that the first step of the gating mechanism can be described with three low-frequency modes only. The movement associated to these modes is clearly an iris-like movement involving both tilt and twist rotation.