Interaction of an Amphiphilic Peptide with a Phospholipid Bilayer Surface by Molecular Dynamics Simulation Study

Abstract

Corticotropin-releasing factor (CRF) is the principal neuroregulator of adrenocorticotropic hormone (ACTH) secretion. Previous experiments have demonstrated that CRF binds avidly to the surface of single egg phosphatidylcholine vesicles and its amphiphilic secondary structure might play an important role in the function. In this study, the interaction of the residues 13–41 in human CRF with the surface of a DOPC bilayer was investigated by molecular dynamics (MD) simulation in order to understand the role of the membrane surface in the formation of the amphiphilic α helix as well as to determine the effects of the peptide on the lipid bilayer. The model used included 60 DOPC molecules, 1 helical peptide (CRF_13–41) on the bilayer surface, and explicit waters of solvation in the lipid polar head group regions, together with constant-volume periodic boundary conditions in three dimensions. The MD simulation was carried out for 510 ps. In addition, CRF_13–41, initially in a helical form, was simulated in vacuo as a control. The results indicate that while it was completely unstable in vacuo, the peptide helical form was generally maintained on the bilayer surface, but with distortions near the terminal ends. The peptide was confined to the bilayer headgroup/water region, similar to that reported from neutron diffraction measurement of tripeptides bound to the phosphatidylcholine bilayer surface (Ref 1). The amphiphilicity of the peptide matched that of the bilayer headgroup environment, with the hydrophilic side oriented toward water and the hydrophobic side making contact with the bilayer hydrocarbon core. These results support the hypothesis that the amphiphilic environment of a membrane surface is important in the induction of peptide amphiphilic α-helical secondary structure. Two major effects of the peptide on the lipids were found: the first CH₂ segment in the lipid chains was significantly disordered and the lipid headgroup distribution was broadened towards the water region.     

ESR and monolayer study of the localization of coenzyme Q10 in artificial membranes

The data obtained from the ESR experiments show a complex, depth dependent effect of CoQ10 on the lipid molecules mobility in the bilayer. These effects depend both on its concentration and the temperature. CoQ10 disturbs not only the hydrophobic core of the membrane but also the region close to the hydrophilic headgroups of phospholipids. Both these effects could be explained by the fact that the high hydrophobicity of CoQ10 causes the molecules to position itself in the interior of the bilayer, but at the same time its water seeking headgroup is located close to the region of the polar headgrops of membrane lipids. The presence of CoQ10 in the hydrophobic core has further implications on the properties of membrane intrinsic domain. Results of monolayer experiments indicate that CoQ10 may form aggregates when mixed with PC molecules in the lipid hydrocarbon chain-length dependent manner. CoQ10 is not fully miscible with DMPC or DPPC but it is well miscible with the long-chain DSPC molecules. Our suggestion is that CoQ10 when present in long-chain phospholipid bilayer, interacts with saturated fatty acyl-chains and adapt the structure which allows such interactions: either parallel to the saturated acyl chains or "pseudo-ring" conformation resembling sterol structure.     

Comparative hydrocarbon utilization by hydrophobic and hydrophilic variants of Pseudomonas aeruginosa

Aims: To investigate hydrocarbon degradation by hydrophobic, hydrophilic and parental strains of Pseudomonas aeruginosa. Methods and Results: Partitioning of hydrocarbon‐degrading P. aeruginosa strain in a solvent/aqueous system yielded hydrophobic and hydrophilic fractions. Exhaustive partitioning of aqueous‐phase cells yielded the hydrophilic variants (L), while sequential fractionation of the hydrophobic phase cells yielded successive fractions exhibiting increasing cell‐surface hydrophobicity (CSH). In hydrocarbon adherence assays (bacterial attachment to hydrocarbon), L had a value of 20%, which increased from 61·7% in first hydrophobic fraction (H₁) to 72·2% in the third (H₃). Crude oil degradation by L was 70%, but increased from 82% in H₁ to 93% in H₃. L variant produced most exopolysaccharides and reduced surface tension from about 73 to 49 mN m^?1. Rhamnolipid production was highest in L, but was not detected in all crude oil cultures. Conclusions: Hydrophobic subpopulations of hydrocarbon‐degrading P. aeruginosa exhibited greater hydrocarbon‐utilizing ability than hydrophilic ones, or the parental strain. Significance and Impact of the Study: Results demonstrate that a population of P. aeruginosa consists of cells with different CSH which affect hydrocarbon utilization. This potentially provides the population with the capacity to utilize different hydrophobic substrates found in petroleum. Judicious selection of such hydrophobic subpopulations can enhance hydrocarbon pollution bioremediation.     

Membrane hydration and structure on a subnanometer scale as seen by high resolution solid state nuclear magnetic resonance: POPC and POPC/C12EO4 model membranes.

The position on a subnanometer scale and the dynamics of structurally important water in model membranes was determined using a combination of proton magic-angle spinning NMR (MAS) with two-dimensional NOESY NMR techniques. Here, we report studies on phosphocholine lipid bilayers that were then modified by the addition of a nonionic surfactant that is shown to dehydrate the lipid. These studies are supplemented by 13C magic-angle spinning NMR investigations to get information on the dynamics of segmental motions of the membrane molecules. It can be shown that the hydrophilic chain of the surfactant is positioned at least partially within the hydrophobic core of the lipid bilayer. With the above NMR approach, we are able to establish molecular contacts between water and the lipid headgroup as well as with certain groups of the hydrocarbon chains and the glycerol backbone. This is possible because high resolution proton and 13C-NMR spectra of multilamellar bilayer membranes are obtained using MAS. A phase-sensitive NOESY must also be applied to distinguish positive and negative cross-peaks in the two-dimensional plot. These studies have high potential to investigate membrane proteins hydration and structural organization in a natural lipid bilayer surrounding.     

Dependence of lipid chain and head group packing of the inverted hexagonal phase on hydration

A model which positions the hydrophobic/hydrophilic boundary in phosphatidylethanolamine lipids at the first CH₂ group in the acyl or alkyl chain is used to calculate the surface area per lipid, the mean chain and head-group dimensions and diameters of the hydrophilic tubes of the inverted hexagonal phase of didodecylphosphatidylethanolamine. The calculated surface areas compare favorably with areas obtained for the lamellar liquid crystal phase of the same lipid using the same boundary. Placement of the boundary within the lipid structure permits a determination of the maximum headgroup packing at hydration levels down to complete dehydration. The headgroup dimensions are consistent with a 5 Å diam void at the center of a hydrophilic tube at zero hydration. The calculated mean fluid chain length is ~2 Å smaller than the mean chain length of the lamellar phase at comparable levels of hydration. Comparison of the calculated mean fluid chain length and distances between hydrophobic boundaries shows that the fluid chains are interdigitated between adjacent tubes, and not interdigitated in the central space between three tubes. At low hydration the chains interdigitate in both spaces. The number of lipids packed around a tube at low hydration is only a function of the headgroup geometry, whereas at high hydration, it is a function of the number of carbon atoms in the chains.     

Effects of n-alkanes on the morphology of lipid bilayers A freeze-fracture and negative stain analysis

The effect of $\text{n-}\text{alkanes}$ on the ultrastructure of lipid bilayers has been investigated using freeze-fracture and negative stain electron microscopy. It has been found that the morphology of bilayers containing the long alkane tetradecane is quite different from bilayers containing the short alkane hexane. The smooth fracture faces of gel and liquid crystalline state bilayers are unmodified by tetradecane. However, hexane dramatically alters the hydrophobic bilayer interior, producing large (20 to 50 nm) mounds and depressions in the fracture faces. The fracture steps in these multilayer preparations containing hexane are variable in thickness and often considerably wider than the corresponding fracture steps in multilayers which contain tetradecane or are solvent-free. Alkanes also modify the structure of the P_β′ or ‘banded’ phase of phosphatidylcholine bilayers. The incorporation of tetradecane removes the banded structure from both the bilayer's hydrophilic surface, as viewed by negative staining, and the bilayer's hydrophobic interior, as viewed by the freeze-fracture technique. These results are consistent with X-ray diffraction data which imply that long alkanes are primarily located between adjacent lipid hydrocarbon chains in each monolayer of the bilayer, while short alkanes can partition into the geometric center of the bilayer between apposing monolayers.     

Constant pressure and temperature molecular dynamics simulation of a fully hydrated liquid crystal phase dipalmitoylphosphatidylcholine bilayer.

We report a constant pressure and temperature molecular dynamics simulation of a fully hydrated liquid crystal (L alpha) phase bilayer of dipalmitoylphosphatidylcholine at 50 degrees C and 28 water molecules/lipid. We have shown that the bilayer is stable throughout the 1550-ps simulation and have demonstrated convergence of the system dimensions. Several important aspects of the bilayer structure have been investigated and compared favorably with experimental results. For example, the average positions of specific carbon atoms along the bilayer normal agree well with neutron diffraction data, and the electron density profile is in accord with x-ray diffraction results. The hydrocarbon chain deuterium order parameters agree reasonably well with NMR results for the middles of the chains, but the simulation predicts too much order at the chain ends. In spite of the deviations in the order parameters, the hydrocarbon chain packing density appears to be essentially correct, inasmuch as the area/lipid and bilayer thickness are in agreement with the most refined experimental estimates. The deuterium order parameters for the glycerol and choline groups, as well as the phosphorus chemical shift anisotropy, are in qualitative agreement with those extracted from NMR measurements.     

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.     

Structure and dynamics of model pore insertion into a membrane

A cylindrical transmembrane molecule is constructed by linking hydrophobic sites selected from a coarse grain model. The resulting hollow tube assembly serves as a representation of a transmembrane channel, pore, or a carbon nanotube. The interactions of a coarse grain di-myristoyl-phosphatidyl-choline hydrated bilayer with both a purely hydrophobic tube and a tube with hydrophilic caps are studied. The hydrophobic tube rotates in the membrane and becomes blocked by lipid tails after a few tens of nanoseconds. The hydrophilic sites of the capped tube stabilize it by anchoring the tube in the lipid headgroup/water interfacial region of each membrane leaflet. The capped tube remains free of lipid tails. The capped tube spontaneously conducts coarse grain water sites; the free-energy profile of this process is calculated using three different methods and is compared to the barrier for water permeation through the lipid bilayer. Spontaneous tube insertion into an undisturbed lipid bilayer is also studied, which we reported briefly in a previous publication. The hydrophobic tube submerges into the membrane core in a carpetlike manner. The capped tube laterally fuses with the closest leaflet, and then, after plunging into the membrane interior, rotates to assume a transbilayer orientation. Two lipids become trapped at the end of the tube as it penetrates the membrane. The hydrophilic headgroups of these lipids associate with the lower tube cap and assist the tube in crossing the interior of the membrane. When the rotation is complete these lipids detach from the tube caps and fuse with the lower leaflet lipids.     

Bilayer conformation of fusion peptide of influenza virus hemagglutinin: a molecular dynamics simulation study

Unraveling the conformation of membrane-bound viral fusion peptides is essential for understanding how those peptides destabilize the bilayer topology of lipids that is important for virus-cell membrane fusion. Here, molecular dynamics (MD) simulations were performed to investigate the conformation of the 20 amino acids long fusion peptide of influenza hemagglutinin of strain X31 bound to a dimyristoyl phosphatidylcholine (DMPC) bilayer. The simulations revealed that the peptide adopts a kinked conformation, in agreement with the NMR structures of a related peptide in detergent micelles. The peptide is located at the amphipathic interface between the headgroups and hydrocarbon chains of the lipid by an energetically favorable arrangement: The hydrophobic side chains of the peptides are embedded into the hydrophobic region and the hydrophilic side chains are in the headgroup region. The N-terminus of the peptide is localized close to the amphipathic interface. The molecular dynamics simulations also revealed that the peptide affects the surrounding bilayer structure. The average hydrophobic thickness of the lipid phase close to the N-terminus is reduced in comparison with the average hydrophobic thickness of a pure dimyristoyl phosphatidylcholine bilayer.     

Preparation and structural studies of cholesterol bilayers

Bilayers consisting, in their hydrophobic core, entirely of cholesterol can be constructed if a hydrophilic molecular anchor is supplied. $\text{O-}\text{Methoxyethoxyethoxyethylcholesterol}$ and cholesterol sulfate form multilayered liposomes in water. With equimolar cholesterol added, cholesterol sulfate, cholesterolphosphocholine, and $\text{O-}\text{methoxyethoxyethoxyethylcholesterol}$ form small unilamellar liposomes on prolonged sonication. The dimensions of cholesterol-cholesterolphosphocholine vesicles are comparable to those of phospholipid vesicles. ¹³C-NMR spectra suggest that the centers of the bilayers are liquid. The permeability of the cholesterol-cholesterolphosphocholine bilayer against glycerol is lower than that of dipalmitoylphosphatidylcholine-cholesterol bilayer; the activation energy of permeation is two times larger, an indication of a higher degree of structural organization in the ‘hydrogen belts’ of the cholesterol-cholesterolphosphocholine bilayer.     

Molecular Dynamics Simulations of Phospholipid Bilayers

Abstract

Molecular dynamics (MD) simulations at 37°C have been performed on three phospholipid bilayer systems composed of the lipids DLPE, DOPE, and DOPC. The model used included 24 explicit lipid molecules and explicit waters of solvation in the polar head group regions, together with constant-pressure periodic boundary conditions in three dimensions. Using this model, a MD simulation samples part of an infinite planar lipid bilayer. The lipid dynamics and packing behavior were characterized. Furthermore, using the results of the simulations, a number of diverse properties including bilayer structural parameters, hydrocarbon chain order parameters, dihedral conformations, electron density profile, hydration per lipid, and water distribution along the bilayer normal were calculated. Many of these properties are available for the three lipid systems chosen, making them well suited for evaluating the model and protocols used in these simulations by direct comparisons with experimental data. The calculated MD behavior, chain disorder, and lipid packing parameter, i.e. the ratio of the effective areas of hydrocarbon tails and head group per lipid (a_t/a_h), correctly predict the aggregation preferences of the three lipids observed experimentally at 37°C, namely: a gel bilayer for DLPE, a hexagonal tube for DOPE, and a liquid crystalline bilayer for DOPC. In addition, the model and conditions used in the MD simulations led to good agreement of the calculated properties of the bilayers with available experimental results, demonstrating the reliability of the simulations. The effects of the cis unsaturation in the hydrocarbon chains of DOPE and DOPC, compared to the fully saturated one in DLPE, as well as the effects of the different polar head groups of PC and PE with the same unsaturated chains on the lipid packing and bilayer structure have been investigated. The results of these studies indicate the ability of MD methods to provide molecular-level insights into the structure and dynamics of lipid assemblies.     

Anchoring of a monotopic membrane protein: the binding of prostaglandin H2 synthase-1 to the surface of a phospholipid bilayer

Prostaglandin H₂ synthases (PGHS-1 and -2) are monotopic peripheral membrane proteins that catalyse the synthesis of prostaglandins in the arachidonate cascade. Picot et al. (1994) proposed that the enzyme is anchored to one leaflet of the bilayer by a membrane anchoring domain consisting of a right-handed spiral of amphipathic helices (residues 73–116) forming a planar motif. Two different computational approaches are used to examine the association of the PGHS-1 membrane anchoring domain with a membrane via the proposed mechanism. The electrostatic contribution to the free energy of solvation is obtained by solving numerically the finite-difference Poisson equation for the protein attached to a membrane represented as a planar slab of low dielectric. The nonpolar cavity formation and van der Waals contributions to the solvation free energy are assumed to be proportional to the water accessible surface area. Based on the optimum position determined from the continuum solvent model, two atomic models of the PGHS-1 anchoring domain associated with an explicit dimyristoylphosphatidylcholine (DMPC) bilayer differing by the thickness of the membrane bilayer were constructed. A total of 2 ns molecular dynamics simulation were performed to study the details of lipid- protein interactions at the microscopic level. In the simulations the lipid hydrocarbon chains interacting with the anchoring domain assume various shapes, suggesting that the plasticity of the membrane is significant. The hydrophobic residues in the membrane side of the helices interact with the hydrophobic membrane core, while the positively charged residues interact with the lipid polar headgroups to stabilize the anchoring of the membrane domain to the upper half of the bilayer. The phosphate headgroup of one DMPC molecule disposed at the center of the spiral formed by helices A, B, C and D interacts strongly with Arg120, a residue on helix D that has previously been identified as being important in the activity of PGHS-1. In the full enzyme structure, this position corresponds to the entrance of a long hydrophobic channel leading to the cyclooxygenase active site. These observations provide insights into the association of the arachidonic acid substrate to the cyclooxygenase active site of PGHS-1. Received: 20 December 1999 / Revised version: 26 March 2000 / Accepted: 26 March 2000     

Structure and interactive properties of highly fluorinated phospholipid bilayers.

Because liposomes containing fluoroalkylated phospholipids are being developed for in vivo drug delivery, the structure and interactive properties of several fluoroalkylated glycerophosphocholines (PCs) were investigated by x-ray diffraction/osmotic stress, dipole potential, and hydrophobic ion binding measurements. The lipids included PCs with highly fluorinated tails on both alkyl chains and PCs with one hydrocarbon chain and one fluoroalkylated chain. Electron density profiles showed high electron density peaks in the center of the bilayer corresponding to the fluorine atoms. The height and width of these high density peaks varied systematically, depending on the number of fluorines and their position on the alkyl chains, and on whether the bilayer was in the gel or liquid crystalline phase. Wide-angle diffraction showed that in both gel and liquid crystalline bilayers the distance between adjacent alkyl chains was greater in fluoroalkylated PCs than in analogous hydrocarbon PCs. For interbilayer separations of less than about 8 A, pressure-distance relations for fluoroalkylated PCs were similar to those previously obtained from PC bilayers with hydrocarbon chains. However, for bilayer separations greater than 8A, the total repulsive pressure depended on whether the fluoroalkylated PC was in a gel or liquid-crystalline phase. We argue that these pressure-distance relations contain contributions from both hydration and entropic repulsive pressures. Dipole potentials ranged from -680 mV for PCs with both chains fluoroalkylated to -180 mV for PCs with one chain fluoroalkylated, compared to +415 mV for egg PC. The change in dipole potential as a function of subphase concentration of tetraphenyl-boron was much larger for egg PC than for fluorinated PC monolayers, indicating that the fluorine atoms modified the binding of this hydrophobic anion. Thus, compared to conventional liposomes, liposomes made from fluoroalkylated PCs have different binding properties, which may be relevant to their use as drug carriers.     

Differences in the interaction of inorganic and organic (hydrophobic) cations with phosphatidylserine membranes

The interaction of phosphatidylserine dispersions with “hydrophobic”, organic cations (acetylcholine, tetraethylammonium ion) is compared with that of simple inorganic cations (Na⁺, Ca²⁺); differences in the hydration properties of the two classes of ions exist in the bulk phase as evident from spin-lattice relaxation time T₁, measurements. It is shown that the reaction products (cation-phospholipid) differ markedly in their physicochemical behaviour. With increasing concentration both classes of ions reduce the ζ-potential of phosphatidylserine surfaces, the monovalent inorganic cations being only slightly more effective than the hydrophobic cations. Inorganic cations cause precipitation of the lipid once the surface charge of the bilayer is reduced to a certain threshold value. This is not the case with the organic cations. The difference is probably associated with the different hydration properties of the resulting complexes. Thus binding of Ca²⁺ causes displacement of water of hydration and formation of an anhydrous, hydrophobic calcium-phosphatidylserine complex which is insoluble in water, whereas the product of binding of the organic cations is hydrated, hydrophilic and water soluble. The above findings are consistent with NMR results which show that the phosphodiester group is involved in the binding of both classes of cations as well as being the site of the primary hydration shell. Besides affecting interbilayer membrane interactions such as those involved in cell adhesion and membrane fusion, the binding of both classes of cation can affect the molecular packing within a bilayer.     

The nature of the hydrophobic binding of small peptides at the bilayer interface: implications for the insertion of transbilayer helices

One method of obtaining useful information about the physical chemistry of peptide/bilayer interactions is to relate thermodynamic parameters of the interactions to structural parameters obtained by diffraction methods. We report here the results of the application of this approach to interactions of hydrophobic tripeptides of the form Ala-X-Ala-O-tert-butyl with lipid bilayers. The thermodynamic constants (delta Gt, delta Ht, and delta St) for the transfer of the tripeptides from water into DMPC vesicles were determined for X = Leu, Phe, and Trp and found to be consistent with those expected for hydrophobic interactions above the phase transition of DMPC. Combining these results with the earlier ones of Jacobs and White [(1986) Biochemistry 25, 2605-2612], the favorable free energies of transfer with different amino acids in the -X- position increase in the order Gly less than Ala less than Leu less than Phe less than Trp in agreement with the Nozaki and Tanford [(1971) J. Biol. Chem. 246, 2211-2217] hydrophobicity scale. Determination of the location of Ala-[2H5]Trp-Ala-O-tert-butyl in oriented DOPC bilayers by neutron diffraction shows that the most hydrophobic peptide of the series is confined to the bilayer headgroup/water region. Refinement of the diffraction measurements shows that only 13% of the tryptophan is associated with the hydrocarbon core. The distribution of the water tends to mirror that of the peptide. Unlike peptide-free bilayers, 5% of the water penetrates the hydrocarbon, which is about 100-fold greater than expected. A quantitative thermodynamic analysis of the interfacial binding of the peptides suggests that (1) the hydrophobic interactions are 60-70% complete upon binding at the bilayer interface, (2) the interface is likely to play an important role in helix formation and insertion, (3) the hydrogen bond status of amino acid side chains is crucial to insertion, and (4) an a priori lack of knowledge of the status of such bonds could limit the precision of hydrophobicity plots. We introduce an interfacial hydrophobicity scale, IFH(h), with a variable hydrogen bond parameter (h) that permits one to consider explicitly hydrogen bonding in transbilayer helix searches.     

Structural basis for carbon dioxide binding by 2-ketopropyl coenzyme M oxidoreductase/carboxylase

The structure of 2-ketopropyl coenzyme M oxidoreductase/carboxylase (2-KPCC) has been determined in a state in which CO₂ is observed providing insights into the mechanism of carboxylation. In the substrate encapsulated state of the enzyme, CO₂ is bound at the base of a narrow hydrophobic substrate access channel. The base of the channel is demarcated by a transition from a hydrophobic to hydrophilic environment where CO₂ is located in position for attack on the carbanion of the ketopropyl group of the substrate to ultimately produce acetoacetate. This binding mode effectively discriminates against H₂O and prevents protonation of the ketopropyl leaving group.

Structured summary

2-KPCCbinds to 2-KPCC by x-ray crystallography (View interaction)     

On the use of deuterium nuclear magnetic resonance as a probe of chain packing in lipid bilayers

The packing of hydrocarbon chains in the bilayers of lamellar (L alpha) phases of soap/water and phospholipid/water mixtures has been studied by deuterium NMR spectroscopy and X-ray diffraction. A universal correlation is shown to exist between the average C-D bond order parameter SCD of hydrocarbon chains and the average area per chain ach, irrespective of the chemical structure of the surfactant (hydrophilic group, number of chains per molecule, and chain length), composition, and temperature. The practical utility of the correlation is illustrated by its application to the characterization of the distribution of various hydrophobic and amphiphilic solutes in bilayers. The distribution of hydrocarbons within a bilayer is shown to depend upon their molecular structure in a manner which highlights the nature of the molecular interactions involved. For example, benzene is shown to be fairly uniformly distributed across the bilayer with an increasing tendency to distribute into the center at high concentrations. In contrast, the more complex hydrocarbon tetradecane preferentially distributes into the center of the bilayer at low concentrations, while at higher concentrations it intercalates between the surfactant chains. Alcohols such as benzyl alcohol, octanol, and decanol all interact similarly with the bilayer in so far as they are pinned to the polar/apolar interface, presumably by involvement of the hydroxyl group in a hydrogen bond. But the response of the surfactant chains to the void volume created in the center of the bilayer is dependent upon the distance of penetration of the alcohol into the bilayer. For benzyl alcohol, the shortest molecule, this void volume is taken up by the disordering of the chains, while for decanol, the longest molecule, it is absorbed by interdigitation of the chains of apposing monolayers. For octanol, the chain interdigitation mechanism is dominant at low concentrations, but there is a transition to chain disordering at high concentrations. Finally, it is shown that the correlation provides a useful test for statistical mechanical models of chain ordering in lipid bilayers.     

Indole localization in lipid membranes revealed by molecular simulation

It is commonly known that the amino acid residue tryptophan and its side-chain analogs, e.g., indole, are strongly attracted to the interfacial region of lipid bilayers. Phenylalanine and its side-chain analogs, e.g., benzene, do not localize in the interface but are distributed throughout the lipid bilayer. We use molecular dynamics to investigate the details of indole and benzene localization and orientation within a POPC bilayer and the factors that lead to their different properties. We identify three sites in the bilayer at which indole is localized: 1), a site in the interface near the glycerol moiety; 2), a weakly bound site in the interface near the choline moiety; and 3), a weakly bound site in the center of the bilayer's hydrocarbon core. Benzene is localized in the same three positions, but the most stable position is the hydrocarbon core followed by the site near the glycerol moiety. Transfer of indole from water to the hydrocarbon core shows a classic hydrophobic effect. In contrast, interfacial binding is strongly enthalpy driven. We use several different sets of partial charges to investigate the factors that contribute to indole's and benzene's orientational and spatial distribution. Our simulations show that a number of electrostatic interactions appear to contribute to localization, including hydrogen bonding to the lipid carbonyl groups, cation-pi interactions, interactions between the indole dipole and the lipid bilayer's strong interfacial electric field, and nonspecific electrostatic stabilization due to a mismatch in the variation of the nonpolar forces and local dielectric with position in the bilayer.     

A theoretical model for the association of amphiphilic transmembrane peptides in lipid bilayers

A theoretical model is proposed for the association of trans-bilayer peptides in lipid bilayers. The model is based on a lattice model for the pure lipid bilayer, which accounts accurately for the most important conformational states of the lipids and their mutual interactions and statistics. Within the lattice formulation the bilayer is formed by two independent monolayers, each represented by a triangular lattice, on which sites the lipid chains are arrayed. The peptides are represented by regular objects, with no internal flexibility, and with a projected area on the bilayer plane corresponding to a hexagon with seven lattice sites. In addition, it is assumed that each peptide surface at the interface with the lipid chains is partially hydrophilic, and therefore interacts with the surrounding lipid matrix via selective anisotropic forces. The peptides would therefore assemble in order to shield their hydrophilic residues from the hydrophobic surroundings. The model describes the self-association of peptides in lipid bilayers via lateral and rotational diffusion, anisotropic lipid-peptide interactions, and peptide-peptide interactions involving the peptide hydrophilic regions. The intent of this model study is to analyse the conditions under which the association of trans-bilayer and partially hydrophilic peptides (or their dispersion in the lipid matrix) is lipid-mediated, and to what extent it is induced by direct interactions between the hydrophilic regions of the peptides. The model properties are calculated by a Monte Carlo computer simulation technique within the canonical ensemble. The results from the model study indicate that direct interactions between the hydrophilic regions of the peptides are necessary to induce peptide association in the lipid bilayer in the fluid phase. Furthermore, peptides within each aggregate are oriented in such a way as to shield their hydrophilic regions from the hydrophobic environment. The average number of peptides present in the aggregates formed depends on the degree of mismatch between the peptide hydrophobic length and the lipid bilayer hydrophobic thickness: The lower the degree of mismatch is the higher this number is. Received: 30 December 1996 / Accepted: 9 May 1997