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.     

Purification and functional reconstitution of N- and C-halves of the MscL channel

MscL is a mechanosensitive channel gated by membrane tension in the lipid bilayer alone. Its structure, known from x-ray crystallography, indicates that it is a homopentamer. Each subunit comprises two transmembrane segments TM1 and TM2 connected by a periplasmic loop. The closed pore is lined by five TM1 helices. We expressed in Escherichia coli and purified two halves of the protein, each containing one of the transmembrane segments. Their electrophysiological activity was studied by the patch-clamp recording upon reconstitution in artificial liposomes. The TM2 moiety had no electrophysiological activity, whereas the TM1 half formed channels, which were not affected by membrane tension and varied in conductance between 50 and 350 pS in 100 mM KCl. Coreconstitution of the two halves of MscL however, yielded mechanosensitive channels having the same conductance as the native MscL (1500 pS), but exhibiting increased sensitivity to pressure. Our results confirm the current view on the functional role of TM1 and TM2 helices in the MscL gating and emphasize the importance of helix-helix interactions for the assembly and functional properties of the channel protein. In addition, the results indicate a crucial role of the periplasmic loop for the channel mechanosensitivity.     

The dynamics of protein-protein interactions between domains of MscL at the cytoplasmic-lipid interface

The bacterial mechanosensitive channel of large conductance, MscL, is one of the best characterized mechanosensitive channels serving as a paradigm for how proteins can sense and transduce mechanical forces. The physiological role of MscL is that of an emergency release valve that opens a large pore upon a sudden drop in the osmolarity of the environment. A crystal structure of a closed state of MscL shows it as a homopentamer, with each subunit consisting of two transmembrane domains (TM). There is consensus that the TM helices move in an iris like manner tilting in the plane of the membrane while gating. An N-terminal amphipathic helix that lies along the cytoplasmic membrane (S1), and the portion of TM2 near the cytoplasmic interface (TM2_ci), are relatively close in the crystal structure, yet predicted to be dynamic upon gating. Here we determine how these two regions interact in the channel complex, and study how these interactions change as the channel opens. We have screened 143 double-cysteine mutants of E. coli MscL for their efficiency in disulfide bridging and generated a map of protein-protein interactions between these two regions. Interesting candidates have been further studied by patch clamp and show differences in channel activity under different redox potentials; the results suggest a model for the dynamics of these two domains during MscL gating.     

The dynamics of protein-protein interactions between domains of MscL at the cytoplasmic-lipid interface

The bacterial mechanosensitive channel of large conductance, MscL, is one of the best characterized mechanosensitive channels serving as a paradigm for how proteins can sense and transduce mechanical forces. The physiological role of MscL is that of an emergency release valve that opens a large pore upon a sudden drop in the osmolarity of the environment. A crystal structure of a closed state of MscL shows it as a homopentamer, with each subunit consisting of two transmembrane domains (TM). There is consensus that the TM helices move in an iris like manner tilting in the plane of the membrane while gating. An N-terminal amphipathic helix that lies along the cytoplasmic membrane (S1), and the portion of TM2 near the cytoplasmic interface (TM2_ci), are relatively close in the crystal structure, yet predicted to be dynamic upon gating. Here we determine how these two regions interact in the channel complex, and study how these interactions change as the channel opens. We have screened 143 double-cysteine mutants of E. coli MscL for their efficiency in disulfide bridging and generated a map of protein-protein interactions between these two regions. Interesting candidates have been further studied by patch clamp and show differences in channel activity under different redox potentials; the results suggest a model for the dynamics of these two domains during MscL gating.     

The Combined Effect of Hydrophobic Mismatch and Bilayer Local Bending on the Regulation of Mechanosensitive Ion Channels

The hydrophobic mismatch between the lipid bilayer and integral membrane proteins has well-defined effect on mechanosensitive (MS) ion channels. Also, membrane local bending is suggested to modulate MS channel activity. Although a number of studies have already shown the significance of each individual factor, the combined effect of these physical factors on MS channel activity have not been investigated. Here using finite element simulation, we study the combined effect of hydrophobic mismatch and local bending on the archetypal mechanosensitive channel MscL. First we show how the local curvature direction impacts on MS channel modulation. In the case of MscL, we show inward (cytoplasmic) bending can more effectively gate the channel compared to outward bending. Then we indicate that in response to a specific local curvature, MscL inserted in a bilayer with the same hydrophobic length is more expanded in the constriction pore region compared to when there is a protein-lipid hydrophobic mismatch. Interestingly in the presence of a negative mismatch (thicker lipids), MscL constriction pore is more expanded than in the presence of positive mismatch (thinner lipids) in response to an identical membrane curvature. These results were confirmed by a parametric energetic calculation provided for MscL gating. These findings have several biophysical consequences for understanding the function of MS channels in response to two major physical stimuli in mechanobiology, namely hydrophobic mismatch and local membrane curvature.     

Physical principles underlying the transduction of bilayer deformation forces during mechanosensitive channel gating

In mechanosensitive (MS) channels, gating is initiated by changes in intra-bilayer pressure profiles originating from bilayer deformation. Here we evaluated two physical mechanisms as triggers of MS channel gating: the energetic cost of protein-bilayer hydrophobic mismatches and the geometric consequences of bilayer intrinsic curvature. Structural changes in the Escherichia coli large MS channel (MscL) were studied under nominally zero transbilayer pressures using both patch clamp and EPR spectroscopic approaches. Changes in membrane intrinsic curvature induced by the external addition of lysophosphatidylcholine (LPC) generated massive spectroscopic changes in the narrow constriction that forms the channel 'gate', trapping the channel in the fully open state. Hydrophobic mismatch alone was unable to open the channel, but decreasing bilayer thickness lowered MscL activation energy, stabilizing a structurally distinct closed channel intermediate. We propose that the mechanism of mechanotransduction in MS channels is defined by both local and global asymmetries in the transbilayer pressure profile at the lipid-protein interface.     

The (beta)gamma subunits of G proteins gate a K(+) channel by pivoted bending of a transmembrane segment

The molecular mechanism of ion channel gating remains unclear. Using approaches such as proline scanning mutagenesis and homology modeling, we localize the gate of the K(+) channels controlled by the (beta)gamma subunits of G proteins at the pore-lining bundle crossing of the second transmembrane (TM2) helices. We show that the flexibility afforded by a highly conserved glycine residue in the middle of TM2 is crucial for channel gating. In contrast, flexibility introduced immediately below the gate disrupts gating. We propose that the force produced by channel-G(beta)gamma interactions is transduced through the rigid region below the helix bundle crossing to bend TM2 at the glycine that serves as a hinge and open the gate.     

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.     

Membrane tension, lipid adaptation, conformational changes, and energetics in MscL gating

This study aims to explore gating mechanisms of mechanosensitive channels in terms of membrane tension, membrane adaptation, protein conformation, and energetics. The large conductance mechanosensitive channel from Mycobacterium tuberculosis (Tb-MscL) is used as a model system; Tb-MscL acts as a safety valve by releasing small osmolytes through the channel opening under extreme hypoosmotic conditions. Based on the assumption that the channel gating involves tilting of the transmembrane (TM) helices, we have performed free energy simulations of Tb-MscL as a function of TM helix tilt angle in a dimyristoylphosphatidylcholine bilayer. Based on the change in system dimensions, TM helix tilting is shown to be essentially equivalent to applying an excess surface tension to the membrane, causing channel expansion, lipid adaptation, and membrane thinning. Such equivalence is further corroborated by the observation that the free energy cost of Tb-MscL channel expansion is comparable to the work done by the excess surface tension. Tb-MscL TM helix tilting results in an expanded water-conducting channel of an outer dimension similar to the proposed fully open MscL structure. The free energy decomposition indicates a possible expansion mechanism in which tilting and expanding of TM2 facilitates the iris-like motion of TM1, producing an expanded Tb-MscL.     

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.     

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.     

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.     

Gating of the ATP-sensitive K+ channel by a pore-lining phenylalanine residue

ATP-sensitive K(+) (K(ATP)) channels are gated by intracellular ATP, proton and phospholipids. The pore-forming Kir6.2 subunit has all essential machineries for channel gating by these ligands. It is known that channel gating involves the inner helix bundle of crossing in which a phenylalanine residue (Phe168) is found in the TM2 at the narrowest region of the ion-conduction pathway in the Kir6.2. Here we present evidence that Phe168-Kir6.2 functions as an ATP- and proton-activated gate via steric hindrance and hydrophobic interactions. Site-specific mutations of Phe168 to a small amino acid resulted in losses of the ATP- and proton-dependent gating, whereas the channel gating was well maintained after mutation to a bulky tryptophan, supporting the steric hindrance effect. The steric hindrance effect, though necessary, was insufficient for the gating, as mutating Phe168 to a bulky hydrophilic residue severely compromised the channel gating. Single-channel kinetics of the F168W mutant resembled the wild-type channel. Small residues increased P(open), and displayed long-lasting closures and long-lasting openings. Kinetic modeling showed that these resulted from stabilization of the channel to open and long-lived closed states, suggesting that a bulky and hydrophobic residue may lower the energy barrier for the switch between channel openings and closures. Thus, it is likely that the Phe168 acts as not only a steric hindrance gate but also potentially a facilitator of gating transitions in the Kir6.2 channel.     

Nanomechanical properties of MscL α helices: A steered molecular dynamics study

Gating of mechanosensitive (MS) channels is driven by a hierarchical cascade of movements and deformations of transmembrane helices in response to bilayer tension. Determining the intrinsic mechanical properties of the individual transmembrane helices is therefore central to understanding the intricacies of the gating mechanism of MS channels. We used a constant-force steered molecular dynamics (SMD) approach to perform unidirectional pulling tests on all the helices of MscL in M. tuberculosis and E. coli homologs. Using this method, we could overcome the issues encountered with the commonly used constant-velocity SMD simulations, such as low mechanical stability of the helix during stretching and high dependency of the elastic properties on the pulling rate. We estimated Young's moduli of the α-helices of MscL to vary between 0.2 and 12.5 GPa with TM2 helix being the stiffest. We also studied the effect of water on the properties of the pore-lining TM1 helix. In the absence of water, this helix exhibited a much stiffer response. By monitoring the number of hydrogen bonds, it appears that water acts like a ‘lubricant’ (softener) during TM1 helix elongation. These data shed light on another physical aspect underlying hydrophobic gating of MS channels, in particular MscL.     

Gating of the ATP-sensitive K channel by a pore-lining phenylalanine residue

ATP-sensitive K⁺ (K_ATP) channels are gated by intracellular ATP, proton and phospholipids. The pore-forming Kir6.2 subunit has all essential machineries for channel gating by these ligands. It is known that channel gating involves the inner helix bundle of crossing in which a phenylalanine residue (Phe168) is found in the TM2 at the narrowest region of the ion-conduction pathway in the Kir6.2. Here we present evidence that Phe168-Kir6.2 functions as an ATP- and proton-activated gate via steric hindrance and hydrophobic interactions. Site-specific mutations of Phe168 to a small amino acid resulted in losses of the ATP- and proton-dependent gating, whereas the channel gating was well maintained after mutation to a bulky tryptophan, supporting the steric hindrance effect. The steric hindrance effect, though necessary, was insufficient for the gating, as mutating Phe168 to a bulky hydrophilic residue severely compromised the channel gating. Single-channel kinetics of the F168W mutant resembled the wild-type channel. Small residues increased P_open, and displayed long-lasting closures and long-lasting openings. Kinetic modeling showed that these resulted from stabilization of the channel to open and long-lived closed states, suggesting that a bulky and hydrophobic residue may lower the energy barrier for the switch between channel openings and closures. Thus, it is likely that the Phe168 acts as not only a steric hindrance gate but also potentially a facilitator of gating transitions in the Kir6.2 channel.     

Energetic and spatial parameters for gating of the bacterial large conductance mechanosensitive channel, MscL

MscL is multimeric protein that forms a large conductance mechanosensitive channel in the inner membrane of Escherichia coli. Since MscL is gated by tension transmitted through the lipid bilayer, we have been able to measure its gating parameters as a function of absolute tension. Using purified MscL reconstituted in liposomes, we recorded single channel currents and varied the pressure gradient (P) to vary the tension (T). The tension was calculated from P and the radius of curvature was obtained using video microscopy of the patch. The probability of being open (Po) has a steep sigmoidal dependence on T, with a midpoint (T1/2) of 11.8 dyn/cm. The maximal slope sensitivity of Po/Pc was 0.63 dyn/cm per e-fold. Assuming a Boltzmann distribution, the energy difference between the closed and fully open states in the unstressed membrane was DeltaE = 18.6 kBT. If the mechanosensitivity arises from tension acting on a change of in-plane area (DeltaA), the free energy, TDeltaA, would correspond to DeltaA = 6.5 nm2. MscL is not a binary channel, but has four conducting states and a closed state. Most transition rates are independent of tension, but the rate-limiting step to opening is the transition between the closed state and the lowest conductance substate. This transition thus involves the greatest DeltaA. When summed over all transitions, the in-plane area change from closed to fully open was 6 nm2, agreeing with the value obtained in the two-state analysis. Assuming a cylindrical channel, the dimensions of the (fully open) pore were comparable to DeltaA. Thus, the tension dependence of channel gating is primarily one of increasing the external channel area to accommodate the pore of the smallest conducting state. The higher conducting states appear to involve conformational changes internal to the channel that don't involve changes in area.     

Voltage-induced gating of the mechanosensitive MscL ion channel reconstituted in a tethered lipid bilayer membrane

The mechanosensitive (MS) ion channel is gated by changes in bilayer deformation. It is functional without the presence of any other proteins and gating of the channel has been successfully achieved using conventional patch clamping techniques where a voltage has been applied together with a pressure over the membrane. Here, we have for the first time analyzed the large conducting (MscL) channel in a supported membrane using only an external electrical field. This was made possible using a newly developed technique utilizing a tethered lipid bilayer membrane (tBLM), which is part of an engineered microelectronic array chip. Single ion channel activity characteristic for MscL was obtained, albeit with lower conductivity. The ion channel was gated using solely a transmembrane potential of 300 mV. Computations demonstrate that this amount of membrane potential induces a membrane tension of 12 dyn/cm, equivalent to that calculated to gate the channel in patch clamp from pressure-induced stretching of the bilayer. These results strengthen the supposition that the MscL ion channel gates in response to stress in the lipid membrane rather than pressure across it. Furthermore, these findings illustrate the possibility of using the MscL as a release valve for engineered membrane devices; one step closer to mimicking the true function of the living cell.     

Modulation of channel activity and gadolinium block of MscL by static magnetic fields

The magnetic field of the Earth has for long been known to influence the behaviour and orientation of a variety of living organisms. Experimental studies of the magnetic sense have, however, been impaired by the lack of a plausible cellular and/or molecular mechanism providing meaningful explanation for detection of magnetic fields by these organisms. Recently, mechanosensitive (MS) ion channels have been implied to play a role in magnetoreception. In this study we have investigated the effect of static magnetic fields (SMFs) of moderate intensity on the activity and gadolinium block of MscL, the bacterial MS channel of large conductance, which has served as a model channel to study the basic physical principles of mechanosensory transduction in living cells. In addition to showing that direct application of the magnetic field decreased the activity of the MscL channel, our study demonstrates for the first time that SMFs can reverse the effect of gadolinium, a well-known blocker of MS channels. The results of our study are consistent with a notion that (1) the effects of SMFs on the MscL channels may result from changes in physical properties of the lipid bilayer due to diamagnetic anisotropy of phospholipid molecules and consequently (2) cooperative superdiamagnetism of phospholipid molecules under influence of SMFs could cause displacement of Gd³⁺ ions from the membrane bilayer and thus remove the MscL channel block.     

Generation and evaluation of a large mutational library from the Escherichia coli mechanosensitive channel of large conductance,MscL: implications for channel gating and evolutionary design

Random mutagenesis of the mechanosensitive channel of large conductance (MscL) from Escherichia coli coupled with a high-throughput functional screen has provided new insights into channel structure and function. Complementary interactions of conserved residues proposed in a computational model for gating have been evaluated, and important functional regions of the channel have been identified. Mutational analysis shows that the proposed S1 helix, despite having several highly conserved residues, can be heavily mutated without significantly altering channel function. The pattern of mutations that make MscL more difficult to gate suggests that MscL senses tension with residues located near the lipid headgroups of the bilayer. The range of phenotypical changes seen has implications for a proposed model for the evolutionary origin of mechanosensitive channels.     

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.