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.     

Experimental evidence for hydrophobic matching and membrane-mediated interactions in lipid bilayers containing gramicidin.

Hydrophobic matching, in which transmembrane proteins cause the surrounding lipid bilayer to adjust its hydrocarbon thickness to match the length of the hydrophobic surface of the protein, is a commonly accepted idea in membrane biophysics. To test this idea, gramicidin (gD) was embedded in 1, 2-dilauroyl-sn-glycero-3-phosphocholine (DLPC) and 1, 2-myristoyl-sn-glycero-3-phosphocholine (DMPC) bilayers at the peptide/lipid molar ratio of 1:10. Circular dichroism (CD) was measured to ensure that the gramicidin was in the beta6.3 helix form. The bilayer thickness (the phosphate-to-phosphate distance, or PtP) was measured by x-ray lamellar diffraction. In the Lalpha phase near full hydration, PtP is 30.8 A for pure DLPC, 32.1 A for the DLPC/gD mixture, 35.3 A for pure DMPC, and 32.7 A for the DMPC/gD mixture. Gramicidin apparently stretches DLPC and thins DMPC toward a common thickness as expected by hydrophobic matching. Concurrently, gramicidin-gramicidin correlations were measured by x-ray in-plane scattering. In the fluid phase, the gramicidin-gramicidin nearest-neighbor separation is 26.8 A in DLPC, but shortens to 23.3 A in DMPC. These experiments confirm the conjecture that when proteins are embedded in a membrane, hydrophobic matching creates a strain field in the lipid bilayer that in turn gives rise to a membrane-mediated attractive potential between proteins.     

Voltage-dependent formation of gramicidin channels in lipid bilayers.

The formation kinetics of gramicidin A channels in lipid bilayer membranes has been characterized as a function of voltage for different solution conditions and membrane composition. The frequency of channel events was measured during the application of voltage ramps and counted in given intervals, a procedure that eliminated the effects of drift in gramicidin concentration. The formation rate was found to increase strongly with voltages up to approximately 50 mV and then to level off slightly. The shape of the voltage dependence was independent of lipid solvent and ramp speed but differed for different ions and different solution concentrations. This suggested an ion occupancy effect on the formation rate that was further supported by the fact that the minimum of the formation rate was shifted toward the equilibrium potential in asymmetric solution concentrations. The effects are explained in terms of a model that contains two contributions to the voltage dependence, a voltage-dependent ion binding to the monomers and a polarization of monomers by the applied electric field and by the occupied ions. The theory is found to give a good fit to experimental data.     

The effects of gramicidin on electroporation of lipid bilayers.

The effects of the channel-forming peptide gramicidin D (gD) on the conductance and electroporation thresholds of planar bilayer lipid membranes, made of the synthetic lipid 1-palmitoyl 2-oleoyl phosphatidylcholine (POPC), was studied. High-amplitude ( approximately 200-900 mV) rectangular voltage pulses of 15 ms duration were used to perturb the bilayers and monitor the transmembrane conductance. Electroporation voltage thresholds were found, and conductance was recorded before and after electroporation. Gramicidin was added to the system in peptide/lipid ratios of 1:10, 000, 1:500, and 1:15. The addition of gD in a ratio of 1:10,000 had no effect on electroporation, but ratios of 1:500 and 1:15 significantly increased the thresholds by 16% (p < 0.0001) and 40% (p < 0.0001), respectively. Membrane conductance before electroporation was measurable only after the addition of gD and increased monotonically as the peptide/lipid ratio increased. The effect of gD on the membrane area expansivity modulus (K) was tested using giant unilamellar vesicles (GUVs). When gD was incorporated into the vesicles in a 1:15 ratio, K increased by 110%, consistent with the increase in thresholds predicted by an electromechanical model. These findings suggest that the presence of membrane proteins may affect the electroporation of lipid bilayers by changing their mechanical properties.     

Small iminium ions block gramicidin channels in lipid bilayers.

Guanidinium and acetamidinium, when added to the bathing solution in concentrations of approximately 0.1M, cause brief blocks in the single channel potassium currents from channels formed in planar lipid bilayers by gramicidin A. Single channel lifetimes are not affected indicating that the channel structure is not modified by the blockers. Guanidinium block durations and interblock times are approximately exponential in distribution. Block frequencies increase with guanidinium concentration whereas block durations are unaffected. Increases in membrane potential cause an increase in block frequency as expected for a positively charged blocker but a decrease in block duration suggesting that the block is relieved when the blocker passes through the channel. At low pH, urea, formamide, and acetamide cause similar blocks suggesting that the protonated species of these molecules also block. Arginine and several amines do not block. This indicates that only iminium ions which are small enough to enter the channel can cause blocks in gramicidin channels.     

Transverse nuclear spin relaxation in phosphatidylcholine bilayers containing gramicidin.

Deuterium nuclear magnetic resonance has been used to study transverse relaxation in samples of 1,2-dimyristoyl-sn-glycero-3-phosphocholine, perdeuterated and specifically deuterated at the alpha position of the chains, containing the polypeptide gramicidin at concentrations of 0, 1, and 4 mol%. For 4 mol% gramicidin, the bilayer is thought to undergo a continuous phase change rather than a phase transition proceeding via two phase coexistence. Information is obtained regarding lipid dynamics in the continuous phase change region of the phase diagram. In the presence of gramicidin, the transverse relaxation time measured by the quadrupole echo technique, T2e, passes through a minimum in the gel phase. The gramicidin concentration dependence of T2e suggests that the polypeptide reduces the temperature sensitivity of the correlation time responsible for the minimum. The polypeptide also increases the sensitivity of the first spectral moment, M1, to the quadrupole echo pulse separation. This behavior is attributed to a polypeptide-induced enhancement of the spread in T2e along the acyl chains. Quadrupole Carr-Purcell-Meiboom-Gill experiments are used to separate contributions to the observed behavior from fast and slow motions.     

Solvent history dependence of gramicidin A conformations in hydrated lipid bilayers.

Reconstituted lipid bilayers of dimyristoylphosphatidylcholine (DMPC) and gramicidin A' have been prepared by cosolubilizing gramicidin and DMPC in one of three organic solvent systems followed by vacuum drying and hydration. The conformational state of gramicidin as characterized by 23Na NMR, circular dichroism, and solid state 15N NMR is dependent upon the cosolubilizing solvent system. In particular, two conformational states are described; a state in which Na+ has minimal interactions with the polypeptide, referred to as a nonchannel state, and a state in which Na+ interacts very strongly with the polypeptide, referred to as the channel state. Both of these conformations are intimately associated with the hydrophobic core of the lipid bilayer. Furthermore, both of these states are stable in the bilayer at neutral pH and at a temperature above the bilayer phase transition temperature. These results with gramicidin suggest that the conformation of membrane proteins may be dictated by the conformation before membrane insertion and may be dependent upon the mechanism by which the insertion is accomplished.     

Small angle x-ray scattering studies of magnetically oriented lipid bilayers.

Magnetically oriented lipid/detergent bilayers are potentially useful for studies of membrane-associated molecules and complexes using x-ray scattering and nuclear magnetic resonance (NMR). To establish whether the system is a reasonable model of a phospholipid bilayer, we have studied the system using x-ray solution scattering to determine the bilayer thickness, interparticle spacing, and orientational parameters for magnetically oriented lipid bilayers. The magnetically orientable samples contain the phospholipid L-alpha-dilauroylphosphatidylcholine (DLPC) and the bile salt analog 3-[(3-cholamidopropyl)dimethylammonio]-2-hydroxy-1-propanesulfonate (CHAPSO) in a 3:1 molar ratio in 70% water (w/v) and are similar to magnetically orientable samples used as NMR media for structural studies of membrane-associated molecules. A bilayer thickness of 30 A was determined for the DLPC/CHAPSO particles, which is the same as the bilayer thickness of pure DLPC vesicles, suggesting that the CHAPSO is not greatly perturbing the lipid bilayer. These data, as well as NMR data on molecules incorporated in the oriented lipid particles, are consistent with the sample consisting of reasonably homogeneous and well dispersed lipid particles. Finally, the orientational energy of the sample suggests that the size of the cooperatively orienting unit in the samples is 2 x 10(7) phospholipid molecules.     

Diffusion in lipid bilayers containing barriers

In epithelial cells, a barrier or tight junction restricts the diffusion of lipid probes from the apical to the basolateral side of the outer membrane bilayer. This phenomenon is studied theoretically with the diffusion equation on planar and spherical surfaces. Two models for the tight junction are considered: a penetrable barrier embedded in a monolayer and an impenetrable obstacle in the outer membrane of a bilayer than must be bypassed by flip-flopping between inner and outer membranes. The rate of passing from one side of the cell to the other is calculated for each of these models under steady state conditions. The results are compared with recent fluorescent photobleaching recovery experiments. The theoretical interpretation indicates that it would be difficult to distinguish experimentally between the flip-flop case and the barrier crossing case. Assuming a flip-flop model, large differences in the magnitude of the flip-flop rates of probes are necessary to explain the experimental results as suggested by Dragsten et al. (Dragsten, P. R., R. Blumenthal, and J. S. Handler, 1981, Nature [Lond.], 294:718--722).     

The alignment of a voltage-sensing peptide in dodecylphosphocholine micelles and in oriented lipid bilayers by nuclear magnetic resonance and molecular modeling.

The S4 segments of voltage-gated sodium channels are important parts of the voltage-sensing elements of these proteins. Furthermore, the addition of the isolated S4 polypeptide to planar lipid bilayers results in stepwise increases of ion conductivity. In order to gain insight into the mechanisms of pore formation by amphipathic peptides, the structure and orientation of the S4 segment of the first internal repeat of the rat brain II sodium channel was investigated in the presence of DPC micelles by multidimensional solution NMR spectroscopy and solid-state NMR spectroscopy on oriented phospholipid bilayers. Both the anisotropic chemical shift observed by proton-decoupled (15)N solid-state NMR spectroscopy and the attenuating effects of DOXYL-stearates on TOCSY crosspeak intensities of micelle-associated S4 indicate that the central alpha-helical portion of this peptide is oriented approximately parallel to the membrane surface. Simulated annealing and molecular dynamics calculations of the peptide in a biphasic tetrachloromethane-water environment indicate that the peptide alpha-helix extends over approximately 12 residues. A less regular structure further toward the C-terminus allows for the hydrophobic residues of this part of the peptide to be positioned in the tetrachloromethane environment. The implications for possible pore-forming mechanisms are discussed.     

Magnetic alignment and orientational order of dipalmitoylphosphatidylcholine bilayers containing palmitoyllysophosphatidylcholine

Mixed bilayers of 1-palmitoyl-sn-glycero-3-phosphocholine (palmitoyllysophosphatidylcholine; PaLPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (dipalmitoyl phosphatidylcholine; DPPC) have been investigated by 2H-NMR and 31P-NMR spectroscopy. Binary phospholipid mixtures were studied in which the acyl chains of one or the other component were perdeuterated. At temperatures below the main order-disorder phase transition, the mixed PaLPC/DPPC bilayers appear to coexist with PaLPC micelles. The micelles disappear at temperatures above the phase transition, where mixed bilayers in the liquid-crystalline state are formed. The orientational order of the alkyl chains of the PaLPC component is essentially identical to that of the DPPC component in the mixed bilayers, both in the low temperature and liquid-crystalline phases. However, the presence of PaLPC perturbs the segmental ordering of DPPC as compared to the pure system. The order is increased in the low-temperature phase, where effective diffusion of the chains about their long axes occurs, but is decreased in the liquid-crystalline phase compared to pure DPPC bilayers. The mixed liquid-crystalline bilayers orient preferentially with their director axes perpendicular to the magnetic field. This alignment is easily observed in 31P- and 2H-NMR spectra, where the intensity of the perpendicular edges of the lineshapes is pronounced. One possible explanation of the magnetic alignment involves alteration of the curvature free energy of the DPPC bilayer due to incorporation of PaLPC in the mixed membranes.     

Material property characteristics for lipid bilayers containing lysolipid.

The apparent area expansion modulus and tensile strength of egg phosphatidylcholine (EPC) membranes are measured in the presence of monooleoylphosphatidylcholine (MOPC). The apparent area expansion modulus decreases from 171 mN m-1 for pure EPC membrane to 82 mN m-1 for a membrane containing 30 mol % MOPC. This significant decrease of the apparent area expansion modulus is attributed to the change of the membrane area due to the tension-dependent exchange of MOPC between the bathing solution and the membrane. Similar to the apparent area expansion modulus, the tensile strength of the membrane decreases with the increase of the molar concentration of MOPC in the membrane. The tensile strength of pure EPC membrane is 9.4 mN m-1 whereas that for a membrane containing 30 mol % MOPC is only 1.8 mN m-1, and for a membrane containing 50 mol % MOPC it is even smaller, on the order of 0.07 mN m-1. The decrease of the tensile strength is coupled with a decrease of the work for membrane breakdown, which changes from 4.3 x 10(-2) kT for pure EPC membrane to 2 x 10(-6) kT for a membrane with 50 mol % MOPC. Overall, these results show that the decrease of the apparent area expansion modulus in the presence of exchangeable molecules is a fundamental property for all membranes and depends on the area occupied by these molecules. The method presented here provides a unique tool for measuring the area occupied by an exchangeable molecule in the bilayer membrane.     

Deuterium NMR of 2HCO-Val1...gramicidin A and 2HCO-Val1-D-Leu2...gramicidin A in oriented DMPC bilayers.

Deuterium NMR is used to study the structure and dynamics of the formyl C-2H bond in selectively deuterated gramicidin molecules. Specifically, the functionally different analogues 2HCO-Val1...gramicidin A and 2HCO-Val1-D-Leu2...gramicidin A are studied by 2H NMR so that any conformational or dynamical differences between the two analogues can be correlated with their difference in lifetime. These analogues are first synthesized, purified, and characterized and then incorporated into oriented bilayers of dimyristoylphosphatidylcholine sandwiched between glass coverslips. Phosphorous NMR line shapes obtained from these samples are consistent with the presence of the bilayer phase and indicate that the disorder exhibited by the lipid matrix is approximately of the same type and degree for both analogues. Deuterium NMR line shapes obtained from these samples indicate that the motional axis of the formyl group of gramicidin is parallel to the coverslip normal, that the distribution of motional axis orientations has a width of 7-9 degrees, and that a similar, major conformational and dynamical state exists for the formyl C-2H bond of both analogues. In this state, if the only motion present is fast axial rotation, then the experimentally derived angle between the formyl C-2H bond and the motional axis is consistent with the presence of a right-handed, single-stranded, beta 6.3 helical dimer but is not consistent with the presence of a left-handed, single-stranded, beta 6.3 helical dimer. However, if fast axial rotation is not the only motion present, then the left-handed, single-stranded, beta 6.3 helical dimer cannot be absolutely excluded as a possibility. Also, a second, minor conformational and dynamical state appears to be present in the spectrum of 2HCO-Val1-D-Leu2...gramicidin A but is not observed in the spectrum of 2HCO-Val1...gramicidin A. This minor conformational and dynamical state may reflect the presence of monomers, while the major conformational and dynamical state may reflect the presence of dimers.     

Proton conduction in gramicidin A and in its dioxolane-linked dimer in different lipid bilayers.

Gramicidin A (gA) molecules were covalently linked with a dioxolane ring. Dioxolane-linked gA dimers formed ion channels, selective for monovalent cations, in planar lipid bilayers. The main goal of this study was to compare the functional single ion channel properties of natural gA and its covalently linked dimer in two different lipid bilayers and HCl concentrations (10-8000 mM). Two ion channels with different gating and conductance properties were identified in bilayers from the product of dimerization reaction. The most commonly observed and most stable gramicidin A dimer is the main object of this study. This gramicidin dimer remained in the open state most of the time, with brief closing flickers (tau(closed) approximately 30 micros). The frequency of closing flickers increased with transmembrane potential, making the mean open time moderately voltage dependent (tau(open) changed approximately 1.43-fold/100 mV). Such gating behavior is markedly different from what is seen in natural gA channels. In PEPC (phosphatidylethanolamine-phosphatidylcholine) bilayers, single-channel current-voltage relationships had an ohmic behavior at low voltages, and a marked sublinearity at relatively higher voltages. This behavior contrasts with what was previously described in GMO (glycerylmonooleate) bilayers. In PEPC bilayers, the linear conductance of single-channel proton currents at different proton concentrations was essentially the same for both natural and gA dimers. g(max) and K(D), obtained from fitting experimental points to a Langmuir adsorption isotherm, were approximately 1500 pS and 300 mM, respectively, for both the natural gA and its dimer. In GMO bilayers, however, proton affinities of gA and the dioxolane-dimer were significantly lower (K(D) of approximately 1 and 1.5 M, respectively), and the g(max) higher (approximately 1750 and 2150 pS, respectively) than in PEPC bilayers. Furthermore, the relationship between single-channel conductance and proton concentration was linear at low bulk concentrations of H+ (0.01-2 M) and saturated at concentrations of more than 3 M. It is concluded that 1) The mobility of protons in gramicidin A channels in different lipid bilayers is remarkably similar to proton mobilities in aqueous solutions. In particular, at high concentrations of HCl, proton mobilities in gramicidin A channel and in solution differ by only 25%. 2) Differences between proton conductances in gramicidin A channels in GMO and PEPC cannot be explained by surface charge effects on PEPC membranes. It is proposed that protonated phospholipids adjacent to the mouth of the pore act as an additional source of protons for conduction through gA channels in relation to GMO bilayers. 3) Some experimental results cannot be reconciled with simple alterations in access resistance to proton flow in gA channels. Said differences could be explained if the structure and/or dynamics of water molecules inside gramicidin A channels is modulated by the lipid environment and by modifications in the structure of gA channels. 4) The dioxolane ring is probably responsible for the closing flickers seen in the dimer channel. However, other factors can also influence closing flickers.     

Calculation of deformation energies and conformations in lipid membranes containing gramicidin channels.

In this paper we calculate surface conformation and deformation free energy associated with the incorporation of gramicidin channels into phospholipid bilayer membranes. Two types of membranes are considered. One is a relatively thin solvent-free membrane. The other is a thicker solvent-containing membrane. We follow the approach used for the thin membrane case by Huang (1986) in that we use smectic liquid crystal theory to evaluate the free energy associated with distorting the membrane to other than a flat configuration. Our approach is different from Huang, however, in two ways. One is that we include a term for surface tension, which Huang did not. The second is that one of our four boundary conditions for solving the fourth-order differential equation describing the free energy of the surface is different from Huang's. The details of the difference are described in the text. Our results confirm that for thin membranes Huang's neglect of surface tension is appropriate. However, the precise geometrical form that we calculate for the surface of the thin membrane in the region of the gramicidin channel is somewhat different from his. For thicker membranes that have to deform to a greater extent to accommodate the channel, we find that the contribution of surface tension to the total energy in the deformed surface is significant. Computed results for the shape of the deformed surface, the total energy in the deformed surface, and the contributions of different components to the total energy, are presented for the two types of membranes considered. These results may be significant for understanding the mechanisms of dimer formation and breakup, and the access resistance for ions entering gramicidin channels.     

Solid-state 15N NMR of oriented lipid bilayer bound gramicidin A'

Highly oriented samples of lipid and gramicidin A' (8:1 molar ratio) have been prepared with the samples extensively hydrated (approximately 70% water v/w). These preparations have been shown to be completely in a bilayer phase with a transition temperature of 28 degrees C, and evidence is presented indicating that the gramicidin is in the channel conformation. An estimate of the disorder in the alignment of the bilayers parallel with the glass plates used to align the bilayers can be made from the asymmetry of the nuclear magnetic resonances (NMR). Such an analysis indicates a maximal range of disorder of +/- 3 degrees. Uniformly 15N-labeled gramicidin has been biosynthesized by Bacillus brevis grown in a media containing 15N-labeled Escherichia coli cells as the only nitrogen source. When prepared with labeled gramicidin, the oriented samples result in high-resolution 15N NMR spectra showing 12 resonances for the 20 nitrogen sites of the polypeptide. The frequency of the three major multiple resonance peaks has been interpreted to yield the approximate orientation of the N-H bonds in the peptide linkages with respect to the magnetic field. These bond orientations are only partially consistent with the extant structural models of gramicidin.     

Calcium-dependent association of annexins with lipid bilayers modifies gramicidin A channel parameters

In order to examine whether calcium-dependent binding of annexin to acidic phospholipids could change the lipid bilayer environment sufficiently to perturb channel-mediated transmembrane ion-transport, gramicidin A channel activity in planar lipid bilayers was investigated in the presence of calcium and annexins II, III or V. The experiments were performed with membranes consisting of phosphatidylcholine, phosphatidylethanolamine and phosphatidylserine in 300 mM KCl solution buffered to pH 7.4 and with either 0.1 or 1 mM calcium added to the solution. Annexin (1 microM) was subsequently applied to the cis side of the membrane. All three annexins (II, III and V) when tested at 1 mM calcium decreased the gramicidin single-channel conductance. Annexins II and III increased the mean lifetime of the channels whereas annexin V seemed to have no influence on the mean lifetime. Since the lifetime of gramicidin A channels is a function of the rate constant for dissociation of the gramicidin dimer, which is dependent on the physical properties of the lipid phase, binding of annexins II and III seems to stabilize the gramicidin channel owing to a change of the bilayer structure.     

Phase and surface properties of lipid bilayers containing neoglycolipids

The physical properties conferred to DPPC bilayers by including neoglycolipids composed by two different trisaccharides: mannose-mannose-mannose (3M) and glucose-mannose-glucose (GMG) attached to a cholesterol (cho) and a distearylglycerol (diC18) lipid moiety by a spacer were evaluated by means of the measurement of the electrokinetic potential and interfacial fluorescent probes. The phase properties measured with diphenylhexatriene (DPH) were correlated with the surface properties measured with merocyanine 540, dansyl, and Laurdan probes. The results show that the surface properties of large unilamellar vesicles depend on the sugar exposure to the water phase and also on the hydrocarbon moiety by which it is anchored to the bilayer. The combination of the cholesterol moiety with the saccharide attenuates the cooperativity decrease induced by the cholesterol moiety without the sugar portion. The neoglycolipid GMG-diC18 promotes opposite effects affecting slightly the cooperativity at the hydrocarbon core of DPPC and displacing the phase transition temperature to higher values. The presence of neoglycolipid with diC18 introduces defects in the packing at the interface of the membrane in the gel state. It is concluded that a relatively low proportion of neoglycolipids affects significantly the interfacial properties of DPPC bilayers in large unilamellar vesicles in the absence of changes at the membrane bulk at 25 degrees C.     

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers

Gramicidin A (gA) is a 15-amino-acid antibiotic peptide with an alternating L-D sequence, which forms (dimeric) bilayer-spanning, monovalent cation channels in biological membranes and synthetic bilayers. We performed molecular dynamics simulations of gA dimers and monomers in all-atom, explicit dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dioleoylphosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers. The variation in acyl chain length among these different phospholipids provides a way to alter gA-bilayer interactions by varying the bilayer hydrophobic thickness, and to determine the influence of hydrophobic mismatch on the structure and dynamics of both gA channels (and monomeric subunits) and the host bilayers. The simulations show that the channel structure varied little with changes in hydrophobic mismatch, and that the lipid bilayer adapts to the bilayer-spanning channel to minimize the exposure of hydrophobic residues. The bilayer thickness, however, did not vary monotonically as a function of radial distance from the channel. In all simulations, there was an initial decrease in thickness within 4–5 Å from the channel, which was followed by an increase in DOPC and POPC or a further decrease in DLPC and DMPC bilayers. The bilayer thickness varied little in the monomer simulations—except one of three independent simulations for DMPC and all three DLPC simulations, where the bilayer thinned to allow a single subunit to form a bilayer-spanning water-permeable pore. The radial dependence of local lipid area and bilayer compressibility is also nonmonotonic in the first shell around gA dimers due to gA-phospholipid interactions and the hydrophobic mismatch. Order parameters, acyl chain dynamics, and diffusion constants also differ between the lipids in the first shell and the bulk. The lipid behaviors in the first shell around gA dimers are more complex than predicted from a simple mismatch model, which has implications for understanding the energetics of membrane protein-lipid interactions.