Mechanosensitivity of gramicidin A channels in semispherical bilayer membranes at constant tension

Mechanoelectrical transduction in gramicidin A channels was studied in macroscopic planar lipid bilayer membranes bulged at constant tension. We found a supralinear increase in the single channel activity, which was proportional to the square of membrane radius but could not be accounted for by the increase in membrane surface area or by recruitment of new channels. When being extrapolated to biological membranes, these observations may suggest that the activity of permeability of ion channels can be influenced simply by changing the shape of the membrane, with or without stretching.     

Theoretical analysis of hydrophobic matching and membrane-mediated interactions in lipid bilayers containing gramicidin.

We present a quantitative analysis of the effects of hydrophobic matching and membrane-mediated protein-protein interactions exhibited by gramicidin embedded in dimyristoylphosphatidylcholine (DMPC) and dilauroylphosphatidylcholine (DLPC) bilayers (Harroun et al., 1999. Biophys. J. 76:937-945). Incorporating gramicidin, at 1:10 peptide/lipid molar ratio, decreases the phosphate-to-phosphate (PtP) peak separation in the DMPC bilayer from 35.3 A without gramicidin to 32.7 A. In contrast, the same molar ratio of gramicidin in DLPC increases the PtP from 30.8 A to 32.1 A. Concurrently, x-ray in-plane scattering showed that the most probable nearest-neighbor separation between gramicidin channels was 26.8 A in DLPC, but reduced to 23.3 A in DMPC. In this paper we review the idea of hydrophobic matching in which the lipid bilayer deforms to match the hydrophobic surface of the embedded proteins. We use a simple elasticity theory, including thickness compression, tension, and splay terms to describe the membrane deformation. The energy of membrane deformation is compared with the energy cost of hydrophobic mismatch. We discuss the boundary conditions between a gramicidin channel and the lipid bilayer. We used a numerical method to solve the problem of membrane deformation profile in the presence of a high density of gramicidin channels and ran computer simulations of 81 gramicidin channels to find the equilibrium distributions of the channels in the plane of the bilayer. The simulations contain four parameters: bilayer thickness compressibility 1/B, bilayer bending rigidity Kc, the channel-bilayer mismatch Do, and the slope of the interface at the lipid-protein boundary s. B, Kc, and Do were experimentally measured; the only free parameter is s. The value of s is determined by the requirement that the theory produces the experimental values of bilayer thinning by gramicidin and the shift in the peak position of the in-plane scattering due to membrane-mediated channel-channel interactions. We show that both hydrophobic matching and membrane-mediated interactions can be understood by the simple elasticity theory.     

Mechanosensitivity of gramicidin A channels in bulged bilayer membranes at constant tension

Mechanoelectrical transduction in gramicidin A channels was studied in macroscopic planar lipid bilayer membranes bulged at constant tension. We found a supralinear increase in the single channel activity that was proportional to the square of membrane radius, but could not be accounted for by the increase in membrane surface area, or by recruitment of new channels. Extrapolated to biological membranes, these observations may suggest that the permeability of ion channels can be influenced simply by changing shape of the membrane, with or without stretching. Published in Russian in Biofizika, 2006, Vol. 51, No. 6, pp. 1014–1018. The text was submitted by the authors in English.     

Influence of ion occupancy and membrane deformation on gramicidin A channel stability in lipid membranes.

The average lifetime of gramicidin A channels in monoolein/decane bilayer membranes was measured. The results support the hypothesis of channel stabilization by ion occupancy. The effects of electric field and salt concentration are consistent with the expected effects on both occupancy and membrane compression. The lifetime in asymmetric solutions with divalent cation blockers on one side of the membrane shows a voltage dependence such that the lifetime decreases for positive voltages applied from the blocking side and increases for negative voltages. This result strongly supports the occupancy hypothesis. The lifetime increases with permeant ion concentration, and at the one molar level it also increases with voltage. The voltage dependence of lifetime for a low concentration of permeant ion depends on the total salt level. The results for these conditions are consistent with the assumption that membrane compression also influences the lifetime, even for the "soft" solvent-containing membrane considered here. It is proposed that the channel nearest neighbor lipids need not be fixed in a plane at the channel end. Using a liquid crystal model it may then be shown that surface tension is the major component of the membrane deformation free energy, which may explain the significant effects of the membrane compression on the lifetime.     

The structure of the gramicidin A transmembrane channel

Gramicidin A is a linear polypeptide antibiotic that facilitates the diffusion of monovalent cations across lipid bilayer membranes by forming channels. It has been proposed that the conducting channel is a dimer which is in equilibrium with nonconducting monomers in the membrane. To directly test this model in several independent ways, we have prepared and purified a series of gramicidin C derivatives. All of these derivatives are fully active analogs of gramicidin A, and each derivative has a useful chromophore esterified to the phenolic hydroxyl of tyrosine #11. Simultaneous conductance and fluorescence measurements on planar lipid bi-layer membranes containing dansyl gramicidin C yielded four conclusions: (1) A plot of the logarithm of the membrane conductance versus the logarithm of the membrane fluorescence had a slope of 2.0 ± 0.3, over a concentration range for which nearly all the gramicidin was monomeric. Hence, the active channel is a dimer of the nonconducting species. (2) In a membrane in which nearly all of the gramicidin was dimeric, the number of channels was approximately equal to the number of dimers. Thus, most dimers are active channels and so it should be feasible to carry out spectroscopic studies of the conformation of the transmembrane channel. (3) The association constant for dimerization is more than 1,000-fold larger in a glycerolester membrane with 26 Å-hydrocarbon thickness than in a 47 Å-glycerolester membrane. The dimerization constant in a 48 Å-phosphatidyl choline membrane was 200 times larger than in a 47 Å-glycerolester membrane, showing that it depends on the type of lipid as well as on the thickness of the hydrocarbon core. (4) We were readily able to detect 10^?14 mole cm^?2 of dansyl gramicidin C in a bilayer membrane, which corresponds to 60 fluorescent molecules per square μm. The fluorescent techniques described here should be sufficiently sensitive for fluorescence studies of reconstituted gates and receptors in planar bilayer membranes. An alternative method of determining the number of molecules of gramicidin in the channel is to measure the fraction of hybrid channels present in a mixture of 2 chemically different gramicidins. The single-channel conductance of p-phenylazo-benzene-sulfonyl ester gramicidin C (PABS gramicidin C) was found to be 0.68 that of gramicidin A. In membranes containing a mixture of these 2 gramicidins, a hybrid channel was evident in addition to 2 pure channels. The hybrid channel conductance was 0.82 that of gramicidin A. Fluorescence energy transfer from dansyl gramicidin C to diethylamino-phenylazobenzene-sulfonyl ester gramicidin C (DPBS gramicidin C), provided an independent way to measure the fraction of hybrid channels on liposomes. For both techniques the fraction of hybrid channels was found to be 2ad where a² and d² were the fractions of the 2 kinds of pure channels. This result strongly supports a dimer channel and the hybrid data excludes the possibility of a tetramer channel. The study of hybrid species by conductance and fluorescence techniques should be generally useful in elucidating the subunit structure of oligomeric assemblies in membranes. The various models which have been proposed for the conformation of the gramicidin transmembrane channel are briefly discussed.     

Stability of an ion channel in lipid bilayers: implicit solvent model calculations with gramicidin

Gramicidin is a helical peptide, 15 residues in length, which dimerizes to form ion-conducting channels in lipid bilayers. Here we report calculations of its free energy of transfer from the aqueous phase into bilayers of different widths. The electrostatic and nonpolar contributions to the desolvation free energy were calculated using implicit solvent models, in which gramicidin was described in atomic detail and the hydrocarbon region of the membrane was described as a slab of hydrophobic medium embedded in water. The free energy penalties from the lipid perturbation and membrane deformation effects, and the entropy loss associated with gramicidin immobilization in the bilayer, were estimated from a statistical thermodynamic model of the bilayer. The calculations were carried out using two classes of experimentally observed conformations: a head-to-head dimer of two single-stranded (SS) beta-helices and a double-stranded (DS) intertwined double helix. The calculations showed that gramicidin is likely to partition into the bilayer in all of these conformations. However, the SS conformation was found to be significantly more stable than the DS in the bilayer, in agreement with most of the experimental data. We tested numerous transmembrane and surface orientations of gramicidin in bilayers of various widths. Our calculations indicate that the most favorable orientation is transmembrane, which is indeed to be expected from a channel-forming peptide. The calculations demonstrate that gramicidin insertion into the membrane is likely to involve a significant deformation of the bilayer to match the hydrophobic width of the peptide (22 A), again in good agreement with experimental data. Interestingly, deformation of the bilayer was induced by all of the gramicidin conformations.     

Lysophospholipids modulate channel function by altering the mechanical properties of lipid bilayers

Lipid metabolites, free fatty acids and lysophospholipids, modify the function of membrane proteins including ion channels. Such alterations can occur through signal transduction pathways, but may also result from "direct" effects of the metabolite on the protein. To investigate possible mechanisms for such direct effects, we examined the alterations of gramicidin channel function by lysophospholipids (LPLs): lysophosphatidylcholine (LPC), lysophosphatidylethanolamine (LPE), lysophosphatidylserine (LPS), and lysophosphatidylinositol (LPI). The experiments were done on planar bilayers formed by diphytanoylphosphatidylcholine in n-decane a system where receptor- mediated effects can be excluded. At aqueous concentrations below the critical micelle concentration (CMC), LPLs can increase the dimerization constant for membrane-bound gramicidin up to 500-fold (at 2 microM). The relative potency increases as a function of the size of the polar head group, but does not seem to vary as a function of head group charge. The increased dimerization constant results primarily from an increase in the rate constant for channel formation, which can increase more than 100-fold (in the presence of LPC and LPI), whereas the channel dissociation rate constant decreases only about fivefold. The LPL effect cannot be ascribed to an increased membrane fluidity, which would give rise to an increased channel dissociation rate constant. The ability of LPC to decrease the channel dissociation rate constant varies as a function of channel length (which is always less than the membrane's equilibrium thickness): as the channel length is decreased, the potency of LPC is increased. LPC has no effect on membrane thickness or the surface tension of monolayers at the air/electrolyte interface. The bilayer-forming glycerolmonooleate does not decrease the channel dissociation rate constant. These results show that LPLs alter gramicidin channel function by altering the membrane deformation energy, and that the changes in deformation energy can be related to the molecular "shape" of the membrane-modifying compounds. Similar alterations in the mechanical properties of biological membranes may form a general mechanism by which one can alter membrane protein function.     

Brief closures of gramicidin A channels in lipid bilayer membranes

Brief closures, so called flickers, gramicidin A channels were observed for glycerol monooleate/n-decane membranes for cesium chloride and hydrochloric acid solutions. The flickers, similar in nature to the flickers observed for physiological channels, were of the order of 1 ms and the interval between flickers was of the order of 50 ms. The flicker-duration and interval between flickers both decrease with voltage. The field dependence of the flickers is consistent with the hypothesis that the membrane forms a dimple when accomodating a dimer in the membrane and that the monomers, on breaking up, are associated over displacements of the order of 2 nm. For similar measurements for glycerol monoleate/hexadecane membranes only rare occurrences of flickers were observed. It is suggested that the flicker phenomenon is governed by the physical and chemical properties of the membrane and the influence of membrane thickness and interfacial free energy is emphasized.     

Molecular dynamics - potential of mean force calculations as a tool for understanding ion permeation and selectivity in narrow channels

Ion channels catalyze the permeation of charged molecules across cell membranes and are essential for many vital physiological functions, including nerve and muscle activity. To understand better the mechanisms underlying ion conduction and valence selectivity of narrow ion channels, we have employed free energy techniques to calculate the potential of mean force (PMF) for ion movement through the prototypical gramicidin A channel. Employing modern all-atom molecular dynamics (MD) force fields with umbrella sampling methods that incorporate one hundred 1-2 ns trajectories, we find that it is possible to achieve semi-quantitative agreement with experimental binding and conductance measurements. We also examine the sensitivity of the MD-PMF results to the choice of MD force field and compare PMFs for potassium, calcium and chloride ions to explore the basis for the valence selectivity of this narrow and uncharged ion channel. A large central barrier is observed for both anions and divalent ions, consistent with lack of experimental conductance. Neither anion or divalent cation is seen to be stabilized inside the channel relative to the bulk electrolyte and each leads to large disruptions to the protein and membrane structure when held deep inside the channel. Weak binding of calcium ions outside the channel corresponds to a free energy well that is too shallow to demonstrate channel blocking. Our findings emphasize the success of the MD-PMF approach and the sensitivity of ion energetics to the choice of biomolecular force field.     

Internal forces, tension and energy density in tethered cellular membranes

We analyze tethered cellular membranes by considering the membrane resultants, tension and densities of two modes of energy, bending and adhesion. These characteristics are determined based on a computational (finite-difference) analysis of membrane shape. We analyze the relative contribution and distribution of the membrane characteristics in four typical zones of the membrane surface. Using an axisymmetric model, we found that the meridional and circumferential components of the resultant are different near the tether body and they converge to the value of membrane tension farther from the tether. At the beginning of the area of membrane detachment from the cytoskeleton, the density of bending energy is on the same order of magnitude as membrane tension (resultant). Away from the tether, the bending energy density quickly decreases and becomes of the same order as that of the adhesion energy in the membrane-cytoskeleton attachment area. In that area, both modes of energy are significantly smaller than the membrane tension. We also consider the effect of the membrane bending modulus on the distribution of the membrane characteristics. An increase in the bending modulus results in changing the length and position on the membrane surface of zone 1 characterized by significant evolution of the resultant components. It also results in shortening zone 2 that covers the rest of the area of membrane detachment. The obtained results can help in a better interpretation of the measurements of membrane mechanical properties as well as in analyses of proteins and channels in curved membranes.     

Modulation of gramicidin channel conformation and organization by hydrophobic mismatch in saturated phosphatidylcholine bilayers

The matching of hydrophobic lengths of integral membrane proteins and the surrounding lipid bilayer is an important factor that influences both structure and function of integral membrane proteins. The ion channel gramicidin is known to be uniquely sensitive to membrane properties such as bilayer thickness and membrane mechanical properties. The functionally important carboxy terminal tryptophan residues of gramicidin display conformation-dependent fluorescence which can be used to monitor gramicidin conformations in membranes [S.S. Rawat, D.A. Kelkar, A. Chattopadhyay, Monitoring gramicidin conformations in membranes: a fluorescence approach, Biophys. J. 87 (2004) 831-843]. We have examined the effect of hydrophobic mismatch on the conformation and organization of gramicidin in saturated phosphatidylcholine bilayers of varying thickness utilizing the intrinsic conformation-dependent tryptophan fluorescence. Our results utilizing steady state and time-resolved fluorescence spectroscopic approaches, in combination with circular dichroism spectroscopy, show that gramicidin remains predominantly in the channel conformation and gramicidin tryptophans are at the membrane interfacial region over a range of mismatch conditions. Interestingly, gramicidin conformation shifts toward non-channel conformations in extremely thick gel phase membranes although it is not excluded from the membrane. In addition, experiments utilizing self quenching of tryptophan fluorescence indicate peptide aggregation in thicker gel phase membranes.     

Modulation of gramicidin channel conformation and organization by hydrophobic mismatch in saturated phosphatidylcholine bilayers

The matching of hydrophobic lengths of integral membrane proteins and the surrounding lipid bilayer is an important factor that influences both structure and function of integral membrane proteins. The ion channel gramicidin is known to be uniquely sensitive to membrane properties such as bilayer thickness and membrane mechanical properties. The functionally important carboxy terminal tryptophan residues of gramicidin display conformation-dependent fluorescence which can be used to monitor gramicidin conformations in membranes [S.S. Rawat, D.A. Kelkar, A. Chattopadhyay, Monitoring gramicidin conformations in membranes: a fluorescence approach, Biophys. J. 87 (2004) 831-843]. We have examined the effect of hydrophobic mismatch on the conformation and organization of gramicidin in saturated phosphatidylcholine bilayers of varying thickness utilizing the intrinsic conformation-dependent tryptophan fluorescence. Our results utilizing steady state and time-resolved fluorescence spectroscopic approaches, in combination with circular dichroism spectroscopy, show that gramicidin remains predominantly in the channel conformation and gramicidin tryptophans are at the membrane interfacial region over a range of mismatch conditions. Interestingly, gramicidin conformation shifts toward non-channel conformations in extremely thick gel phase membranes although it is not excluded from the membrane. In addition, experiments utilizing self quenching of tryptophan fluorescence indicate peptide aggregation in thicker gel phase membranes.     

The dimeric nature of the gramicidin A transmembrane channel: Conductance and fluorescence energy transfer studies of hybrid channels

Gramicidin A, a linear peptide antibiotic, makes membranes permeable to alkali cations and hydrogen ions by forming transmembrane channels. We report here conductance and fluorescence energy transfer studies of channels containing two kinds of gramicidin. These studies of hybrid channels were designed to determine the number of molecules in a channel. The gramicidins studied were gramicidin A, dansyl gramicidin C, the p-phenylazobenzene sulfonyl derivative of gramicidin C (PABS⁴ gramicidin C), and the 4-(diethylamino)-phenylazobenzene-4-sulfonyl chloride derivative of gramicidin C (DPBS gramicidin C). The dansyl, PABS and DPBS groups were linked to the hydroxyl group of tyrosine 11 in gramicidin C. The single-channel conductance of PABS gramicidin C in planar bilayer membranes is 0.68 that of gramicidin A. Membranes containing both PABS gramicidin C and gramicidin A exhibit three kinds of channels: a pure gramicidin A, a pure PABS gramicidin C channel, and a hybrid channel with an intermediate conductance (0.82 that of gramicidin A). The dependence of the frequencies of these three kinds of channels on the mole fractions of gramicidin A and PABS gramicidin C in the membrane-forming solution fits a dimer model. Fluorescence energy transfer was used as a complementary means of ascertaining the frequency of hybrid channels. Dansyl gramicidin C was the fluorescent energy donor and DPBS gramicidin C was the energy acceptor. The efficiency of energy transfer between these chromophores in hybrid channels in liposomes was 75%. The relative quantum yield of the dansyl fluorescence was measured as a function of the mole fraction of DPBS gramicidin C. These fluorescence studies, like the single-channel conductance measurements, showed that there are two molecules of gramicidin in a channel. The study of hybrid species by conductance and fluorescence techniques should be generally useful in elucidating the subunit structure of oligomeric assemblies in membranes.     

Structure of gramicidin A.

Gramicidin A, a hydrophobic linear polypeptide, forms channels in phospholipid membranes that are specific for monovalent cations. Nuclear Magnetic Resonance (NMR) spectroscopy provided the first direct physical evidence that the channel conformation in membranes is an amino terminal-to-amino terminal helical dimer, and circular dichroism (CD) spectroscopy has shown the sensitivity of its conformation to different environments and the structural consequences of ion binding. The three-dimensional structure of a gramicidin/cesium complex has been determined by x-ray diffraction of single crystals using single wavelength anomalous scattering for phasing. The left-handed double helix in this crystal form corresponds to one of the intermediates in the process of folding and insertion into membranes. Co-crystals of gramicidin and lipid that appear to have gramicidin in their membrane channel conformation have also been formed and are presently under investigation. Hence, we have used a combination of spectroscopic and diffraction techniques to examine the conformation and functionally-related structural features of gramicidin A.     

Molecular dynamics simulations of membrane proteins

Membrane proteins control the traffic across cell membranes and thereby play an essential role in cell function from transport of various solutes to immune response via molecular recognition. Because it is very difficult to determine the structures of membrane proteins experimentally, computational methods have been increasingly used to study their structure and function. Here we focus on two classes of membrane proteins—ion channels and transporters—which are responsible for the generation of action potentials in nerves, muscles, and other excitable cells. We describe how computational methods have been used to construct models for these proteins and to study the transport mechanism. The main computational tool is the molecular dynamics (MD) simulation, which can be used for everything from refinement of protein structures to free energy calculations of transport processes. We illustrate with specific examples from gramicidin and potassium channels and aspartate transporters how the function of these membrane proteins can be investigated using MD simulations.     

Entropic tension in crowded membranes

Unlike their model membrane counterparts, biological membranes are richly decorated with a heterogeneous assembly of membrane proteins. These proteins are so tightly packed that their excluded area interactions can alter the free energy landscape controlling the conformational transitions suffered by such proteins. For membrane channels, this effect can alter the critical membrane tension at which they undergo a transition from a closed to an open state, and therefore influence protein function in vivo. Despite their obvious importance, crowding phenomena in membranes are much less well studied than in the cytoplasm. Using statistical mechanics results for hard disk liquids, we show that crowding induces an entropic tension in the membrane, which influences transitions that alter the projected area and circumference of a membrane protein. As a specific case study in this effect, we consider the impact of crowding on the gating properties of bacterial mechanosensitive membrane channels, which are thought to confer osmoprotection when these cells are subjected to osmotic shock. We find that crowding can alter the gating energies by more than [Formula: see text] in physiological conditions, a substantial fraction of the total gating energies in some cases. Given the ubiquity of membrane crowding, the nonspecific nature of excluded volume interactions, and the fact that the function of many membrane proteins involve significant conformational changes, this specific case study highlights a general aspect in the function of membrane proteins.     

Deformation free energy of bilayer membrane and its effect on gramicidin channel lifetime.

The deformation free energy of a lipid bilayer is presented based on the principle of a continuum theory. For small deformations, the free energy consists of a layer-compression term, a splay-distortion term, and a surface-tension term, equivalent to the elastic free energy of a two-layer smectic liquid crystal with surface tension. Minimization of the free energy leads to a differential equation that, with boundary conditions, determines the elastic deformation of a bilayer membrane. When a dimeric gramicidin channel is formed in a membrane of thickness greater than the length of the channel, the membrane deformation reduces the stability of the channel. Previously this effect was studied by comparing the variation of channel lifetime with the surface tension of bilayers (Elliott, J. R., D. Needham, J. P. Dilger, and D. A. Hayden, 1983, Biochim. Biophys. Acta, 735:95-103). The tension was assumed to pull a dimer for a distance z before the channel loses ion conductivity. To account for the data, z was found to be 18 A. With the deformation free energy, the data can be accounted for with z less than or approximately to 1 A, which is consistent with the breaking of hydrogen bonds in a dimer dissociation. Increasing the strength of lipid-protein interactions is not the only consequence of the complete free energy compared with the previous discussions. It also changes the shape of membrane deformation around an embedded channel from convex to concave, and increases the range of deformation from less than 10 A to greater than 20 A. Clearly these will be important factors in the general considerations of lipid-protein interactions and membrane-mediated interactions between proteins. In addition, thermal fluctuations of a membrane are calculated; in particular, we calculate the relations between the intrinsic thickness and the experimentally measured values. The experimental parameters of monoolein-squalene membranes are used for quantitative analyses.