From the * Department of Biology, University of Maryland, College Park, Maryland 20742; Department of Physiology and Biophysical
Sciences, State University of New York, Buffalo, New York 14214; and § Laboratory of Molecular Biology and Department of Genetics,
University of Wisconsin, Madison, Wisconsin 53706
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
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 E = 18.6 kBT. If the mechanosensitivity arises from tension acting on a change of in-plane area (
A), the free energy, T
A, would correspond to
A = 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
A. 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
A.
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.
![]() |
INTRODUCTION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Mechanosensitive (MS)1 channels or stretch-sensitive
channels, discovered in chick skeletal muscle cells (Guharay and Sachs, 1984; Brehm et al., 1984
) are likely
candidates for the role of primary mechanoreceptors
in unicellular and multicellular organisms. In the hair
cells responsible for hearing and balance in vertebrates, MS channels are implicated in generation of
primary potentials associated with the hair bundle displacement (Corey and Hudspeth, 1983
). In nonsensory cells, MS channels have been shown to mediate mechanical stress-induced changes in membrane permeability to monovalent ions and Ca2+ (Sachs and Morris, 1998
), potentially triggering cascades of second
messenger signaling. In bacteria that live in a rapidly changing environment, mechanosensitive channels
mediate permeation of small osmolytes from the cytoplasm through the periplasm to the extracellular space,
potentially permitting a rapid regulation of turgor pressure (Saimi et al., 1993
; Berrier et al., 1996
; Blount et al.,
1998
).
MS channels are identified by their activities, not sequence similarities. Given the variety of observed ionic
selectivities and gating properties, there is little evidence that they form a family by sequence similarity
(Martinac, 1993). In eucaryotes, MS channel activity
appears to require coupling to the cytoskeleton and/or
the extracellular matrix, probably for the transmission of force (Sachs and Morris, 1998
). In contrast, bacterial
MscL is functional in lipid bilayers (Sukharev et al.,
1994b
).
Although some MS channels have been characterized biophysically, very little is known about their molecular structure. An extensive genetic dissection of the
nematode Caenorhabditis elegans revealed ~20 genes involved in touch sensation (Driscoll and Chalfie, 1993;
Tavernakis and Driscoll, 1998
). Two of them are possible candidates for MS channel subunits, featuring close
sequence similarity to the ENaC, the amiloride-sensitive Na+ channel (Rossier et al., 1994
). Other touch
genes code for cytoskeletal or extracellular matrix components, presumably those elements necessary for the
efficient transmission of forces to the channels. Unfortunately, there are only two descriptions of electrophysiological evidence for the mechanosensitive role of any
of these genes (Awayda et al., 1995
; Kizer et al., 1997
).
As an alternative to the complexity of metazoa, MS
channels in microorganisms are biochemically and genetically accessible (Martinac et al., 1992; Martinac,
1993
; Sukharev et al., 1993
, 1994a
, 1997
; Le Dain et al.,
1998
). Patch-clamp studies on native and reconstituted
Escherichia coli membranes revealed three types of
mechanosensitive channel (Msc) activities: MscL, MscS,
and MscM (L, S, and M stand for large, small, and
"mini" conductances, respectively; Berrier et al., 1996
).
Using a variety of chromatographic techniques followed by reconstitution and patch recording of channels in liposomes, MscL, the most conductive of MS
channels in E. coli was identified as an ~17-kD protein,
and the corresponding mscL gene was then cloned
(Sukharev et al., 1994a
).
Biochemical studies have shown that MscL resides in
the inner membrane of E. coli (Blount et al., 1996a).
Each MscL subunit is a 15-kD protein with two putative
transmembrane domains and a high
-helical content
(Arkin et al., 1998
). Whereas the functional channel
complex was proposed to be a homohexamer (Blount
et al., 1996b
), a recent crystallographic study indicates
a pentameric structure for the closed state of the channel (Chang et al., 1998
) MscL can be activated by 70-
180-mmHg pressure gradients across patches of bacterial spheroplasts, and purified MscL channels reconstituted into phospholipid liposomes produced similar
currents. This indicates that MscL can be gated directly
by tension transmitted via the lipid bilayer alone. In
vivo, the channel may be opened by osmotic gradients
of 200-700 mOsm (Blount et al., 1998
). MscL-like channels are found in several groups of eubacteria
(Sukharev et al., 1997
; Moe et al., 1998
), and there is
increasing evidence that MscL plays the role of a "safety
valve" in prokaryotes, releasing small osmolytes, there-
by reducing the turgor pressure and the chance of cell lysis.
In the present work, we quantitatively evaluate the
energetic parameters for MscL gating in reconstituted
liposomes using a similar approach as described by Opsahl and Webb (1994a). We used high-resolution video
microscopy to measure the curvature of patches at different activating pressures, permitting us to calculate
the absolute tension. Kinetic and thermodynamic analysis of the channel allows us to calculate the free energy
differences between states and their tension dependence (Hamill and McBride, 1994
). Close examination
of the single channel conductance has identified at
least four open conductance classes. Analysis of transition rate constants, using a simple linear Markov
model, has permitted identification of the tension dependence of each rate constant between the closed and
open conductance classes. This work represents the
first calibration of a cloned, biological mechanosensitive ion channel.
![]() |
METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
MscL Isolation and Reconstitution
The procedure for MscL purification using a 6His-tag has been
described previously (Blount et al., 1996b; Sukharev et al., 1996
).
Briefly, a tag of six sequential histidines was added to the COOH
terminus of MscL by a two-step PCR amplification and the extended gene was cloned into the pB10a vector (Sukharev et al.,
1994a
). The PB104 cells expressing MscL-6His were French-pressed and the total membrane fraction was isolated. The MscL-6His protein was extracted from membranes by solubilization in
3%
-octylglucoside and purified in one step using a Ni-NTA column (QIAGEN Inc.) as described (Sukharev et al., 1996
). Azolectin (Soybean lecithin, type II; Sigma Chemical Co.), a lipid component of proteoliposomes, was partially purified from oxidized
and lyso forms using chloroform/water separation. Briefly, 500 mg
of azolectin beads were washed five times with 5-7 ml of acetone
(electron microscopy grade), with gentle swirling after each
change. After the last wash, the residual acetone was removed by
vacuum (20 min) and the beads were dissolved in 5 ml of chloroform. The solution was placed in a thick-walled glass tube with a
Teflon stopper, overlaid with 5 ml of water, and shaken vigorously. The cream-colored mixture was separated by 2 h centrifugation in a swinging-bucket rotor at 13,000 rpm at 15°C. The
lower chloroform portion was retrieved by a long-needle glass syringe, placed in an airtight Teflon-capped vial, and could be
stored at
20°C under nitrogen for ~1 mo.
MscL-6His was reconstituted into azolectin liposomes by dialyzing the -octylglucoside-solubilized mixture of a protein-to-lipid ratio of 1:500 to 1:2,000 (Sukharev et al., 1994b
). In special cases
when we needed to record single MscL currents for kinetic analysis, the protein-to-lipid ratio was lowered to ~1:10,000. Proteoliposomes were subjected to a dehydration-rehydration cycle on
glass slides (Sukharev et al., 1994b
) and the resultant multilayer aggregates were placed in the recording buffer (see below) for 30-60 min before the patch-clamp experiment. Large, and apparently unilamellar, blisters formed on the surface of multilayer
liposomes were examined as excised inside-out patches.
Channel Recording
Borosilicate glass pipettes with 1-2-µm bore diameter were used
in all experiments. The pipette pulling protocol was adjusted to
form long-tipped pipettes with an almost cylindrical, 10-20-µm- thick by 100-µm-long region before the tip. This narrow taper was useful for observation of liposome patches that tend to creep up the pipette under pressure gradients. To compensate for the headstage tilt and make the observed part of the pipette nearly parallel to the focal plane, pipette tips were bent as described previously (Sokabe et al., 1991). All recordings were performed in a
symmetrical buffer containing 200 M KCl, 40 mM MgCl2, and 10 mM HEPES, pH 7.2. Pressure gradients were delivered by either a
pneumatic screw-driven syringe and monitored by an electronic pressure transducer or with a laboratory-built hydraulic pressure servo (Sachs, 1987
; Sokabe et al., 1991
). Patch currents were recorded at hyperpolarizing voltages (+20 mV in the pipette) using an Axopatch 200 (Axon Instruments) and stored on a PCM tape
(48 kHz sampling rate; Instrutech Corp.). The second data channel on the recorder was allocated to recording pressure.
To determine Po as a function of pressure, the data were analyzed using PCLAMP6. Po was calculated as the mean patch conductance GP divided by n · GMscL, where GMscL is the maximal conductance of a single MscL channel (typically 3.5-3.7 nS) and n is
the number of channels in the patch. n was estimated either by
measuring the current at saturating pressures or, with less precision, by using the average number of channels per patch in a
given proteoliposome preparation determined in separate experiments. When n was >10, the patch conductance, GP, was corrected to account for the series pipette resistance (RS 1.5-2 M
)
by the equation, GP = I/(V
IRS), where V and I are the transmembrane voltage and current, respectively. We made >200
attempts to record complete activation curves for MscL in different settings and, of these, fifteen were considered extensive
enough to warrant analysis. They compose the data presented below. Experiments were usually terminated by lysis of the patch.
Kinetic Analysis
For the multistate analysis, we used the QuB program suite
(). To determine the rate constants between states, the digitized data was first idealized using SKM, a Hidden Markov algorithm (Feng et al., 1996; Qin et al., 1996
). The
"events list" outputs from SKM were then grouped as a global collection of data sets over defined tension. This tension series was
fit to a kinetic model using MIL, a maximum likelihood interval
analysis program that permits data to be fit across independent
variables and corrects for missed events (Qin et al., 1995
). Since
we determined there were five states (four subconducting and
one shut state), the number of possible kinetic connections was
extremely large (728 models). We were able to exhaustively
search all models and connectivities using the program
MSEARCH, which employs MIL to rank all models, based on the
likelihood calculated for each model. For all nonlooping models,
the linear sequential model gave the greatest likelihood, and
therefore was used for the analysis: C1-S2-S3-S4-O5, where C,
S, and O refer to the closed, subconductance, and open states, respectively.
To determine how the single channel MscL conductance depends on the bath conductivity, we recorded currents from reconstituted channels in baths containing 10 mM HEPES, pH 7.2, and 40 mM MgCl2, with KCl concentrations varied between 0.1 and 2 M. The specific conductivity of each buffer was measured directly using a YSI 34 Conductance Meter equipped with a dip-type glass cell with platinized platinum-iridium electrodes (1.0 cm cell constant; Yellow Springs Instrument Co.). Single-channel currents were recorded at +20 mV (pipette voltage), and the unitary current of the fully open state was determined using FETCHAN's `Measure' function (PCLAMP suite; Axon Instruments).
Patch Imaging and Determination of Curvature
Patches were imaged using an inverted microscope (Axiovert;
Carl Zeiss, Inc.) equipped for Differential Interference Contrast with a CCD camera as initially described by Sokabe et al. (1991). The digitized images were analyzed with an algorithm written in JAVA by Akinlaja (1997)
that solved for the radius of curvature.
To avoid confusion within equations, T will reference tension, while temperature T will only appear as a term with the Boltzmann constant; i.e., kBT.
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Previous patch-clamp experiments revealed that activation of MscL requires pressures close to the lytic tension of the unmodified bilayer (Martinac, 1993; Sachs
and Morris, 1998
). This proximity made it difficult to
obtain many data sets with saturating responses. Typically, activation in spheroplast patches required 70-200
mmHg, while liposome patches needed 40-150 mmHg.
Patches formed with large-diameter pipettes activated
at lower pressures, suggesting that the actual parameter
that drives MscL gating is tension, not pressure. Prolonged exposure of patches to high pressures often
caused lysis. This imposed strict requirements on the
stability of the bilayer in which the channels were reconstituted. A partially purified crude phospholipid
fraction from soybean (azolectin) gave us acceptable
stability and reproducibility of results. Under a small
pressure gradient, liposome patches, lacking cellular components, are spherical caps and are large enough
to be viewed by conventional optical microscopy. Experiments were considered successful if we were able to
obtain enough points to fit the partial activation curve
with Po > 0.3 with images of the patch of sufficient quality to calculate the curvature.
Analysis of the Data as a Two-State System
As the data below will show, MscL is not a binary channel, but has multiple conducting states. However, some of the essential features of the gating process can be gained from considering the simpler two-state model, a model that would correspond to the data viewed at low bandwidth. This analysis is based on setting a threshold for being open at half the fully open state amplitude. Previous data indicated that mean current MscL pressure-activation curves can be well fit with single component Boltzmann distribution. Our initial experiments were designed to determine the midpoint (T1/2) and the maximum slope of this distribution, the two key parameters of the dose-response curve. Only one parameter, T1/2, actually requires the measurement of patch curvature. An accurate measure of the slope can be obtained from Po(P) curves measured "blindly" (i.e., without geometric measurements), provided that the radius of curvature is independent of the pressure (see below).
Fig. 1 shows a typical trace from a patch containing ~100 MscL channels in response to a stepwise increase of pressure gradient P (bottom). The MscL current activates at ~40 mmHg and increases with P in a nonlinear manner. Note that in the beginning of the trace the variance of the current (amplitude of fluctuations around the mean level) increases with each step, reaches a maximum, and finally decreases during the last pressure step before the patch ruptures. The maximum of the current variance indicates the point at which Po = 0.5; thus the half-maximal pressure P1/2 is between 48 and 52 mmHg for this particular recording (Fig. 1, D). The Rs-corrected dose-response curve derived from this trace is shown in Fig. 1, inset.
|
|
Fig. 2 shows the geometry of the same patch at pressure gradients of 0, 20, 44, and 52 mmHg. In the absence of suction, the patch is essentially flat, subjected
to the "resting" tension arising from membrane adhesion to the glass surface (Fig. 2 A; Opsahl and Webb,
1994b). At low suction, the patch appears as a spherical
cap, with a progressively increasing curvature (Fig. 2
B). At higher suction, the curvature saturates because
the membrane is nearly inextensible (see calculation
below). At higher suctions, the patch may creep up the
pipette (Fig. 2 D). This creep does not affect the geometry significantly as long as the pipette taper is small. As
seen in Fig. 1, MscL channels are only active for P > 40 mmHg. This is the range of pressures where the patch
curvature reaches saturation so that one measurement
of the curvature was adequate to calculate the tension
for all pressures where MscL is active. The procedure of
fitting the patch curvature is depicted in Fig. 3.
|
We have been able to obtain partial activation curves from four independent patches, in which we were also able to measure the radius of curvature, r. The membrane tension was calculated for every pressure according to Laplace's law, T = p × r/2, and the Po data were plotted against T as shown in Fig. 4. The patches were fit as a group to a single Boltzmann function with T1/2 = 11.8 ± 0.8 dyn/cm and a maximal slope of 0.61 ± 0.17 dyn/cm per e-fold change of Po/Pc.
|
More precise and statistically reliable estimates for the slope of Po(T) curves were obtained from the analysis of additional Po(P) experiments done without measurement of the patch curvature, but aligned to T1/2 = 11.81 dyn/cm. Given that the patch curvature does not change significantly in the range of pressures where MscL is active (compare Fig. 3), the two scales, T and P, are equivalent with T1/2 = r × P1/2/2. Measuring Po(P) curves alone is much easier than the simultaneous geometric measurement. 11 complete Po(P) curves obtained from three independent MscL preparations were fit independently with a Boltzmann curve. P1/2 was determined individually for each curve, and then rescaled to the same midpoint, T1/2 = 11.81 dyn/cm. As shown in Fig. 5, the curves display a remarkable consistency and, when fitted together, indicate the slope factor of 0.63 ± 0.08 dyn/cm.
|
If interpreted in the framework of a two-state Boltzmann model with the change of in-plane area being
the dominant energy term, TA, we have:
![]() |
(1) |
Fitting to Eq. 1 gave E = 18.61 kBT (46.3 kJ/mol) for
the free energy of the closed-to-open transition in the
absence of stress, and
A = 6.52 nm2 for the in-plane
change in area between closed and open (kBT is Boltzmann's constant × the absolute temperature = 4.04 × 10
14 erg at room temperature).
In a more specific representation of the two-state model, we can write the forward and backward rates as:
![]() |
(2) |
![]() |
(3) |
where k0 is a scaling factor incorporating T1/2 and a and
b are functions of the energies between the energy wells
and the top of the barrier separating the states (Fig. 6).
For example, a = Ac/kBT, where
Ac is the increase of
in-plane area between the closed state and the top of
the energy barrier. The units of a are inverse tension
(cm2/erg = cm/dyne). Since the probability of being
open, Po = k12/(k12 + k21), the Boltzmann equation assumes the form Po/Pc = exp[(a
b) T].
|
Fig. 7 represents typical single-channel traces illustrating the MscL kinetics at different tensions (A) and
semi-logarithmic tension dependencies of rate constants for opening and closing transitions (kon and koff)
derived from these traces according to a two-state
model (B). The dependence of ln(Po/Pc) for the same
patch illustrates the rate at which the equilibrium occupancies of open and closed states change with tension.
The tension dependence of the rate of opening (kon) is
about equal to that of ln(Po/Pc). The tension dependence of the rate of closing is four times less. This suggests that the shortening of the closed state makes a larger contribution to the change of Po than the lengthening of the open state. From the tension dependence
of the rate constants, we estimate that Ac Ab is
4.42 nm2 (equivalent to a disk of radius 1.19 nm) and Ao
Ab
1.37 nm2, (equivalent to a disk of radius 0.7 nm).
The total area change, based on transition rates between closed and open, is therefore 5.8 nm2, as compared with 6.5 nm2 calculated from Po. As shown below,
the channel is not a two-state system and therefore a
discrepancy may be expected.
|
The Multi-State Analysis
The above analysis was done as though the channel
were a two-state system with rate constants calculated
from dwell times generated according to a half-amplitude detection criterion (Sachs et al., 1982; Colquhoun
and Sigworth, 1983
). Closer examination of the data,
however, shows multiple conducting states (Fig. 8). The
conductances illustrated in Fig. 8 are shown in Table I, along with estimates of the pore radius calculated assuming a cylindrical pore 4.2-nm long (Cruickshank et
al., 1997
).
|
|
While the minimal reaction scheme for so many
states would normally be difficult to determine reliably,
exploration of the likelihood of all possible interconnections (Feng et al., 1996; see METHODS) indicated
that the simple sequential model (shown below) was
the best fit for all nonlooping models (Table II).
|
|
The rate constants for all patches are plotted as a
function of tension and are shown in Fig. 9. The
straight lines are nonlinear regressions to a simple exponential of k = k0 exp(T) with the residuals
weighted by the inverse of the variance provided by MIL. The parameters of the regression are shown in Table III. Most striking is that only k12 has a significant
positive slope; all other rates are insensitive or have
negative slopes. In the simple free energy model presented above, the tension sensitivity
can be identified
with
A/kBT. One needs to bear in mind, however, that
the sign of
will change with the direction of the reaction. A transition from a well to a barrier will produce
the opposite change in area than moving from the barrier to the same well.
|
|
To summarize the kinetics contained in Fig. 9 and Table III, we have plotted in Fig. 10 the energy profiles of the states and the changes of in-plane area. Fig. 10 shows again that the rate limiting step to opening channel is k12, for which the energy barrier is ~38 kBT. At zero tension, the energy difference between the closed and any of the conducting states is >18 kBT, accounting for the fact that the channel is almost never open at rest. All states of conductance >S1 have about the same energy and are insensitive to tension. Applying tension lowers the energy of all the conducting states (energy wells) and barriers equally. The dimensional changes of the channel are shown in Fig. 10, bottom, where we plotted a running sum of the area changes accompanying each transition. As expected from the tension sensitivity, the area changes associated with changing states are small except for the transition from closed to S1. The total change of in-plane area between the closed and open states is 5-6 nm2, as estimated from the two-state analysis.
|
How well does the above analysis fit reasonable physical models of the channel? There are two types of information we can use to estimate the physical structure:
the pore diameter and the channel protein properties.
The channel conductance gives information about the
aqueous pore, and the protein properties give estimates of the wall thickness. We measured the MscL full conductance state as a function of bath conductivity. As
shown in Fig. 11, the fully open single channel conductance was linear with the conductivity up to 2 M KCl.
The absence of saturation and the lack of anion/cation
selectivity is consistent with the representation of the
open MscL as a wide aqueous pore. The conductance data were fit with the Hille equation describing conductance of a cylindrical pore (Hille, 1992), and the results are shown in Table IV. The first column corresponds to the assumed length of the channel. The calculated cross-sectional area suggests that the open
MscL channel has a diameter of 2.7-3.6 nm, which is in good agreement with the data reported recently by
Cruickshank et al. (1997)
. The value for in-plane area
expansion during gating estimated from the two-state
analysis (above),
A = 6.5 nm2, is in the same range estimated for the pore cross-section, 5.8 > A > 9.7 nm2
(see Table IV).
|
|
Having an estimate of the open channel pore diameter, we now need to estimate the wall thickness to obtain the outer diameter where the bilayer tension is applied. The channel is a pentamer, with each monomer
containing two helical transmembrane segments (Blount et al., 1996b; Arkin et al., 1998
; Cheng et al.,
1999). Given that the fully open channel will have a
pore diameter of ~4 nm, there are barely enough
helices to coat an aqueous pore of the estimated diameter
if the helices make a close-contact circle around the
pore. Using 1 nm as an estimate for helix diameter, the
outer diameter of the open channel will be ~6 nm. Alternatively, adding the change in radius of ~0.5 nm
(calculated from the change in area seen by gating) to
the crystal radius of the closed channel (2.5 nm; Cheng
et al., 1999) predicts an outer diameter of open channel of 5.5 nm.
The structure of the lower conductance states, however, is less well constrained since the in-plane area changes associated with gating do not match the changes in conductance. If internal pore reorganization is involved in setting the conductance of the closed and substates, we would predict each substate to exhibit a distinct dependency on the ionic composition of the bath.
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
MscL, the first mechanosensitive channel for which
some structural characteristics are available, now has
quantitative parameters that can be used to judge the
effects of mutagenesis and pharmacology. The ability to
use a fluid bilayer for reconstitution eliminated the
need to deal with the complex mechanics of heterogeneous biological membranes. The low bending moment of bilayers assured the membranes were spherical
so that the pressure gradient could be simply translated into tension according to Laplace's law. Since bilayers
have substantial resistance to area changes, there was
little change in curvature once the membrane assumed
a spherical shape. The image-fitting algorithm used to
determine patch curvature was limited by the optical
resolution of the microscope (~0.35 µm), such that
changes in the radius of curvature could not be reliably
detected at the higher pressures. The expected changes
can be predicted from an equation describing the
stretch of an elastic membrane attached to the pipette
walls (Sokabe et al., 1993):
![]() |
(4) |
where KA is the area elastic constant of the membrane and a is the radius of the pipette where it meets the patch. Using Laplace's law, we can then solve for the radius of curvature r as a function of pressure. Inserting a few constants makes the behavior clear. Consider a soft membrane with KA = 100 dyn/cm and a patch with a = 2 µm. For P = 50 mmHg (1 mmHg = 1.33 × 103 dyn/ cm2), r = 3.33 µm and for 100 mmHg, r = 2.76 µm, a change of 20%. For a stiffer membrane with KA = 500 dyn/cm, the same pressures give r = 5.44 and 4.38 µm, a change of 24%. Thus, in this pressure range, the tension is nearly proportional to the pressure.
The set of data in Fig. 4 shows that the tension to activate MscL is near the lytic strength of the bilayer. The
Po(T) curves have a midpoint of 11.8 dyn/cm, which
exceeds the critical tension for mechanical breakdown
of many phospholipid bilayers (Bloom et al., 1991).
This is also close to the lytic tension for bacterial membranes as suggested by frequent lysis of spheroplast
patches at tensions that activate MscL. This proximity is
consistent with the proposed "safety valve" function for
MscL as a pore that can dissipate osmotic gradients
when the membrane stability is in danger. In a cell of
1-µm diameter, the osmotic gradient required to create
a tension of 12 dyn/cm is 20 mOsm, corresponding to a pressure of 360 mmHg. Soil bacteria may experience
much larger osmotic stresses during rainfall. Knockouts
of MscL have not shown osmotic fragility, however, and
this is probably caused by the presence of other more
sensitive mechanical channels, including MscS and
MscM (Martinac, 1993
; Berrier et al., 1996
). MscL may become active and release larger osmolytes only at extremely strong downshocks, >200 mOsm, as indicated
by osmotically-induced ATP efflux (Sukharev, unpublished observations).
The nonsaturable character of MscL conductance
(Fig. 11) and the absence of anion/cation selectivity
(Sukharev et al., 1993b) strongly suggest a large water-filled pore. The specific conductivity of electrolyte inside the pore would then be similar to its bulk value
and, therefore, the "macroscopic" equations deriving
the channel conductance from its geometry and conductivity are reasonable (Hille, 1992
). The length of
the channel pore was estimated to be
4 nm from sieving data (Cruickshank et al., 1997
).
The energetics of MscL gating suggests that this
channel undergoes large changes in dimensions, accounting for its steep dose-response curve. MS channel
gating is, by definition, a function of force (Corey and
Hudspeth, 1983; Sachs and Lecar, 1991
; Sachs and Morris, 1998
). The simplest model for planar mechanical free
energy is based on the notion that if the open channel
occupies a larger area than the closed channel, then
tension favors opening. This model was used in the energy diagram of Fig. 6, although it should be pointed
out that while the figure is drawn with finite widths for
the wells and barriers, the mathematics assumes the
barrier is a delta function and the wells are boxes.
The dose-response curves in Fig. 5 were well fit by
Eq. 1 with parameters E = 18.6 kBT (46.3 kJ/mol) and
A = 6.5 nm2. Both numbers are large at the molecular
scale and their combination makes the MscL dose-
response curve steep. The large value of
E predicts that
the open probability at low tension (T < 1 dyn/cm) would be ~10
8. Indeed, with no mechanical stimulus,
the impedance of a patch with >200 MscLs is >1 G
closed channels are not leaky and have almost no spontaneous activity. If MscL opened spontaneously, it
would dissipate the proton gradient within ~1 µs and
disrupt cell energetics.
Using the slope of the tension-sensitivity curve, we can
compare MscL gating energetics with other channels.
Alamethicin (Opsahl and Webb, 1994a) and the mechanically sensitive cation-selective channel found in
chick skeletal muscle (Sokabe et al., 1991
) have similar
sensitivities on the order of 3-8 dyn/cm per e-fold
change in Po. The yeast (Gustin, 1991
) and bacterial
MS channels described here have much steeper slope
sensitivities (0.05 and 0.63, respectively), indicative of
large area changes between closed and open states. The midpoints of the activation curves were quite different, however, with the yeast MS channel reaching
half max with only 0.7 dyn/cm tension (Gustin, 1991
)
as compared with MscL, which required 11.8 dyn/cm.
The setpoints and sensitivities of the MS channels so far
examined appeared to be tuned to respond to membrane tension in very different ways, perhaps reflecting
the different physiological roles.
The similarity between the size of the MscL pore and
the change of in-plane area Aoc is striking, and the basic relationship is maintained in the multistate analysis,
where the C-S1, instead of C-O, is rate limiting. These
dimensional changes strongly suggest that the opening
of the pore constitutes the major part of the entire protein complex expansion. When the rate constants of
Table I are taken explicitly in terms of an Eyring model,
differences in k0s reflect differences in the entropy of
activation. In turn, the entropy is a measure of deformability, with narrow energy wells representing stiff conformations and wide energy wells representing soft
conformations (Sachs and Lecar, 1991
). The rate constant k12 has the largest entropy change of all the rates,
suggesting that the closed state is quite flexible relative
to the excited state (the peak of the barrier between 1 and 2). If we examine the change of stiffness between
states by calculating the entropy differences,
S, the
only significant change of stiffness is between the
closed state and S1. Thus, the closed state seems to be
the only one sufficiently soft to allow the available tension to do significant work on the channel.
Do these results agree with structural information?
Each MscL subunit (15 kD, 136 amino acids) spans the
membrane twice with both termini intracellular. The
two transmembrane domains are most likely helical
and oriented normal to the membrane as judged by circular dichroism (Arkin et al., 1998
). Although initial
cross-linking experiments had suggested that the functional MscL complex is a homohexamer (Blount et al.,
1996b
; Sukharev et al., 1997
), more recent and extensive cross-linking studies and the crystal structure of the
closed channel indicate a pentameric stoichiometry
(Cheng et al., 1999). Images of negatively stained two-dimensional crystals of tag-purified MscLs revealed
hexagonal lattices of doughnut-shaped particles (Le
Dain et al., 1998
) that were recognized as hexamers. The
resolution of unsymmetrized projection maps achieved
in this work does, however, seem to be insufficient to
draw an unambiguous conclusion on the number of subunits in the complex. As discussed by Cruickshank et al.
(1997)
, twelve transmembrane helices would be just sufficient to line a pore
4 nm in diameter. The 30 amino acid periplasmic loop might be a part of the pore lining
of the open channel, but is clearly within the cytoplasm
in the closed channel (Cheng et al., 1998).
In examining the energetic model as illustrated in
Fig. 10, the question arises as to why the increase in
conductance from S1 to O5 (presumably an increase
in cross-sectional area of the channel) is not correlated
with a significant increase in the tension sensitivitywe
know that the conductance, and presumably pore diameter, is increasing. There are two types of explanations: (a) the in-plane area really doesn't increase very
much, or (b) the free energy available from tension described as T
A is incomplete. Explanations for a could
involve increases in conductance from shortening the
pore, rather than increasing the cross-section. Alternatively, there may be more complicated conformational
changes to the pore interior that affect ion transport.
Explanations for b might involve components of free
energy such as the line tension, a term proportional to
the perimeter of the channel that has the opposite sign
to T
A. Line tension plays a role in the stability of
pores in lipid bilayers (Zhelev and Needham, 1993
).
Or, we could have stress-induced changes in channel
(water?) entropy that would affect the free energy
(Sachs and Lecar, 1991
). Currently, we cannot separate
these various components, but the relative simplicity of
internal pore reorganization as a method to alter conductance (Fig. 12) has appeal from studies on mechanically insensitive channels.
|
To account for the similarity of A for the rate-limiting step and the pore diameter, the closed conformation would have to exclude most of the water from its
interior (~15-30 nm3), presumably becoming some
sort of compact, ion-impermeable, bundle. If the helices remain normal to the membrane, the pore cannot
be closed by steric interactions unless the channel
forms a close packed trimer. If the helices twist about
the pore axis, the effects of pore diameter are minor as
long as the pore diameter is comparable to the membrane thickness. If the helices rotate axially, forcing hydrophobic faces toward the pore, the channel might
close using ordered water as the "gate" (Sachs and
Feng, 1993
), but this seems unlikely given the large diameter of the structure. If the channel closed by using
the extracellular domain as a gate, then the agreement
of the pore cross-section with the
Aco must be assumed coincidental.
Regardless of details, there are large dimensional
changes involved in opening MscL; changes much
larger than proposed for other channels. Our analysis
has focussed on state models because of the clear finite
residence times in the different conducting states.
However, the large and rapid dimensional changes associated with gating probably also involve inertial components that affect the current rise times. Since the
propagation velocity of shear waves in lipids are on the
order of 10 nm/ms (Debregeas et al., 1998), there may
be useful physical details in the form of the transition
currents, particularly those between C and S1.
Recently, Gu et al. (1998) proposed an electromechanical model of MscL gating that involves the NH2-terminal domains as gates, pivoting under stress from a
position parallel to the membrane to one normal to the
plasmalemma. The model attempts to calculate the
electrostatic force between specific charged residues located on the NH2- and COOH-terminal domains and
those on the membrane-spanning helices. Membrane
tension causes the extracellular end of the helices to tilt
inward, changing the distances between the charges, lowering the force that causes the NH2-terminal domains to swing to a normal position with respect to the
bilayer. However, the balance of electrostatic forces between charged residues was calculated without accounting for electrolyte screening, which dictates that in 0.2 M
salt there will be almost no interaction between two unitary charges positioned 1-nm apart. Our experimental data show there is no significant change in the gating pattern in the range of salt concentration between
0.05 and 1 M. This practically excludes the role of long-range electrostatic interactions in MscL gating. MscL is
also weakly voltage dependent, which is contrary to the model where an external field must strongly influence
the distribution between open and closed states.
The Gu et al. (1998) model also imposes strict constraints on the length and charge of the NH2-terminal
domain. This region must be about the pore radius in
length, and six of them must occlude the pore completely, as we know that closed MscL is absolutely nonleaky. This is difficult to satisfy and also contradicts the work of Blount et al. (1996c)
and Hase et al. (1997)
,
whose data show that the removal of 3, substitution of
8, or addition of 20 new residues to this domain doesn't
significantly change the channel gating.
The model assumes the closed channel is already in a fully patent configuration; i.e. the 4-nm pore is present through all of the open and most of the closed states. Although the model predicts a small increase in area during gating, the ~6-nm2 area increase we calculate (which accounts for the steepness of the dose-response curve) is much larger than predicted. Even if the channel operates through a combination of area changes that pull the gates to a lower energy state followed by the NH2-terminal swinging gates, our evidence for large area changes does not seem compatible with their model. However, the subconductance states could result from individual terminal domains partially interfering with permeation.
Considerations of Our Results in the Context of the Recent Three-Dimensional Structure of MscL
As this paper was in the process of review, Dr. Doug Rees
(California Institute of Technology, Pasadena, CA) kindly
shared with us a preprint of the full x-ray structure of an
MscL homolog (Chang et al., 1998) and we felt that it was
useful to make a first-order comparison with our results.
The structure looks like two barrels in series
a wide one
in the transmembrane portion and a narrow one in the
intracellular compartment. The channel is a pentamer
with each subunit having two alpha helical transmembrane domains tilted at ~28° and an intracellular helical
domain. The intracellular domains form a 3.5-nm continuation of the pore. The result is a channel ~8.5 nm in
length with a diameter that varies from ~1.8 to 0.2 nm,
the latter representing the gate in the closed channel.
This putative gating region consists of a ring of hydrophobic residues located near the intracellular depth of
the bilayer. The outer diameter of the transmembrane
portion, where tension is applied, is ~5 nm.
The most striking contrast with our results (and
other published results) comes from the expected
channel conductance. Knowing the open channel conductance, and modeling it as a cylindrical pore, we can
calculate that a channel 4 nm in length must have a diameter >3 nm to have a conductance of ~3 nS (Table
IV). If the pore were opened to its maximal 1.8-nm diameter along its entire 8.5-nm length, the predicted
conductance would be ~0.5 nS instead of the calculated 3.2 nS. If, upon opening, the intracellular pore
domain were assumed to unfold completely, the transmembrane length would be ~5 nm and the conductance 0.85 nS. Since the simple cylindrical model
(Hille, 1992) assumes no interaction of ions with the
channel, the ions are point charges, and there are no image forces, the conductance estimate should be a
maximum. Furthermore, since the channel is nonselective between anions and cations and its conductance is
exactly proportional to the solution conductance up to
2 M KCl, we cannot invoke local fixed charges as concentrators of ions to increase conductance. It would appear that opening of the channel must involve major alterations in conformation that both shorten and widen
the channel. The constraints are clear if we suppose
tension causes the channel to splay into a cone with the
narrow end of the pore remaining at the observed 1.8-nm diameter (presumably with the narrow end extracellular). As with a cylindrical pore, the conductance of
a tapered pore consists of two convergence resistances
and the pore resistance given by Rpore =
l/
r1 × r2,
where the r is the radii at each end,
is the solution resistivity, and l is the pore length. If the intracellular portion of the channel were folded out of the way against
the bilayer so the pore was only 4-nm long, we would
still have to expand the internal diameter to 100 nm to
get the observed 3-nS conductance. Constraining even
one end of the channel to 1.8-nm diameter strongly
limits the possible pore conductance. However, it is
perhaps not surprising that the pore dimensions of the
closed channel differ greatly from the dimensions predicted for the open pore.
Concerning the structural origin of the mechanical sensitivity, the tension sensitivity can be explained by a modest increase in the outer diameter of the transmembrane domain from ~5.0 to 5.5 nm. If each of the 10 transmembrane helices were 12-13 nm in diameter and arranged in a close-packed ring normal to the membrane, the outer diameter would be 5-6 nm, in the range necessary to account for the mechanical sensitivity. If tension untwists the transmembrane helices to form a set of barrel staves perpendicular to the membrane, the cytoplasmic domains may peel away, shortening the length of the pore. The substates we observed may reflect such movements of the cytoplasmic helices. These domains, which are outside the bilayer, should not be strongly driven by membrane tension, making them compatible with the lack of tension sensitivity of the higher conductance substates.
There are methodological questions to be resolved between the crystallography and the electrophysiology. The crystals were formed in solutions at pH 3.7 and were stabilized with glutaraldehyde as well as heavy metal compounds, including Gd+3, which is a known blocker of the channel. These conditions may create structures different from the native state and physiological tests need to be made on channels treated this way. Nonetheless, it is a thrill to have a real structure to examine and we eagerly await a structure of the open channel.
In summary, we have performed the first calibration of a biological mechanosensitive ion channel. These measurements place strong constraints on kinetic and structural models of MscL and related channels.
![]() |
FOOTNOTES |
---|
Address correspondence to Wade J. Sigurdson, Ph.D., Dept. Physiology and Biophysics, 320 Cary Hall, SUNY at Buffalo, Buffalo, NY 14214. Fax: 716-829-2028; E-mail: wjs{at}buffalo.edu
Original version received 17 June 1998 and accepted version received 3 February 1999.
S.I. Sukharev and W.J. Sigurdson contributed equally to this work and should be considered co-first authors.We thank the United States Army Research Office, National Aeronautics and Space Administration, and the National Institutes of Health for grant support.
![]() |
Abbreviation used in this paper |
---|
MS, mechanosensitive.
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
1. | Akinlaja, J. 1997. Electromechanical rupture of cell-attached patches. Physics. Thesis/Dissertation. SUNY at Buffalo, Buffalo, NY. 1-56. |
2. | Arkin, I.T., S. Sukharev, P. Blount, C. Kung, and A.T. Brunger. 1998. Helicity, membrane incorporation, orientation and thermal stability of the large conductance mechano-sensitive ion channel from E. coli. Biochim. Biophys. Acta. 1369: 131-140 [Medline]. |
3. |
Awayda, M.S.,
I.I. Ismailov,
B.K. Berdiev, and
D.J. Benos.
1995.
A
cloned renal epithelial Na+ channel protein displays stretch activation in planar lipid bilayers.
Am. J. Physiol.
268:
C1450-C1459
|
4. | Berrier, C., M. Besnard, B. Ajouz, A. Coulombe, and A. Ghazi. 1996. Multiple mechanosensitive ion channels from Escherichia coli, activated at different thresholds of applied pressure. J. Membr. Biol. 151: 175-187 [Medline]. |
5. | Bloom, M., E. Evans, and O.G. Mouritsen. 1991. Physical properties of the fluid lipid-bilayer component of cell membranes: a perspective. Q. Rev. Biophys. 24: 293-397 [Medline]. |
6. |
Blount, P.,
M.J. Schroeder, and
C. Kung.
1998.
Mutations in a bacterial mechanosensitive channel change the cellular response to
osmotic stress.
J. Biol. Chem.
272:
32150-32157
|
7. | Blount, P., S. Sukharev, M. Schroeder, S. Nagle, and C. Kung. 1996a. Mutations that change gating properties of a mechanosensitive channel in E. coli. Biophys. J. 70: A366 . (Abstr.) . |
8. | Blount, P., S.I. Sukharev, P.C. Moe, S.K. Nagle, and C. Kung. 1996b. Towards an understanding of the structural and functional properties of MscL, a mechanosensitive channel in bacteria. Biol Cell. 87: 1-8 [Medline]. |
9. | Blount, P., S.I. Sukharev, P.C. Moe, M.J. Schroeder, R.H. Guy, and C. Kung. 1996c. Membrane topology and multimeric structure of a mechanosensitive channel protein of Escherichia coli. EMBO (Eur. Mol. Biol. Organ.) J. 15: 4798-4805 [Abstract]. |
10. | Brehm, P., R. Kullberg, and F. Moody-Corbett. 1984. Properties of non-junctional acetylcholine receptor channels on innervated muscle of Xenopus Laevis. J. Physiol. (Lond.). 350: 631-648 [Abstract]. |
11. |
Chang, G.,
R.H. Spencer,
A.T. Lee,
M.T. Barclay, and
D.C. Rees.
1998.
Structure of the MscL homolog from Mycobacterium tuberculosis: a gated mechanosensitive ion channel.
Science
282:
2220-2226
|
12. | Colquhoun, D., and F.J. Sigworth. 1983. Fitting and statistical analysis of single-channel records. In Single-Channel Recording. B. Sakmann and E. Neher, editors. Plenum Publishing Corp., New York. 191. |
13. | Corey, D.P., and A.J. Hudspeth. 1983. Kinetics of the receptor current in bullfrog saccular hair cells. J. Neurosci. 3: 962-976 [Abstract]. |
14. | Cruickshank, C.C., R.F. Minchin, A.C. Le Dain, and B. Martinac. 1997. Estimation of the pore size of the large-conductance mechanosensitive ion channel of Escherichia coli. Biophys. J. 73: 1925-1931 [Abstract]. |
15. |
Debregeas, D.,
P.G. de Gennes, and
F. Brochard-Wyart.
1998.
The
life and death of "bare" viscous bubbles.
Science.
279:
1704-1707
|
16. | Driscoll, M., and M. Chalfie. 1993. The mec-4 gene is a member of a family of Caenorhabditis elegans genes that can mutate to induce neuronal degeneration. Nature. 349: 588-593 . |
17. | Feng, Q., A. Auerbach, and F. Sachs. 1996. Estimating single channel kinetic parameters from idealized patch-clamp data containing missed events. Biophys. J. 70: 264-280 [Abstract]. |
18. |
Gu, L.,
W. Liu, and
B. Martinac.
1998.
Electromechanical coupling
model of gating the large mechanosensitive ion channel (MscL)
of Escherichia coli by mechanical force.
Biophys J.
74:
2889-2902
|
19. | Guharay, F., and F. Sachs. 1984. Stretch-activated single ion channel currents in tissue-cultured embryonic chick skeletal muscle. J. Physiol. (Lond.). 352: 685-701 [Abstract]. |
20. | Gustin, M.C.. 1991. Single-channel mechanosensitive currents. Science. 253: 800 . |
21. | Hamill, O.P., and D.W. McBride Jr.. 1994. The cloning of a mechano-gated membrane ion channel. Trends Neurosci. 17: 439-443 [Medline]. |
22. | Hase, C.C., A.C. Le Dain, and B. Martinac. 1997. Molecular dissection of the large mechanosensitive ion channel (MscL) of E. coli: mutants with altered channel gating and pressure sensitivity. J. Membr. Biol. 157: 17-25 [Medline]. |
23. | Hille, B. 1992. Ionic Channels of Excitable Membranes. Sinauer Associates Inc., Sunderland, MA. 185-188. |
24. |
Kizer, N.,
X.L. Guo, and
K. Hruska.
1997.
Reconstitution of stretch-activated cation channels by expression of the alpha-subunit of
the epithelial sodium channel cloned from osteoblasts.
Proc. Natl.
Acad. Sci. USA.
94:
1013-1018
|
25. |
Le Dain, A.C.,
N. Saint,
A. Kloda,
A. Ghazi, and
B. Martinac.
1998.
Mechanosensitive ion channels of the archaeon
Haloferax volcanii. J. Biol. Chem.
273:
12116-12119
.
|
26. | Martinac, B. 1993. Mechanosensitive ion channels: biophysics and physiology. In Thermodynamics of Cell Surface Receptors. M. Jackson, editor. CRC Press, Boca Raton, FL. 327-351. |
27. | Martinac, B., A.H. Delcour, M. Buechner, J. Adler, and C. Kung. 1992. Mechanosensitive ion channels in bacteria. In Comparative Aspects of Mechanoreceptors. F. Ito, editor. Springer-Verlag, Berlin. 3-18. |
28. | Moe, P.C., P. Blount, and C. Kung. 1998. Functional and structural conservation in the mechanosensitive channel MscL implicates elements crucial for mechanosensation. Mol. Microbiol. 28: 583 [Medline]. |
29. | Opsahl, L.R., and W.W. Webb. 1994a. Lipid-glass adhesion in giga-sealed patch-clamped membranes. Biophys. J. 66: 75-79 [Abstract]. |
30. | Opsahl, L.R., and W.W. Webb. 1994b. Transduction of membrane tension by the ion channel alamethicin. Biophys. J. 66: 71-74 [Abstract]. |
31. | Qin, F., A. Auerbach, and F. Sachs. 1995. Maximum likelihood estimation of kinetic parameters for single channels from dwell-time sequences. Biophys. J. 68: A79 . (Abstr.) . |
32. | Qin, F., A. Auerbach, and F. Sachs. 1996. Idealization of single channel currents using the segmental K-means method. Biophys. J. 70: A227 . (Abstr.) . |
33. | Rossier, B.C., C.M. Canessa, L. Schild, and J.D. Horisberger. 1994. Epithelial sodium channels. Curr. Opin. Nephrol. Hypertens. 3: 487-496 [Medline]. |
34. | Sachs, F.. 1987. Baroreceptor mechanisms at the cellular level. Fed. Proc. 46: 12-16 [Medline]. |
35. | Sachs, F., and Q. Feng. 1993. Gated, ion-selective channels observed without membranes: novel properties of the gigaseal. Biophys. J. 65: 1101-1107 [Abstract]. |
36. | Sachs, F., and H. Lecar. 1991. Stochastic models for mechanical transduction. Biophys. J. 59: 1143-1145 [Medline]. |
37. | Sachs, F., and C. Morris. 1998. Mechanosensitive ion channels in non specialized cells. In Reviews of Physiology and Biochemistry and Pharmacology. M.P. Blaustein, R. Greger, H. Grunicke, R. Jahn, L.M. Mendell, A. Miyajima, D. Pette, G. Schultz, and M. Schweiger, editors. Springer-Verlag, Berlin. 1-78. |
38. | Sachs, F., J. Neil, and N. Barkakati. 1982. The automated analysis of data from single ionic channels. Pflügers Arch. 395: 331-340 [Medline]. |
39. | Saimi, Y., B. Martinac, A.H. Delcour, P.V. Minorsky, M.C. Gustin, M.R. Culbertson, J. Adler, and C. Kung. 1993. Patch clamp studies of microbial ion channels. Methods Enzymol. 207: 681-691 . |
40. |
Sokabe, M.,
F. Sachs, and
Z. Jing.
1991.
Quantitative video microscopy of patch clamped membranes![]() |
41. | Sokabe, M., W.S. Sigurdson, and F. Sachs. 1993. Effect of excision and cytochalasin on the viscoelastic properties of patch clamped membranes in heart and skeletal muscle. Biophys. J. 61: A513 . (Abstr.) . |
42. | Sukharev, S., P. Blount, M. Schroeder, and C. Kung. 1996. Multimeric structure of bacterial mechanosensitive channel MscL. Biophys. J. 70: A366 . (Abstr.) . |
43. | Sukharev, S.I., P. Blount, B. Martinac, F.R. Blattner, and C. Kung. 1994a. A large conductance mechanosensitive channel in E. coli encoded by MscL alone. Nature. 368: 265-268 [Medline]. |
44. | Sukharev, S.I., P. Blount, B. Martinac, and C. Kung. 1997. Mechanosensitive channels of Escherichia coli: the MscL gene, protein and activities. Annu. Rev. Physiol. 59: 633-657 [Medline]. |
45. | Sukharev, S.I., B. Martinac, P. Blount, and C. Kung. 1994b. Functional reconstitution as an assay for biochemical isolation of channel proteins: application to the molecular identification of a bacterial mechanosensitive channel. Methods: A Companion to Methods in Enzymology. 6:51-59. |
46. | Sukharev, S.I., B. Martinac, and C. Kung. 1993. Reconstitution of two distinct types of mechanosensitive channels from the E. Coli envelope. Biophys. J. 61: A93 . (Abstr.) . |
47. | Tavernakis, N., and M. Driscoll. 1998. Molecular modeling of mechanotransduction in the nematode Caenorhabditis elegans. Annu. Rev. Physiol. 59: 659-689 [Medline]. |
48. | Zhang, B., W.J. Sigurdson, and F. Sachs. 1997. Storing analog data in a video record. J. Neurosci. Methods. 76: 151-155 [Medline]. |
49. | Zhelev, D.V., and D. Needham. 1993. Tension stabilized pores in giant vesicles: determination of pore size and pore line tension. Biochim. Biophys. Acta. 1147: 89-104 [Medline]. |