The effect of local bending on gating of MscL using a representative volume element and finite element simulation

Many physiological processes such as cell division, endocytosis and exocytosis cause severe local curvature of the cell membrane. Local curvature has been shown experimentally to modulate numerous mechanosensitive (MS) ion channels. In order to quantify the effects of local curvature we introduced a coarse grain representative volume element for the bacterial mechanosensitive ion channel of large conductance (MscL) using continuum elasticity. Our model is designed to be consistent with the channel conformation in the closed and open states to capture its major continuum rheological behavior in response to the local membrane curvature. Herein we show that change in the local curvature of the lipid bilayer can modulate MscL activity considerably by changing both bilayer thickness and lateral pressure profile. Intriguingly, although bending in any direction results in almost the same free-energy cost, inward (cytoplasmic) bending favors channel opening, whereas outward (periplasmic) bending facilitates closing of the narrowest part of the MscL pore. This quantitative study using MscL as a model channel may have wide reaching consequences for the effect of local curvature on the physiological function of other types of prokaryotic and eukaryotic membrane proteins.     

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.     

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.     

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.     

Membrane-protein interactions in mechanosensitive channels

In this article, we examine the mechanical role of the lipid bilayer in ion channel conformation and function with specific reference to the case of the mechanosensitive channel of large conductance (MscL). In a recent article we argued that mechanotransduction very naturally arises from lipid-protein interactions by invoking a simple analytic model of the MscL channel and the surrounding lipid bilayer. In this article, we focus on improving and expanding this analytic framework for studying lipid-protein interactions with special attention to MscL. Our goal is to generate simple scaling relations which can be used to provide qualitative understanding of the role of membrane mechanics in protein function and to quantitatively interpret experimental results. For the MscL channel, we find that the free energies induced by lipid-protein interaction are of the same order as the measured free energy differences between conductance states. We therefore conclude that the mechanics of the bilayer plays an essential role in determining the conformation and function of the channel. Finally, we compare the predictions of our model to experimental results from the recent investigations of the MscL channel by a variety of investigators and suggest a suite of new experiments.     

A finite element framework for studying the mechanical response of macromolecules: application to the gating of the mechanosensitive channel MscL

The gating pathways of mechanosensitive channels of large conductance (MscL) in two bacteria (Mycobacterium tuberculosis and Escherichia coli) are studied using the finite element method. The phenomenological model treats transmembrane helices as elastic rods and the lipid membrane as an elastic sheet of finite thickness; the model is inspired by the crystal structure of MscL. The interactions between various continuum components are derived from molecular-mechanics energy calculations using the CHARMM all-atom force field. Both bacterial MscLs open fully upon in-plane tension in the membrane and the variation of pore diameter with membrane tension is found to be essentially linear. The estimated gating tension is close to the experimental value. The structural variations along the gating pathway are consistent with previous analyses based on structural models with experimental constraints and biased atomistic molecular-dynamics simulations. Upon membrane bending, neither MscL opens substantially, although there is notable and nonmonotonic variation in the pore radius. This emphasizes that the gating behavior of MscL depends critically on the form of the mechanical perturbation and reinforces the idea that the crucial gating parameter is lateral tension in the membrane rather than the curvature of the membrane. Compared to popular all-atom-based techniques such as targeted or steered molecular-dynamics simulations, the finite element method-based continuum-mechanics framework offers a unique alternative to bridge detailed intermolecular interactions and biological processes occurring at large spatial scales and long timescales. 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.     

Novel splice variants of amphiphysin I are expressed in retina

The MscL channel is a mechanosensitive channel which is gated by membrane stress or tension. Here, we describe a series of simulations which apply simulated mechanical stress to a molecular model of the MscL channel using two methods - direct force application to the transmembrane segments, and anisotropic pressure coupling. In the latter simulations, pressures less than that equivalent to a bilayer tension of 12 dyn/cm did not cause the channel to open, while pressures in excess of this value resulted in the channel opening. These results are in approximate agreement with experimental findings.     

Bilayer-mediated clustering and functional interaction of MscL channels

Mechanosensitive channels allow bacteria to respond to osmotic stress by opening a nanometer-sized pore in the cellular membrane. Although the underlying mechanism has been thoroughly studied on the basis of individual channels, the behavior of channel ensembles has yet to be elucidated. This work reveals that mechanosensitive channels of large conductance (MscL) exhibit a tendency to spatially cluster, and demonstrates the functional relevance of clustering. We evaluated the spatial distribution of channels in a lipid bilayer using patch-clamp electrophysiology, fluorescence and atomic force microscopy, and neutron scattering and reflection techniques, coupled with mathematical modeling of the mechanics of a membrane crowded with proteins. The results indicate that MscL forms clusters under a wide range of conditions. MscL is closely packed within each cluster but is still active and mechanosensitive. However, the channel activity is modulated by the presence of neighboring proteins, indicating membrane-mediated protein-protein interactions. Collectively, these results suggest that MscL self-assembly into channel clusters plays an osmoregulatory functional role in the membrane.     

Force Transduction and Lipid Binding in MscL: A Continuum-Molecular Approach

The bacterial mechanosensitive channel MscL, a small protein mainly activated by membrane tension, is a central model system to study the transduction of mechanical stimuli into chemical signals. Mutagenic studies suggest that MscL gating strongly depends on both intra-protein and interfacial lipid-protein interactions. However, there is a gap between this detailed chemical information and current mechanical models of MscL gating. Here, we investigate the MscL bilayer-protein interface through molecular dynamics simulations, and take a combined continuum-molecular approach to connect chemistry and mechanics. We quantify the effect of membrane tension on the forces acting on the surface of the channel, and identify interactions that may be critical in the force transduction between the membrane and MscL. We find that the local stress distribution on the protein surface is largely asymmetric, particularly under tension, with the cytoplasmic side showing significantly larger and more localized forces, which pull the protein radially outward. The molecular interactions that mediate this behavior arise from hydrogen bonds between the electronegative oxygens in the lipid headgroup and a cluster of positively charged lysine residues on the amphipathic S1 domain and the C-terminal end of the second trans-membrane helix. We take advantage of this strong interaction (estimated to be 10–13 kT per lipid) to actuate the channel (by applying forces on protein-bound lipids) and explore its sensitivity to the pulling magnitude and direction. We conclude by highlighting the simple motif that confers MscL with strong anchoring to the bilayer, and its presence in various integral membrane proteins including the human mechanosensitive channel K2P1 and bovine rhodopsin.     

Curvature Generation and Pressure Profile Modulation in Membrane by Lysolipids: Insights from Coarse-Grained Simulations

Although many membrane additives are known to modulate the activities of membrane proteins via perturbing the properties of lipid membrane, the underlying mechanism is often not precisely understood. In this study, we investigate the impact of asymmetrically incorporating single-tailed lysophosphatidylcholine (LPC) into a membrane bilayer using coarse-grained molecular dynamics simulations. Using a simple computational protocol designed to approximately mimic a micropipette setting, we show that asymmetric incorporation of LPC can lead to significant curvature in a bilayer. Detailed analysis of geometrical and mechanical properties (pressure profile) of the resulting mound structure indicates that the degree of pressure profile perturbation is determined not by the local curvature per se but by the packing of lipid headgroups (i.e., area-per-lipid). The findings help provide a concrete basis for understanding the activation mechanism of mechanosensitive channels by asymmetric incorporation of LPC into membrane patches in patch-clamp experiments. The calculated local pressure profiles are valuable to the construction of realistic membrane models for the analysis of mechanosensation in a continuum mechanics framework.     

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.     

In vitro synthesis and oligomerization of the mechanosensitive channel of large conductance, MscL, into a functional ion channel

Elucidation of high-resolution structures of integral membrane proteins is drastically lagging behind that of cytoplasmic proteins. In vitro synthesis and insertion of membrane proteins into synthetic membranes could circumvent bottlenecks associated with the overexpression of membrane proteins, producing sufficient membrane-inserted, correctly folded protein for structural studies. Using the mechanosensitive channel of large conductance, MscL, as a model protein we show that in vitro synthesized MscL inserts into YidC-containing proteoliposomes and oligomerizes to form a homopentamer. Using planar membrane bilayers, we show that MscL forms functional ion channels capable of ion transport. These data demonstrate that membrane insertion of MscL is YidC mediated, whereas subsequent oligomerization into a functional homopentamer is a spontaneous event.     

Correlating a protein structure with function of a bacterial mechanosensitive channel

MscL, a mechanosensitive channel found in many bacteria, protects cells from hypotonic shock by reducing intracellular pressure through release of cytoplasmic osmolytes. First isolated from Escherichia coli, this protein has served as a model for how a protein senses and responds to membrane tension. Recently the structure of a functionally uncharacterized MscL homologue from Mycobacterium tuberculosis was solved by x-ray diffraction to a resolution of 3.5 A. Here we demonstrate that the protein forms a functional MscL-like mechanosensitive channel in E. coli membranes and azolectin proteoliposomes. Furthermore, we show that M. tuberculosis MscL crystals, when re-solubilized and reconstituted, yield wild-type channel currents in patch clamp, demonstrating that the protein does not irreversibly change conformation upon crystallization. Finally, we apply functional clues acquired from the E. coli MscL to the M. tuberculosis channel and show a mechanistic correlation between these channels. However, the inability of the M. tuberculosis channel to gate at physiological membrane tensions, demonstrated by in vivo E. coli expression and in vitro reconstitution, suggests that the membrane environment or other additional factors influence the gating of this channel.     

Site-directed spin-labeling analysis of reconstituted Mscl in the closed state

The mechanosensitive channel from Escherichia coli (Eco-MscL) responds to membrane lateral tension by opening a large, water-filled pore that serves as an osmotic safety valve. In an attempt to understand the structural dynamics of MscL in the closed state and under physiological conditions, we have performed a systematic site-directed spin labeling study of this channel reconstituted in a membrane bilayer. Structural information was derived from an analysis of probe mobility, residue accessibility to O(2) or NiEdda and overall intersubunit proximity. For the majority of the residues studied, mobility and accessibility data showed a remarkable agreement with the Mycobacterium tuberculosis crystal structure, clearly identifying residues facing the large water-filled vestibule at the extracellular face of the molecule, the narrowest point along the permeation pathway (residues 21-26 of Eco-MscL), and the lipid-exposed residues in the peripheral transmembrane segments (TM2). Overall, the present dataset demonstrates that the transmembrane regions of the MscL crystal structure (obtained in detergent and at low pH) are, in general, an accurate representation of its structure in a membrane bilayer under physiological conditions. However, significant differences between the EPR data and the crystal structure were found toward the COOH-terminal end of TM2.     

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.     

Multifunctional,Micropipette-based Method for Incorporation And Stimulation of Bacterial Mechanosensitive Ion Channels in Droplet Interface Bilayers

MscL, a large conductance mechanosensitive channel (MSC), is a ubiquitous osmolyte release valve that helps bacteria survive abrupt hypo-osmotic shocks. It has been discovered and rigorously studied using the patch-clamp technique for almost three decades. Its basic role of translating tension applied to the cell membrane into permeability response makes it a strong candidate to function as a mechanoelectrical transducer in artificial membrane-based biomolecular devices. Serving as building blocks to such devices, droplet interface bilayers (DIBs) can be used as a new platform for the incorporation and stimulation of MscL channels. Here, we describe a micropipette-based method to form DIBs and measure the activity of the incorporated MscL channels. This method consists of lipid-encased aqueous droplets anchored to the tips of two opposing (coaxially positioned) borosilicate glass micropipettes. When droplets are brought into contact, a lipid bilayer interface is formed. This technique offers control over the chemical composition and the size of each droplet, as well as the dimensions of the bilayer interface. Having one of the micropipettes attached to a harmonic piezoelectric actuator provides the ability to deliver a desired oscillatory stimulus. Through analysis of the shapes of the droplets during deformation, the tension created at the interface can be estimated. Using this technique, the first activity of MscL channels in a DIB system is reported. Besides MS channels, activities of other types of channels can be studied using this method, proving the multi-functionality of this platform. The method presented here enables the measurement of fundamental membrane properties, provides a greater control over the formation of symmetric and asymmetric membranes, and is an alternative way to stimulate and study mechanosensitive channels.     

Molecular dynamics study of MscL interactions with a curved lipid bilayer

Mechanosensitivity is a ubiquitous sensory mechanism found in living organisms. The simplest known mechanotransducing mechanism is found in bacteria in the form of the mechanosensitive membrane channel of large conductance, MscL. This channel has been studied extensively using a variety of methods at a functional and structural level. The channel is gated by membrane tension in the lipid bilayer alone. It serves as a safety valve protecting bacterial cells against hypoosmotic shock. MscL of Escherichia coli embedded in bilayers composed of asymmetric amounts of single-tailed and double-tailed lipids has been shown to gate spontaneously, even in the absence of membrane tension. To gain insight into the effect of the lipid membrane composition and geometry on MscL structure, a fully solvated, all-atom model of MscL in a stress-free curved bilayer composed of double- and single-tailed lipids was studied using a 9.5-ns molecular dynamics simulation. The bilayer was modeled as a domed structure accommodating the asymmetric composition of the monolayers. During the course of the simulation a spontaneous restructuring of the periplasmic loops occurred, leading to interactions between one of the loops and phospholipid headgroups. Previous experimental studies of the role of the loops agree with the observation that opening starts with a restructuring of the periplasmic loop, suggesting an effect of the curved bilayer. Because of limited resources, only one simulation of the large system was performed. However, the results obtained suggest that through the geometry and composition of the bilayer the protein structure can be affected even on short timescales.     

Thermodynamics of mechanosensitivity

Mechanosensitivity of ion channels is conventionally interpreted as being driven by a change of their in-plane cross-sectional area A(msc). This, however, does not include any factors relating to membrane stiffness, thickness, spontaneous curvature or changes in channel shape, length or stiffness. Because the open probability of a channel is sensitive to all these factors, we constructed a general thermodynamic formalism. These equations provide the basis for the analysis of the behaviour of mechanosensitive channels in lipids of different geometric and chemical properties such as the hydrophobic mismatch at the boundary between the protein and lipid or the effects of changes in the bilayer intrinsic curvature caused by the adsorption of amphipaths. This model predicts that the midpoint gamma(1/2) and the slope(1/2) of the gating curve are generally not independent. Using this relationship, we have predicted the line tension at the channel/lipid border of MscL as approximately 10 pN, and found it to be much less than the line tension of aqueous pores in pure lipid membranes. The MscL channel appears quite well matched to its lipid environment. Using gramicidin as a model system, we have explained its observed conversion from stretch-activated to stretch-inactivated gating as a function of bilayer thickness and composition.     

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_BT 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_BT on the gating energy of MscL. This contribution arises largely from the interfacial tension between hydrophobic and hydrophilic regions in a lipid bilayer.     

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.