Interactions of fatty acids with phosphatidylethanolamine membranes: X-ray diffraction and molecular dynamics studies.

An experimental and theoretical study on 1,2-dielaidoyl-sn-glycero-3-phosphoethanolamine (DEPE) membranes containing fatty acids (FAs) was performed by means of X-ray diffraction analysis and molecular dynamics (MD) simulations. The study was aimed at understanding the interactions of several structurally related FAs with biomembranes, which is necessary for further rational lipid drug design in membrane-lipid therapy. The main effect of FAs was to promote the formation of a H(II) phase, despite a stabilization of the coexisting L(alpha) + H(II) phases. Derivatives of OA exhibited a specific density profile in the direction perpendicular to the bilayer that reflects differences in the relative localization of the carboxylate group within the polar region of the membrane as well as in the degree of membrane penetration of the FA acyl chain. Hydroxyl and methyl substituents at carbon-2 in the FA acyl chain were identified as effective modulators of the position of carboxylate group in the lipid bilayer. Our data highlight the specific potential of each FA in modulating the membrane structure properties.

It is a well-known fact that dietary fat uptake infl uences the fatty acid (FA) composition of cellular membranes. For example, the Mediterranean diet characterized by a high consumption of olive oil, which contains a large proportion of oleic acid ( ‫ف‬ 80%), was associated with increases in the levels of this FA in various plasma membrane lipid species in rat and human cells ( 1,2 ). Such a change in FA composition has a substantial impact on structural membrane properties. Indeed, for various diseases, including hypertension, it has been demonstrated that changes in biophysical membrane parameters are closely associated to the pathology of these disorders ( 3,4 ). In this context, we recently demonstrated that hypertensive patients on a Mediterranean-style diet containing virgin olive oil showed a reduction of elevated blood pressure, which correlated with changes in membrane lipid composition and which was linked with plasma membrane structural properties, such as an increase in phosphatidylethanolamine (PE) lipid class and a higher propensity to form a nonlamellar H II structure of reconstituted plasma membranes ( 5 ). This is an important fi nding, because it has been shown that PE lipids are able to modulate peripheral protein-membrane interactions via their H II phase propensity ( 6 ). Moreover,

Model membrane preparation
Multilamellar lipid vesicles (MLV; 15% w/w) were prepared in 10 mM HEPES, 100 mM NaCl, 1 mM EDTA, pH 7.4 (HEPES buffer). Lipid mixtures were hydrated in the presence or absence of FAs at the specifi ed molar ratio, homogenized with a pestle-type minihomogenizer (Sigma) and vortexed until a suspension was obtained. The suspensions were then submitted to fi ve heating and cooling cycles (heating to 70°C and cooling to 4°C) and equilibrated prior to data acquisition, as previously reported ( 13 ).

X-ray diffraction
Small-and wide-angle synchrotron radiation X-ray scattering data (SAXS and WAXS) were collected simultaneously on the Soft Condensed Matter beamline A2 of Hasylab at the storage ring DORIS III of the Deutsches Elektronen Synchrotron (DESY) using standard procedures. Samples were heated from 27°C to 75°C and then cooled to 27°C at a scan rate of 1°C/min. Data collection was performed as described before ( 13 ).

Molecular dynamics simulations
All computer simulations were performed using the GRO-MACS 3.3.2 package ( 21 ). Molecular dynamics (MD) simulations were performed at 323 K (above the L ␤ -to-L ␣ phase transition and below the temperature of L ␣ and H II coexistence) using periodic boundary conditions and a time step of 4 fs. The temperature was kept constant using separate thermostats for lipids, water, and ions with a time constant of 0.1 ps. The pressure in the three coordinate directions was kept at 0.1 MPa by independent Berendsen barostats using a time constant of 1.0 ps. Lennard-Jones interactions were computed with a cutoff of 10 Å, and electrostatic interactions were treated with the Particle Mesh Ewald (PME) procedure.
In a fi rst step, a box containing 216 DEPE lipids arranged in a bilayer, more than 10,000 water molecules, 28 sodium, and 28 chloride ions was constructed as described in a previous study ( 22 ) and simulated for 200 ns. In a second step, boxes for DEPE:FA systems were constructed by replacing 5 lipids on each leafl et by FAs and removing 6 additional lipids. This procedure allows simulating systems with the desired DEPE:FA molecular ratio of 20:1. Similar simulations working with phospholipid/FA systems have recently been reported (23)(24)(25). To preserve electroneutrality, 10 water molecules were replaced by 10 additional Na + ions. Each system was energy minimized and subsequently, subjected to 200 ns MD simulations at 323 K following the protocol previously described in various reports ( 22,26,27 ). For each trajectory, the fi rst 100 ns were considered as equilibration period and were not included in the analysis.
In the absence of specifi c force fi eld parameters for DEPE lipids available, a combination of DPPC ( 28 ) and POPE parameters (taken from Tieleman's group [http://moose.bio.ucalgary.ca/]) alteration of membrane FA composition is able to infl uence pivotal elements of important signaling transduction pathways, which fi nally control vital physiological functions, such as cardiovascular parameters ( 4,7 ). Thus, if cell functions are altered when membrane lipid composition changes, then modulation of biophysical membrane properties by administration of appropriate FAs might be a strategy to revert disease status. This novel pharmacological approach has recently been introduced as membrane-lipid therapy ( 8 ). An important step to a rational design of new lipid molecules with pharmacological activity is a comprehensive study of the differential effects of closely related FAs on the structural properties of phospholipid membranes.
The asymmetric distribution of PE phospholipids and the relatively high proportion of this lipid in the cytosolic leafl et of the mammalian plasma membrane ( ‫ف‬ 30%) indicates that PE should have a great functional relevance. In fact, PE lipids have an intrinsic propensity to segregate forming nonlamellar H II structures in vitro ( 9 ), and some special features were attributed to H II -prone phospholipids, such as the functional control of membrane proteins and the structural organization inside cells (10)(11)(12). On the other hand, FAs incorporated in phospholipids or present in the membrane in their unesterifi ed form are able to modulate the structure of model PE membranes by interacting specifi cally with the latter ( 13,14 ). This is also true for synthetic fatty acids. 2-Hydroxyoleic acid (2OHOA), a synthetic derivative of the naturally occurring oleic acid, induced the formation of H II structures in model membranes ( 15 ). In vivo, this observation was paralleled by an alteration of the membrane lipid composition which was related to the observed increase of G proteins in the plasma membrane of 2OHOA-treated animals ( 12,16 ). Moreover, 2OHOA showed a strongly enhanced antihypertensive action in this animal model of hypertension ( 17,18 ), thus demonstrating the viability of this new pharmacological approach.
In the present study, we explore the effect of the structural parameters acyl chain length, double bond position, and stereochemical properties of series of FAs ( Fig.  1 ) on the structural properties of PE model membranes, using X-ray diffraction and MD simulations. Our study comprises a series of both natural and synthetic FAs, which are described below. OA and SA (18:0) are the major lipid components of animal tissues. Palmitoleic acid ( cis-16:1n7, PoA) is the delta-9 desaturase product of palmitic acid (16:0) and a minor component of most lipid classes, generally found in higher concentrations in the liver. Cis -vaccenic acid ( cis -18:1n7, cVcc) is the direct elongation product of palmitoleic acid and an essential fatty acyl-component of cardiolipin. Palmitelaidic ( trans -16:1n7, PeA) and EA are dietary trans-FAs, which have received growing attention because of their negative infl uence on cardiovascular health ( 19 ). Four additional synthetic FAs are included to elucidate the effect of chemical modifi cations on membrane structural parameters: 2-methyloleic, 2-hydroxyoleic, 2-methylstearic, and 2-hydroxystearic FAs. rameters ( Fig. 2 and Table 1 ). X-ray heating scans of DEPE:FA (20:1, mol:mol) mixtures showed a L ␤ phase up to 38°C-40°C with a repeat distance value of 6.5-6.7 nm similar to the value shown for DEPE alone. A single L ␣ phase was present in the range from 38-39°C to 51-57°C, with a repeat distance value of 5.4-5.5 nm at 40°C, which decreased linearly with temperature (compressibility factor ~ Ϫ 0.013 nm/°C). Above these temperatures, the L ␣ phase developed into a binary system, where L ␣ and H II phases coexisted in a range that depended on each FA. Finally, a single H II phase appeared with a repeat distance value of ‫ف‬ 6.4 nm at 72°C and a constant compressibility factor in the range from Ϫ 0.020 to Ϫ 0.025 nm/°C. The main effect of FAs on membrane structure was to induce the formation of a H II phase at a lower temperature than DEPE alone, despite a stabilization of the coexisting L ␣ + H II phases. X-ray measurements did not highlight the differences in the interaction between DEPE and FAs, a detail that was further investigated by MD simulations.

Molecular dynamics study
Analysis of the MD trajectories of the different DEPE:FA systems permits the characterization of the molecular arrangement of each FA in the DEPE bilayer at the atomistic was used. The area per DEPE lipid is not known experimentally, although it is reasonable to expect to be similar to DPPE, for which it has been determined at two different temperatures ( 29 ). Preliminary simulations (data not shown) employing these parameters revealed that DEPE bilayers were half way to the gel phase, in agreement with recent simulations of DPPE membranes ( 30,31 ). This suggested the necessity to introduce a correction for PE lipids. Scaling Lennard-Jones parameters for the acyl carbons by a factor of 3/4 provided an area compatible with all the observations above. FA topologies were generated by PRODRG 2.5 ␤ server ( 32 ), whereas the parameters were taken by analogy to the lipids force fi eld ( 28 ). All FAs were modeled as negatively charged according to the results of our previous report ( 15 ). The apparent pKa of FAs is typically about 7 when incorporated in the bilayer of a phospholipid vesicle, although the specifi c value varies depending upon the environment ( 33,34 ).

X-ray diffraction study
We addressed the effect of FAs on the structural properties of model DEPE membranes. The sequence of diffraction patterns collected showed defi ned phase transitions that allowed unequivocal characterization of the L ␤ , L ␣ and H II phases, and their respective lattice pa- level, as well as the identifi cation of specifi c interactions. According to the system design (see "Materials and Methods"), the two leafl ets in the DEPE bilayer can be considered as equivalent, and therefore, the analysis presented in this report is an average of the two layers. Whenever possible, the analysis of the systems is displayed in three groups corresponding to the type of FA present: (1) OA series -OA, 2OHOA and 2MOA-; (2) SA series -SA, 2OHSA and 2MSA-; and (3) the remaining FAs not substituted at carbon -2 of the acyl chain (PeA, PoA, EA, and cVcc).
FA positioning in the membrane. Fig. 3 displays snapshots of equilibrated DEPE:FA systems showing the disposition of the FAs in the bilayer for each of the three members the OA series. The average position of each type of FA within a DEPE layer was monitored using the atom density distribution along the axis perpendicular to the bilayer plane ( Fig. 4A ). Decomposition of the total atom density distribution into its components allows us to assign the distinct features of each FA. The density of the DEPE molecules exhibits a maximum in the region 0.5-2.5 nm relative to the bilayer center with a peak at 2.1 nm that corresponds to the head groups (phosphate and ethanolamine). The density of water suggests that molecules penetrate up to 1.75 nm with an interfacial region up to 3 nm that overlaps with the density provided by the polar head groups ( 35 ).
All FAs occupy a region between the bilayer center and up to 3.0 nm with differences depending on their specifi c structural features. In the OA series, the atom density distribution of 2OHOA is shifted toward the water phase when compared with OA, and in the opposite direction in the case of the 2MOA. 2OHOA and (to a lesser extent) 2MOA systems exhibit a peak at 2.0-2.1 nm (at the same position as the maximum in the profi le showed by DEPE alone), which originates from the additional contribution of the hydroxyl/methyl group at position C2. A similar trend is found within the SA series. The remaining FAs with 18 carbons (EA and cVcc) showed similar profi les to OA and SA. On the contrary, FAs with shorter chains (PoA and PeA) exhibited a much smaller density in the region close to the bilayer center, indicating that 16 carbons are not enough to go across the whole DEPE layer. Fig. 4B displays density distributions for FA carboxyl carbons along the axis perpendicular to the bilayer plane. Table 2 shows the (time and system) average positions, which provides a measure of the relative distances with the ethanolamine moieties. It can be seen that there is a small difference both in the amplitude of the peak and also in puted from the distance between the planes defi ned by the phosphate groups in each layer) and the average area per lipid ( Table 3 ). The different DEPE:FA systems exhibit average thickness ranging 4.90-4.92 ± 0.02 nm, slightly smaller than the value computed without FAs of 4.93 ± 0.02 nm. Comparatively, theoretical estimates obtained from simulations are ‫ف‬ 0.5-0.8 nm smaller than the experimental data of the lamellar repeat distance (d L ␣ ) in excess of water obtained by X-ray diffraction experiments (see Table 1 ). Such difference can be attributable to the water layer between lamellar repeats. The area per lipid was estimated as described elsewhere ( 35 )  Effect of the fatty acid on lipid order. The lipid bilayer order was assessed by measuring its average thickness (com-   ). The angular coeffi cient of the interplanar distance dependence on the temperature ( ␦ d/ ␦ T) < 0 indicates a compression process. 2MOA, 2-methyloleic acid; 2MSA, 2-methylstearic acid; 2OHOA, 2-hydroxyoleic acid; 2OHSA, 2-hydroxystearic acid; cVcc, cis -vaccenic acid; DEPE, 1,2-dielaidoyl -sn -glycero -3-phosphoethanolamine; EA, elaidic acid; FA, fatty acid; H II , inverted hexagonal phase; L ␤ , gel lamellar phase, L ␣ , liquid-crystalline lamellar phase; OA, oleic acid; PeA, palmitelaidic acid; PoA, palmitoleic acid; SA, stearic acid.
a The compressibility factor of the L ␣ phase was linear in the single-or two-phase regions. Single H II phases also had linear compressibility factor. The temperature range where the L ␣ phase was observed is shown in ⌬ T L ␣ .
b The parenthesis indicates the temperature limit of the L ␣ phase in a two-phase region. Values on the left correspond to L ␤ + L ␣ and on the right to L ␣ + H II temperature range coexistence. c d L ␣ at 40°C and d HII at 72°C. The correlation between critical phosphorus (P) and nitrogen (N) atoms in DEPE molecules was analyzed by computing radial density distribution function (rdf) for all possible atom-atom combinations. As shown in Fig. 5 , all DEPE:FA systems exhibit an increase in the value of the maximum of the N-N function with regard to DEPE simulation without FAs, due to the fact that the negative charge of the FAs allows ethanolamine nitrogens to be closer to each other. This property indicates that FAs alter the lateral order of DEPE bilayers, which affects the hexagonal H II phase transition. Furthermore, it seems reasonable to hypothesize that a larger effect could be observed for larger DEPE:FA ratios. Conversely, the computed rdf functions for P-P and P-N pairs exhibit a peak centered at 0.53 nm and at 0.40 nm, respectively, in all systems analyzed independently of the presence of different FAs (data not shown), indicating no differences from one system to another. . On the other hand, an average area of 0.517 nm 2 for free-FA system was found, in reasonable agreement with the available experimental measurements ( 29 ). Compared with recent simulations of DOPC/ FA systems ( 25 ), the present results indicate that the stronger interactions between PE groups give a smaller decrease in the area per lipid in presence of FAs. Thus, the incorporation of FAs to the lipid bilayer has a small effect on the L ␣ phase order in DEPE:FA systems at 20:1 molar ratio, in agreement with the experimental data reported.   Among OA series, the number of interactions decrease as 2OHOA > 2MOA > OA, whereas in SA series, it follows the order 2MSA > 2OHSA > SA. The results suggest that a bulkier FA head (a C2 substitution in the acyl chain) enhances the number of interactions, especially for the SA series. Moreover, the specifi c trends seem to be related to the different ability of each series to interact with the aliphatic chains of DEPE lipids, which result in a different relative positioning as discussed below. The presence of FA molecules within DEPE lipids reduced the number of DEPE-DEPE H-bonds per molecule from 2.0 to 1.8, considering a maximum interatomic distance of 0.4 nm. The number of DEPE-water H-bonds also decreased from 5.7 up to 4.7 interactions per phospholipid, considering the same cutoff owing to the different localization of the FA carboxyl group in the polar region of the bilayer. The number of H-bonds was 1-2 units smaller than in DOPC/ FA systems ( 25 ) due to the existence of H-bonds between PE groups. The presence of the FA hydroxyl group gives a small increase in the number of direct interactions with DEPE lipids (10%) but allows the possibility to form watermediated interactions, providing an additional increase in the number of hydrogen bonds with water molecules.
Electrostatic features of the lipid membrane. Table 3 displays the computed electrostatic potentials (measured from the bilayer center to the water interface) for all simulated systems. The values were obtained by integrating charge densities twice according to the one-dimensional Poisson equation ( 36 ). A DEPE bilayer without FAs exhibits a potential of Ϫ 795 mV, which originates mainly from the macroscopical charge distribution created from the sum of individual phosphorus-nitrogen (P-N) dipoles. It is reason-nm ( Fig. 6A ). This pattern suggests the existence of at least two groups of interactions, one closer than the other. Remarkably, the area under the peak is smaller for 2OHOA than for the remaining FAs of the OA series, indicating that the hydroxyl group diminishes the magnitude of the correlation.
The C(COO) FA -P DEPE rdf ( Fig. 6B ) exhibits a rough profi le with a maximum at about 0.6 nm for all the systems. Since the carboxyl and phosphate groups are both negatively charged, this effect must be considered a consequence of the strong correlation with the neighboring N DEPE atom. However, it can be seen that DEPE:2OHOA system exhibits a tail at short distances, refl ecting the existence of additional interactions between the 2-hydroxyl groups with DEPE phosphate oxygens. The same effect is observed for the DEPE:2OHSA system (data not shown). Regarding the C(COO) FA -C(O) DEPE rdf ( Fig. 6C ), there is a maximum correlation at 0.45 nm both for the sn-1 and sn-2 chains. For all systems except 2OHOA and 2OHSA, the correlation is larger for the sn-1 chain than for the sn-2. The larger difference is observed for 2MOA and 2MSA, pointing to the preference of FAs to interact with DEPE ethanolamine groups from the side of the sn-1 chain, and more importantly, that this can be modulated by the substitution at position C2.
An analysis of hydrogen bond interactions between FA carboxyl and DEPE ethanolamine groups was also undertaken ( Table 4 ). For 2OHOA and 2OHSA, the hydrogen binding hydroxyl group was also considered. We identifi ed a large network of H-bond and charge interactions in the interfacial region of the membrane involving DEPE lipids, FA, water molecules, and ions. They provide strong interactions that are probably responsible for stabilizing the bilayer in the L ␣ phase. Each FA forms roughly two H-bonds with DEPE donors/acceptors, with the specifi c numbers ranging between 1.8 and 2.5 depending on the FA type, considering a cutoff of 0.3 nm. Excluding C2-substituted FAs, the number of interactions decreases as SA > OA > PoA > EA > PeA > cVcc, indicating a dependence on Bilayer thickness was computed from the average distance between phosphate planes (std. dev. < ± 0.02). The area per lipid (SD < ± 0.007) was computed as in reference (35). The variation of the electrostatic potential was calculated taking DEPE electrostatic potential as reference (795 mv). 2MOA, 2-methyloleic acid; 2MSA, 2-methylstearic acid; 2OHOA, 2-hydroxyoleic acid; 2OHSA, 2-hydroxystearic acid; cVcc, cisvaccenic acid; DEPE, 1,2-dielaidoyl -sn -glycero -3-phosphoethanolamine; EA, elaidic acid; FA, fatty acid; OA, oleic acid; PeA, palmitelaidic acid; PoA, palmitoleic acid; SA, stearic acid.  Table 2 , together with the relative distances to phosphorus and nitrogen DEPE atoms. The results shown are qualitatively similar to the remaining systems studied. It can be seen that the density of DEPE phosphorus atoms is maximal at 2.4 nm whereas the density for nitrogens exhibits the maximum at 2.3 nm. This generates a dipole whose (time and system) average orientation is 6.5° measured from the bilayer plane and facing toward the bilayer in a system without FAs. This value is in good agreement with previous results of PE bilayers ( 30 ) and differs from those of PCs, where able to expect various contributions to this potential including the P-N dipole, but also the degree of water polarization and the presence of bound ions. In DEPE:FA systems there is an additional contribution due to the FA carboxyl group. The total electrostatic potential computed for DEPE:FA systems are the same or smaller than that computed for DEPE alone (with a maximum decrease of 74 mV). The specifi c trends that can be observed in the table are not fully understood, although they indicate the existence of distinct charge distributions of the different system components within the interfacial region of the membrane that tend to compensate each other. Figure 7 displays the number density distributions for atoms of the functional groups in DEPE:2OHOA and also bound to FAs, although the number of both FAs and sodium ions are two small for a quantitative analysis. Remarkably, the binding of ions to the DEPE bilayer without FAs is negligible. Finally, FAs provide a charge density with a single minimum at ‫ف‬ 2.2 nm (see Table 2 ), close to the maxima distance shown by the solvent, sodium ions, and DEPE lipids.

DISCUSSION
The strategy by which the physiology of a cell can be modulated by drugs that interact with membrane lipids or lipid structures has recently been named 'membrane lipid therapy' ( 8 ). The modulation of the structural and functional properties of cell membranes is a novel pharmacological approach whose molecular basis has not been considered previously but which may be valuable in a variety of human diseases ( 38 ). In this context, elucidating the specifi c interactions of FAs with biomembranes represents the fi rst step toward a design of lipids with potential pharmacological activity. Our study provides a comparative insight into the structural features that alter the molecular location and organization of FAs within DEPE membranes and are able to affect their structural properties.
X-ray diffraction experiments reveal that the main effect of FAs (structurally related to OA) on DEPE membranes is to promote the formation of H II , a nonlamellar phase. In addition, FAs permitted L ␣ and H II structures to coexist in a wider range of temperatures. The repeat distance of L ␣ phases did not appear to be sensitive to the presence of structurally related FAs, differing only in the stereochemistry of the double bond, its position in the acyl chain or the chain length (e.g., EA, PoA, PeA, cVcc). When considering these experimental data in conjunction with data from previous studies ( 13,14 ), we conclude that the ability of the FAs to promote H II phase formation in DEPE membranes follows the order 2MOA > OA > 2OHOA, EA, PoA, PeA, cVcc, 2MSA > SA > 2OHSA. To explain the different abilities of FAs to modulate the membrane thermotropic properties, it is necessary to evaluate some structural features, including the degree of FA ionization, the localization of FAs within the membrane, the interactions between FAs and DEPE lipids and water molecules, and possible ordering effects of FAs on bilayers in the L ␣ phase. The apparent pKa of FAs is about 7 when incorporated in the bilayer of a phospholipid vesicle. Thus, at physiological pH, FAs could be partially ionized, although the pKa varies depending upon the environment ( 33,34 ). In previous reports we analyzed the effect of having a net negative charge in 2OHOA and OA within a DEPE matrix, concluding that the extent of FA ionization affected the surface charge density and modulated the number of H-bond interactions in the polar region of the membrane ( 15,39 ). In the present report, the remaining structural features are examined by means of a set of MD trajectories of various DEPE:FA systems simulated in the L ␣ phase.
For all FAs, the analysis of MD simulations revealed small or negligible effects both in membrane thickness the P-N angle is about 15° in the opposite direction, that is, pointing toward the water face ( 22 ). DEPE:FA systems exhibit almost the same value for the P-N angle as the DEPE bilayer with the small increase of 0.8-2.6° that can be attributed to a shift in the N DEPE density toward lower distances. This suggests that the extensive network of H-bonds between DEPE molecules allows only a small response of the phospholipid head-group to the negatively charge densities provided by FAs. These results indicate that PE membranes are much poorer charge sensors than PC bilayers ( 22,37 ).
DEPE carbonyl oxygens are located in a much deeper position (2.0-2.1 nm) with the group at the sn-2 chain always at a smaller distance than that at sn-1. The maximum carboxyl carbon density value corresponds to 2.15 and 2.25 nm for OA and 2OHOA, respectively. Interestingly, the carboxyl groups in the OA system lie close to carbonyl oxygens and overlap with half of the nitrogen density. On the contrary, in the 2OHOA system they are located between carbonyl and phosphate oxygens, almost at the same position than the maximum in the nitrogen density. Moreover, the width of the phosphorus density is slightly smaller for 2OHOA than for OA. The total DEPE charge density shows a rough profi le with three maxima and two minima in the OA system and a fl atter profi le with only one minimum in the case of 2OHOA. It is known from previous studies that the lipid charge distributions in a bilayer are almost compensated by the remaining system components in the interfacial region ( 22 ), in this case, water molecules, FAs, and ions. Indeed, water molecules orient cooperatively generating a bimodal charge density, with a maximum at ‫ف‬ 2.25 nm and a minimum at ‫ف‬ 2.95 nm. Below 3.3 nm, the density of Na + ions is larger than the density of Cl Ϫ , and therefore, they provide an additional contribution to the positive charge density. The Na + density ranges 0.06-0.08 ± 0.02 ions bound per lipid. Some of these ions of the membrane, whereas a methyl group (2MOA) moves them toward the inner part. On the other hand, the presence of a hydroxyl group provides only a 10% increase in the number of direct interactions with DEPE lipids and gives an additional increase in the number of H-bonds with water molecules. The structural data here presented indicated that the main effect of the C2-substitution is the modulation of the FA density distributions rather than the extent of the H-bond network within the interfacial region of the membrane. Considering the macroscopic "molecular shape model" ( 9 ) and assuming a similar molecular geometry for all DEPE:FA systems studied, the primary reason for the 2OHOA, OA, and 2MOA differential effect on the thermotropic behavior of DEPE membranes could be related to the FA effect on the headgroup area at the lipid-water interface, because the acyl chain packing would be similar for all systems. The shift of 2OHOA toward the outer side of the membrane will be associated with a higher optimum headgroup area that could account for the appearance of the L ␣ -to-H II phase transition at higher temperatures and the stabilization of the L ␣ phase, as observed experimentally.
The present study demonstrates that changes in the chemical structure of FAs affect both FA-DEPE and DEPE-DEPE interactions. These are important results, because and in the area per DEPE molecule with the presence of FAs in a 20:1 molar ratio, in agreement with the preceding X-ray diffraction experiments. On the contrary, the average density distributions showed the existence of different profi les both for the density of carboxyl carbons and for the acyl chain FA molecule in the axis perpendicular to the membrane plane. These differences are found both in the relative localization and amplitude of the distribution of the carboxyl group within the polar region of the membrane and in the degree of membrane penetration of the FA acyl chain. The chemical group at position C2 strongly modulates the former, whereas the nature of the acyl chain affects mainly the latter. The number of H-bonds, either with DEPE lipids or water molecules, was sensitive to the localization of the FA carboxyl group in the DEPE membrane. C2-substituted FAs exhibit the largest number of DEPE-FA H-bonds as the C2-group provides a larger head group area and, therefore, a more accessible carboxyl group.
Due to the importance of 2OHOA as a hypotensive drug ( 15 ), characterizing the differences between DEPE:OA, DEPE:2OHOA, and DEPE:2MOA systems at a molecular level is of large interest. Considering the profi le exhibited by OA as a reference, a hydroxyl group at the C2 position (2OHOA) shifts FA atomic densities toward the outer side not all structural changes in the FA molecules were equally effective in the infl uence on the membrane structural properties. When trying to extrapolate these fi ndings on biological systems, the results of in vivo studies are indicating that the above-mentioned extent of FA-membrane interactions have a correlation to their physiological effi cacy. For example, a recent study in rats has demonstrated that saturated SA or trans -monounsatured EA were less effective in reducing blood pressure than the monounsaturated OA ( 17 ). In turn, 2OHOA showed a highly increased antihypertensive action in genetically hypertensive rats ( 18 ). Moreover, reconstituted plasma membranes from these rats showed a clear shift of their biophysical properties toward a higher propensity to form H II phases, paralleled by a signifi cant increase in the level of unsaturated FAs ( 16 ). Although these results do not fi rmly establish the molecular mechanism of action of 2OHOA, they indicate that the observed alteration of biophysical membrane parameters in vitro also occurs in vivo and should have pharmacological relevance. In addition, a similar observation regarding body weight homeostasis could be made because neither SA nor EA were able to reduce body weight in rats, while OA showed a slight and 2OHOA a strong lipolytic effect on adipose tissue mass ( 40 ).
Summing up, both dietary and synthetic FAs are fi nally incorporated, at least partially, into cell membranes either as free FAs or in the form of phospholipids. FAs able to affect biophysical properties of cell membranes, in turn will also alter localization and/or function of membrane protein involved in the regulation of cellular processes. For the above reasons, FAs or other lipids could be a tool to modulate pathophysiological conditions via cell membrane properties, a pharmacological target that has not been explored yet.