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

Sensing bilayer tension: bacterial mechanosensitive channels and their gating mechanisms

Mechanosensitive channels sense and respond to changes in bilayer tension. In many respects, this is a unique property: the changes in membrane tension gate the channel, leading to the transient formation of open non-selective pores. Pore diameter is also high for the bacterial channels studied, MscS and MscL. Consequently, in cells, gating has severe consequences for energetics and homoeostasis, since membrane depolarization and modification of cytoplasmic ionic composition is an immediate consequence. Protection against disruption of cellular integrity, which is the function of the major channels, provides a strong evolutionary rationale for possession of such disruptive channels. The elegant crystal structures for these channels has opened the way to detailed investigations that combine molecular genetics with electrophysiology and studies of cellular behaviour. In the present article, the focus is primarily on the structure of MscS, the small mechanosensitive channel. The description of the structure is accompanied by discussion of the major sites of channel–lipid interaction and reasoned, but limited, speculation on the potential mechanisms of tension sensing leading to gating.

Adaptive MscS gating in the osmotic permeability response in E. coli: the question of time

Microorganisms adapt to osmotic downshifts by releasing small osmolytes through mechanosensitive (MS) channels. We want to understand how the small mechanosensitive channel's (MscS) activation and inactivation, both driven by membrane tension, optimize survival in varying hypoosmotic shock situations. By measuring light scattering with a stopped-flow device, we estimate bacterial swelling time as 30-50 ms. A partial solute equilibration follows within 150-200 ms, during which optical responses from cells with WT MscS deviate from those lacking MS channels. MscS opening rates estimated in patch clamp show the channels readily respond to tensions below the lytic limit with a time course faster than 20 ms and close promptly upon tension release. To address the role of the tension-insensitive inactivated state in vivo, we applied short, long, and two-step osmotic shock protocols to WT, noninactivating G113A, and fast-inactivating D62N mutants. WT and G113A showed a comparable survival in short 1 min 800 mOsm downshock experiments, but G113A was at a disadvantage under a long 60 min shock. Preshocking cells carrying WT MscS for 15 s to 15 min with a 200 mOsm downshift did not sensitize them to the final 500 mOsm drop in osmolarity of the second step. However, these two-step shocks induced death in D62N more than just a one-step 700 mOsm downshift. We conclude MscS is able to activate and exude osmolytes faster than lytic pressure builds inside the cell under abrupt shock. During prolonged shocks, gradual inactivation prevents continuous channel activity and assists recovery. Slow kinetics of inactivation in WT MscS ensures that mild shocks do not inactivate the entire population, leaving some protection should conditions worsen.

Protein shape change has a major effect on the gating energy of a mechanosensitive channel

Increasing experimental evidence has shown that membrane protein functionality depends on molecular composition of cell membranes. However, the origin of this dependence is not fully understood. It is reasonable to assume that specific lipid-protein interactions are important, yet more generic effects due to mechanical properties of lipid bilayers likely play a significant role too. Previously it has been demonstrated using models for elastic properties of membranes and lateral pressure profiles of lipid bilayers that the mechanical properties of a lipid bilayer can contribute as much as ∼10 k(B)T to the free energy difference associated with a change in protein conformational state. Here, we extend those previous approaches to a more realistic model for a large mechanosensitive channel (MscL). We use molecular dynamics together with the MARTINI model to simulate the open and closed states of MscL embedded in a DOPC bilayer. We introduce a procedure to calculate the mechanical energy change in the channel gating using a three-dimensional pressure distribution inside a membrane, computed from the molecular dynamics simulations. We decompose the mechanical energy to terms associated with area dilation and shape contribution. Our results highlight that the lateral pressure profile of a lipid bilayer together with the shape change in gating can induce a contribution of ∼30 k(B)T on the gating energy of MscL. This contribution arises largely from the interfacial tension between hydrophobic and hydrophilic regions in a lipid bilayer.

Effects of GsMTx4 on bacterial mechanosensitive channels in inside-out patches from giant spheroplasts

GsMTx4 is a 34-residue peptide isolated from the tarantula Grammostola spatulata folded into an inhibitory cysteine knot and it selectively affects gating of some mechanosensitive channels. Here we report the effects of cytoplasmic GsMTx4 on the two bacterial channels, MscS and MscL, in giant Escherichia coli spheroplasts. In excised inside-out patches, GsMTx4 sensitized both channels to tension by increasing the opening rate and decreasing the closing rate. With ascending and descending pressure ramps, GsMTx4 increased the gating hysteresis for MscS, a consequence of slower gating kinetics. Quantitative kinetic analysis of the primary C↔O transition showed that the hysteresis is a result of the decreased closing rate. The gating barrier location relative to the open state energy well was unaffected by GsMTx4. A reconstructed energy profile suggests that the peptide prestresses the resting state of MscS, lowering the net barrier to opening and stabilizes the open conformation by ∼8 kT. In excised patches, both MscL and MscS exhibit reversible adaptation, a process separable from inactivation for MscS. GsMTx4 decreased the rate of reversible adaptation for both channels and the MscS recovery rate from the inactivation. These measurements support a mechanism where GsMTx4 binds to the lipid interface of the channel, increasing the local stress that is sensed by the channels and stabilizing the expanded conformations.

An improved open-channel structure of MscL determined from FRET confocal microscopy and simulation

Mechanosensitive channels act as molecular transducers of mechanical force exerted on the membrane of living cells by opening in response to membrane bilayer deformations occurring in physiological processes such as touch, hearing, blood pressure regulation, and osmoregulation. Here, we determine the likely structure of the open state of the mechanosensitive channel of large conductance using a combination of patch clamp, fluorescence resonance energy transfer (FRET) spectroscopy, data from previous electron paramagnetic resonance experiments, and molecular and Brownian dynamics simulations. We show that structural rearrangements of the protein can be measured in similar conditions as patch clamp recordings while controlling the state of the pore in its natural lipid environment by modifying the lateral pressure distribution via the lipid bilayer. Transition to the open state is less dramatic than previously proposed, while the N terminus remains anchored at the surface of the membrane where it can either guide the tilt of or directly translate membrane tension to the conformation of the pore-lining helix. Combining FRET data obtained in physiological conditions with simulations is likely to be of great value for studying conformational changes in a range of multimeric membrane proteins.

Hydration properties of mechanosensitive channel pores define the energetics of gating

Opening of ion channels directly by tension in the surrounding membrane appears to be the most ancient and simple mechanism of gating. Bacterial mechanosensitive channels MscL and MscS are the best-studied tension-gated nanopores, yet the key physical factors that define their gating are still hotly debated. Here we present estimations, simulations and experimental results showing that hydration of the pore might be one of the major parameters defining the thermodynamics and kinetics of mechanosensitive channel gating. We associate closing of channel pores with complete dehydration of the hydrophobic gate (occlusion by 'vapor lock') and formation of two water-vapor interfaces above and below the constriction. The opening path is the expansion of these interfaces, ultimately leading to wetting of the hydrophobic pore, which does not appear to be the exact reverse of the closing path, thus producing hysteresis. We discuss specifically the role of polar groups (glycines) buried in narrow closed conformations but exposed in the open states that change the wetting characteristics of the pore lining and stabilize conductive states of the channels.

Release of content through mechano-sensitive gates in pressurized liposomes

Mechano-sensitive channels are ubiquitous membrane proteins that activate in response to increasing tension in the lipid membrane. They facilitate a sudden, nonselective release of solutes and water that safeguards the integrity of the cell in hypo- or hyper-osmotic shock conditions. We have simulated the rapid release of content from a pressurized liposome through a particular mechano-sensitive protein channel, MscL, embedded in the liposomal membrane. We show that a single channel is able to relax the liposome, stressed to the point of bursting, in a matter of microseconds. We map the full activation-deactivation cycle of MscL in near-atomic detail and are able to quantify the rapid decrease in liposomal stress as a result of channel activation. This provides a computational tool that opens the way to contribute to the rational design of functional nano-containers.

Adaptive behavior of bacterial mechanosensitive channels is coupled to membrane mechanics

Mechanosensitive channel of small conductance (MscS), a tension-driven osmolyte release valve residing in the inner membrane of Escherichia coli, exhibits a complex adaptive behavior, whereas its functional counterpart, mechanosensitive channel of large conductance (MscL), was generally considered nonadaptive. In this study, we show that both channels exhibit similar adaptation in excised patches, a process that is completely separable from inactivation prominent only in MscS. When a membrane patch is held under constant pressure, adaptation of both channels is manifested as a reversible current decline. Their dose–response curves recorded with 1–10-s ramps of pressure are shifted toward higher tension relative to the curves measured with series of pulses, indicating decreased tension sensitivity. Prolonged exposure of excised patches to subthreshold tensions further shifts activation curves for both MscS and MscL toward higher tension with similar magnitude and time course. Whole spheroplast MscS recordings performed with simultaneous imaging reveal activation curves with a midpoint tension of 7.8 mN/m and the slope corresponding to ∼15-nm² in-plane expansion. Inactivation was retained in whole spheroplast mode, but no adaptation was observed. Similarly, whole spheroplast recordings of MscL (V23T mutant) indicated no adaptation, which was present in excised patches. MscS activities tried in spheroplast-attached mode showed no adaptation when the spheroplasts were intact, but permeabilized spheroplasts showed delayed adaptation, suggesting that the presence of membrane breaks or edges causes adaptation. We interpret this in the framework of the mechanics of the bilayer couple linking adaptation of channels in excised patches to the relaxation of the inner leaflet that is not in contact with the glass pipette. Relaxation of one leaflet results in asymmetric redistribution of tension in the bilayer that is less favorable for channel opening.

YbdG in Escherichia coli is a threshold-setting mechanosensitive channel with MscM activity

We describe a mechanosensitive (MS) channel that has mechanosensitive channel of miniconductance (MscM) activity, and displays unique properties with respect to gating. Mechanosensitive channels respond to membrane tension, are ubiquitous from bacteria to man, and exhibit a great diversity in structure and function. These channels protect Bacteria and Archaea against hypoosmotic shock and are critical determinants of shape in chloroplasts. Given the dominant roles played in bacteria by the mechanosensitive channel of small conductance (MscS) and the mechanosensitive channel of large conductance (MscL), the role of the multiple MS channel homologs observed in most organisms remains obscure. Here we demonstrate that a MscS homolog, YbdG, extends the range of hypoosmotic shock that Escherichia coli cells can survive, but its expression level is insufficient to protect against severe shocks. Overexpression of the YbdG protein provides complete protection. Transcription and translation of the ybdG gene are enhanced by osmotic stress consistent with a role for the protein in survival of hypoosmotic shock. Measurement of the conductance of the native channel by standard patch clamp methods was not possible. However, a fully functional YbdG mutant channel, V229A, exhibits a conductance in membrane patches consistent with MscM activity. We find that MscM activities arise from more than one gene product because ybdG deletion mutants still exhibit an occasional MscM-like conductance. We propose that ybdG encodes a low-abundance MscM-type MS channel, which in cells relieves low levels of membrane tension, obviating the need to activate the major MS channels, MscS and MscL.

Identification and experimental verification of a novel family of bacterial cyclic nucleotide-gated (bCNG) ion channels

Studies of bacterial ion channels have provided significant insights into the structure-function relationships of mechanosensitive and voltage-gated ion channels. However, to date, very few bacterial channels that respond to small molecules have been identified, cloned, and characterized. Here, we use bioinformatics to identify a novel family of bacterial cyclic nucleotide-gated (bCNG) ion channels containing a channel domain related by sequence homology to the mechanosensitive channel of small conductance (MscS). In this initial report, we clone selected members of this channel family, use electrophysiological measurements to verify their ability to directly gate in response to cyclic nucleotides, and use osmotic downshock to demonstrate their lack of mechanosensitivity. In addition to providing insight into bacterial physiology, these channels will provide researchers with a useful model system to investigate the role of ligand-gated ion channels (LGICs) in the signaling processes of higher organisms. The identification of these channels provides a foundation for structural and functional studies of LGICs that would be difficult to perform on mammalian channels. Moreover, the discovery of bCNG channels implies that bacteria have cyclic nucleotide-gated and cyclic nucleotide-modulated ion channels, which are analogous to the ion channels involved in eukaryotic secondary messenger signaling pathways.

Mechanosensitivity of ion channels based on protein-lipid interactions

Ion channels form a group of membrane proteins that pass ions through a pore beyond the energy barrier of the lipid bilayer. The structure of the transmembrane segment of membrane proteins is influenced by the charges and the hydrophobicity of the surrounding lipids and the pressure on its surface. A mechanosensitive channel is specifically designed to change its conformation in response to changes in the membrane pressure (tension). However, mechanosensitive channels are not the only group that is sensitive to the physical environment of the membrane: voltage-gated channels are also amenable to the lipid environment. In this article, we review the structure and gating mechanisms of the mechanosensitive channels and voltage-gated channels and discuss how their functions are affected by the physical properties of the lipid bilayer.

Mechanosensitivity of ion channels based on protein-lipid interactions

Ion channels form a group of membrane proteins that pass ions through a pore beyond the energy barrier of the lipid bilayer. The structure of the transmembrane segment of membrane proteins is influenced by the charges and the hydrophobicity of the surrounding lipids and the pressure on its surface. A mechanosensitive channel is specifically designed to change its conformation in response to changes in the membrane pressure (tension). However, mechanosensitive channels are not the only group that is sensitive to the physical environment of the membrane: voltage-gated channels are also amenable to the lipid environment. In this article, we review the structure and gating mechanisms of the mechanosensitive channels and voltage-gated channels and discuss how their functions are affected by the physical properties of the lipid bilayer.

Gadolinium ions block mechanosensitive channels by altering the packing and lateral pressure of anionic lipids

Effects of polyvalent ions on the lateral packing of phospholipids have been known for decades, but the physiological consequences have not been systematically studied. Gd(3+) is a relatively nonspecific agent that blocks mechano-gated channels with a variable affinity. In this study, we show that the large mechanosensitive channel MscL of Escherichia coli is effectively blocked by Gd(3+) only when reconstituted with negatively charged phospholipids (e.g., PS). Taking this lead, we studied effects of Gd(3+) on monolayers and unilamellar vesicles made of natural brain PS, DMPS, and its mixtures with DMPC. In monolayer experiments, we found that muM Gd(3+) present in the subphase leads to approximately 8% lateral compaction of brain PS (at 35 mN/m). Gd(3+) more strongly shrinks and rigidifies DMPS films causing a spontaneous liquid expanded-to-compact transition to the limiting 40 A(2)/mol. Pressure-area isotherms of uncharged DMPC were unaffected by Gd(3+), and neutralization of DMPS surface by low pH did not produce strong compaction. Upshifts of surface potential isotherms of DMPS monolayers reflected changes in the diffuse double layer due to neutralization of headgroup charges by Gd(3+), whereas the increased packing density produced up to a 200 mV change in the interfacial dipole potential. The slopes of surface potential versus reciprocal area predicted that Gd(3+) induced a modest ( approximately 18%) increase in the magnitude of the individual lipid dipoles in DMPS. Isothermal titration calorimetry indicated that binding of Gd(3+) to DMPS liposomes in the gel state is endothermic, whereas binding to liquid crystalline liposomes produces heat consistent with the isothermal liquid-to-gel phase transition induced by the ion. Both titration curves suggested a K(b) of approximately 10(6) M(-1). We conclude that anionic phospholipids serve as high-affinity receptors for Gd(3+) ions, and the ion-induced compaction generates a lateral pressure increase estimated as tens of mN/m. This pressure can "squeeze" the channel and shift the equilibrium toward the closed state.

Mechanosensitive closed-closed transitions in large membrane proteins: osmoprotection and tension damping

Multiconformation membrane proteins are mechanosensitive (MS) if their conformations displace different bilayer areas. Might MS closed-closed transitions serve as tension buffers, that is, as membrane "spandex"? While bilayer expansion is effectively instantaneous, transitions of bilayer-embedded MS proteins are stochastic (thermally activated) so spandex kinetics would be critical. Here we model generic two-state (contracted/expanded) stochastic spandexes inspired by known bacterial osmovalves (MscL, MscS) then suggest experimental approaches to test for spandex-like behaviors in these proteins. Modeling shows: 1), spandex kinetics depend on the transition state location along an area reaction coordinate; 2), increasing membrane concentration of a spandex right-shifts its midpoint (= tension-Boltzmann); 3), spandexes with midpoints below the activating tension of an osmovalve could optimize osmovalve deployment (required: large midpoint, barrier near the expanded state); 4), spandexes could damp bilayer tension excursions (required: midpoint at target tension, and for speed, barrier halfway between the contracted and expanded states; the larger the spandex Delta-area, the more precise the maintenance of target tension; higher spandex concentrations damp larger amplitude strain fluctuations). One spandex species could not excel as both first line of defense for osmovalve partners and tension damper. Possible interactions among MS closed-closed and closed-open transitions are discussed for MscS- and MscL-like proteins.

Cloning and functional expression of an MscL ortholog from Rhizobium etli: characterization of a mechanosensitive channel

Rhizobium etli is equipped with several systems to handle both hyper- and hypo-osmotic stress. For adaptation to hypo-osmotic stress, R. etli possesses a single gene with clear homology to MscS, four MscS-like channels and one ortholog of MscL (ReMscL, identity approximately 44% compared to Escherichia coli MscL). We subcloned and expressed the ReMscL channel ortholog from R. etli in E. coli to examine its activity by patch clamp in giant spheroplasts and characterized it at the single-channel level. We obtained evidence that ReMscL prevents the lysis of E. coli null mutant log-phase cells upon a rapid, osmotic downshock and identified a slight pH dependence for ReMscL activation. Here, we describe the facilitation of ReMscL activation by arachidonic acid (AA) and a reversible inhibitory effect of Gd(3+). The results obtained in these experiments suggest a stabilizing effect of micromolar AA and traces of Gd(3+) ions in the partially expanded conformation of the protein. Finally, we discuss a possible correlation between the number of gene paralogs for MS channels and the habitats of several microorganisms. Taken together, our data show that ReMscL may play an important role in free-living rhizobacteria during hypo-osmotic shock in the rhizosphere.

Phase-transition-induced protein redistribution in lipid bilayers

We report an atomic force microscopy study on the lateral spatial redistribution of an integral membrane protein reconstituted in supported lipid bilayers (SLBs) subjected to a thermally induced phase transition. KcsA proteins were reconstituted in proteoliposomes of POPE/POPG (3:1, mol/mol), and SLBs, including the proteins, were then obtained by the vesicle fusion technique on mica. By decreasing the temperature, the lipid bilayer passed from a liquid disordered (l(d)) phase in which the proteins are homogeneously distributed to a coexistence of solid ordered (s(o)) and l(d) domains with the proteins preferentially distributed in the l(d) domains. The inhomogeneous distribution eventually led to protein clustering. The obtained results are discussed in light of the role that the lipid/protein interaction can have in determining the function of integral membrane proteins.

Structural changes in the cytoplasmic domain of the mechanosensitive channel MscS during opening

The bacterial mechanosensitive channel MscS forms a homoheptamer of subunits composed of a transmembrane (TM) domain and a large cytoplasmic (CP) domain. Recent studies suggest that a lateral expansion of the TM domain, structural change in the CP domain, and TM-CP interactions are essential to open the channel. However, it has not been examined whether the CP domain undergoes structural changes during channel opening. The aim of this study was to estimate structural changes in the CP domain during channel opening using fluorescence resonance energy transfer (FRET) spectroscopy. To monitor changes in the horizontal diameter of the CP domain, four point mutants (A132C, F178C, L246C, and R259C), all of which had channel activity, were created and labeled with Alexa488 and Alexa568 for FRET analysis. The FRET efficiency of these mutants decreased when lysophosphatidylcholine was applied to open the channel, suggesting that the CP domain swells up when the channel opens. The degree of the decease in FRET efficiency after lysophosphatidylcholine treatment was smaller in the D62N/F178C mutant, which was deficient in the TM-CP interactions, than in the F178C mutant. These findings provide the first, to our knowledge, experimental evidence that the CP domain swells up during channel opening, and the swelling is mediated by the TM-CP interactions.

Structure of a tetrameric MscL in an expanded intermediate state

The ability of cells to sense and respond to mechanical force underlies diverse processes such as touch and hearing in animals, gravitropism in plants, and bacterial osmoregulation1, 2. In bacteria, mechanosensation is mediated by the mechanosensitive channels of large (MscL), small (MscS), potassium-dependent (MscK), and mini (MscM) conductances. These channels act as “emergency relief valves” protecting bacteria from lysis upon acute osmotic downshock3. Among them, MscL has been intensively studied since the original identification and characterization 15 years ago by Kung and co-workers4. MscL is reversibly and directly gated by changes in membrane tension. In the open state, MscL forms a nonselective 3 nS-conductance channel which gates at tensions close to the lytic limit of the bacterial membrane. An earlier crystal structure at 3.5 Å resolution of a pentameric MscL from Mycobacterium tuberculosis (TbMscL) represents a closed-state or nonconducting conformation5, 6. MscL has a complex gating behaviour; it exhibits several intermediates between the closed and open states, including one putative nonconductive expanded state and at least three sub-conducting states7. Although our understanding of the closed5, 6 and open8-10 states of MscL has been increasing, little is known about the structures of the intermediate states despite their importance in elucidating the complete gating process of MscL. Here we present the crystal structure of a truncation mutant (Δ95-120) of MscL from Staphylococcus aureus (SaMscL-CΔ26) at 3.8 Å resolution. Strikingly, SaMscL-CΔ26 forms a tetrameric channel with both transmembrane (TM) helices tilted away from the membrane normal at angles close to that inferred for the open state9, likely corresponding to a nonconductive but partially expanded intermediate state.

Emerging roles for lipids in shaping membrane-protein function

Studies of membrane proteins have revealed a direct link between the lipid environment and the structure and function of some of these proteins. Although some of these effects involve specific chemical interactions between lipids and protein residues, many can be understood in terms of protein-induced perturbations to the membrane shape. The free-energy cost of such perturbations can be estimated quantitatively, and measurements of channel gating in model systems of membrane proteins with their lipid partners are now confirming predictions of simple models.

A simple automated stimulator of mechanically induced arrhythmias in the isolated rat heart

Transient stretching of the ventricle can trigger arrhythmias and evoke ventricular fibrillation, especially when the stimulation occurs in the vulnerable period. To explore the sensitivity of small hearts we used a commercial pressure servo to study the kinetic relationship of left ventricular pressure to excitability and arrhythmias in the rat heart. Stimulation protocols were readily composed on the computer and programmed to vary the stimulus amplitude and timing relative to pacing. The pressure-induced premature ventricular excitations were similar to those observed in larger hearts, but the convenience of using small hearts allows the use of inexpensive transgenic animals to explore the molecular basis of transduction.

Conserved motifs in mechanosensitive channels MscL and MscS

Mechanosensitive (MS) channels play a major role in protecting bacterial cells against hypo-osmotic shock. To understand their function, it is important to identify the conserved motifs using sequence analysis methods. In this study, the sequence conservation was investigated by an in silico analysis to generate sequence logos. We have identified new conserved motifs in the domains TM1, TM2 and the cytoplasmic helix from 231 homologs of MS channel of large conductance (MscL). In addition, we have identified new motifs for the TM3 and the cytoplasmic carboxy-terminal domain from 309 homologs of MS channel of small conductance (MscS). We found that the conservation in MscL homologs is high for TM1 and TM2 in the three domains of life. The conservation in MscS homologs is high only for TM3 in Bacteria and Archaea.

State-stabilizing Interactions in Bacterial Mechanosensitive Channel Gating and Adaptation

We outline several principles that we believe define the gating of two bacterial mechanosensitive channels, MscL and MscS. Serving as turgor regulators in bacteria and other walled cells, these molecules are tangible models for studying conformational transitions in membrane proteins driven directly by membrane tension. MscL, a compact pentamer, reversibly opens a gigantic 30-A pore at near-lytic tensions. MscS, a heptameric complex, exhibits transient activation of a smaller pore at moderate tensions, thereby entering a tension-insensitive inactivated state. By comparing the structures and predicted transitions in these channels, we concluded that opening is commonly achieved through tilting and outward motion of the pore-lining helices, which is kinetically limited by hydration of the pore. The intricate adaptive behavior in MscS appears to depend on specific interhelical associations and the flexibility of the pore-lining helices. We discuss physical factors that may direct the transitions and stabilize main functional states in these channels.

NUMERICAL SIMULATION OF NANOINDENTATION AND PATCH CLAMP EXPERIMENTS ON MECHANOSENSITIVE CHANNELS OF LARGE CONDUCTANCE IN ESCHERICHIA COLI

A hierarchical simulation framework that integrates information from all-atom simulations into a finite element model at the continuum level is established to study the mechanical response of a mechanosensitive channel of large conductance (MscL) in bacteria Escherichia Coli (E.coli) embedded in a vesicle formed by the dipalmitoylphosphatidycholine (DPPC) lipid bilayer. Sufficient structural details of the protein are built into the continuum model, with key parameters and material properties derived from molecular mechanics simulations. The multi-scale framework is used to analyze the gating of MscL when the lipid vesicle is subjective to nanoindentation and patch clamp experiments, and the detailed structural transitions of the protein are obtained explicitly as a function of external load; it is currently impossible to derive such information based solely on all-atom simulations. The gating pathways of E.coli-MscL qualitatively agree with results from previous patch clamp experiments. The gating mechanisms under complex indentation-induced deformation are also predicted. This versatile hierarchical multi-scale framework may be further extended to study the mechanical behaviors of cells and biomolecules, as well as to guide and stimulate biomechanics experiments.

Voltage-dependent K+ channel gating and voltage sensor toxin sensitivity depend on the mechanical state of the lipid membrane

Voltage-dependent K(+) (Kv) channels underlie action potentials through gating conformational changes that are driven by membrane voltage. In this study of the paddle chimera Kv channel, we demonstrate that the rate of channel opening, the voltage dependence of the open probability, and the maximum achievable open probability depend on the lipid membrane environment. The activity of the voltage sensor toxin VsTx1, which interferes with voltage-dependent gating by partitioning into the membrane and binding to the channel, also depends on the membrane. Membrane environmental factors that influence channel function are divisible into two general categories: lipid compositional and mechanical state. The mechanical state can have a surprisingly large effect on the function of a voltage-dependent K(+) channel, including its pharmacological interaction with voltage sensor toxins. The dependence of VSTx1 activity on the mechanical state of the membrane leads us to hypothesize that voltage sensor toxins exert their effect by perturbing the interaction forces that exist between the channel and the membrane.

Membrane mechanics as a probe of ion-channel gating mechanisms

The details of conformational changes undergone by transmembrane ion channels in response to stimuli, such as electric fields and membrane tension, remain controversial. We approach this problem by considering how the conformational changes impose deformations in the lipid bilayer. We focus on the role of bilayer deformations in the context of voltage-gated channels because we hypothesize that such deformations are relevant in this case as well as for channels that are explicitly mechanosensitive. As a result of protein conformational changes, we predict that the lipid bilayer suffers deformations with a characteristic free-energy scale of 10 k{B}T . This free energy is comparable to the voltage-dependent part of the total gating energy, and we argue that these deformations could play an important role in the overall free-energy budget of gating. As a result, channel activity will depend upon mechanical membrane parameters such as tension and leaflet thickness. We further argue that the membrane deformation around any channel can be divided into three generic classes of deformation that exhibit different mechanosensitive properties. Finally, we provide the theoretical framework that relates conformational changes during gating to tension and leaflet thickness dependence in the critical gating voltage. This line of investigation suggests experiments that could discern the dominant deformation imposed upon the membrane as a result of channel gating, thus providing clues as to the channel deformation induced by the stimulus.

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° 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 (∼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.

The structure of an open form of an E. coli mechanosensitive channel at 3.45 A resolution

How ion channels are gated to regulate ion flux in and out of cells is the subject of intense interest. The Escherichia coli mechanosensitive channel, MscS, opens to allow rapid ion efflux, relieving the turgor pressure that would otherwise destroy the cell. We present a 3.45 angstrom-resolution structure for the MscS channel in an open conformation. This structure has a pore diameter of approximately 13 angstroms created by substantial rotational rearrangement of the three transmembrane helices. The structure suggests a molecular mechanism that underlies MscS gating and its decay of conductivity during prolonged activation. Support for this mechanism is provided by single-channel analysis of mutants with altered gating characteristics.

Mechanosensitive channels: insights from continuum-based simulations

Mechanotransduction plays an important role in regulating cell functions and it is an active topic of research in biophysics. Despite recent advances in experimental and numerical techniques, the intrinsic multiscale nature imposes tremendous challenges for revealing the working mechanisms of mechanosensitive channels. Recently, a continuum-mechanics-based hierarchical modeling and simulation framework has been established and applied to study the mechanical responses and gating behaviors of a prototypical mechanosensitive channel, the mechanosensitive channel of large conductance (MscL) in bacteria Escherichia coli (E. coli), from which several putative gating mechanisms have been tested and new insights are deduced. This article reviews these latest findings using the continuum mechanics framework and suggests possible improvements for future simulation studies. This computationally efficient and versatile continuum-mechanics-based protocol is poised to make contributions to the study of a variety of mechanobiology problems.

Atom-scale molecular interactions in lipid raft mixtures

We review the relationship between molecular interactions and the properties of lipid environments. A specific focus is given on bilayers which contain sphingomyelin (SM) and sterols due to their essential role for the formation of lipid rafts. The discussion is based on recent atom-scale molecular dynamics simulations, complemented by extensive comparison to experimental data. The discussion is divided into four sections. The first part investigates the properties of one-component SM bilayers and compares them to bilayers with phosphatidylcholine (PC), the focus being on a detailed analysis of the hydrogen bonding network in the two bilayers. The second part deals with binary mixtures of sterols with either SM or PC. The results show how the membrane properties may vary substantially depending on the sterol and SM type available, the membrane order and interdigitation being just two of the many examples of this issue. The third part concentrates on the specificity of intermolecular interactions in three-component mixtures of SM, PC and cholesterol (CHOL) under conditions where the concentrations of SM and CHOL are dilute with respect to that of PC. The results show how SM and CHOL favor one another, thus acting as nucleation sites for the formation of highly ordered nanosized domains. Finally, the fourth part discusses the large-scale properties of raft-like membrane environments and compares them to the properties of non-raft membranes. The differences turn out to be substantial. As a particularly intriguing example of this, the lateral pressure profiles of raft-like and non-raft systems indicate that the lipid composition of membrane domains may have a major impact on membrane protein activation.

Bacterial mechanosensitive channels: Experiment and theory

Since their discovery in Escherichia coli some 20 years ago, studies of bacterial mechanosensitive (MS) ion channels have been at the forefront of the MS channel research field. Two major events greatly advanced the research on bacterial MS channels: (i) cloning of MscL and MscS, the MS channels of Large and Small conductance, and (ii) solving their 3D crystal structure. These events enabled further experimental studies employing EPR and FRET spectroscopy in addition to patch clamp and molecular biological techniques that have successfully been used in characterization of the structure and function of bacterial MS channels. In parallel with the experimental studies computational modelling has been applied to elucidate the molecular dynamics of MscL and MscS, which has significantly contributed to our understanding of basic physical principles of the mechanosensory transduction in living organisms.

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 mechanisms of mechanosensitive channels of large conductance, II: systematic study of conformational transitions

Part II of this study is based on the continuum mechanics-based molecular dynamics-decorated finite element method (MDeFEM) framework established in Part I. In Part II, the gating pathways of Escherichia coli-MscL channels under various basic deformation modes are simulated. Upon equibiaxial tension (which is verified to be the most effective mode for gating), the MDeFEM results agree well with both experiments and all-atom simulations in literature, as well as the analytical continuum models and elastic network models developed in Part I. Different levels of model sophistication and effects of structural motifs are explored in detail, where the importance of mechanical roles of transmembrane helices, cytoplasmic helices, and loops are discussed. The conformation transitions under complex membrane deformations are predicted, including bending, torsion, cooperativity, patch clamp, and indentation. Compared to atom-based molecular dynamics simulations and elastic network models, the MDeFEM framework is unusually well-suited for simulating complex deformations at large length scales. The versatile hierarchical framework can be further applied to simulate the gating transition of other mechanosensitive channels and other biological processes where mechanical perturbation is important.

Energetics of hydrophobic matching in lipid-protein interactions

Lipid chain length modulates the activity of transmembrane proteins by mismatch between the hydrophobic span of the protein and that of the lipid membrane. Relative binding affinities of lipids with different chain lengths are used to estimate the excess free energy of lipid-protein interaction that arises from hydrophobic mismatch. For a wide range of integral proteins and peptides, the energy cost is much less than the elastic penalty of fully stretching or compressing the lipid chains to achieve complete hydrophobic matching. The chain length dependences of the free energies of lipid association are described by a model that combines elastic chain extension with a free energy term that depends linearly on the extent of residual mismatch. The excess free energy densities involved lie in the region of 0.5-2.0 k(B)T x nm(-2). Values of this size could arise from exposure of hydrophobic groups to polar portions of the lipid or protein, but not directly to water, or alternatively from changes in tilt of the transmembrane helices that are energetically comparable to those activating mechanosensitive channels. The influence of hydrophobic mismatch on dimerization of transmembrane helices and their transfer between lipid vesicles, and on shifts in chain-melting transitions of lipid bilayers by incorporated proteins, is analyzed by using the same thermodynamic model. Segmental order parameters confirm that elastic lipid chain distortions are insufficient to compensate fully for the mismatch, but the dependence on chain length with tryptophan-anchored peptides requires that the free energy density of hydrophobic mismatch should increase with increasing extent of mismatch.

Mechanosensitive membrane channels in action

The tension-driven gating process of MscL from Mycobacterium tuberculosis, Tb-MscL, has been addressed at near-atomic detail using coarse-grained molecular dynamics simulations. To perform the simulations, a novel coarse-grained peptide model based on a thermodynamic parameterization of the amino-acid side chains has been applied. Both the wild-type Tb-MscL and its gain-of-function mutant V21D embedded in a solvated lipid bilayer have been studied. To mimic hypoosmotic shock conditions, simulations were performed at increasing levels of membrane tension approaching the rupture threshold of the lipid bilayer. Both the wild-type and the mutant channel are found to undergo significant conformational changes in accordance with an irislike expansion mechanism, reaching a conducting state on a microsecond timescale. The most pronounced expansion of the pore has been observed for the V21D mutant, which is consistent with the experimentally shown gain-of-function phenotype of the V21D mutant.

Gating mechanisms of mechanosensitive channels of large conductance, I: a continuum mechanics-based hierarchical framework

A hierarchical simulation framework that integrates information from molecular dynamics (MD) simulations into a continuum model is established to study the mechanical response of mechanosensitive channel of large-conductance (MscL) using the finite element method (FEM). The proposed MD-decorated FEM (MDeFEM) approach is used to explore the detailed gating mechanisms of the MscL in Escherichia coli embedded in a palmitoyloleoylphosphatidylethanolamine lipid bilayer. In Part I of this study, the framework of MDeFEM is established. The transmembrane and cytoplasmic helices are taken to be elastic rods, the loops are modeled as springs, and the lipid bilayer is approximated by a three-layer sheet. The mechanical properties of the continuum components, as well as their interactions, are derived from molecular simulations based on atomic force fields. In addition, analytical closed-form continuum model and elastic network model are established to complement the MDeFEM approach and to capture the most essential features of gating. In Part II of this study, the detailed gating mechanisms of E. coli-MscL under various types of loading are presented and compared with experiments, structural model, and all-atom simulations, as well as the analytical models established in Part I. It is envisioned that such a hierarchical multiscale framework will find great value in the study of a variety of biological processes involving complex mechanical deformations such as muscle contraction and mechanotransduction.

Anionic phospholipids affect the rate and extent of flux through the mechanosensitive channel of large conductance MscL

The mechanosensitive channel of large conductance MscL from Escherichia coli has been reconstituted into sealed vesicles, and the effects of lipid structure on the flux of the fluorescent molecule calcein through the open channel have been studied. The channel was opened by reaction of the G22C mutant of MscL with the reagent [2-(triethylammonium)ethyl]methanethiosulfonate (MTSET) which introduces five positive charges within the pore constriction. Flux through the channel was small when the lipid was phosphatidylcholine, but addition of the anionic lipids phosphatidylglycerol, phosphatidic acid, or cardiolipin up to 50 mol % resulted in increases in the amplitudes and rates of release of calcein. Similar effects were seen when either wild-type MscL or the G22C mutant was opened by osmotic pressure difference; rates of release of calcein were very slow in the absence of anionic lipid but increased with increasing concentrations of phosphatidylglycerol to 50 mol %. The observed partial release of trapped calcein following activation of MscL was attributed to the formation of a long-lived subconductance state of MscL following channel opening. Effects of anionic lipid were attributed to an increase in the rate of the transition from closed to fully open state and to a decrease in the rate of the transition from the fully open state to the subconductance state. Higher concentrations of anionic lipid led to a decrease in the rate and amplitude of release of calcein, possibly due to a decreased rate of flux through the open channel. In mixtures with anionic lipids, phosphatidylethanolamine resulted in lower rates and amplitude of release than phosphatidylcholine.

Mechanical and biochemical modeling of cortical oscillations in spreading cells

Actomyosin-based cortical contractility is a common feature of eukaryotic cells and is involved in cell motility, cell division, and apoptosis. In nonmuscle cells, oscillations in contractility are induced by microtubule depolymerization during cell spreading. We developed an ordinary differential equation model to describe this behavior. The computational model includes 36 parameters. The values for all but two of the model parameters were taken from experimental measurements found in the literature. Using these values, we demonstrate that the model generates oscillatory behavior consistent with current experimental observations. The rhythmic behavior occurs because of the antagonistic effects of calcium-induced contractility and stretch-activated calcium channels. The model makes several experimentally testable predictions: 1), buffering intracellular calcium increases the period and decreases the amplitude of cortical oscillations; 2), increasing the number or activity of stretch activated channels leads to an increase in period and amplitude of cortical oscillations; 3), inhibiting Ca(2+) pump activity increases the period and amplitude of oscillations; and 4), a threshold exists for the calcium concentration below which oscillations cease.

Gating-associated conformational changes in the mechanosensitive channel MscL

Bacterial cells avoid lysis in response to hypoosmotic shock through the opening of the mechanosensitive channel MscL. Upon channel opening, MscL is thought to expand in the plane of the membrane and form a large pore with an estimated diameter of 3-4 nm. Here, we set out to analyze the closed and open structure of cell-free MscL. To this end, we characterized the function and structure of wild-type MscL and a mutant form of the protein (G22N MscL) that spontaneously adopts an open substate. Patch-clamp analysis of MscL that had been reconstituted into liposomes revealed that wild-type MscL was activated only by mechanical stimuli, whereas G22N MscL displayed spontaneous opening to the open substate. In accord with these results, Ca(2+) influx into G22N MscL-containing liposomes occurred in the absence of mechanical stimulation. The electrophoretic migration of chemically cross-linked G22N MscL was slower than that of cross-linked wild-type MscL, suggesting that G22N MscL is in an expanded form. Finally, electron microscopy using low-angle rotary shadowing revealed the presence of a pore at the center of G22N MscL. No pore could be detected in wild-type MscL. However, wild-type MscL possessed a protrusion at one end, which was absent in G22N MscL. The deletion of carboxyl-terminal 27 residues resulted in the loss of protrusion and proper multimerization. The structures of wild-type and G22N MscL reveal that the opening of MscL is accompanied by the dissociation of a carboxyl-terminal protrusion and pore formation.

Actin stress fibers transmit and focus force to activate mechanosensitive channels

Mechanosensitive (MS) channels are expressed in various cells in a wide range of phylogenetic lineages from bacteria to humans. Understanding the molecular and biophysical mechanisms of their activation is an important research pursuit. It is controversial whether eukaryotic MS channels need accessory proteins – typically cytoskeletal structures – for activation, because MS channel activities are modulated by pharmacological treatments that affect the cytoskeleton. Here we demonstrate that direct mechanical stimulation (stretching) of an actin stress fiber using optical tweezers can activate MS channels in cultured human umbilical vein endothelial cells (HUVECs). Furthermore, by using high-speed total internal reflection microscopy, we visualized spots of Ca2+ influx across individual MS channels distributed near focal adhesions in the basal surface of HUVECs. This study provides the first direct evidence that the cytoskeleton works as a force-transmitting and force-focusing molecular device to activate MS channels in eukaryotic cells.

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 A 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.

Nanomagnetic actuation of receptor-mediated signal transduction

Complex cell behaviours are triggered by chemical ligands that bind to membrane receptors and alter intracellular signal transduction. However, future biosensors, medical devices and other microtechnologies that incorporate living cells as system components will require actuation mechanisms that are much more rapid, robust, non-invasive and easily integrated with solid-state interfaces. Here we describe a magnetic nanotechnology that activates a biochemical signalling mechanism normally switched on by binding of multivalent chemical ligands. Superparamagnetic 30-nm beads, coated with monovalent ligands and bound to transmembrane receptors, magnetize when exposed to magnetic fields, and aggregate owing to bead-bead attraction in the plane of the membrane. Associated clustering of the bound receptors acts as a nanomagnetic cellular switch that directly transduces magnetic inputs into physiological cellular outputs, with rapid system responsiveness and non-invasive dynamic control. This technique may represent a new actuator mechanism for cell-based microtechnologies and man-machine interfaces.