Simulation studies of protein-induced bilayer deformations, and lipid-induced protein tilting, on a mesoscopic model for lipid bilayers with embedded proteins

Biological membranes are complex and highly cooperative structures. To relate biomembrane structure to their biological function it is often necessary to consider simpler systems. Lipid bilayers composed of one or two lipid species, and with embedded proteins, provide a model system for biological membranes. Here we present a mesoscopic model for lipid bilayers with embedded proteins, which we have studied with the help of the dissipative particle dynamics simulation technique. Because hydrophobic matching is believed to be one of the main physical mechanisms regulating lipid-protein interactions in membranes, we considered proteins of different hydrophobic length (as well as different sizes). We studied the cooperative behavior of the lipid-protein system at mesoscopic time- and lengthscales. In particular, we correlated in a systematic way the protein-induced bilayer perturbation, and the lipid-induced protein tilt, with the hydrophobic mismatch (positive and negative) between the protein hydrophobic length and the pure lipid bilayer hydrophobic thickness. The protein-induced bilayer perturbation was quantified in terms of a coherence length, xi(P), of the lipid bilayer hydrophobic thickness profile around the protein. The dependence on temperature of xi(P), and the protein tilt-angle, were studied above the main-transition temperature of the pure system, i.e., in the fluid phase. We found that xi(P) depends on mismatch, i.e., the higher the mismatch is, the longer xi(P) becomes, at least for positive values of mismatch; a dependence on the protein size appears as well. In the case of large model proteins experiencing extreme mismatch conditions, in the region next to the so-called lipid annulus, there appears an undershooting (or overshooting) region where the bilayer hydrophobic thickness is locally lower (or higher) than in the unperturbed bilayer, depending on whether the protein hydrophobic length is longer (or shorter) than the pure lipid bilayer hydrophobic thickness. Proteins may tilt when embedded in a too-thin bilayer. Our simulation data suggest that, when the embedded protein has a small size, the main mechanism to compensate for a large hydrophobic mismatch is the tilt, whereas large proteins react to negative mismatch by causing an increase of the hydrophobic thickness of the nearby bilayer. Furthermore, for the case of small, peptidelike proteins, we found the same type of functional dependence of the protein tilt-angle on mismatch, as was recently detected by fluorescence spectroscopy measurements.     

Consequences of hydrophobic mismatch between lipids and melibiose permease on melibiose transport

The structural and functional consequences of a mismatch between the hydrophobic thickness d(P) of a transmembrane protein and that d(L) of the supporting lipid bilayer were investigated using melibiose permease (MelB) from Escherichia coli reconstituted in a set of bis saturated and monounsaturated phosphatidylcholine species differing in acyl-chain length. Influence of MelB on the midpoint gel-to-liquid-phase transition temperature, T(m), of the saturated lipids was investigated through fluorescence polarization experiments, with 1,6-diphenyl-1,3,5-hexatriene as the probe, for varying protein/lipid molar ratio. Diagrams in temperature versus MelB concentration showed positive or negative shifts in T(m) with the short-chain lipids DiC12:0-PC and DiC14:0-PC or the long-chain lipids DiC16:0-PC and DiC18:0-PC, respectively. Theoretical analysis of the data yielded a d(L) value of 3.0 +/- 0.1 nm for the protein, similar to the 3.02 nm estimated from hydropathy profiles. Influence of the acyl chain length on the carrier activity of MelB was investigated in the liquid phase, using the monounsaturated PCs. Binding of the sugar to the transporter showed no dependence on the acyl chain length. In contrast, counterflow and Deltapsi-driven experiments revealed strong dependence of melibiose transport on the lipid acyl chain length. Similar bell-shaped transport versus acyl chain length profiles were obtained, optimal activity being supported by diC16:1-PC. On account of a d(P) value of 2.65 nm for the lipid and of various local constraints which would all tend to elongate the acyl chains in contact with the protein, one can conclude that maximal activity was obtained when the hydrophobic thickness of the bilayer matched that of the protein.     

Models of protein lateral arrangements in lipid bilayer membranes. Application to electron spin resonance studies of cytochrome c oxidase

We consider the situation of integral membrane proteins in a lipid bilayer matrix where the size of the polar group of the protein is important in determining the lateral packing of the proteins. We represent the cross-section of the protein hydrophobic core as a hexagon moving on a lattice, and represent the projection of the polar group onto the plane of the bilayer as a shape, parts of which overlap the hexagon. Lattice sites represent lipid molecules. We calculate the fraction of lipid molecules which are adjacent to the hydrophobic core of at least one protein. We use this data to consider the "motion restricted" spectrum observed in electron spin resonance (ESR) probe studies, and compute the dependence of the "motion restricted" fraction upon protein concentration. The resulting curves can be used to analyse ESR data in order to deduce the size and shape of the proteins' polar segment. We have used the range of models examined to study the dependence upon protein concentration of the particular case of the "motion restricted" spectrum of a spin-labelled lipid freely diffusing or, alternatively, covalently bound to cytochrome c oxidase. We find that our calculations are in accord with a model where approximately 60 lipid molecules can fit around an isolated such protein in both halves of the bilayer, and where the polar segment is substantially anisotropic and extends laterally beyond the limits of the hydrophobic core. The latter is in accord with what is known about the structure of cytochrome c oxidase. We indicate further measurements that should be performed in order to establish more definitively the dependence of the "motion restricted" component upon protein concentration, giving the lipid protein ratios at which they should be performed, and we make predictions concerning the results. Finally we argue for a particular unified way of plotting experimental data.     

Interaction of a peptide model of a hydrophobic transmembrane alpha-helical segment of a membrane protein with phosphatidylcholine bilayers: differential scanning calorimetric and FTIR spectroscopic studies.

High-sensitivity differential scanning calorimetry (DSC) and Fourier transform infrared (FTIR) spectroscopy were used to study the interaction of a synthetic model hydrophobic peptide, Lys2-Gly-Leu24-Lys2-Ala-amide, and members of the homologous series of n-saturated diacylphosphatidylcholines. In the low range of peptide mole fractions, the DSC thermograms exhibited by the lipid/peptide mixtures are resolvable into two components. One of these components is fairly narrow, highly cooperative, and exhibits properties which are similar to but not identical with those of the pure lipid. In addition, the fractional contribution of this component to the total enthalpy change, the peak transition temperature, and cooperativity decrease with an increase in peptide concentration, more or less independently of acyl chain length. The other component is very broad and predominates in the high range of peptide concentration. These two components have been assigned to the chain-melting phase transitions of populations of bulk lipid and peptide-associated lipid, respectively. Moreover, when the mean hydrophobic thickness of the PC bilayer is less than the peptide hydrophobic length, the peptide-associated lipid melts at higher temperatures than does the bulk lipid and vice versa. In addition, the chain-melting enthalpy of the broad endotherm does not decrease to zero even at high peptide concentrations, suggesting that this peptide reduces but do not abolish the cooperative gel/liquid-crystalline phase transition of the lipids with which it is in contact. Our DSC results indicate that the width of the phase transition observed at high peptide concentration is inversely but discontinuously related to hydrocarbon chain length and that gel phase immiscibility occurs when the hydrophobic thickness of the bilayer greatly exceeds the hydrophobic length of the peptide. The FTIR spectroscopic data indicate that the peptide forms a very stable alpha-helix under all of our experimental conditions but that small distortions of its alpha-helical conformation are induced in response to any mismatch between peptide hydrophobic length and bilayer hydrophobic thickness. These results also indicate that the peptide alters the conformational disposition of the acyl chains in contact with it and that the resultant conformational changes in the lipid hydrocarbon chains tend to minimize the extent of mismatch of peptide hydrophobic length and bilayer hydrophobic thickness.     

Energetics of hydrophobic matching in lipid-protein interactions

Lipid chain length modulates the activity of transmembrane proteins by mismatch between the hydrophobic span of the protein and that of the lipid membrane. Relative binding affinities of lipids with different chain lengths are used to estimate the excess free energy of lipid-protein interaction that arises from hydrophobic mismatch. For a wide range of integral proteins and peptides, the energy cost is much less than the elastic penalty of fully stretching or compressing the lipid chains to achieve complete hydrophobic matching. The chain length dependences of the free energies of lipid association are described by a model that combines elastic chain extension with a free energy term that depends linearly on the extent of residual mismatch. The excess free energy densities involved lie in the region of 0.5-2.0 k_BT.nm⁻². Values of this size could arise from exposure of hydrophobic groups to polar portions of the lipid or protein, but not directly to water, or alternatively from changes in tilt of the transmembrane helices that are energetically comparable to those activating mechanosensitive channels. The influence of hydrophobic mismatch on dimerization of transmembrane helices and their transfer between lipid vesicles, and on shifts in chain-melting transitions of lipid bilayers by incorporated proteins, is analyzed by using the same thermodynamic model. Segmental order parameters confirm that elastic lipid chain distortions are insufficient to compensate fully for the mismatch, but the dependence on chain length with tryptophan-anchored peptides requires that the free energy density of hydrophobic mismatch should increase with increasing extent of mismatch.     

Molecular sorting of lipids by bacteriorhodopsin in dilauroylphosphatidylcholine/distearoylphosphatidylcholine lipid bilayers.

A combined experimental and theoretical study is performed on binary dilauroylphosphatidylcholine/distearoylphosphatidylcholine (DLPC/DSPC) lipid bilayer membranes incorporating bacteriorhodopsin (BR). The system is designed to investigate the possibility that BR, via a hydrophobic matching principle related to the difference in lipid bilayer hydrophobic thickness and protein hydrophobic length, can perform molecular sorting of the lipids at the lipid-protein interface, leading to lipid specificity/selectivity that is controlled solely by physical factors. The study takes advantage of the strongly nonideal mixing behavior of the DLPC/DSPC mixture and the fact that the average lipid acyl-chain length is strongly dependent on temperature, particularly in the main phase transition region. The experiments are based on fluorescence energy transfer techniques using specifically designed lipid analogs that can probe the lipid-protein interface. The theoretical calculations exploit a microscopic molecular interaction model that embodies the hydrophobic matching as a key parameter. At low temperatures, in the gel-gel coexistence region, experimental and theoretical data consistently indicate that BR is associated with the short-chain lipid DLPC. At moderate temperatures, in the fluid-gel coexistence region, BR remains in the fluid phase, which is mainly composed of short-chain lipid DLPC, but is enriched at the interface between the fluid and gel domains. At high temperatures, in the fluid phase, BR stays in the mixed lipid phase, and the theoretical data suggest a preference of the protein for the long-chain DSPC molecules at the expense of the short-chain DLPC molecules. The combined results of the experiments and the calculations provide evidence that a molecular sorting principle is active because of hydrophobic matching and that BR exhibits physical lipid selectivity. The results are discussed in the general context of membrane organization and compartmentalization and in terms of nanometer-scale lipid-domain formation.     

Interaction of myelin basic protein with different ionization states of phosphatidic acid and phosphatidylserine

Myelin basic protein (BP) has a perturbing effect on some lipids, causing, among other effects, a decrease in the temperature and enthalpy of the phase transition. This is believed to be a result of penetration of some hydrophobic residues of the protein partway into the lipid bilayer. Variations in the perturbing effect of BP on different acidic lipids has been attributed to the ability of the lipids to participate in intermolecular hydrogen bonding which inhibits penetration of the protein. Participation in intermolecular hydrogen bonding depends on the ionization state of the lipid as well as the type of lipid. In order to further test the dependence of the degree of penetration of BP on the intermolecular hydrogen bonding properties of lipids, the effect of BP on the phase transition of lipids in different ionization states was studied using differential scanning calorimetry. Dipalmitoylphosphatidic acid (DPPA) and dimyristoylphosphatidylserine (DMPS) were studied at different pH-values from 4 to 9.5. The results were compared to data obtained earlier with phosphatidylglycerol (PG), which is in the same ionization state at pH-values above 4, in order to distinguish the effects of pH on the protein from effects on the lipids. The perturbing effect of BP on PG increases with increase in pH. This is probably a result of the increasing hydrophobicity of the protein as the histidines become deprotonated, which allows greater penetration of the protein into the bilayer. In contrast, the effect on DPPA was greatest at low pH, where the state of ionization of the lipid is less than 1 and protein binding utilizes all of the hydrogen bond accepting sites (P-O-) on the lipid. BP had no perturbing effect on DPPA at higher pH where the state of ionization is between 1 and 1.5, and hydrogen bond accepting and donating sites (P-OH) are still available even after binding of the protein. Thus hydrogen bonding occurs at high pH and penetration of hydrophobic residues of the protein into DPPA is inhibited. BP had a large perturbing effect on DMPS at all pH values above 4 suggesting that lipid intermolecular hydrogen bonding does not occur in the presence of the protein and its hydrophobic residues consequently can penetrate into the bilayer. The protein may inhibit hydrogen bonding by binding electrostatically to the anionic hydrogen bond accepting group of PS.     

Effect of alterations in lipid packing order by hydrophobic solutes on the association state of protein assemblies in model membranes

The effects of two hydrophobic solutes which perturb lipid packing order, permethrin and allethrin, on the aggregated state of a lipid membrane-incorporated protein, bacteriorhodopsin (BR), have been determined by resonance energy transfer measurements. As temperature is increased from well below the main gel-fluid phase transition temperature (Tc) of the lipid, patches of aggregated BR dissociate into monomers, a few degrees below the Tc (M.P. Heyn, A. Blume, M. Rehorek and N.A. Dencher, Biochemistry 20 (1981) 7109; M.P. Heyn, R.J. Cherry and N.A. Dencher, Biochemistry 20 (1981) 840). Permethrin and allethrin were found to cause a decrease in the temperature of BR disaggregation which was associated with a decrease in the Tc of the lipid. In gel phase dipalmitoylphosphatidylcholine at 25 degrees C, the pertubing effects of permethrin on lipid packing order were associated with a decrease in the average patch radius from 123 to 33 A. It is concluded that perturbation of lipid packing order by small hydrophobic molecules may alter the stability of protein assemblies in membranes.     

Calorimetric investigation of polymyxin binding to phosphatidic acid bilayers

The cooperative binding process between the antibiotic peptide polymyxin-B and negatively-charged phosphatidic acid bilayers was investigated by differential thermal analysis and completed by fluorescence polarization measurements. The sigmoidal binding curves were analyzed in terms of the interaction energy within a domain formed by polymyxin and phosphatidic acid molecules. The formation of such a heterogeneous domain structure was favoured by high concentration of external monovalent ions. The cooperativity of the binding increased while a charge-induced decrease in the phase transition temperature of the pure lipid phase was observed with increasing ion concentration at a given pH. The reduced lateral coupling within the lipid bilayer in the presence of salt ions, as demonstrated by an increase in the lipid phase transition enthalpy, was considered to facilitate the cooperative domain formation. Moreover, an increase in the cooperativity of the polymyxin binding could be observed if phosphatidic acids of smaller chain length and thus of a lowered phase transition temperature were used. By the use of chemically-modified polymyxin we were able to demonstrate the effect of electrostatic and hydrophobic interaction. Acetylated polymyxin with a reduced positive charge was used to demonstrate the pure hydrophobic effect of polymyxin binding leading to a decrease in the phosphatidic acid phase transition temperature by about 20°C. The cooperativity of the binding was strongly reduced. Cleavage of the hydrophobic polymyxin tail yielded a colistinnonapeptide which caused an electrostatically-induced increase in the phosphatidic acid phase transition temperature. With unmodified polymyxin we observed the combined effects of electrostatic as well as hydrophobic interaction making this model system interesting for the understanding of lipid-protein interactions. Evidence is presented that the formation of the polymyxin-phosphatidic acid complex is a lateral phase separation phenomenon.     

Interaction of poly(l-lysines) with negatively charged membranes: an FT-IR and DSC study

The influence of the binding of poly(l-lysine) (PLL) to negatively charged membranes containing phosphatidylglycerols (PG) was studied by DSC and FT-IR spectroscopy. We found a general increase in the main transition temperature as well as increase in hydrophobic order of the membrane upon PLL binding. Furthermore we observed stronger binding of hydration water to the lipid head groups after PLL binding. The secondary structure of the PLL after binding was studied by FT-IR spectroscopy. We found that PLL binds in an α-helical conformation to negatively charged DPPG membranes or membranes with DPPG-rich domains. Moreover we proved that PLL binding induces domain formation in the gel state of mixed DPPC/DPPG or DMPC/DPPG membranes as well as lipid remixing in the liquid–crystalline state. We studied these effects as a function of PLL chain length and found a significant dependence of the secondary structure, phase transition temperature and domain formation capacity on PLL chain length and also a correlation between the peptide secondary structure and the phase transition temperature of the membrane. We present a system in which the membrane phase transition triggers a highly cooperative secondary structure transition of the membrane-bound peptide from α-helix to random coil. Dedicated to Prof. K. Arnold on the occasion of his 65th birthday.     

The influence of protein-lipid interactions on the order-disorder conformational transitions of the hydrocarbon chain.

The phases of simple systems involving one type of protein (lysozyme or cytochrome c) and one type of lipid (phosphatidic acid) have been characterized by X-ray crystallography, chemical analysis and spin-labeling technique as a function of temperature. They are of the lamellar type with alternative protein monolayers and lipid bilayers. According to the pH, two types of lamellar phases are obtained, one where the lipid-protein interactions are mainly hydrophobic, the other where they are electrostatic. In both cases, a phase transition occurs as temperature is lowered, between a high temperature phase, where all the lipids are in the liquid-like state, and another phase where some lipid chains are rigid. In the case of the phases with electrostatic interaction, it is shown that the onset of the order-disorder transition is shifted towards low temperature as compared with the homologous lipid-water phase and that the protein content of the phase decreases as the ratio of the liquid to rigid hydrocarbon chains decreases. This leads us to suggest that in the systems studied in this work the proteins interact only with lipid in the liquid-like state. In the case of the phases with hydrophobic interaction, it is shown that the extent of hydrophobic interaction between protein and lipid increases as the unsaturation of the hydrocarbon chains increases. The onset of the order-disorder transition shows a greater shift towards low temperature than the one observed in the case of the phase with electrostatic interaction.     

Modulation of Na,K-ATPase and Na-ATPase activity by phospholipids and cholesterol. I. Steady-state kinetics

The effects of phospholipid acyl chain length (n(c)), degree of acyl chain saturation, and cholesterol on Na,K-ATPase reconstituted into liposomes of defined lipid composition are described. The optimal acyl chain length of monounsaturated phosphatidylcholine in the absence of cholesterol was found to be 22 but decreased to 18 in the presence of 40 mol % cholesterol. This indicates that the hydrophobic matching of the lipid bilayer and the transmembrane hydrophobic core of the membrane protein is a crucial parameter in supporting optimal Na,K-ATPase activity. In addition, the increased bilayer order induced by both cholesterol and saturated phospholipids could be important for the conformational mobility of the Na,K-ATPase changing the distribution of conformations. Lipid fluidity was important for several parameters of reconstitution, e.g., the amount of protein inserted and the orientation in the liposomes. The temperature dependence of the Na,K-ATPase as well of the Na-ATPase reactions depends both on phospholipid acyl chain length and on cholesterol. Cholesterol increased significantly both the enthalpy of activation and entropy of activation for Na,K-ATPase activity and Na-ATPase activity of Na,K-ATPase reconstituted with monounsaturated phospholipids. In the presence of cholesterol the free energy of activation was minimum at a lipid acyl chain length of 18, the same that supported maximum turnover. In the case of ATPase reconstituted without cholesterol, the minimum free energy of activation and the maximum turnover both shifted to longer acyl chain lengths of about 22.     

Calorimetric investigation of polymyxin binding to phosphatidic acid bilayers

The cooperative binding process between the antibiotic peptide polymyxin-B and negatively-charged phosphatidic acid bilayers was investigated by differential thermal analysis and completed by fluorescence polarization measurements. The sigmoidal binding curves were analyzed in terms of the interaction energy within a domain formed by polymyxin and phosphatidic acid molecules. The formation of such a heterogeneous domain structure was favoured by high concentration of external monovalent ions. The cooperativity of the binding increased while a charge-induced decrease in the phase transition temperature of the pure lipid phase was observed with increasing ion concentration at a given pH. The reduced lateral coupling within the lipid bilayer in the presence of salt ions, as demonstrated by an increase in the lipid phase transition enthalpy, was considered to facilitate the cooperative domain formation. Moreover, an increase in the cooperativity of the polymyxin binding could be observed if phosphatidic acids of smaller chain length and thus of a lowered phase transition temperature were used. By the use of chemically-modified polymyxin we were able to demonstrate the effect of electrostatic and hydrophobic interaction. Acetylated polymyxin with a reduced positive charge was used to demonstrate the pure hydrophobic effect of polymyxin binding leading to a decrease in the phosphatidic acid phase transition temperature by about 20 degrees C. The cooperativity of the binding was strongly reduced. Cleavage of the hydrophobic polymyxin tail yielded a colistinnonapeptide which caused an electrostatically-induced increase in the phosphatidic acid phase transition temperature. With unmodified polymyxin we observed the combined effects of electrostatic as well as hydrophobic interaction making this model system interesting for the understanding of lipid-protein interactions. Evidence is presented that the formation of the polymyxin-phosphatidic acid complex is a lateral phase separation phenomenon.     

Elastic Interactions of Photosynthetic Reaction Center Proteins Affecting Phase Transitions and Protein Distributions

Reaction-center proteins of Rhodopseudomonas Sphaeroides reconstituted into phosphatidylcholine vesicles shift and broaden the fluid-gel transition of the lipid bilayer. The amount of broadening and temperature shift of the transition depend both on protein concentration and on lipid chain length. In particular, the direction of the transition temperature shift is very sensitive to lipid chain length. Electron micrographs show homogeneous protein distribution on the fluid surface whereas the solid phase contains protein aggregates the type depending on chain length. The results can qualitatively be understood in the framework of a mattress model of lipid/protein interactions in membranes.     

Peptide models of the helical hydrophobic transmembrane segments of membrane proteins: interactions of acetyl-K2-(LA)12-K2-amide with phosphatidylethanolamine bilayer membranes

High-sensitivity differential scanning calorimetry (DSC) and Fourier transform infrared (FTIR) spectroscopy were used to study the interaction of a synthetic alpha-helical hydrophobic transmembrane peptide, acetyl-Lys(2)-(Leu-Ala)(12)-Lys(2)-amide [(LA)(12)], and members of a homologous series of n-saturated diacylphosphatidylethanolamines (PEs). In the lower range of peptide mole fractions, the DSC endotherms exhibited by the lipid/peptide mixtures consist of two components. The temperature and cooperativity of the sharper, higher temperature component are very similar to those of pure PE bilayers and are almost unaffected by variations in the protein/lipid ratio. However, the fractional contribution of this component to the total enthalpy changes decreases with increases in peptide concentration, and this component completely disappears at higher protein mole fractions. The other component, which is less cooperative and occurs at a lower temperature, predominates at higher protein concentrations. These two components of the DSC endotherm have been assigned to the chain-melting phase transitions of peptide-nonassociated and peptide-associated PE molecules, respectively. Although the temperature at which the peptide-associated PE molecules melt is progressively decreased by increases in (LA)(12) concentration, the magnitude of this downward shift is progressively greater as the length of the PE hydrocarbon chain decreases. As well, mixtures of (LA)(12) with the longer chain PEs exhibit unusual biomodal enthalpy variations, suggesting peptide immiscibility in thicker gel state bilayers. Moreover, the enthalpy of the chain-melting transition of the peptide-associated PE does not decrease to zero even at high peptide concentrations, indicating that (LA)(12) attenuates but does not abolish the cooperative gel/liquid-crystalline phase transition of the lipids with which it is in contact. Our FTIR spectroscopic data indicate that (LA)(12) remains in a predominantly alpha-helical conformation in liquid-crystalline PE bilayers of various hydrophobic thickness but that the helical conformation is altered in gel-state PE bilayers generally, probably due to peptide lateral aggregation. These data also suggest that (LA)(12) significantly disorders the hydrocarbon chains of adjacent PE molecules in both the gel and liquid-crystalline states, relatively independently of lipid hydrocarbon chain length. Many aspects of PE/(LA)(12) interactions exhibit a different dependence on the hydrophobic thickness of the host bilayer than was observed in our previous study of (LA)(12)-phosphatidylcholine (PC) model membranes [Zhang et al. (1995) Biochemistry 34, 2362-2371]. The differing effects of (LA)(12) incorporation on PE and PC bilayers is ascribed primarily to the much stronger lipid polar headgroup interactions characteristic of the former system. Finally, the considerable differences observed in the behavior of (LA)(12) and the related polyleucine-based peptide P(24) in both PC and PE bilayers indicate that the structure of the hydrophobic core of alpha-helical transmembrane peptides can affect their conformational plasticity and state of aggregation and thus the nature of their interactions with different phospholipid bilayers.     

Lateral diffusion of gramicidin S, M-13 coat protein and glycophorin in bilayers of saturated phospholipids. Mean field and Monte Carlo studies

We have developed a general model that relates the lateral diffusion coefficient of one isolated large intrinsic molecule (mol. wt. greater than or approximately 1000) in a phosphatidylcholine bilayer to the static lipid hydrocarbon chain order. We have studied how protein lateral diffusion can depend upon protein-lipid interactions but have not investigated possible non-specific contributions from gel-state lattice defects. The model has been used in Monte Carlo simulations or in mean-field approximations to study the lateral diffusion coefficients of Gramicidin S, the M-13 coat protein and glycophorin in dimyristoyl- and dipalmitoylphosphatidylcholine (DMPC and DPPC) bilayers as functions of temperature. Our calculated lateral diffusion coefficients for Gramicidin S and the M-13 coat protein are in good agreement with what has been observed and suggest that Gramicidin S is in a dimeric form in DMPC bilayers. In the case of glycophorin we find that the 'ice breaker' effect can be understood as a consequence of perturbation of the lipid polar region around the protein. In order to understand this effect is necessary that the protein hydrophilic section perturb the polar regions of at least approx. 24 lipid molecules, in good agreement with the numbers of 29-30 measured using 31P-NMR. Because of lipid-lipid interactions this effect extends itself out to four or five lipid layers away from the protein so that the hydrocarbon chains of between approx. 74 and approx. 108 lipid molecules are more disordered in the gel phase, so contributing less to the transition enthalpy, in agreement with the numbers of 80-100 deduced from differential scanning calorimetry (DSC). An understanding of the abrupt change in the diffusion coefficient at a temperature below the main bilayer transition temperature requires an additional mechanism. We propose that this change may be a consequence of a 'coupling-uncoupling' transition involving the protein hydrophilic section and the lipid polar regions, which may be triggered by the lipid bilayer pretransition. Our calculation of the average number of gauche bonds per lipid chain as a function of temperature and distance away from an isolated polypeptide or integral protein shows the extent of statically disordered lipid around such molecules. The range of this disorder depends upon temperature, particularly near the main transition.     

Stability and phase separation in mixed monopolar lipid/bolalipid layers

The phase stability of a fluid lipid layer that is a mixture of conventional monopolar lipids and C20 bipolar bolalipids was studied using a mean field theory that explicitly includes molecular details and configurational properties of the lipid molecules. The effect of changing the fraction of bolalipids, as well as the length of the hydrocarbon chain of the monopolar lipids, was probed. A phase separation between two liquid lipid phases was found when a mismatch exists in the optimal hydrophobic thicknesses of the pure bolalipid and monopolar lipid layers. The lipid mixture phase separates into a thin bolalipid-rich layer and a thicker monopolar-rich layer. The thin membrane phase is mainly composed of transmembrane bolalipid molecules whose polar heads are positioned at opposite membrane-water interfaces. In the monopolar lipid-rich phase, bolalipids are the minor component and most of them assume a looping configuration where both headgroups are present at the same membrane-water interface. For mixed layers that form a single lipid phase across all bolalipid concentrations, the hairpin-transmembrane ratio strongly depends on the hydrocarbon chain length of the monopolar lipid and the bolalipid concentration. The C-D bond order parameters of the different species have been calculated. Our findings suggest that the concentration-dependent phase transition should be experimentally observable by measuring of the order parameters through quadrupolar splitting experiments. The driving force for the phase separation in the monopolar lipid/bolalipid mixture is the packing mismatch between hydrophobic regions of the monopolar lipid hydrocarbon chains and the membrane-spanning bolalipid chains. The results from the molecular theory may be useful in the design of stable lipid layers for integral membrane protein sensing.     

The Hydrophobic Temperature Dependence of Amino Acids Directly Calculated from Protein Structures

The hydrophobic effect is the main driving force in protein folding. One can estimate the relative strength of this hydrophobic effect for each amino acid by mining a large set of experimentally determined protein structures. However, the hydrophobic force is known to be strongly temperature dependent. This temperature dependence is thought to explain the denaturation of proteins at low temperatures. Here we investigate if it is possible to extract this temperature dependence directly from a large set of protein structures determined at different temperatures. Using NMR structures filtered for sequence identity, we were able to extract hydrophobicity propensities for all amino acids at five different temperature ranges (spanning 265-340 K). These propensities show that the hydrophobicity becomes weaker at lower temperatures, in line with current theory. Alternatively, one can conclude that the temperature dependence of the hydrophobic effect has a measurable influence on protein structures. Moreover, this work provides a method for probing the individual temperature dependence of the different amino acid types, which is difficult to obtain by direct experiment.     

The influence of protein-lipid interactions on the order-disorder conformational transitions of the hydrocarbon chain

The phases of simple systems involving one type of protein (lysozyme or cytochrome $\text{c}$) and one type of lipid (phosphatidic acid) have been characterized by X-ray crystallography, chemical analysis and spin-labeling technique as a function of temperature. They are of the lamellar type with alternative protein monolayers and lipid bilayers. According to the pH, two types of lamellar phases are obtained, one where the lipid-protein interactions are mainly hydrophobic, the other where they are electrostatic. In both cases, a phase transition occurs as temperature is lowered, between a high temperature phase, where all the lipids are in the liquid-like state, and another phase where some lipid chains are rigid. In the case of the phases with electrostatic interaction, it is shown that the onset of the order-disorder transition is shifted towards low temperature as compared with the homologous lipid-water phase and that the protein content of the phase decreases as the ratio of the liquid to rigid hydrocarbon chains decreases. This leads us to suggest that in the systems studied in this work the proteins interact only with lipid in the liquid-like state. In the case of the phases with hydrophobic interaction, it is shown that the extent of hydrophobic interaction between protein and lipid increases as the unsaturation of the hydrocarbon chains increases. The onset of the order-disorder transition shows a greater shift towards low temperture than the one observed in the case of the phase with electrostatic interaction.