Open access peer-reviewed chapter

Application of Nuclear Magnetic Resonance Spectroscopy (NMR) to Study the Properties of Liposomes

By Anna Timoszyk

Submitted: November 24th 2016Reviewed: March 13th 2017Published: October 25th 2017

DOI: 10.5772/intechopen.68522

Downloaded: 906

Abstract

The liposomes are well‐known lipid aggregates. The lipid composition and size of the liposomes can be controlled. The method of preparation, lipid composition, temperature, and pH have an influence on the liposome size and bilayer structure. The physicochemical properties of liposomes allow them to various applications. Nuclear magnetic resonance (NMR) is one of the methods used to study liposome properties. The abilities of the method are the high sensitivity and high resolution. Moreover, it provides information about dynamics and structure of molecules. 1H and 31P NMR are most convenient methods to study liposomes, because liposomes are typically formed from phospholipids. Additionally, two‐dimensional NMR spectroscopy reveals information about the nature of intermolecular and intramolecular interactions (scalar and dipole‐dipole interactions) that makes easier to interpret the structure of molecules. The chapter aims to introduce the NMR phenomenon, interactions between spins in magnetic field, dynamics of molecules and physical parameters of NMR spectra, and the necessary information for analyzing and interpreting high‐resolution NMR spectra. It also aims to show how various changes in the bilayer structure or dynamics of lipid molecules are visible in the NMR spectra.

Keywords

  • liposomes
  • nuclear magnetic resonance
  • dynamics
  • half‐width of signal
  • signal splitting
  • 1H NMR
  • 31P NMR
  • 2D NMR

1. Introduction

Nuclear magnetic resonance (NMR) spectroscopy provides information about the structure and dynamics at the molecular level. The knowledge about the peculiarity of NMR phenomenon and about the physical parameters of NMR spectra would help in analyzing obtained results.

1.1. The magnetic properties and magnetic moment of nuclei

Studies of the atomic structure of discrete spectra have shown that, similar to electrons, nuclei have momentum, which is called nuclear spin. Nuclear spin is directly connected to a given magnetic moment. The values of these magnetic moments can vary. The magnetic moment μis proportional to momentum Κ(spin) [1]:

μ=γΚE1

where γ is the gyromagnetic ratio, characteristic for a given nucleus (Table 1)

NucleiNuclear spinGyromagnetic ratio [T‐1.s‐1]Resonance frequency in B0 = 14.092 T [MHz]
1H1/22.6752 × 108600
13C1/26.7266 × 107150.9
15N1/2−2.7108 × 10760.8
17O5/2−3.6267 × 10881.4
19F1/22.5167 × 108564.5
25Mg5/2−1.6371 × 10736.7
31P1/21.0829 × 108243.9
33S3/22.0518 × 10746
35Cl3/22.6213 × 10758.8
39K3/21.2482 × 10728

Table 1.

The values of nuclear spin (Κ), gyromagnetic ratio (γ), and resonance frequency (ω0) [1, 3].

1.2. Resonance condition

A given nucleus with its magnetic moment and spin will precess in intense static magnetic field B0 with frequency ω0:

ω0=2πν.E2

The value of the frequency is proportional to γ and to the B0 field [1]:

ω0=γB0.E3

The frequency of precession in the B0 field is known as Larmor frequency. In precession, the nuclear dipoles move around a cone in the B0 field at frequency ω0. The precessing nuclei transverse magnetic field B1, which is a linearly oscillating magnetic field along the x‐axis, and they then submit to the combined action of both the B0 and B1 fields [1]. If B1 field oscillates at the Larmor frequency, the nuclear magnetic resonance phenomenon will be observed, according to resonance condition as follows:

ω1=ω0=γB0.E4

1.3. Nuclear magnetic relaxation

The nuclear magnetic relaxation phenomenon is connected to the interaction between the magnetization vector M(equilibrium magnetization) and the B0 field. The Mvector is the sum of the magnetic moments of the given nuclei:

M=i=1nμi.E5

The value of magnetization is equal to zero when the nuclei are not in the B0 field. Then, the μmoments of the nuclei are oriented chaotically in accordance with statistical distribution. In the B0 field, the μvectors are ordered, which results in thermodynamic equilibrium. The time in which the equilibrium is set depends on the type of sample and the temperature [2]. In the equilibrium state, more μmoments are oriented parallel to the B0 field than are antiparallel to it. It is in agreement with the normal Boltzmann equilibrium between the spin states. Thus, the Mvector is also oriented along the direction of the B0 field (Figure 1).

Figure 1.

Equilibrium magnetization M → in B0 field.

The μmoments are not coherent in precession; thus, there is no gain in the transverse magnetization Mxy. When the B1 field interacts with the precessing spins at the Larmor frequency, the Mxyvector gains (Figure 2) [2].

Figure 2.

The M → vector emerged at angle α with respect to the z‐axis and the transverse M x y → .

After switching off the B1 field at time τ, the Mvector will emerge at angle α = γB1 τ with respect to the z‐axis. Then, the Mvector begins to precess around the direction of B0 field, which is now active [1]. Thus, the transverse Mxywill lose at time T2 (transverse relaxation time or spin‐spin relaxation time) (Figure 3).

Figure 3.

The loss of transverse M x y → vector at time T2.

1.3.1. Spin‐lattice relaxation process

The orientation of the μin magnetic field B0 depends on the value of the interaction energy between μand B0 field [3]:

E=μB0=mγB0,E6

where m is the magnetic quantum number.

In the equilibrium, the nuclei at every energetic level are in accordance with the normal Boltzmann equilibrium law [3]:

N1N2=expΔERT1+γB0kTL,E7

where N1 is the nucleus population on the energy state, which is less than ΔE energy in comparison to the N2 population, TL is the sample temperature, and =h/2πwhere h is Planck constant.

The B1 field induces the transition between quantum states, which gives rise to the NMR signal. The equalization of Boltzmann populations (extinction of NMR signal) in the quantum states, known as saturation phenomenon, is accompanied by an increase in the spin temperature TS. At the same moment, the sample temperature TL does not change. The condition in which TS > TLmeans that the equilibrium state has not been archived; however, after switching off the B1 field, the TS temperature will equalize with TL (lattice) temperature. The spin existing in the excited quantum state cannot return to the ground state spontaneously. Thus, forced emission is one way that the spins lose energy [4]. This type of emission is possible only throughout spin‐lattice interactions. The transitions between the spin states can be enforced only via the local magnetic field Bloc at the Larmor frequency. Fluctuations of Bloc field are generated by the thermal motion of the atoms and molecules from which a network is formed [1, 2]. Each magnetic moment that participates in random rotational or translational Brownian motion causes a fluctuation of the Bloc field, which causes a Fourier's frequency spectrum for those fields. For the spin ½ nuclei 1H, 13C, and 31P, the dominant mechanism of relaxation is a dipole‐dipole interaction, which is a result of the interaction between μand Bloc that is generated by neighboring magnetic moments. Thus, for the spin ½ nuclei, the dominant mechanism of relaxation is the interaction between the gradient of the electric field generated at the location of the observed nucleus by its electrical surroundings and nuclear electric quadrupole moment [46].

The simplest way to describe dipole‐dipole interaction is by using the system of two spins. If the spin I (μI)is near the second spin S (μS), then the Bloc field created by the nucleus S at the position of nucleus I is equal to [13]

Bloc=±μSRIS(3cos2θ1),E8

where θ is the angle between the B0 field and the RISvector, and RISvector is the distance between spin I and spin S (Figure 4).

Figure 4.

The system of two spins in magnetic field.

If spins S and I are from the same molecule, then the RISis constant, and fluctuations in the Bloc will be associated with random changes in the angle. This mechanism of relaxation is known as rotational or intermolecular relaxation. If spin I and spin S are not from the same molecule, then the Bloc field fluctuations will depend on changes in both angle and the length of the RIS[2]. This is called as translational or intramolecular mechanism of relaxation, which primarily relates to a liquid. In fact, in the spin I position, the Bloc field is produced by more than one spin S; however, if a system has more than two spins, the calculation of relaxation time rate is limited and the additivity of the effects is assumed. Then, the rate of the spin‐lattice relaxation time T1 can be expressed as the sum of the probabilities of the transitions between spin energy levels [6]:

1T1=MM'WMM'(EMEM')2MEM2,E9

where WMM'is the probability of the transition between energy level EM and level EM'.

In the case in which the two protons will be taken into consideration, for example, in an H2O molecule, then the rate of T1 can be summarized as follows [2]:

1T1=2(W1+W2),E10

where w1 is the probability of a single spin‐flip transition, and w2 is the probability of a double spin‐flip transition.

In a system of AX spins, for example, in a CH group, the cross‐relaxation effect occurs. In such a case, it is necessary to use the proton‐decoupling method to prevent the splitting of signals (multiplets) caused by the spin‐spin couplings. Then, the rate of T1 can be expressed by the formula [2]:

1T1=W0+2W1+W2,E11

where w0, w1, and w2 represent the probabilities of the two‐spin system (Figure 5).

Figure 5.

The unlimited transitions in two spin system.

After determining the transition probabilities and the perturbation Hamiltonian, it is possible to write the T1 relaxation time for the homonuclear spins (II) and heteronuclear spins (IS) as follows [5]:

1T1=310γI42RII6[τc1+ωI2τc2+4τc1+4ωI2τc2]E12

and

1T1=110γI2γS22RIS6[τc1+(ωSωI)2τc2+3τc1+ωI2τc2+6τc1+(ωS+ωI)2τc2],E13

where τc is the correlation time.

The time T1 is very sensitive to the length of RISbetween near spins, as directly proportional to its sixth power [7].

The high‐temperature approximation is fulfilled, ω02τc21, for liquids, thus the rate of T1 is independent of the frequency and Eqs. (12) and (13) are simplified to the following forms [2]:

1T1=32γI42RII6τcE14

and

1T1=32γI2γS22RII6τc.E15

As a result of the high‐temperature approximation, the rate of T1 increases as the temperature increases.

In general, the total rate of T1 for protons is the sum of the rates of relaxation times caused by intermolecular and intramolecular factors [2]:

1T1=(1T1)intra+(1T1)inter,E16

where the rate of relaxation time, denoted as intramolecular (1/T1)intra, fulfills relation (14).

The calculation of the rate of relaxation time caused by the intermolecular factors was presented in reference [8]. The (1/T1)inter is directly proportional to the population of spins N1, and it is inversely proportional to the translational diffusion coefficient D:

(1T1)inter=1730πγI4N1aD,E17

where a is the closest possible distance between two spins belonging to two molecules (a is usually equal to the particle diameter).

Thus, D can be obtained directly while measuring viscosity based on the Stokes‐Einstein relationship:

D=6πηa.E18

The τc obtained by measuring the rates of T1 for each chemical group in the molecule should have the same value when the molecule is subjected to isotropic rotation. In fact, various correlation time values can be observed, and this proves the existence of internal motions in the molecule [46].

1.3.2. The Lipari‐Szabo model‐free approach

In liquids, molecules are subject to rotational motion around the symmetry axis, and individual groups of molecules demonstrate internal movement. A correlation function for this complex motion of the RISvector was derived in reference [9]. The model assumed that the total rate of the rotational movements is equal to the sum of the rates of the entire molecule (overall rotation—isotropic rotation) and internal rotations:

1τi=1τiall+1τiinter,E19

where τiallis the time of isotropic rotation of the entire molecule, and τiinteris the time of the internal isotropic rotation.

According to this model, the rate of T1 is dependent on the complex motion coefficient C1[9]:

C1=[134(sin4β+sin22β)],E20

where β is the angle between the axis of rotation and the vector RIS.

The coefficient C1 can take any value between 0 and 1 (0> C1> 1). When C1 = 1, there is no internal movement, when C1 <1, the vector RISis subject to the movement [9].

1.3.3. Spin‐spin relaxation process

The spin‐spin relaxation process is associated with losing the phase coherence via the nuclear spin system, which leads to the loss of the Mxyvector. The effectiveness of the process depends on the rate of the molecular reorientations. In viscous liquids and solids, molecular reorientation is slow, it either takes a few microseconds or it does not occur at all [2]. The strength of magnetic field comes from the μmoments, decreases as the distance between the spins increases. Thus, only the nearest nuclei have a significant contribution to the Blocfield. Therefore, from the position of the individual spins, the strength of B0 field in which they exist differs in the range of Bloc, and it may increase by several Gauss. The consequence of this phenomenon is that the resonance frequency of the spins will also be different [1]. Spins can transfer absorbed energy to other spins that are located at a lower state of energy. The rate of T2 time determines a spin’s lifetime at a given energy state. This phenomenon is associated with the broadening of quantum states, which is explained via the Heisenberg uncertainty principle [3]:

Δεte2.E21

The shorter is the spin’s lifetime, the greater the broadening of the quantum states. Therefore, the transverse relaxation time is related to the width of the resonance signal [1, 3]. The dependency of the half‐width of the resonance signal 1/2 and the relaxation time T2 is

Δν12=12πT2.E22

The values of the T2 relaxation times in liquids are similar to the values of the T1 times and they are relatively long; however, T2 cannot be greater than T1[1].

1.4. Spin‐spin couplings

The interaction between μmoments, located in the magnetic field B0, is known as spin‐spin coupling J (scalar coupling). Each nucleus interacts with every other nucleus through their valence electrons, and the value of the coupling J depends on the gyromagnetic ratio and the distance between the coupling nuclei (number of bonds). The consequence of the interaction (attraction or repulsion of the moments ì) is the splitting of the quantum states (discrete spectrum). Thus, a given multiplet is found on the spectrum and the distance between the split signals is equal to J. In general, if the proton n is coupled with the other protons, then (n + 1) signals are obtained in the multiplet [1]. Thus, the methyl protons (CH3) coupled with the two protons of the CH2 group cause the splitting of the resonance signal on (2 + 1) peaks. In multiplet, the intensity of the signals can differ, and this can be determined using Pascal's triangle. These rules go into effect only if the differences between the chemical shift ranges and the value of spin coupling J are sufficiently large [3].

The spin‐spin coupling J of spin Y and spin X can provide a relaxation mechanism for spin Y if spin X undergoes relaxation with time T1. Spin Y is subjected to a fluctuating field due to the rapid spin X reorientation [1].

1.4.1. Amplifying the signals via polarization transfer

In NMR, the polarization transfer method, or the spin population transfer, is used to amplify the weak resonance signals, for example, in the 13C spectra. The intensity of the signal is directly proportional to the difference between the N1 and N2 spin populations at energy levels. The N1/N2 ratio fulfils the normal Boltzmann equilibrium between the spin states. The greater the B0 field, the greater the difference between the spin populations, and this depends on the ratio of the γ of the spins [10]. The nuclei 1H and 31P have large values, hence, the resonance signals are easier to observe than a signal from 13C nuclei. The use of selective pulses of B1 field might increase signal intensity of the spins in coupled systems. In the case of a system in which a sensitive nuclei A (proton) is coupled with an insensitive nuclei X (13C) in 13C spectra, signals with a coupling value of J(C, H) = 209 Hz (e.g., for 13CHCl3) are obtained [10].

Amplifying the parameter signal p may be represented as follows [10]:

p=1+γAγXp=1γAγX.E23

1.4.2. Nuclear Overhauser effect (NOE)

The nuclear Overhauser effect is a double‐coupled nuclei resonance, which results in a change in the signal intensity (typically an increase). Intramolecular and intermolecular factors have an impact on the value of signal amplification [11]. The theoretical maximum value of the signal amplification, in the case of coupling of 1H and 13C nuclei and dipole‐dipole relaxation, is 2.989. In fact, the value can change from 1 (which means no gain) to a maximum theoretical value [10]. Because dipole‐dipole relaxation is the major relaxation pathway for protons and the carbons that are directly bonded with protons, the intensity of resonance signals corresponding to them will be more than the intensity from other nuclei. Thus, the spectrum will show signals with various intensities [3]. Therefore, the proton‐decoupled resonance technique is used for 13C nuclei. The application of the proton‐decoupling method and the NOE effect makes it possible to amplify the signal by up to 200% [10].

The NOE effect can also be observed in the 1H NMR spectroscopy, when two protons, Hα and Hβ, are not directly bonded, but they are sufficiently close to each other. If the signal from Hα is gained, then the intensity of the signal corresponding to Hβ will also increase about 45% [10]. Two mechanisms are responsible for the amplification of the signal, which causes the transfer of polarization:

  • dipole‐dipole interaction through space and

  • chemical exchange (in the two spins system AX, the nucleus A is polarized and then the polarization is transferred from A to X with the exchange constant k).

1.5. Detection of NMR spectra and NMR spectra parameters

After the discovery of the phenomenon of NMR in solids, it was found that this phenomenon can also be observed by treating the sample with a sequence of short pulses of B1 field at resonance frequency [12]. After that, the sample can induce alternating voltage with frequency ω0 in the coil. An increase in the coil‐alternating current is observed as a signal of free induction decay (FID).

The pulse duration is given by the following formula [1]:

ti=θγB1.E24

The selected pulse rotates Mvector by a given angle θ. The most commonly used pulses are the ones that rotate the Mvector about 90o (π/2) or 180° (π) (Figure 6). The duration time of the pulses ranges from 1 to 100 μs [3].

Figure 6.

The B1 field pulses rotate the M → vector about (B) 90° (π/2) and (C) 180° (π); (A) the equilibrium state; the evaluation of the M → vector and return to the equilibrium state in (D) the case (B) and (E) the case (C).

The shape of the FID signal is a fading oscillation curve as a function of time. This function is archived in the acquisition time of approximately 1 s. Then, using Fourier transformation, the periodic changes in time are converted to the frequency spectrum. The advantage of using this method is that it enables the fast recording of the spectrum as well as accumulation and averaging, which increases the signal‐to‐noise ratio [1].

1.5.1. Chemical shift

The chemical shift is the most common parameter analyzed in the 1H NMR spectra. This is due to the fact that the distribution of electrons in the molecule is varied, so the different chemical groups of the same molecule have a different screening constant [1, 3]. Consequently, the same nuclei require a different B0 field to achieve a resonance condition at a predetermined frequency ω0.In this situation, the resonance condition must take the following form [1]:

Bef(ofsample)=Bst(ofstandard)E25

Thus, the effective field acting on the nucleus must satisfy the following equation [1]:

Bef=B0(1κ),E26

where κ is the screening constant that is a measure of the density of the electron cloud around the nucleus.

Substituting formula (26) to (25), we get

B0(1κ)=Bst(1κst).E27

Assuming, that (1κ)0, because κ is the order of 10‐6–10‐2, so κ1, we then get

BstB0Bstκstκ=σ,E28

where σ is the chemical shift.

However, in the given B0 field the frequency changes ω0, and then σ can be presented in the following form [3]:

σ=ννstνst106[ppm].E29

As the physical parameter used to analyze 1H NMR spectra, σ is defined as the difference between the positions of the sample and standard signals [1]:

σ=(ννst)106[ppm].E30

The measurement range for hydrogen nuclei is about 15 ppm, which is about 3–10% of the spectroscopic range of other magnetic nuclei.

1.5.2. Presence of paramagnetic ions as a factor affecting a value of σ

During the interaction between metals, such as Eu3+ or Pr3+, and molecules containing oxygen or nitrogen atoms, metal ions increase their coordination number and form unstable associations [1]. This changes the chemical environment of the protons as well as the dipolar interactions between the unpaired electrons of the metal ions and the protons. This causes a change in the chemical shifts of hydrogen. The value of σ depends on the distance between the proton and the paramagnetic ion, and the value decreases as the distance increases [1].

Paramagnetic ions are used to distinguish the signals assigned to the choline groups of phospholipids. The most frequently used ions are those from the lanthanide group [13]. The concentration of paramagnetic ions added to the external environment of liposomes varies and is unique to each ion. If the concentration is too great, it could broaden all NMR spectra due to the dominant paramagnetic interactions [14].

1.5.3. Intensity of the signal

The intensity of the signals in the 1H NMR spectra is a measure of adsorption, and it is proportional to the number of protons that induced the signal [1]:

I=(1+γ2B12T1T2)12.E31

The above formula takes into consideration the relaxation times of different protons in the molecule, and the strength of B1 field. The physicochemical analysis of the NMR spectrum assumed that the area under the obtained signal is proportional to the number of protons inducing the signals [3]. Therefore, based on the measured area under the signal, the relative number of the protons in the molecule can be determined. Sometimes, the determined relative number of the protons in the molecule is not equal to the number of the protons resulting from the molecular formula. This usually occurs when oxygen, nitrogen, or sulfur are present in a molecular structure; in such a case, the proton may be exchanged for deuterium from the deuterated solvent [3].

1.5.4. Chemical shift anisotropy (CSA)

The 31P lineshape is directly related to the CSA tensor and to the orientation of lipid molecule (relative to the B0 field) [15]. The value of CSA depends on the phosphate group motion and the temperature. The 31P spectra exhibit a characteristic narrow peak (σ—high‐field maximum; isotropic part) and a low‐field shoulder (σ—anisotropic part). The CSA can be calculated using the fallowing formula [16]:

Δσ=3(σ+σ),E32

where σand σare the values of 31P shielding of the lipid molecules, oriented parallel or perpendicular relative to the magnetic field.

The value of CSA for multilamellar vesicles (MLVs) is about 40–50 ppm, and depends from the size of liposome. For small unilamellar vesicles (SUVs), the CSA value may decrease to 10 ppm [15].

1.6. Two‐dimensional (2D) NMR spectra

The two‐dimensional NMR experiments are most often used to determine the third and fourth structure of macromolecules.

1.6.1. Two‐dimensional correlation spectroscopy (COSY) and total correlation spectroscopy (TOCSY)

The 2D homonuclear (H,H)‐COSY experiment showed that the spectra with 1H chemical shifts along the axes were correlated with each other [10, 17]. The pulse sequence of the COSY experiment is shown in Figure 7A. In the COSY spectra, the diagonal and cross peaks are visible and they always differ by 90o.

Figure 7.

The pulse sequence of (A) the COSY experiment and (B) the TOCSY experiment; t1 is a time between pulses; t2 is an acquisition time; π and π/2 are pulses; spin‐lock is a time of the spin states mixing.

When the system of two spins corrects the phase, the absorption signal appears as the cross peak, while the dispersion signal (incorrect phases) appears as the diagonal peak. The cross peaks in the horizontal and vertical directions are the absorption signals with positive and negative amplitudes [10]. When a proton is coupled with more than one proton, the diagonal peak is seen in the corner and it occupies more than one square. This simple rule of spectrum analysis makes the COSY technique an ideal tool for evaluating the 1H spectra. The COSY spectra reveal information about the scalar coupling of protons within a few bonds [10, 17]. One disadvantage of the COSY method is that the cross peaks overlap with the diagonal signal when the chemical shift differences between the coupled nuclei are too small [10]. The 2D homonuclear (H,H)‐TOCSY experiment is similar to the COSY experiment, but it exhibits the peaks from all scalar‐coupled protons from the spin system [10, 17]. The pulse sequence of the TOCSY experiment is shown in Figure 7B. The method used to analyze the TOCSY spectra is similar to the method used to analyze the COSY spectra.

1.6.2. Rotating frame nuclear Overhauser enhancement spectroscopy (ROESY)

The nuclear Overhauser enhancement spectroscopy (NOESY) and ROESY techniques reveal information about the dipole‐dipole‐coupled protons and the transfer of polarization (the cross‐polarization effect) [10, 17]. Usually, the ROESY experiment is used when the studied molecules are large. Then, the dipole‐dipole relaxation is less effective because the rotation of the molecule is slow (long τc time), which extends the T1 time [10]. Thus, the ROESY method is more convenient, because it ensures that the NOE effect will be positive [18, 19]. The pulse sequence of the ROESY experiment is shown in Figure 8.

Figure 8.

The pulse sequence of the ROESY experiment; t1 is a time between pulses; t2 is an acquisition time; π and π/2 are pulses; Δ is internal fixed time.

Both processes, that is, dipole‐dipole interaction through space and cross‐polarization, influence the relaxation pathway. The mechanism of magnetization transfer is schematically shown in Figure 9.

Figure 9.

The evaluation of the magnetization vector MA (black arrow) belonging to A protons and magnetization vector MX (gray arrow) belonging to X protons. Spins A and X are close enough to interact through space. The diagrams from A to F show evaluation of the MA and MX vectors after two gradient pulses (π/2). The small contreating arrows superposed at the ends of the vectors represent the portion of the magnetization that is transferred from the nuclei of the other sort by cross-polarization during mixing time ∆.

Dipole‐dipole interactions strongly depend on the distance between the coupled protons. Thus, the cross peaks differ in size. The ability to detect the interactions through space and to obtain general information about the distance between interacting protons provides valuable data about the stereochemistry of a molecule [10].

2. NMR studies of liposomes

The NMR spectra of liposomes differ from typical solution‐NMR spectra. It is due to the specificity of lipid aggregates. The size of liposomes, hydration level, and the packing regime of lipid molecules in such structures have an influence on the NMR spectra.

2.1. Analysis of the 1H NMR spectra

MLVs are not suitable for measurements using 1H NMR method due to the large number of lipid bilayers and the large size of the liposomes. In this case, the signals in the 1H NMR spectra are drastically broadened, and it is impossible to analyze any of the physical parameters of the spectrum. Even the signals on the spectrum of SUV/LUV are slightly broadened, and the broadening is not great enough to preclude spectrum analysis. This is due to the fact that liposomes are a kind of lipid aggregate, which distinguishes them from any substance which is soluble in water. From the physical point of view, this kind of sample is neither a homogeneous liquid nor a solid. The diffusion of water molecules and various ions through the lipid bilayer is the basic phenomenon that can be examined by the 1H NMR spectra. The observation of this process is possible, thanks to the use of paramagnetic ions. It changes the chemical environment and the screening constant of each proton from the outer layer, which changes the σ values. Consequently, the choline group signal is split (δ) into two signals and it is assigned to the protons from the outer and inner layers of the liposome (Figure 10).

Figure 10.

The characteristic 1H spectra of egg‐lecithin SUV (A) before and (B) after addition of 5 mM Pr3+ions.

The type and concentration of used lanthanide ions is crucial, for example, if it is too great, the concentration of Eu3+ ions may exhibit as broadened 1H spectra or they may even destroy the membrane structure. This effect is associated with the properties of the Eu3+ ions, which interact to the same extent with the hydrophilic and hydrophobic part of the liposome membrane. Moreover, the signal corresponded to water is broadened, which means that the Eu3+ ions also interact with the water molecules from the hydration shell of the liposome [14]. Thus, Pr3+ ions are most often used, which can split the choline signals within a few ppms without the effect of broadening the signals. The preferred concentration of Pr3+ ions ranges from 4 to 7 mM [20].

The split signals showed the different intensities. Since the area under a signal is directly proportional to the number of protons that induce the signal, a more intense signal is assigned to the protons from the outer layer and the lower intensity signal is assigned to the protons from the inner layer. This phenomenon is related to the asymmetric distribution of the lipid molecules in each layer. The splitting of the choline group signal creates new possibilities for research. The ratio of the area under the signal corresponds to the outer layer and the area under the signal corresponds to the inner layer (Io/Ii); this provides information about the size of liposomes (Figure 11) [13, 20].

Figure 11.

The measure of the area under the signals (integral of the signal) assigned to the protons from the choline groups.

The splitting of the choline signal also offers the possibility to observe the ion transport through the membrane, since every change in the chemical environment of the protons from the outer layer and the inner layer will be visible as a change in the σ of the signals. The value of the δ increases as the difference between the chemical environment of the protons from the inner and the outer layer increases [20, 21]. It is very important to maintain a constant temperature and pH for the sample during this kind of experiment because the binding of metal cations is dependent on both of these parameters [22]. The effect is clearly seen in the 1H spectra of PE/PS/PC (phosphatidylethanolamine/phosphatidylserine/phosphatidylcholine) SUV (Figure 12).

Figure 12.

The time‐dependent changes of 1H resonance signals assigned to choline groups of PE/PS/PC SUVs after addition of 5.0 mM Pr3+ ions in the presence of Ca2+/PS molar ratio of 2.0.

The splitting of the choline groups is caused by 5 mM of the Pr3+ ions. Because the Pr3+ is associated only with the outer layer of a liposome, it is easy to observe the changes caused by the diffusion of Ca2+ions in the distance between the splitting signals and the intensity. The value of δ decreased as the concentration of the Ca2+ ions increased inside the liposome [21]. The intensity of the signals assigned to the choline groups changed because the fusion process occurred in the SUVs. Ca2+ ions are a well‐known fusogenic reagent. The fusion process of vesicles caused the increase in their size. Thus, the difference in the intensity of the signals assigned to the choline groups decreased. The fusion process also had an impact on the increase in the concentration of the Ca2+ ions in the liposome [21]. In a similar way, it is possible to conduct the experiment with compounds that are adsorbed on the surface of liposomes. Studying the area under each signal may also provide information about the processes that occur in the hydrophilic part of the lipid bilayer.

The half‐width of the signal is another parameter that can be analyzed in the 1H NMR spectra (Figure 13). The Δν1/2 of the signal is closely related to the dynamics of the chemical groups. The slower the movement, the greater the Δν1/2 of the signal.

Figure 13.

The Δν1/2 of the signals from 1H spectrum of egg‐lecithin SUVs.

The processes that occur on the surface of the liposomes can cause the choline head to either become rigid or become more fluid, thereby slowing down/speeding up the rotational motion, which results in an increase/decrease in the Δν1/2 of the signal. Additionally, cholesterol, antioxidants, and drugs contained in the liposomal membrane also may increase/decrease the fluidity of the membrane. For instance, the presence of azithromycin molecules increased the fluidity of 1,2‐dipalmitoyl‐sn‐glycero‐3‐phosphocholine (DPPC) liposomes decrease of Δν1/2 below the main phase transition temperature. However, amphotericin B rigidified the hydrophilic part (increase of Δν1/2) of phospholipid bilayer, which increased the fluidity (decrease of Δν1/2) of the hydrophobic core of PC membrane [16]. The opposite effect can be observed in the 1H spectra of PC liposomes in the presence of polysialic acid (Figure 14). The interaction of well‐hydrated and anionic polysialic acid with the PC headgroups increased the fluidity of the hydrophilic part of the membrane and rigidified the hydrophobic core of the PC bilayer [23, 24].

Figure 14.

The effect of polysialic acid on changes of the PC SUVs membrane fluidity (the Δν1/2 changes). (A) The 1H NMR spectra of PC SUVs and (B) in the presence of polysialic acid.

The effect manifested as a decrease in the Δν1/2 of the signal assigned to the choline groups from the outer layer of the liposome and as an increase in the Δν1/2 of the signal assigned to the choline groups from the inner layer and to the –(CH2)n and –CH3 groups from the fatty acid chains. In fact, the observed effect is connected to the restricted motion of lipid molecules in the bilayer structure. The increase of the fluidity of the headgroups is connected to the increase of their rotational motion, their reorientation, and better hydration [17, 24]. The motion of the headgroups is restricted by the hydrophobic interactions, hydrogen bonding, and the intermolecular force between the lipid molecules. Thus, the presence of polysialic acid indicates that the strength of the PO4 and N+(CH3)3 interactions is weakened. The increase of hydration (unrestricted motion) of headgroups and their reorientation causes the hydrocarbon chains to be more exposed to water molecules. Thus, the membrane polar‐apolar interface is more hydrophobic [25].

Studies of liposomes using 1H‐NMR can also be conducted with various physical parameters, including temperature. The influence of temperature on the 1H spectra manifests as an increase in the σ values of all the signals (Figure 15).

Figure 15.

The influence of temperature on the 1H resonance signals shifting. The 1H NMR spectra of PC SUVs at (A) room temperature and (B) 308 K.

The resonance signals shift toward the direction of the lower magnetic field. This effect is typical for lipid bilayers [26]. The increase in temperature has an impact on the increase in the fluidity of the membrane. It manifests as a decrease of the Δν1/2 of the signals. The 1H NMR temperature‐dependent studies may be used to analyze the properties of temperature‐sensitive liposomes. The concept of using temperature‐sensitive liposomes as drug carriers in local hyperthermia is based on the increase in their therapeutic effect, the ability to reduce drug toxicity for normal cells, and the increase in the permeability of the lipid bilayer at the proper temperature [26, 27]. The release rate of a drug depends on the temperature changes and the serum compounds (lipoproteins); thus, liposomes should be stable in serum and they should release drugs slowly under a proper temperature [28]. The 1H spectra of PC and PC/octadecylamine liposomes (positively charged LUV) showed a narrowing of the resonance signals (decrease in the Δν1/2) assigned to –N+(CH3)3, ‐(CH2)n, and –CH3(Figure 16). The largest changes in the Δν1/2 were observed for the signal corresponding to the fatty acid chains, –(CH2)ngroups [26]. The effect was observed in temperatures ranging from 5 to 50°C.

Figure 16.

The 1H NMR spectra of PC/octadecylamine liposomes at (A) room temperature and (B) 318 K.

Studies on changes of the splitting and intensity of signals assigned to choline groups revealed that the size of the liposomes increases [26]. In fact, the size of PC liposomes changed from 20–30 nm to 1 μm. Thus, as the temperature increases, the size of the PC liposome also increases. Additionally, when the temperature ranges from 30 to 40°C, the structure of the liposome membrane is damaged; this makes the PC liposomes completely unstable and useless as thermosensitive drug carriers [26]. Quite the opposite effect was observed in the case of the PC/octadecylamine liposomes. The temperature had a smaller influence on the liposome size changes (from 20 to 60 nm), and the size was changed slowly. The PC/octadecylamine liposomes also seemed to be stable at temperatures ranging from 40 to 50°C [26]. At higher temperatures, the temperature‐sensitive liposomes may aggregate or fuse, which makes it possible to transfer the drug to cells by fusion or via an endocytosis process in hyperthermia [26, 29].

2.2. Analysis of the 31P NMR spectra

Liposomes can be composed of one or more kinds of phospholipids. Their molecular structure contains phosphorus atoms, which makes the 31P NMR method extremely useful for studying them. Both MLVs and LUVs/SUVs can be studied using that method. In MLVs, the signal assigned to the phosphate groups is drastically broadened, and its shape is not like the Lorentzian function, unlike the SUVs/LUVs. The 31P NMR spectra are mainly used to study the thermotropic properties of liposomes. The various phases change the structure of membrane, and there are specific transition temperatures for each lipid.

The 31P spectra of PC/Ch MLVs exhibited changes in the lineshape in the temperature range of 10–40°C. The characteristic 31P lineshape of the phase was observed before the phase transition. The monitoring of 31P spectra at short temperature intervals led to observing the intermediate lineshapes between those characteristics for the and HII phases [21]. Substances that are either added to liposome membrane or are found in the liposomal environment can increase/decrease the phase transition temperature. For instance, this can be seen in the effect of the change of phase transition temperature caused by various drugs. The 31P spectra of DPPC MLVs used in various concentrations of Piracetam showed that an additional narrow signal was assigned to the drug [30]. As the concentration of Piracetam increases, the intensity of the signal assigned to that drug also increases. The temperature studies showed that the main phase transition temperature of DPPC MLVs decreased in the presence of Piracetam. The results suggest that hydrophilic Piracetam molecules are associated with the hydrophilic part of the liposome membrane, which increases the fluidity of the membrane. Thus, the temperature of the main phase transition decreases [30].

The 31P NMR method reveals information about the mobility of phosphate groups and about local order. Various substances impact the dynamics of the hydrophilic and hydrophobic part of the lipid bilayer. The number of narrow peaks in the 31P spectra depends on the number of phospholipid types used to form the liposomes [31]. The effect is due to differences in the chemical surroundings of the phosphate groups in each type of phospholipid molecule. In the 31P spectra, the splitting between the signals depends more on the averaged CSA motions than on an isotropic chemical shift [31]. It influences the value of chemical shift, and it may be difficult to assign the signals.

As previously mentioned, the effect of azithromycin on DPPC MLVs [16] was studied. Temperatures ranging between 35 and 45°C did not change the position of the narrow peak and CSA value, but the 31P lineshape in the low‐field shoulder was changed, and the presence of azithromycin in the external environment of the liposome decreased the CSA value. In fact, above 40°C, only a narrow peak stays in the spectrum because the CSA value is averaged to zero. The azithromycin caused an increase in the fluidity of the DPPC membrane below the temperature of the main phase transition [16].

When hydrophobic β‐carotene is added to a lipid membrane, changes in DPPC membrane fluidity can also be observed. In temperatures above the main phase transition, β‐carotene increases the fluidity of the DPPC MLV membrane, and in temperatures below the main phase transition, it decreases the fluidity of that membrane [32].

The opposite effect may be observed for PC MLVs in the presence of polysialic acid (Figure 17). In the temperature range of 10–30°C, the 31P spectra show a narrowing of the isotropic part and broadening of the anisotropic part (increase in the CSA value) [23]. The increase of well‐hydrated polysialic acid in the membrane increased the fluidity of the headgroups; this resulted in a decrease in the hydrophobic core fluidity.

Figure 17.

The effect of polysialic acid on 31P NMR spectra of PC MLVs at room temperature. The 31P spectra (A) before and (B) after addition of polysialic acid.

The 31P NMR spectra also are used to study the fusion process that occurs between liposomes. These types of experiments use fusogenic factors. The 31P NMR is the best method for examining the fusion process because it has been proven that the fusion of two vesicles is accompanied by a transient structure, that is, the inverse hexagonal phase [29, 33].

The 31P lineshape of PE/PS/PC MLVs showed changes with an increased molar ratio of Ca2+ ions (Figure 18). In fact, the Ca2+ ions are a well‐known fusogenic factor. The characteristic lineshape for the HII phase was obtained when the molar ratio of Ca2+/PS was 2.0 [21]. It means that when the molar ratio of Ca2+/PS is 2.0, the fusion process occurs.

Figure 18.

The 31P NMR spectra of PE/PS/PC liposomes; (A) the characteristic 31P lineshape of MLVs in the Lα and HII phase, (B) the time‐dependent changes of LUVs after addition of 5.0 mM Pr3+ ions in the presence of Ca2+/PS molar ratio of 2.0.

The monitoring of changes in the 31P spectra of PE/PS/PC SUVs after the addition of Ca2+ and Pr3+ ions showed that the signals were assigned to the choline groups from the inner and outer layers of the membrane. The obtained results revealed the decrease of δ‐value and the intensity equalization of the signals corresponded to the choline groups [21]. The Pr3+ ions are only associated with the outer layer of a liposome, while the Ca2+ ions are not (they can diffuse through the membrane). Thus, the obtained results suggest that during the transient phase (inversed micelle) the fusion process, the lipid molecule transition from the outer to the inner layer and the size of the liposome, increases [21]. These results also confirmed the topological model of the fusion.

2.3. Two‐dimensional NMR spectroscopy

2D NMR spectroscopy is the most convenient technique for studying scalar and dipole‐dipole couplings. The results obtained using the method are used to study the structure and stereochemistry of molecules [10].

2.3.1. COSY and TOCSY spectra

In the COSY spectra of the PC/octadecylamine SUVs, the signals from the protons coupled within a few of the chemical bonds are visible. The diagonal peaks corresponded to each proton cross‐correlated with every other proton from spin system [17, 26]. Figure 19 depicts the method used to analyze the COSY spectra.

Figure 19.

The COSY spectra of PC/octadecylamine SUVs. The cross peaks of scalar-coupled protons from the PC and octadecylamine molecules, and between protons from the PC/octadecylamine and water molecules (unsigned) are depicted in gray.

On comparing the TOCSY and COSY spectra of PC/octadecylamine SUVs, it is possible to determine the differences between them (Figure 20). In the TOCSY spectra, all the scalar‐coupled protons in the PC molecule can be seen [17]. Both spectra also exhibit the cross peaks between the fatty acid chain groups from the hydrophobic core of the membrane and the water molecules. This is characteristic of well‐hydrated membranes in the liquid crystalline phases [17, 26].

Figure 20.

The TOCSY spectra of PC/octadecylamine SUVs. Only new cross peaks are depicted in gray.

2.3.2. ROESY spectra

The ROESY spectra of the PC/octadecylamine SUVs reveal information about the protons coupled through space. Now, the diagonal peaks exhibit dipole‐dipole interactions. The observed interactions may occur within one molecule or between neighboring PC molecules. The results also show the interactions between the hydrophilic part (headgroups) and the hydrophobic part (fatty acid chains) of the membrane (Figure 21) [17, 26]. The size of the cross‐peak is directly proportional to the distance between the coupled protons. This dependency is clearly visible on the ROESY spectra of the PC/octadecylamine liposomes.

Figure 21.

The ROESY spectra of PC/octadecylamine SUVs. (A) The full assigned 1H spectra of the PC/octadecylamine liposomes; (B) The ROESY spectra of the PC/octadecylamine SUVs with assigned dipole-dipole interactions (in gray).

The 2D NMR technique is most often used to study the interactions between lipid molecules and substances added to the liposome membrane or to the external environment of liposome. To examine the nature of the interaction (either scalar coupling or through space), COSY/TOCSY or NOESY/ROESY experiments should be used, respectively. An example of this could be the ROESY spectra of PC/octadecylamine SUVs in the presence of sialic acid [17]. When comparing the ROESY spectra of sialic acid and PC/octadecylamine liposomes with PC/octadecylamine SUVs in the presence of sialic acid, it is easy to identify the new cross peaks due to interactions of the protons from the PC/octadecylamine and sialic acid molecules. While the interaction of sialic acid with the PC membrane has a considerable impact on membrane fluidity, the ROESY results only showed one cross peak between sialic acid and the PC molecules [17]. It is interesting to note that the dipole‐dipole interaction occurs between the protons from the acyl group of sialic acid and the (‐CH2)n groups of the PC fatty chains. This result suggests that the other functional groups of sialic acid are well hydrated and, perhaps, the hydrophilic part of the membrane and sialic acid molecule interacts between their hydration shells. The obtained result also explains the strong influence that sialic acid has on decreasing the fluidity of the hydrophobic core of the membrane [17, 26].

© 2017 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Anna Timoszyk (October 25th 2017). Application of Nuclear Magnetic Resonance Spectroscopy (NMR) to Study the Properties of Liposomes, Liposomes, Angel Catala, IntechOpen, DOI: 10.5772/intechopen.68522. Available from:

chapter statistics

906total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Liposomes Used as a Vaccine Adjuvant-Delivery System

By Ning Wang, Tingni Wu and Ting Wang

Related Book

First chapter

Lipid Peroxidation: Chemical Mechanism, Biological Implications and Analytical Determination

By Marisa Repetto, Jimena Semprine and Alberto Boveris

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us