Dmitry V. Leonova,
Sergei A. Dzuba*ab and
Nikolay V. Surovtsevac
aDepartment of Physics, Novosibirsk State University, Novosibirsk, 630090, Russia. E-mail: dzuba@kinetics.nsc.ru
bVoevodsky Institute of Chemical Kinetics and Combustion, Russian Academy of Sciences, Novosibirsk, 630090, Russia
cInstitute of Automation and Electrometry, Russian Academy of Sciences, Novosibirsk, 630090, Russia
First published on 25th October 2019
A lipid bilayer containing a ternary mixture of low- and high-melting lipids and cholesterol (Chol) can give rise to domain formation, referred to as lipid rafts. Low-frequency Raman spectroscopy at reduced temperatures allows detection of normal membrane mechanical vibrations. In this work, Raman spectra were obtained in the spectral range between 5 and 90 cm−1 for bilayers prepared from dioleoyl-glycero-phosphocholine (DOPC), dipalmitoyl-glycero-phosphocholine (DPPC) and Chol. A narrow peak detected between 13 and 16 cm−1 was attributed to the vibrational eigenmode of a lipid monolayer (a leaflet). For the equimolar DOPC/DPPC ratio, the Chol concentration dependence for the peak position, width and amplitude may be divided into three distinct ranges: below 9 mol%, the intermediate range between 9 mol% and 38 mol%, and above 38 mol%. In the intermediate range the peak position attains its minimum, and the peak width drops approximately by a factor of two as compared with the Chol-free bilayers. Meanwhile, this range is known for raft formation in a fluid state. The obtained results may be interpreted as evidence that bilayer structures in the raft-containing fluid state may be frozen at low temperatures. The drop of peak width indicates that at the spatial scale of the experiment (∼2.5 nm) the intermolecular bilayer structure with raft formation becomes more homogeneous and more cohesive.
Vibrational frequencies of specific functional groups obtained in Raman scattering experiment15–20 are sensitive to the local environment of lipids and proteins, providing so information on the lipid order. Recently, low-frequency Raman scattering spectra23,24 were found to be capable for detection of THz acoustic vibrations of the bilayer as a whole. The narrow peak between 9 and 18 cm−1 detected below 250 K was attributed to the eigenmode of normal mechanical vibrations of a lipid monolayer (a leaflet). This peak cannot be observed above 250 K because of mode overdamping by fast relaxation motions. The peak position and width were found to depend on the lipid content; the width decrease with increase of the Chol content observed in binary lipid/Chol bilayers was attributed to increase of the bilayer homogeneity in the spatial scale of the sound wave half-length (that is the monolayer thickness, i.e. ∼2.5 nm).24
So low-frequency Raman scattering provides principally new information on the bilayer normal mechanical vibrations, which may appear useful for the elucidation of the bilayer nanoscale structure. Although these vibrations can be detected only at reduced temperatures, one may suppose that the bilayer nanoscale structure may be frozen with temperature decrease, – like it is assumed in different experimental approaches employing structure cryogenic trapping.
Note that the in-plane mechanical vibrations of the membranes were studied employing inelastic X-ray scattering,25 inelastic neutron scattering26 and molecular dynamics (MD) simulations.27,28 The data obtained have demonstrated that lipid bilayers exhibit collective short-wavelength dynamics showing correlated density fluctuations on the sub-ps time scale; in the presence of Chol these motions are stiffening.
Here, we apply low-frequency Raman scattering to the ternary lipid mixture of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), 1,2-dipoleoyl-sn-glycero-3-phosphocholine (DPPC) and Chol. This system is known for coexistence of liquid disordered (Ld) and gel (g) phases, liquid ordered (Lo), Ld and g phases, Lo and Ld phases,1–4,10,14 which takes place for certain lipid and Chol compositions, – see Fig. 1 where phase diagram retrieved from the literature data2 is presented. These phases are characterized by different molecular ordering and form domain structure with lipid composition differed from the mean composition; the Lo domains coexisting with other phases are lipid rafts. As formation of these domains may influence the lipid layer vibrational eigenmodes, the goal of this work is measurement of low-frequency Raman spectra of the bilayers of different lipid/Chol composition.
![]() | ||
Fig. 1 Phase diagram of the DOPC/DPPC/Chol ternary system at 18 °C retrieved from literature data.2 The vertical line corresponds to the compositions studied in this work. The grey triangle with vertexes A, B and C corresponds to the compositions where the Ld, Lo and g phases are present; see text for meaning of the A′ and B′ points. |
![]() | (1) |
Spectra in Fig. 2 contain several peaks, which correspond to vibrational eigenmodes of the lipid bilayer.23,24 Except of the intensive central peak which can be tentatively ascribed to the relaxation contribution, one can see three resolved peaks: near 14–17 cm−1, near 30 cm−1, and near 50–60 cm−1.
The well-defined peak between 13 cm−1 and 16 cm−1 corresponds to the first vibrational eigenmode of the lipid monolayer.23,24 In the continuum description of the elastic response of the layer, its frequency is described by the relation:
![]() | (2) |
The insert to Fig. 2 shows that the experimental Raman spectra can be decomposed onto the sum of four Lorentzians, taking also into account the contribution from the central peak. The dashed line here corresponds to the central peak; the dotted line reflects contribution from acoustic-like excitations. From this decomposition the position and width of the 13–16 cm−1 peak were evaluated with the precision of ±0.2 cm−1.
Note that one may suppose that upon freezing the macroscopic phase separation may occur, with DOPC, DPPC and Chol ingredients forming distinct different phases. This possibility however may be certainly ruled out by comparison with low-frequency Raman spectra for bilayers of pure DPPC lipids:23 the latter contains below 120 K a well-distinguished peak at 9.3 cm−1 (see Fig. 1 therein), which is not observed in our case of the DOPC/DPPC/Chol mixtures – see Fig. 2. (This peak was ascribed in ref. 23 to normal vibrations of the DPPC bilayer as a whole). So the absence of 9.3 cm−1 peak in Fig. 2 unambiguously evidences that under our experimental conditions phase separation into macroscopic pure DOPC, DPPC and Chol phases does not take place. Also, one more argument for the absence of phase separation is worth to be added: the Raman spectra obtained for ratio of the peaks corresponding to C–H2 symmetric and antisymmetric vibrations in DPPC lipids show that in pure DPPC and in DPPC/Chol bilayers this ratio is remarkably different31 (see Fig. 6a therein) which would not occur if DPPC and Chol are separated into different phases.
In this work, we explore Chol concentration dependence for the compositions with equimolar DPPC and DOPC quantities. The found Raman peak positions and widths at different Chol concentrations are shown in Fig. 3. The data were obtained at two temperatures, 100 K and 170 K. One can see that for both temperatures the Chol concentration dependences are non-monotonic. For the peak position the dependences contain two kinks separating three ranges: below 9 mol%, between 9 mol% and 38 mol%, and above 38 mol%. The concentration dependence of the peak width contains only one kink near ∼ 9–10 mol% Chol, which separates two scenarios for this dependence: sharp decrease from 10 cm−1 to 5 cm−1 and much slower decrease at higher Chol concentration.
![]() | ||
Fig. 3 The peak position (A) and the peak width (B) versus Chol concentration at 100 K (black squares) and 170 K (red triangles). The solid straight-line segments are drawn to guide the eye, and the dotted lines are predictions of eqn (3) (see text). |
Chol dependence of the peak amplitude at 100 K is shown in Fig. 4. (Data at 170 K displayed high scattering due to the higher contribution of the central peak). From Fig. 4 one can see an indication to two kinks near 10 mol% and near 40 mol% Chol concentration.
![]() | ||
Fig. 4 The peak amplitude versus Chol concentration at 100 K. Three solid straight-line segments are drawn to guide the eye, and the dashed lines are predictions of eqn (3) (see text). |
The concentration dependence of the peak width (Fig. 3B) demonstrates a sharp decrease with the Chol concentration increase – almost by a factor of 2 – when approaching the Chol concentration of 9 mol%. Since the peak narrowing implies lower matter defectiveness for the sound propagation, this observation indicates on higher homogeneity of the bilayer above Chol concentration of 9 mol%. And this increase of homogeneity occurs just for the Chol concentration range where rafts are formed at 18 °C – see Fig. 1.
At first look, the conclusion about the increasing membrane homogeneity is paradoxical. Indeed, coexistence of different phases (Lo, Ld, g) just implies system heterogeneity. However the heterogeneity-homogeneity concept depends on the spatial scale in which the system is considered. The spatial scale of the described low-frequency Raman experiment corresponds to the half of the wavelength that is the monolayer thickness, d ∼2.5 nm. So we may admit that for Chol concentrations larger than 9 mol%, the bilayer becomes more homogeneous in the nanometer spatial scale. This increased homogeneity implies more compact intermolecular packing and mechanically more cohesive intermolecular structures. So, lipid rafts may homogenize bilayer structure at the nanoscale distance range.32
Data in Fig. 4 show that the peak amplitude becomes somewhat larger for the Chol concentration range between ∼10 mol% and 40 mol%. As this range is similar to that for the raft formation at 18 °C (see Fig. 1), this coincidence again shows correlation between mechanical vibrations of frozen bilayers and the raft formation in the fluid state. According to the previously suggested interpretation,23,24 the peak increase reflects a reduction of the interleaflet coupling in the bilayers. So the peak increase in Fig. 4 may imply that raft formation reduces this coupling.
One may describe properties of the compositions inside the coexisting range employing those for the individual domains. For this purpose, the composition may be considered as a sum of the constituting domains those weights can be found from the tie lines in the phase diagram and the lever rule.1,2,33–35 For the linear borders encompassing the Lo + Lo + g phase coexistence triangle (see Fig. 1), the Raman spectrum f(ν) at any point inside the coexisting range may be presented as a linear superposition of the spectra for the A, B and C compositions corresponding to the three vertexes in Fig. 1:
f(ν) = afA(ν) + bfB(ν) + cfC(ν), | (3) |
The fA(ν), fB(ν) and fC(ν) spectra were obtained as Lorentzian approximation of the monolayer eigenmodes for the corresponding bilayer composition – see Fig. 5 (left). Then, denoting the DOPC/DPPC/Chol composition for any point O within the triangle as xO, yO, zO in the phase diagram depicted in Fig. 1 (xO + yO + zO = 1), we write
![]() | (4) |
![]() | ||
Fig. 5 Left: The circles show experimental Raman spectra for the DOPC/DPPC/Chol bilayers taken at 100 K for the compositions A, B, and C of the triangle vertexes shown in Fig. 1. Data are shifted upwards for convenience. The solid lines drawn through the experimental points are fits by a sum of Loretzians (cf. Fig. 2). The insert shows the Lorentzians describing the peak between 13 and 16 cm−1, for the A, B and C compositions. Right: the peak weights a, b and c found from eqn (4), as a function of Chol concentration. |
For the Chol concentration above 23 mol% and below 38 mol% only two phases, Lo and Ld, are present (see Fig. 1). The tie-line analysis performed for the DOPC/DPPC/Chol bilayers34,35 has shown that the tie-line endpoints A′ and B′ for the Lo + Ld phase coexistence are located like it is shown in Fig. 1. In our consideration we neglected the difference between the points A′ and A and between the points B′ and B in the phase diagram in Fig. 1, and took c = 0 in eqn (4).
The peak weights a, b and c obtained using eqn (4) in the case of the compositions studied here, xO = yO = (1 − zO)/2, are shown in Fig. 5 (right) as a function of Chol concentration (zO). Then for any point O position the Raman spectra were simulated using eqn (3). The resulting f(ν)Raman line was found to be fitted fairly well by a single Lorentzian (data not shown), and for this Lorentzian the peak position, width, and amplitude were determined. These predictions for the peak position and width are given in Fig. 3 by dotted lines as a function of Chol content. Also, in Fig. 4 the predicted peak amplitudes are also shown by the dashed line. One can see in Fig. 3 and 4 a good agreement with the experiment.
Raman spectra were studied for the compositions with the equimolar DOPC/DPPC ratio and Chol concentration varied. The found peak position possessed two kinks at the Chol concentration dependence, located at 9 mol% and 38 mol% Chol concentrations. These values are close to the borders encompassing at fluid state the Lo + Ld + g and Lo + Ld phase coexistence (the raft formation). This similarity allows to suggest that peculiarities of the bilayer structure at fluid state become frozen when temperature decreases below 170 K.
The found peak width decreases with Chol concentration increase and attains a plateau-like behavior above 9 mol% of Chol concentration, with the magnitude approximately twice smaller than that for the Chol-free bilayer. This evidences a more homogeneous (in the spatial scale of the experiment, ∼2.5 nm) and more cohesive intermolecular bilayer structure in the Chol concentration range in which rafts are formed at fluid state.
The peak amplitude was found to increase between 9 mol% and 38 mol% Chol concentration borders, which implies that in this concentration range the interleaflet coupling in the bilayer becomes smaller.
This journal is © The Royal Society of Chemistry 2019 |