Structural models of the MscL gating mechanism.

Three-dimensional structural models of the mechanosensitive channel of large conductance, MscL, from the bacteria Mycobacterium tuberculosis and Escherichia coli were developed for closed, intermediate, and open conformations. The modeling began with the crystal structure of M. tuberculosis MscL, a homopentamer with two transmembrane alpha-helices, M1 and M2, per subunit. The first 12 N-terminal residues, not resolved in the crystal structure, were modeled as an amphipathic alpha-helix, called S1. A bundle of five parallel S1 helices are postulated to form a cytoplasmic gate. As membrane tension induces expansion, the tilts of M1 and M2 are postulated to increase as they move away from the axis of the pore. Substantial expansion is postulated to occur before the increased stress in the S1 to M1 linkers pulls the S1 bundle apart. During the opening transition, the S1 helices and C-terminus amphipathic alpha-helices, S3, are postulated to dock parallel to the membrane surface on the perimeter of the complex. The proposed gating mechanism reveals critical spatial relationships between the expandable transmembrane barrel formed by M1 and M2, the gate formed by S1 helices, and "strings" that link S1s to M1s. These models are consistent with numerous experimental results and modeling criteria.     

Multimerization of the ligand binding domains of cyclic nucleotide-gated channels

Cyclic nucleotide-gated (CNG) channels comprise four subunits and are activated by the direct binding of cyclic nucleotide to an intracellular domain on each subunit. This ligand binding domain is thought to contain a beta roll followed by two alpha helices, designated the B and C helices. To examine the quaternary structure of CNG channels and how it changes during ion channel gating, we introduced single cysteines along the C helix of each subunit in an otherwise cysteineless channel. We found that cysteines on the C helices could form intersubunit disulfide bonds, even between diagonal subunits. Disulfide bond formation occurred primarily in closed channels and inhibited channel opening. These data suggest that the C helices from all four channel subunits are in close proximity in the closed state and move apart during channel opening.     

Kir6.2 Channel Gating by Intracellular Protons: Subunit Stoichiometry for Ligand Binding and Channel Gating

The adenosine triphosphate-sensitive K⁺ (K_ATP) channels are gated by several metabolites, whereas the gating mechanism remains unclear. Kir6.2, a pore-forming subunit of the K_ATP channels, has all machineries for ligand binding and channel gating. In Kir6.2, His175 is the protonation site and Thr71 and Cys166 are involved in channel gating. Here, we show how individual subunits act in proton binding and channel gating by selectively disrupting functional subunits using these residues. All homomeric dimers and tetramers showed pH sensitivity similar to the monomeric channels. Concatenated construction of wild type with disrupted subunits revealed that none of these residues had a dominant-negative effect on the proton-dependent channel gating. Subunit action in proton binding was almost identical to that for channel gating involving Cys166, suggesting a one-to-one coupling from the C terminus to the M2 helix. This was significantly different from the effect of T71Y heteromultimers, suggesting distinct contributions of M1 and M2 helices to channel gating. Subunits underwent concerted rather than independent action. Two wild-type subunits appeared to act as a functional dimer in both cis and trans configurations. The understanding of K_ATP channel gating by intracellular pH has a profound impact on cellular responses to metabolic stress as a significant drop in intracellular pH is more frequently seen under a number of physiological and pathophysiological conditions than a sole decrease in intracellular ATP levels. Runping Wang, Junda Su contributed equally to this work.     

Gating of MscL studied by steered molecular dynamics

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

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.     

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

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

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.     

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

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

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

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

An NMDA Receptor Gating Mechanism Developed from MD Simulations Reveals Molecular Details Underlying Subunit-Specific Contributions

N-methyl-D-aspartate (NMDA) receptors are obligate heterotetrameric ligand-gated ion channels that play critical roles in learning and memory. Here, using targeted molecular dynamics simulations, we developed an atomistic model for the gating of the GluN1/GluN2A NMDA receptor. Upon agonist binding, lobe closure of the ligand-binding domain produced outward pulling of the M3-D2 linkers, leading to outward movements of the C-termini of the pore-lining M3 helices and opening of the channel. The GluN2A subunits, similar to the distal (B/D) subunits in the homotetrameric GluA2 α-amino-3-hydroxy-5-methyl-4-isoxazoleproprionate receptor, had greater M3 outward movements and thus contributed more to channel gating than the GluN1 subunits. Our gating model is validated by functional studies, including cysteine modification data indicating wider accessibility to the GluN1 M3 helices than to the GluN2A M3 helices from the lumen of the open channel, and reveals why the Lurcher mutation in GluN1 has a stronger ability in maintaining channel opening than the counterpart in GluN2A. The resulting structural model for the open state provides an explanation for the Ca²⁺ permeability of NMDA receptors, and the structural differences between the closed and open states form the basis for drug design.     

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

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

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.     

Capping transmembrane helices of MscL with aromatic residues changes channel response to membrane stretch

Tyrosines and tryptophans that anchor both ends of the helices to membrane interfaces in many transmembrane proteins are not common in MscL and homologous mechanosensitive channels. This characteristic absence of two aromatic "belts" may be critical for MscL function as the opening transition is predicted to be associated with a strong helical reorientation. A single tyrosine (Y75) on the extracellular side of the M2 helix of pentameric EcoMscL is absent in TbMscL, which instead has a single tyrosine (Y87) on the cytoplasmic side of M2. Moving the tyrosine of EcoMscL to the intracellular side (Y75F/F93Y) or capping the TM2 helix on both sides (F93Y/W) slows the kinetics of gating and increases the threshold for activation, leading to a partial loss-of-function in osmotic shock survival assays. Increasing the distance between the caps (L98W, L102Y/W) partially restores channel function presumably by loosening restraints for tilting. Capping the TM2 helix with a charged residue (Y75E) causes a right shift of the activation curve ("stiff" phenotype) and abolishes function. Introducing a "cap" into the TM1 helix (I41W) decreases the activation threshold and shortens the mean open time but unexpectedly leads to a complete loss-of-function in vivo. The data are consistent with the view that restraining helical positions in MscL by introducing specific protein-lipid interactions at membrane interfaces compromises MscL function. Subtle differences in osmotic shock survival are more evident at low levels of mutant protein expression. We observed a correlation between the right shift of tension activation threshold and the loss-of-function channel phenotype, with a few exceptions that point to other parameters of gating that may define the osmotic rescuing ability in vivo.     

Gating-induced Conformational Rearrangement of the γ-Aminobutyric Acid Type A Receptor β-α Subunit Interface in the Membrane-spanning Domain

GABA(A) receptors mediate fast inhibitory synaptic transmission. The transmembrane ion channel is lined by a ring of five α helices, M2 segments, one from each subunit. An outer ring of helices comprising the alternating M1, M3, and M4 segments from each subunit surrounds the inner ring and forms the interface with the lipid bilayer. The structural rearrangements that follow agonist binding and culminate in opening of the ion pore remain incompletely characterized. Propofol and other intravenous general anesthetics bind at the βM3-αM1 subunit interface. We sought to determine whether this region undergoes conformational changes during GABA activation. We measured the reaction rate of p-chloromercuribenzenesulfonate (pCMBS) with cysteines substituted in the GABA(A) receptor α1M1 and β2M3 segments. In the presence of GABA, the pCMBS reaction rate increased significantly in a cluster of residues in the extracellular third of the α1M1 segment facing the β2M3 segment. Mutation of the β2M2 segment 19' position, R269Q, altered the pCMBS reaction rate with several α1M1 Cys, some only in the resting state and others only in the GABA-activated state. Thus, β2R269 is charged in both states. GABA activation induced disulfide bond formation between β2R269C and α1I228C. The experiments demonstrate that α1M1 moves in relationship to β2M2R269 during gating. Thus, channel gating does not involve rigid body movements of the entire transmembrane domain. Channel gating causes changes in the relative position of transmembrane segments both within a single subunit and relative to the neighboring subunits.     

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

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

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

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

Critical protein domains and amino acid residues for gating the KIR6.2 channel by intracellular ATP

K(ATP) channels couple intermediary metabolism to cellular excitability. Such a property relies on the inherent ATP-sensing mechanism known to be located in the Kir6 subunit. However, the molecular basis for the ATP sensitivity remains unclear. Here we showed evidence for protein domains and amino acid residues essential for the channel gating by intracellular ATP. Chimerical channels were constructed using protein domains of Kir6.2 and Kir1.1, expressed in HEK293 cells, and studied in inside-out patches. The N and C termini, although important, were inadequate for channel gating by intracellular ATP. Full ATP sensitivity also required M1 and M2 helices. Cytosolic portions of the M1 and M2 sequences were crucial, in which six amino acid residues were identified, i.e., Thr76, Met77, Ala161, Iso162, Leu164, and Cys166. Site-specific mutation of any of them reduced the ATP sensitivity. Construction of these residues together with the N/C termini produced ATP sensitivity identical to the wild-type channels. The requirement for specific membrane helices suggests that the Kir6.2 gating by ATP is not shared by even two closest relatives in the K(+) channel family, although the general gating mechanisms involving membrane helices appear to be conserved in all K(+) channels.     

Determinant Role of Membrane Helices in K<Subscript>ATP</Subscript> Channel Gating

The ATP-sensitive K⁺ (K_ATP) channels couple chemical signals to cellular activity, in which the control of channel opening and closure (i.e., channel gating) is crucial. Transmembrane helices play an important role in channel gating. Here we report that the gating of Kir6.2, the core subunit of pancreatic and cardiac K_ATP channels, can be switched by manipulating the interaction between two residues located in transmembrane domains (TM) 1 and 2 of the channel protein. The Kir6.2 channel is gated by ATP and proton, which inhibit and activate the channel, respectively. The channel gating involves two residues, namely, Thr71 and Cys166, located at the interface of the TM1 and TM2. Creation of electrostatic attraction between these sites reverses the channel gating, which makes the ATP an activator and proton an inhibitor of the channel. Electrostatic repulsion with two acidic residues retains or even enhances the wild-type channel gating. A similar switch of the pH-dependent channel gating was observed in the Kir2.1 channel, which is normally pH- insensitive. Thus, the manner in which the TM1 and TM2 helices interact appears to determine whether the channels are open or closed following ligand binding.^*These authors contributed equally to this work.     

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.