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

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

From membrane tension to channel gating: A principal energy transfer mechanism for mechanosensitive channels

Mechanosensitive (MS) channels are evolutionarily conserved membrane proteins that play essential roles in multiple cellular processes, including sensing mechanical forces and regulating osmotic pressure. Bacterial MscL and MscS are two prototypes of MS channels. Numerous structural studies, in combination with biochemical and cellular data, provide valuable insights into the mechanism of energy transfer from membrane tension to gating of the channel. We discuss these data in a unified two‐state model of thermodynamics. In addition, we propose a lipid diffusion‐mediated mechanism to explain the adaptation phenomenon of MscS.     

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.     

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.     

On the conformation of the COOH-terminal domain of the large mechanosensitive channel MscL

COOH-terminal (S3) domains are conserved within the MscL family of bacterial mechanosensitive channels, but their function remains unclear. The X-ray structure of MscL from Mycobacterium tuberculosis (TbMscL) revealed cytoplasmic domains forming a pentameric bundle (Chang, G., R.H. Spencer, A.T. Lee, M.T. Barclay, and D.C. Rees. 1998. SCIENCE: 282:2220-2226). The helices, however, have an unusual orientation in which hydrophobic sidechains face outside while charged residues face inside, possibly due to specific crystallization conditions. Based on the structure of pentameric cartilage protein, we modeled the COOH-terminal region of E. coli MscL to better satisfy the hydrophobicity criteria, with sidechains of conserved aliphatic residues all inside the bundle. Molecular dynamic simulations predicted higher stability for this conformation compared with one modeled after the crystal structure of TbMscL, and suggested distances for disulfide trapping experiments. The single cysteine mutants L121C and I125C formed dimers under ambient conditions and more so in the presence of an oxidant. The double-cysteine mutants, L121C/L122C and L128C/L129C, often cross-link into tetrameric and pentameric structures, consistent with the new model. Patch-clamp examination of these double mutants under moderately oxidizing or reducing conditions indicated that the bundle cross-linking neither prevents the channel from opening nor changes thermodynamic parameters of gating. Destabilization of the bundle by replacing conservative leucines with small polar residues, or complete removal of COOH-terminal domain (Delta110-136 mutation), increased the occupancy of subconducting states but did not change gating parameters substantially. The Delta110-136 truncation mutant was functional in in vivo osmotic shock assays; however, the amount of ATP released into the shock medium was considerably larger than in controls. The data strongly suggest that in contrast to previous gating models (Sukharev, S., M. Betanzos, C.S. Chiang, and H.R. Guy. 2001a. NATURE: 409:720-724.), S3 domains are stably associated in both closed and open conformations. The bundle-like assembly of cytoplasmic helices provides stability to the open conformation, and may function as a size-exclusion filter at the cytoplasmic entrance to the MscL pore, preventing loss of essential metabolites.     

Toward elucidating the heat activation mechanism of the TRPV1 channel gating by molecular dynamics simulation

As a key cellular sensor, the TRPV1 cation channel undergoes a gating transition from a closed state to an open state in response to various physical and chemical stimuli including noxious heat. Despite years of study, the heat activation mechanism of TRPV1 gating remains enigmatic at the molecular level. Toward elucidating the structural and energetic basis of TRPV1 gating, we have performed molecular dynamics (MD) simulations (with cumulative simulation time of 3 μs), starting from the high‐resolution closed and open structures of TRPV1 solved by cryo‐electron microscopy. In the closed‐state simulations at 30°C, we observed a stably closed channel constricted at the lower gate (near residue I679), while the upper gate (near residues G643 and M644) is dynamic and undergoes flickery opening/closing. In the open‐state simulations at 60°C, we found higher conformational variation consistent with a large entropy increase required for the heat activation, and both the lower and upper gates are dynamic with transient opening/closing. Through ensemble‐based structural analyses of the closed state versus the open state, we revealed pronounced closed‐to‐open conformational changes involving the membrane proximal domain (MPD) linker, the outer pore, and the TRP helix, which are accompanied by breaking/forming of a network of closed/open‐state specific hydrogen bonds. By comparing the closed‐state simulations at 30°C and 60°C, we observed heat‐activated conformational changes in the MPD linker, the outer pore, and the TRP helix that resemble the closed‐to‐open conformational changes, along with partial formation of the open‐state specific hydrogen bonds. Some of the residues involved in the above key hydrogen bonds were validated by previous mutational studies. Taken together, our MD simulations have offered rich structural and dynamic details beyond the static structures of TRPV1, and promising targets for future mutagenesis and functional studies of the TRPV1 channel. Proteins 2016; 84:1938–1949. © 2016 Wiley Periodicals, Inc.     

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.     

Structure and function of the bacterial mechanosensitive channel of large conductance.

Mechanosensation in bacteria involves transducing membrane stress into an electrochemical response. In Escherichia coli and other bacteria, this function is carried out by a number of proteins including MscL, the mechanosensitive channel of large conductance. MscL is the best characterized of all mechanosensitive channels. It has been the subject of numerous structural and functional investigations. The explosion in experimental data on MscL recently culminated in the solution of the three-dimensional structure of the MscL homologue from Mycobacterium tuberculosis. In this review, much of these data are united and interpreted in terms of the newly published M. tuberculosis MscL crystal structure.     

The "dashpot" mechanism of stretch-dependent gating in MscS

The crystal structure of the small conductance mechanosensitive channel (MscS) has been an invaluable tool in the search for the gating mechanism, however many functional aspects of the channel remain unsettled. Here we characterized the gating of MscS in Escherichia coli spheroplasts in a triple mutant (mscL-, mscS-, mscK-) background. We used a pressure clamp apparatus along with software developed in-lab to generate dose-response curves directly from two-channel recordings of current and pressure. In contrast to previous publications, we found that MscS exhibits essentially voltage-independent activation by tension, but at the same time strong voltage-dependent inactivation under depolarizing conditions. The MscS activation curves obtained under saturating ramps of pressure, at different voltages, gave estimates for the energy, area, and gating charge for the closed-to-open transition as 24 kT, 18 nm2, and +0.8, respectively. The character of activation and inactivation was similar in both K+ and Na+ buffers. Perhaps the most salient and intriguing property of MscS gating was a strong dependence on the rate of pressure application. Patches subjected to various pressure ramps from 2.7 to 240 mmHg/s revealed a midpoint of activation almost independent of rate. However, the resultant channel activity was dramatically lower when pressure was applied slowly, especially at depolarizing pipette voltages. It appears that MscS prefers to respond in full to abrupt stimuli but manages to ignore those applied slowly, as if the gate were connected to the tension-transmitting element via a velocity-sensitive "dashpot." With slower ramps, channels inactivate during the passage through a narrow region of pressures below the activation midpoint. This property of "dumping" a slowly applied force may be important in environmental situations where rehydration of cells occurs gradually and release of osmolytes is not desirable. MscS often enters the inactivated state through subconducting states favored by depolarizing voltage. The inactivation rate increases exponentially with depolarization. Based on these results we propose a kinetic scheme and gating mechanism to account for the observed phenomenology in the framework of available structural information.     

Fundamental gating mechanism of nicotinic receptor channel revealed by mutation causing a congenital myasthenic syndrome

We describe the genetic and kinetic defects in a congenital myasthenic syndrome due to the mutation epsilonA411P in the amphipathic helix of the acetylcholine receptor (AChR) epsilon subunit. Myasthenic patients from three unrelated families are either homozygous for epsilonA411P or are heterozygous and harbor a null mutation in the second epsilon allele, indicating that epsilonA411P is recessive. We expressed human AChRs containing wild-type or A411P epsilon subunits in 293HEK cells, recorded single channel currents at high bandwidth, and determined microscopic rate constants for individual channels using hidden Markov modeling. For individual wild-type and mutant channels, each rate constant distributes as a Gaussian function, but the spread in the distributions for channel opening and closing rate constants is greatly expanded by epsilonA411P. Prolines engineered into positions flanking residue 411 of the epsilon subunit greatly increase the range of activation kinetics similar to epsilonA411P, whereas prolines engineered into positions equivalent to epsilonA411 in beta and delta subunits are without effect. Thus, the amphipathic helix of the epsilon subunit stabilizes the channel, minimizing the number and range of kinetic modes accessible to individual AChRs. The findings suggest that analogous stabilizing structures are present in other ion channels, and possibly allosteric proteins in general, and that they evolved to maintain uniformity of activation episodes. The findings further suggest that the fundamental gating mechanism of the AChR channel can be explained by a corrugated energy landscape superimposed on a steeply sloped energy well.     

Ion channel gating: insights via molecular simulations

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

Conversion of a mechanosensitive channel protein from a membrane-embedded to a water-soluble form by covalent modification with amphiphiles

Covalent modification of integral membrane proteins with amphiphiles may provide a general approach to the conversion of membrane proteins into water-soluble forms for biophysical and high-resolution structural studies. To test this approach, we mutated four surface residues of the pentameric Mycobacterium tuberculosis mechanosensitive channel of large conductance (MscL) to cysteine residues as anchors for amphiphile attachment. A series of modified ion channels with four amphiphile groups attached per channel subunit was prepared. One construct showed the highest water solubility to a concentration of up to 4mg/ml in the absence of detergent. This analog also formed native-like, alpha-helical homo-pentamers in the absence of detergent as judged by circular dichroism spectroscopy, size-exclusion chromatography and various light-scattering techniques. Proteins with longer, or shorter polymers attached, or proteins modified exclusively with polar cysteine-reactive small molecules, exhibited reduced to no solubility and higher-order aggregation. Electron microscopy revealed a homogeneous population of particles consistent with a pentameric channel. Solubilization of membrane proteins by covalent attachment of amphiphiles results in homogeneous particles that may prove useful for crystallization, solution NMR spectroscopy, and electron microscopy.     

Essential function of the N‐termini tails of the proteasome for the gating mechanism revealed by molecular dynamics simulations

Proteasome is involved in the degradation of proteins. Proteasome activators bind to the proteasome core particle (CP) and facilitate opening a gate of the CP, where Tyr8 and Asp9 in the N‐termini tails of the CP form the ordered open gate. In a double mutant (Tyr8Gly/Asp9Gly), the N‐termini tails are disordered and the stabilized open‐gate conformation cannot be formed. To understand the gating mechanism of the CP for the translocation of the substrate, four different molecular dynamics simulations were carried out: ordered‐ and Tyr8Gly/Asp9Gly disordered‐gate models of the CP complexed with an ATP‐independent PA26 and ordered‐ and disordered‐gate models of the CP complexed with an ATP‐dependent PAN‐like activator. The free‐energies of the translocation of a polypeptide substrate moving through the gate were estimated. In the ordered‐gate models, the substrate in the activator was more stable than that in the CP. The conformational entropy of the N‐termini tails of the CP was larger when the substrate was in the activator than in the CP. In the disordered‐gate models, the substrate in the activator was more destabilized than in the ordered‐gate models. The mutated N‐termini tails became randomized and their increased conformational entropy could no longer increase further even when the substrate was in the activator, meaning the randomized N‐termini tails had lost the ability to stabilize the substrate in the activator. Thus, it was concluded that the dynamics of the N‐termini tails entropically play a key role in the translocation of the substrate. Proteins 2014; 82:1985–1999. © 2014 Wiley Periodicals, Inc.     

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.     

Investigating lipid composition effects on the mechanosensitive channel of large conductance (MscL) using molecular dynamics simulations

Previous experimental work has shown that the functional properties of the mechanosensitive channel of large conductance (MscL) are affected by variations in lipid composition. Here, we utilize molecular dynamics simulations of Mycobacterium tuberculosis MscL to investigate such lipid composition effects on a molecular level. In particular, two sets of simulations were performed. In the first, trajectories using lipids with different headgroups (phosphatidylcholine and phosphatidylethanolamine) were compared. Protein-lipid interactions were clearly altered by the headgroup changes, leading to conformational differences in the C-terminal region of M. tuberculosis MscL. In the second set of simulations, lipid tails were gradually shortened, thinning the membrane over a molecular dynamics trajectory. These simulations showed evidence of hydrophobic matching between MscL and the lipid membrane, as previously proposed. For all simulations, protein-lipid interaction energies in the second transmembrane region were correlated to mutagenic data, emphasizing the importance of lipid interactions for proper MscL function.     

Homology modeling and molecular dynamics simulations of the alpha1 glycine receptor reveals different states of the channel

Homology modeling is used to build initial models of the transmembrane domain of the human alpha1 glycine receptor (GlyR) based on the most recently published refined structure of nAChR (PDB ID: 2BG9). Six preliminary GlyR models are constructed using two different approaches. In one approach, five different homopentamers are built by symmetric assembly of alpha1 GlyR subunits using only one of the five unique chains of nAChR as a template. In a second approach, each nAChR subunit serves as a template for an alpha1 GlyR subunit. All six initial GlyR constructs are then embedded into a hydrated POPC lipid bilayer and subjected to molecular dynamics simulation for at least six nanoseconds. Each model is stable throughout the simulation, and the final models fall into three distinct categories. Homopentameric GlyR bundles using a single alpha nAChR subunit as a template appear to be in an open conformation. Under an applied external potential, permeation of Cl(-) ions is observed within several ns in a channel built on an alpha chain. Model channels built on non-alpha chains have a constriction either near the intracellular mouth or more centrally located in the pore domain, both of which may be narrow enough to close the channel and whose locations correspond to putative gates observed in nicotinicoid receptors. The differences between these three general models suggest that channel closure may be effected by either rotation or tangential tilting of TM2.     

Computer simulation of the KvAP voltage-gated potassium channel: steered molecular dynamics of the voltage sensor

The recent crystal structures of the voltage-gated potassium channel KvAP and its isolated voltage-sensing 'paddle' (composed of segments S1-S4) challenge existing models of voltage gating and raise a number of questions about the structure of the physiologically relevant state. We investigate a possible gating mechanism based on the crystal structures in a 10 ns steered molecular dynamics simulation of KvAP in a membrane-mimetic octane layer. The structure of the full KvAP protein has been modified by restraining the S2-S4 domain to the conformation of the isolated high-resolution paddle structure. After an initial relaxation, the paddle tips are pulled through the membrane from the intracellular to the extracellular side, corresponding to a putative change from closed to open. We describe the effect of this large-scale motion on the central pore domain, which remains largely unchanged, on the protein hydrogen-bonding network and on solvent. We analyze the motion of the S3b-S4 portion of the protein and propose a possible coupling mechanism between the paddle motion and the opening of the channel. Interactions between the arginine residues in S4, solvent and chloride ions are likely to play a role in the gating charge.     

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.     

Expression of Arabidopsis MCA1 enhanced mechanosensitive channel activity in the Xenopus laevis oocyte plasma membrane

Higher plants sense and respond to osmotic and mechanical stresses such as turgor, touch, flexure and gravity. Mechanosensitive (MS) channels, directly activated by tension in the cell membrane and cytoskeleton, are supposed to be involved in the cell volume regulation under hypotonic conditions and the sensing of these mechanical stresses based on electrophysiological and pharmacological studies. However, limited progress has been achieved in the molecular identification of plant MS channels. Here, we show that MCA1 (mid1-complementing activity 1; a putative mechanosensitive Ca²⁺-permeable channel in Arabidopsis thaliana) increased MS channel activity in the plasma membrane of Xenopus laevis oocytes. The functional and kinetic properties of MCA1 were examined by using a Xenopus laevis oocytes expression system, which showed that MCA1-dependent MS cation currents were activated by hypo-osmotic shock or by membrane stretch produced by pipette suction. Single-channel analyses suggest that MCA1 encodes a possible MS channel with a conductance of 34 pS.