# A Small-Angle Neutron Scattering, Calorimetry and Densitometry Study to Detect Phase Boundaries and Nanoscale Domain Structure in a Binary Lipid Mixture

Natalie Krzyzanowski, Lionel Porcar, Ursula Perez-Salas

► To cite this version:

Natalie Krzyzanowski, Lionel Porcar, Ursula Perez-Salas. A Small-Angle Neutron Scattering, Calorimetry and Densitometry Study to Detect Phase Boundaries and Nanoscale Domain Structure in a Binary Lipid Mixture. Membranes, 2023, 13 (3), pp.323. 10.3390/membranes13030323. hal-05313439

HAL Id: hal-05313439

https://hal.science/hal-05313439v1

Submitted on 22 Oct 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# A Small-Angle Neutron Scattering, Calorimetry and Densitometry Study to Detect Phase Boundaries and Nanoscale Domain Structure in a Binary Lipid Mixture

Natalie Krzyzanowski <fnref n="1" />,<fnref n="2" />, Lionel Porcar<fnref n="2" /> and Ursula Perez-Salas <fnref n="1" />,* ³

**Abstract:** Techniques that can probe nanometer length scales, such as small-angle neutron scattering (SANS), have become increasingly popular to detect phase separation in membranes. But to extract the phase composition and domain structure from the SANS traces, complementary information is needed. Here, we present a SANS, calorimetry and densitometry study of a mixture of two saturated lipids that exhibits solidus-liquidus phase coexistence: 1,2-dipalmitoyl-d62-*sn*-glycero-3-phosphocholine (dDPPC, tail-deuterated DPPC) and 1,2-dilauroyl-*sn*-glycero-3-phosphocholine (DLPC). With calorimetry, we investigated the phase diagram for this system and found that the boundary traces for both multilamellar vesicles (MLVs) as well as 50 nm unilamellar vesicles overlap. Because the solidus boundary was mostly inaccessible by calorimetry, we investigated it by both SANS and molecular volume measurements for a 1:1 dDPPC:DLPC lipid mixture. From the temperature behavior of the molecular volume for the 1:1 dDPPC:DLPC mixture, as well as the individual molecular volume of each lipid species, we inferred that the liquidus phase consists of only fluid-state lipids while the solidus phase consists of lipids that are in gel-like states. Using this solidus-liquidus phase model, the SANS data were analyzed with an unrestricted shape model analysis software: MONSA. The resulting fits show irregular domains with dendrite-like features as those previously observed on giant unilamellar vesicles (GUVs). The surface pair correlation function describes a characteristic domain size for the minority phase that decreases with temperature, a behavior found to be consistent with a concomitant decrease in membrane mismatch between the liquidus and solidus phases.

**Keywords:** lipid phase separation; binary lipid phase diagram; dipalmitoylphosphocholine; DPPC; dilauroylphosphocholine; DLPC; unilamellar vesicles; model membrane; small-angle neutron scattering; contrast matching; calorimetry; densitometry; MONSA; bead vesicle model

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Krzyzanowski, N.; Porcar, L.; Perez-Salas, U. A Small-Angle Neutron Scattering, Calorimetry and Densitometry Study to Detect Phase Boundaries and Nanoscale Domain Structure in a Binary Lipid Mixture.

*Membranes* **2023**, *13*, 323.
https://doi.org/10.3390/membranes13030323

Academic Editor: Thais A. Enoki

Received: 15 February 2023

Revised: 6 March 2023

Accepted: 7 March 2023

Published: 10 March 2023

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

## 1. Introduction

The major constituents of biological membranes are lipids. Because lipids are amphiphilic, they spontaneously self-assemble into the essential feature of biological membranes: a lipid bilayer. Lipid bilayers thus serve as models for biological membranes and have been the subject of intense research for over half a century. Early on, it was recognized that the lipid composition in biological membranes is immensely diverse—a critical feature for function [1]. This led to the study of relatively simple lipid mixtures, such as binary lipid mixtures. Although simple, even binary mixtures produce somewhat sophisticated compositional phase diagrams. The intent of understanding the behavior of simple lipid mixtures was that they could suggest possible features of biological membranes as well [2]. Indeed, a thermodynamic framework led to a deeper understanding of these mixtures, from molecular cooperativity to non-ideal mixing characteristics [3–7], which became the precursor of what later became the lipid raft hypothesis [8–10].

<footnote n="1">¹ Department of Physics, University of Illinois at Chicago, Chicago, IL 60608, USA</footnote>

<footnote n="2">² Large Scale Structures Group, Institut Laue-Langevin, CEDEX 9, 38042 Grenoble, France</footnote>

<footnote>* Correspondence: ursulaps@uic.edu</footnote>
The lipid raft hypothesis proposes that rather than having biological membrane constituents be randomly mixed, the cell membrane partitions lipids, sterols and proteins into laterally segregated regions. Although some cell systems have been found to exhibit observable micron-scale domains [11,12], the majority of cell membranes do not, which led to the idea that lipid rafts were still present but were mostly small (nanometers), dynamic and possibly transient [13]. Thus, techniques that study membranes at the nanoscale have been heavily pursued. Small-angle neutron scattering (SANS) stands out as a technique to probe this length scale because even though it is considered to be a low-resolution structure technique, with the use of selective deuteration, it has exquisite sensitivity to the emergence of nanoscale heterogeneities in membranes [14] and, as a result, has helped advance the field both using membrane mimics [15] as well as living cells [16]. In addition, analysis of the structures formed by phase separation in model membranes has given insights into, for example, the role of membrane mismatch between domains and the surrounding membrane environment and the size of the domains [17].

Analysis of SANS data from membranes exhibiting domains has mostly been conducted using analytical models with predefined shapes, such as stripes and circular domains [18–20]. More recently, Bolmatov et al. developed a promising scattering model that uses a phenomenological free energy approach and whose fitting parameters could in principle identify domain configurations having diverse shapes; however, its implementation is still ongoing [21]. Indeed, as shown by Heberle et al. [15], the structures of phase-separated membranes can be quite complex. One example is the binary mixture of two saturated phosphatidylcholine lipids having different tail lengths: 16- and 12-length carbons (di-pamitoyl phosphatidylcholine DPPC and di-lauroyl phosphatidylcholine DLPC, respectively), which does not display round domains or simple stripes [22,23]. Indeed, limiting the analysis of SANS data from membranes exhibiting complex structures to these analytical shapes can produce inconsistencies with other aspects of these membranes. For example, Anghel et al. [18], who studied a deuterated analog for this system (dDPPC:DLPC, where dDPPC differs from DPPC by having deuterated tails) by SANS, obtained fits that suggested compositions that conflicted with the known phase diagram for DPPC:DLPC reported by Van Dijck et al., obtained by calorimetry [3].

In this work, we explore the analysis of SANS spectra from 1:1 dDPPC:DLPC vesicles using an alternative strategy: bead vesicles. The roughly 1200 beads forming the vesicle provide for local variations in composition and collectively are not constrained to any particular shape. Using this approach, the resulting membrane structures or patterns found were similar to those observed in large unilamellar vesicles by microscopy [22,23]. The resulting irregular domains in the bead vesicles could be characterized by a surface bead pair correlation function, from which a temperature-dependent characteristic length was obtained.

The resulting phase diagram for the binary mixture of dDPPC and DLPC, which combined calorimetry, densitometry and SANS analysis, produced the elusive solidus boundary for the system. Together, the solidus and liquidus boundaries describe a non-ideal system that follows the thermodynamics of regular solution theory [3]. Interestingly, the combined approach of calorimetry, densitometry and SANS revealed details of the behavior of these lipids in each phase—particularly the volume of DLPC in the solidus phase—which had not been reported before.

## 2. Materials and Methods

### 2.1. Materials

1,2-dipalmitoyl-*sn*-glycero-3-phosphocholine (hDPPC), 1,2-dipalmitoyl-d62-*sn*-glycero-3-phosphocholine (dDPPC) and 1,2-dilauroyl-*sn*-glycero-3-phosphocholine (DLPC) were purchased from Avanti Polar Lipids (Alabaster, AL, USA) as lyophilized powders and stored at -80 °C until further use.
## 2.2. Differential Scanning Calorimetry (DSC)

Using a micro differential scanning calorimeter, µDSC III from SETARAM Instrumentation, we obtained the phase diagram for dDPPC:DLPC. DSC data were taken for both 50 nm unilamellar vesicles (SUVs) and multilamellar vesicles (MLVs). dDPPC and DLPC in powder form were weighed into a glass vial and thoroughly mixed using chloroform (vortexed vigorously until the solution was clear). The removal of this chloroform was performed first with a stream of N₂ gas using a nitrogen evaporation system (N-EVAP) until a thin in-appearance dry lipid film remained on the vial’s walls. The vials were then placed in vacuum at 60 °C overnight to remove any residual chloroform. Once dry, the sample was hydrated with 1 ml of H₂O. MLVs were prepared by leaving the solution on a shaker overnight mixing in a 47 °C oven until the resulting mixture became a homogenous suspension. The SUVs were prepared, subsequently, by extrusion. The extrusion was performed with a modified Avanti Polar Lipids extruder system that included heating via an Anova water bath (see Supplementary Information, Section S1). A syringe pump mechanically manipulated the syringes of the extruder (New Era Pump Systems Inc., Farmingdale, NY, USA). The final extrusion consisted of 41 passes through a polycarbonate filter with a 50 nm average pore size. The extrusion of the vesicles was performed at 50 °C, which is well above the melting temperature of the lipid system. The speed of extrusion can be controlled on the syringe pump, which allows extruding of even relatively high-concentration solutions.

To obtain the phase diagram of dDPPC:DLPC, thermograms for several dDPPC:DLPC ratios were collected. Keeping DLPC constant (15 mg/mL), the amount of dDPPC was varied from 0 mol% to 90 mol%. The 100 mol% dDPPC was run at 15 mg/mL. SUVs were obtained from the MLV solutions after their calorimetry data were collected. Thermogram scans were obtained at 0.2 °C/min. SetSoft 2000 provided by SETARAM was used to collect DSC data and exported as a text file to use with the program Mathematica™ for baseline subtraction and determination of onset and completion temperatures (details in Section S3 of the Supplementary Information).

## 2.3. Densitometry

Densitometry data were obtained on samples having a total 10 mg/mL lipid concentration using the Density and Sound Velocity Meter: DMA 5000 M from Anton-Paar (Ashland, Virginia) at the Partnership for Soft Condensed Matter (PSCM) in Grenoble, France. Three samples were measured by cooling: dDPPC, DLPC and 1:1 dDPPC:DLPC. The lipid samples consisted of multilamellar vesicles (MLVs) in H₂O, which had been incubated and gently rocked overnight at 40 °C. The samples were loaded at room temperature, and densitometry data were taken from 50 °C down to 5 °C in steps of 1 °C.

## 2.4. Small-Angle Neutron Scattering

SANS measurements were taken on 50 nm SUV vesicle solutions using a 1:1 ratio of dDPPC:DLPC. The measurements were performed on the D22 instrument at Institut Laue-Langevin (Grenoble, France). Instrument configuration covered a Q range of 0.003 ≤ Q ≤ 0.5 Å⁻¹. The wavelength used was 6 Å with a wavelength resolution of Δλ/λ = 10%. Data were collected with a 2-D detector and reduced using GRASP, the reduction package provided by the ILL or the reduction package provided by NCNR in IGOR [24].

# 3. Results

SANS is a powerful technique to obtain the structural information of particles that range from a few to tens of nanometers in size [25]. This is because the scattered intensity, I(Q), is directly related to the particle’s shape, size and composition. For a dilute, monodisperse particle solution, I(Q) is given by:

$$I(Q) - I_{\text{incoh}} = vV_{\text{particle}} (\text{SLD}_{\text{particle}} - \text{SLD}_{\text{solvent}})^2 P(Q)_{\text{particle}} \quad (1)$$
where $P(Q)$ is the form factor of the particles and contains information about their shape and size. $v$ is the volume fraction of particles, and $V_{particle}$ is the volume of the particle. The scattering length density (SLD) is a measure of the interaction of neutrons with the atomic and isotopic species in the particles' and solvent's chemical composition. As shown in Equation (1), to be able to extract information about the particles' shape and size through $P(Q)$, contrast is key, i.e., $(SLD_{particle} - SLD_{solvent})^2 > 0$.

Vesicles whose membranes are composed of a 1:1 molar mixture of DPPC and DLPC will phase separate and form lipid domains—solidus and liquidus phases—with distinct lipid compositions given by the phase diagram [3]. The approach to achieve contrast between these domains and be measured by SANS is to substitute one of the lipids with its deuterated version. This approach does not change the chemical identity of the lipid, but it does change its SLD significantly. For example, at room temperature, the tail SLD of DPPC (in its fully hydrogenated form) is $-0.4 \times 10^{-6} \AA^{-2}$, while for dDPPC (with all 62 tail hydrogens substituted with deuteriums), it is $7.5 \times 10^{-6} \AA^{-2}$. Hence, when a mixture consisting of a deuterated lipid species (dDPPC) and a hydrogenated lipid species (DLPC) phase separates, the phase-separated regions acquire different SLDs because of the difference in lipid composition in each phase. To drastically highlight these differences in SLD between phases in the membrane, the solvent's SLD is chosen to match the membrane's mean SLD, typically obtained when the membrane's lipids are fully mixed—in a single phase—above the miscibility transition [14]. At this contrast match condition (when $SLD_{vesicle} = SLD_{solvent}$), the scattering becomes flat (equal to the background) and featureless (as shown by Equation (1)). However, when the temperature is lowered below the miscibility temperature, phase separation emerges and, because the SLDs of the domains do not match each other or that of the solvent, scattering emerges too. Figure 1 shows the scattering from 1:1 dDPPC:DLPC 50 nm in diameter vesicles in their contrast-matched (CM) solvent (see Section S2 in the Supplemental Information for details on obtaining the CM condition for the system) above and below the miscibility temperature. At 40.6 °C, which is above the miscibility transition for this mixture, the scattering signal is flat, confirming that the membrane's SLD has been matched to the solvent's SLD. Upon lowering the temperature, below the miscibility transition, phase separation into solidus and liquidus domains ensues, and their presence is clearly captured by a significant increase in the scattering intensity [14]. This emergent scattering, however, is not due to the vesicle form factor as it is contrast-matched to the solvent but due to the characteristics of the domains in the vesicles. Extracting structure and composition information from these coexisting phases in vesicles, beyond just detecting their presence, requires that the collected scattering signal captures the peak resulting from the size of the domains present in the vesicles. Using 50 nm vesicles allows the capture of this peak [17,21], as it is clearly apparent at 24 °C and 22 °C. As temperature decreases, we observe the scattering signal evolve. Initially, the scattering peak increases in intensity and shifts to higher Q values (from 24 °C to 22 °C), and then low Q scattering increases while the scattering peak lowers in intensity while still shifting to higher Q values. At the lowest temperatures, the scattering curves have significant low Q scattering while the peak has become broad and shoulder-like in the high Q. This evolution in the scattering spectra suggests that the membranes' domains are not only changing in lipid composition, size and number but that the lipids' SLD is changing as well. Therefore, in order to quantitatively analyze the SANS data presented in Figure 1, we need to constrain the analysis with the thermodynamic parameters of the system. The phase diagram provides the lipid composition of each phase and fraction of each phase [26] while the volumes of the lipids provide the SLD of each phase [27].
<figure><img src="image_5.png" /><figcaption>**Figure 1.** SANS spectra for 20 mg/mL, 50 nm in diameter small unilamellar vesicles of a 1:1 mixture of tail-deuterated DPPC (dDPPC) and DLPC at different temperatures. The solvent matches the mean SLD of the vesicles at high temperature (40.6 °C) when no domains are present. At lower temperatures, excess scattering appears due to the emergence of solidus and liquidus domains, which unequally partition deuterated and non-deuterated lipids, creating contrast between regions within the membrane as well as the solvent.</figcaption></figure>

The phase diagram for DPPC:DLPC has been previously reported [3,28]; however, because of our use of the tail-deuterated DPPC (dDPPC), we revise it. Calorimetry is a technique that is commonly used to generate the phase diagram from heat capacity traces. Figure 2A shows the specific heat capacity (cooling) traces for dDPPC:DLPC compositions varying from 100 mol% dDPPC to 100 mol% DLPC. The phase diagram's liquidus and solidus boundaries are determined by the onset and completion temperatures of phase transitions (see Supplementary Information Section S3 for details; Figure S3.1C,F). The onset and completion temperatures of the phase transitions detected in the heat capacity trace of a given mixture are marked with dots in Figure 2A and displayed as a composition versus temperature phase diagram in Figure 2B. Since the SANS experiments were performed on 50 nm in diameter small unilamellar vesicles (SUVs), Figure 2B also shows the position of the onset and completion temperatures for both multilamellar vesicles (MLVs) and 50 nm SUVs (see Supplementary Information S3 showing the heat capacity traces for MLVs during heating and for SUVs during cooling). As observed in Figure 2B, the boundaries of the phase diagram remain the same for both MLVs and SUVs, demonstrating that the high membrane curvature in SUVs does not affect the phase diagram boundaries of the system.

The phase diagram for dDPPC:DLPC, shown in Figure 2B, indicates that the system exhibiting a single phase is at a high temperature, but when the temperature drops and crosses the upper miscibility boundary, the system separates into a solidus and a liquidus phase. The molecular composition of each phase as well as the overall fraction of each phase at any given temperature (through the lever rule—see more details in Section S4 of the Supplementary Information) is extracted from the upper miscibility or liquidus boundary and the lower miscibility or solidus boundary of the phase diagram. The calorimetry phase diagram shown in Figure 2B describes a well-defined liquidus boundary for all dDPPC compositions; however, the solidus boundary appears to have a gap or discontinuity
between $\chi_{dDPPC} = 0.6$ and $\chi_{dDPPC} = 0.7$. Indeed, a similar behavior for the solidus boundary was reported for the DPPC:DLPC system by Van Dijck et al. [3]. As we will show, by combining this calorimetry information with the required thermodynamic volumes of the lipids in the solidus and liquidus phases in the analysis of the SANS data, we will be able to retrieve the elusive solidus boundary for the system.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** (A) Cooling calorimetry traces for MLVs for mixtures of dDPPC:DLPC. Large dots mark the onset and the completion temperatures of the transition. The dots trace the liquidus and the solidus boundaries. (B) The composition, $\chi_{dDPPC}$, vs. temperature, T, plot corresponding to the phase diagram for the dDPPC:DLPC system. The solid dots correspond to MLVs, blue for cooling and red for heating. The open blue circles correspond to 50 nm SUVs during cooling. Error bars correspond to at least two calorimetry scans (see Figure S3.1 and Table S3.1 in Section S3 of the Supplementary Information for details).</figcaption></figure>

In order to obtain the thermodynamic molecular volumes of the lipids, we used densitometry. Our densitometry measurements consisted of measuring the density of a mixture of lipids in water, which is given by:

$$ \rho_{\text{mixture}}(T) = \frac{w_{\text{lipids}} + w_{\text{water}}}{V_{\text{lipids}} + V_{\text{water}}} \qquad (2a) $$

where $w_{lipids}$ and $w_{water}$ are the weights and $V_{lipids}$ and $V_{water}$ are the volumes of lipids and water, respectively. From Equation (2a), we can solve for the molecular volume of the lipids, $V_L$:

$$ V_L(T) = \frac{M_w \times 10^{24}}{N_A} \frac{1}{w_{\text{lipids}}} \left( \frac{w_{\text{lipids}} + w_{\text{water}}}{\rho_{\text{mixture}}(T)} - \frac{w_{\text{water}}}{\rho_{\text{water}}(T)} \right) \qquad (2b) $$

Here, $N_A$ is Avogadro's number and $M_w$ is the molecular weight of the lipid system. The factor $10^{24}$ is a unit conversion factor from cm$^3$ to Å$^3$, since the density was obtained in g·cm$^{-3}$. Because we performed the density measurements in H$_2$O, we also measured $\rho_{H2O}$ separately (see Figure S5.1 in Section S5 of the Supplementary Information). Figure 3 shows a plot of the molecular volume, $V_L$, for DLPC, dDPPC and the 1:1 dDPPC:DLPC mixture as a function of temperature in the range between 5 °C and 50 °C. We observe the expected sharp molecular volume change in dDPPC between the fluid and the gel state at the known $T_m$ (melting or transition temperature) value of 37 °C [29]. Below and above 37 °C, the molecular volume change in dDPPC is linear and slowly varies with temperature. DLPC's molecular volume, which is in the fluid phase between 5 °C and 50 °C, is also observed to be linear and slowly varying with temperature. Interestingly, away from the $T_m$, the change in the molecular volumes for both dDPPC and DLPC is very similar, both displaying similar slopes (see Table S5.2 in Section S5 in the Supplementary Information). In addition, as shown in Figure S5.2 in Section S5 of the Supplementary
Information, we find that these molecular volumes are consistent with those previously reported for dDPPC [30–33] and DLPC [33,34].

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Molecular volumes, $V_L$, for dDPPC, DLPC and the 1:1 mixture of dDPPC and DLPC obtained by densitometry (red circles) and by inversion of calorimetric traces (blue) according to Equation (3). The red arrows indicate corresponding $V_L$ extrapolations: see Table S5.2 in Section S5 of the Supplementary Information for the resulting linear, temperature-dependent molecular volumes for dDPPC and DLPC above and below $T_m$. The black trace shows the 1:1 molar additive mean volume obtained from the densitometry traces for dDPPC and DLPC. Error bars for $V_L$, $\delta V_L$, were calculated using the error in weighing (±0.0001 g) and the error in the density (±0.0001 g/mL) and are shown in Figure S5.3 in Section S5 of the Supplementary Information; for dDPPC and DLPC, the error bars are approximately the size of the symbol ($\delta V_L \approx 3–4$ Å), while for 1:1 dDPPC and DLPC, $\delta V_L \approx 7$ Å.</figcaption></figure>

For the 1:1 dDPPC:DLPC mixture, we also observe that the molecular volume behaves linearly above the miscibility temperature ($T > 27$ °C) with a slope consistent with those of dDPPC and DLPC in the fluid phase. However, for temperatures below its miscibility temperature ($\approx 27$ °C), the volume decreases more rapidly and nonlinearly. In this region, the system exhibits the coexistence of solidus and liquidus phases. As a reference, the molecular volume average of dDPPC and DLPC in the 5 °C to 50 °C range is shown with a black trace in Figure 3. We find that above 37 °C, this molecular volume average trace for dDPPC and DLPC exactly coincides with the molecular volume behavior of the 1:1 mixture. This linear behavior, however, extends to $\approx 27$ °C and thus differs from the molecular volume average of dDPPC and DLPC, which displays a sharp drop at 37 °C. Starting at $\approx 27$ °C, the molecular volume for the 1:1 mixture decreases smoothly but nonlinearly as the temperature drops. Notably, below 18 °C, the molecular volume average of dDPPC and DLPC becomes larger than the measured molecular volume for the mixture. This behavior revealed that a fraction of DLPC must be in the gel state. Hence, we inferred that the solidus phase must only comprise lipids in the gel state while the liquidus phase must only consist of lipids in the fluid state. Indeed, from the molecular volume behavior of the 1:1 mixture, we find that dDPPC remains fluid well below its $T_m$, down to $\approx 27$ °C, and will remain fluid while in the liquidus phase, while DLPC, which is fluid above freezing temperatures, will adopt a gel-state volume in the solidus phase up to $\approx 27$ °C. However, what are the fluid molecular volume of dDPPC below 37 °C and the gel-state volume of DLPC above –5 °C? We assume that the fluid molecular volume of dDPPC below 37 °C is a linear extrapolation of its fluid molecular, and similarly, the gel-state volume of DLPC above –5 °C must be a linear extrapolation of the gel-state molecular volume of DLPC [27]. Support for this assumption comes from the observation that the extrapolated fluid molecular volume of
dDPPC averaged with the fluid molecular volume of DLPC reproduced the molecular volume behavior of the 1:1 mixture above 27 °C.

In order to obtain the gel-phase molecular volume for DLPC, which could not be obtained by the densitometry approach because the lowest measurable temperature was 5 °C, we instead used the calorimetry traces for which lower temperatures were recorded. Indeed, as demonstrated by Ebel et al. [35], molecular volumes can be derived from heat capacity traces as follows:

$$V_L(T) = - \frac{M_w \times 10^{24}}{N_A} \int_T^{T_{max}} (\gamma c_p(T) + \alpha) dT + V_L(T_{max}) \quad T < T_{max} \quad (3)$$

where $\alpha$ and $\gamma$ are constants; $\alpha$ determines the linear behavior of the molecular volume away from the phase transition, while $\gamma$ is a pre-factor controlling the magnitude of the volume change at the transition.

Since the molecular volumes obtained from calorimetry had to coincide with those obtained by densitometry, the values for $\gamma$ and $\alpha$ were varied until the inverted calorimetry curve overlapped with that obtained from densitometry. In all cases, $T_{max}$ was chosen to be well above the miscibility temperature. Figure 3 shows inverted calorimetry curves for dDPPC, DLPC and 1:1 dDPPC:DLPC, using Equation (3). The values for $\gamma$ and $\alpha$ in each case are in Table S5.1 in Section S5 of the Supplementary Information. As shown in the table, for dDPPC, $\gamma$ was found to be $8.5 \pm 0.5 \times 10^{-4} \text{ cm}^3 \cdot \text{J}^{-1}$, which is within the reported value by Ebel et al. [35] ($8.5 \pm 0.8 \times 10^{-4} \text{ cm}^3 \cdot \text{J}^{-1}$). Ebel et al. also reported that this value of $\gamma$ was independent of the lipid type; however, for the 1:1 dDPPC:DLPC mixture and for DLPC, we found that $\gamma$ was slightly larger ($13.7 \pm 0.6 \times 10^{-4} \text{ cm}^3 \cdot \text{J}^{-1}$ and $15.2 \pm 0.8 \times 10^{-4} \text{ cm}^3 \cdot \text{J}^{-1}$, respectively).

Because the molecular volumes obtained for DLPC from densitometry did not reach the region where DLPC transitions to the gel phase, the value of $\gamma$ for DLPC had to be obtained following a self-consistency condition: at -5 °C, where DLPC and dDPPC are in the gel phase, the molecular volume for the 1:1 dDPPC:DLPC mixture had to be equal to the mean molecular volume of dDPPC and DLPC. Since molecular volumes from the inverted calorimetry curves for dDPPC and the 1:1 dDPPC:DLPC mixture were known at -5 °C, the molecular volume of DLPC in the gel phase was obtained and $\gamma$ was determined. Figure 3 shows the inverted calorimetry curve for DLPC resulting from the condition that all lipids are in the gel phase at -5 °C. In the figure, the extrapolated molecular volumes for dDPPC in the fluid phase below 37 °C and DLPC in the gel state above -5 °C are marked with red arrows. Table S5.2 in Section S5 of the Supplementary Information provides the linear equations that describe both the gel and fluid molecular volumes of dDPPC and DLPC for all temperatures.

With this information, it is now possible to obtain the temperature-dependent molecular SLDs for dDPPC and DLPC:

$$\text{SLD}_{\text{dDPPC}_j} = \frac{b_{\text{dDPPC}}}{V_{\text{dDPPC}_j}} \quad (4a)$$

$$\text{SLD}_{\text{DLPC}_j} = \frac{b_{\text{DLPC}}}{V_{\text{DLPC}_j}} \quad (4b)$$

The index j indicates gel (g) or fluid (f). b is the scattering length, which only depends on the atomic and isotopic make-up of, in this case, the lipid indicated in the subscript. The scattering lengths used in this work are listed in Table S2.1 in Section S2 of the Supplementary Information.

Fitting of the SANS data requires the SLD of the liquidus and solidus domains given by:

$$\text{SLD}_s = v_{\text{dDPPC}_s} \text{SLD}_{\text{dDPPC}_g} + (1 - v_{\text{dDPPC}_s}) \text{SLD}_{\text{DLPC}_g} \quad (5a)$$

$$\text{SLD}_l = v_{\text{dDPPC}_l} \text{SLD}_{\text{dDPPC}_f} + (1 - v_{\text{dDPPC}_l}) \text{SLD}_{\text{DLPC}_f} \quad (5b)$$
For all scattering data in the two-phase coexistence region ($T \le 24$ °C), the beads were allowed to take either a liquidus or a solidus SLD value: SLDs or SLDl (Equation (5)). The value for SLDl was directly obtained from the well-defined liquidus boundary of the phase diagram (Figure 2B and Equation (S4.17) in Section S4 of the Supplementary Information); however, because the phase diagram shows a discontinuity in the solidus boundary in the temperature range probed by SANS, the value of SLDs was varied by varying the ratio of dDPPC to DLPC in this phase (i.e., by varying the $f_s$ value in Equation (S4.18) in Section S4 of the Supplementary Information). As shown in Figures S8.1 and S8.2 in Section S8 of the Supplementary Information, we found that MONSA obtains several equivalent excellent fits, as assessed by $\chi^2$—the mean deviation between the fit and the data—for most $f_s$ values probed. Hence, $\chi^2$ could not inform us of the loci of the solidus boundary directly. In order to extract this information, we included two additional constraints; the first required that the volume fraction of the liquidus obtained from the fit coincided with the expected liquidus volume fraction for that given $f_s$ value (Equation (S4.16) in Section S4 of the Supplementary Information). The second constraint required that the molecular volume obtained from the fit (Equation (S4.19) in Section S4 of the Supplementary Information) coincided exactly with the molecular volume measured by densitometry/calorimetry (Figure 3). Plots of the difference between the fit and the expected liquidus volume fractions as a function of the $f_s$ value for the lower temperature data, which included 9 °C, 11.9 °C and 14.9 °C, produced only one possible solution for each temperature, as shown in Figure S8.1 in Section S8 of the Supplementary Information. Further, as shown in Figure 4B, these values of $f_s$ satisfied the molecular volume obtained by calorimetry/densitometry.

Although we successfully modeled the solidus as being composed of only gel-phase lipids in the low-temperature regime, we found that this model failed for the higher-temperature data, i.e., 18 °C, 22 °C and 24 °C. Instead, we modeled the solidus as composed of dDPPC in the gel phase and DLPC in some intermediate volume between the fluid state and the gel state using a sliding parameter $\alpha$:

$$V_{DLP C_s} = \alpha V_{DLP C_g} + (1 - \alpha) V_{DLP C_f} \quad (9)$$

where $\alpha = 1$ corresponds to the gel state and $\alpha = 0$ the fluid state.

DLPC is known to have a broad melting transition occurring over nearly 10 degrees, as shown in Figure S3.1C,F in Section S3 of the Supporting Information, and therefore, as the system approaches ≈ 27 °C, the molecular volume of DLPC could take a value between its fluid and gel state. Using this approach, we found that for the 18 °C data, MONSA fits satisfied the expected liquidus volume fraction condition for multiple values of $\alpha$, including $\alpha = 1$, which corresponds to the volume of DLPC to its gel state (see Figure S8.2A in Section S8 of the Supplementary Information). However, only for one value of $\alpha$ ($\alpha = 0.75$) did the molecular volume obtained (using Equation (S4.19) in Section S4 of the Supporting Information) coincide with the measured value (as shown in Figure S8.3A in Section S8 of the Supplementary Information). In the case of 22 °C and 24 °C, setting the volume of DLPC to its gel-state value in the solidus gave no solution (see Figure S8.2B in Section S8 of the Supplementary Information). For these temperatures, the most consistent solution suggested a DLPC volume closer to its fluid state, with $\alpha = 0.3$ (see Figures S8.2B and S8.3A in Section S8 of the Supplementary Information). Alternative models that varied the volume of both lipids did not produce fits that satisfied all constraints either. Indeed, this suggests that this behavior is due to the unique melting profile of DLPC (see Figure S3.1C,F in Section S3 of the Supplementary Information). In Section S8 of the Supplemental Information, Figure S8.3A shows the family of molecular volumes obtained, highlighting the effect of a sliding value of $\alpha$ for the higher-temperature data and, in Figure S8.3B, their concomitant position marking the solidus boundary in the phase diagram. Figure 4B shows the molecular volume for the 1:1 mixture obtained by densitometry and calorimetry with the corresponding best-fit molecular volumes obtained from the SANS analysis, and in Figure 4C are their loci in the phase diagram. Figure 4A shows the corresponding best fits against the scattering data.
SANS is a technique that detects nanoscale structures, and using real-space bead vesicle models allowed us to extract this information. Shown in Figure 5A are real-space bead vesicles whose structures produced the scattering curves shown in Figure 4A. Interestingly, the features observed are reminiscent of those observed in giant unilamellar vesicles made from a similar mixture, as shown in Figure S9.1 in Section S9 of the Supplementary Information. Due to the bead grid nature of the model, we were able to analyze the domain characteristics using the pair correlation function on the surface of the vesicle shell: $g(R_{\text{arc length}})$ [39]. As shown in Figure 5B, we obtained an oscillatory behavior for the bead density variation of the liquidus phase in the membrane relative to the mean density. The wavelength, $\lambda_{\text{domain}}$, of the domains was obtained by fitting the pair correlation function to an exponentially attenuated cosine function [40]. The wavelength of the oscillation in the pair correlation function is seen to increase with temperature, as highlighted by the arrow. The wavelength of the oscillation corresponds to a characteristic length—in this case, a characteristic domain size—as shown in Figure 5C.

The correlation peak positions in the data (highlighted with black dots in Figure S10.1 in Section S10 of the Supplementary Information) also reflect a characteristic length through the relation $2\pi/Q_{\text{peak}}$. We find that both correlation lengths nearly overlap for $T \le 18$ °C. For higher temperatures, the oscillation behavior of the pair correlation function is less obvious, and the fits show a more drastic increase in the characteristic length than the one derived from the position of the correlation peaks in the data as shown in Figure 5C.

## 4. Discussion

As shown above, the data and analysis presented here had the necessary combination of techniques to unveil the nature of the two phases in the coexistence region for the DPPC:DLPC system. With calorimetry, the phase diagram shows that this system forms a homogenous liquid phase at high temperatures. At lower temperatures, a region of coexistence between a liquidus phase and a solidus phase emerges. The upper boundary of the coexistence region, also referred to as the liquidus boundary, was found to be well defined; that is, the boundary shows a smooth and continuous variation in the composition of the liquidus phase as a function of temperature for all compositions ($0 \le \chi_{\text{DPPC}} \le 1$), as has been reported by Van Dijck et al. [3]. The lower or solidus boundary for $0 \le \chi_{\text{DPPC}} \le 0.6$ was found to have a very slight variation in composition with temperature but not a horizontal isothermal behavior as reported by Van Dijck et al. [3]. For $0.6 > \chi_{\text{DPPC}} \ge 0.7$, the calorimetry-derived phase diagram displays a gap or discontinuity in the solidus boundary. This behavior of the solidus boundary is also observed in the binary system DSPC:DMPC (distearoyl phosphocholine (DSPC) and dimyristoyl phosphatidylcholine (DMPC)) reported originally by Mabrey et al. [41] and others [42,43] and more recently by Losada-Perez et al. [44]. Because DSPC: DMPC and DPPC:DLPC are mixtures of saturated lipids that differ in tail length by four CH₂ groups, it is not surprising to find that their phase diagrams have the same features.

Continuity between the lower branch ($\chi_{\text{dDPPC}} < 0.6$) and the upper branch ($\chi_{\text{dDPPC}} > 0.7$) of the solidus boundary has been sought theoretically [3] and experimentally using nano-probes (Foster Resonance Energy Transfer) [28] as well as by fluorescence confocal microscopy using giant unilamellar vesicles [23], and they suggest that the lower branch ($\chi_{\text{dDPPC}} < 0.6$) of the solidus boundary continues smoothly—without a sharp turn where $0.6 < \chi_{\text{dDPPC}} < 0.7$—to $\chi_{\text{dDPPC}} = 1$, thus not coinciding with the upper branch described by the calorimetry result. Interestingly, consistency between molecular volumes and SANS fits resulted in a solidus boundary that connects the lower and upper solidus branches (see Figure 4C). Knoll and collaborators, in their exploration of the DSPC:DMPC phase diagram using a similar combination of approaches, calorimetry and densitometry [27] and SANS [45], found a similar discontinuity region suggested by densitometry and had a near-constant composition solidus boundary connecting the lower and upper solidus branches as well. SANS, however, was applied differently than was presented here. In their case, they performed solvent contrast variation experiments (varying D₂O/H₂O ratios) for two mixtures with
different initial compositions [45] from which they were able to deduce the solidus and liquidus boundaries for the system, albeit for only a small number of data points. Using this approach, they additionally observed that in the solidus region, below the solidus boundary, the solidus phase was composed of two solidus phases. From the molecular volume behavior for the 1:1 mixture, we could deduce that this also occurs in the DPPC:DLPC system. Although most phase diagrams with peritectic characteristics are shown with an isothermal (horizontal) solidus boundary branch, we found that in the DPPC:DLPC system, the lower branch of the solidus boundary ($\chi_{dDPPC} < 0.6$) is not horizontal but instead shows a gradually changing composition with temperature. Lozada-Perez et al. found similar behavior in the DSPC:DMPC system [44].

From the molecular volumes measured for the individual lipids as well as for the 1:1 mixture, we deduced that both DLPC and DPPC behaved gel-like in the solidus phase and liquid-like in the liquidus phase of the coexistence region. The combination of the molecular volume behavior for the 1:1 mixture and the SANS fits provided details on the nature of each phase; the liquidus phase was found to always consist of fluid-phase lipids while the solidus phase was found to be composed of gel-phase lipids only at low temperatures, and at higher temperatures (at and above 18 °C), DLPC was found to take intermediate volume states between its gel and fluid phase. To the best of our knowledge, this behavior had not been reported earlier.

The binary mixture DPPC:DLPC, which is not an ideal mixture, is described by regular solution theory. We follow the thermodynamic approach of Van Dijk et al. [3] to obtain the phenomenological thermodynamic curves that reflect the boundaries of the phase diagram and are represented in Equation (10).

$$T_{EGV}(\chi_{dDPPC}) = \frac{\frac{\Delta H_{DLPC}}{R}(1 - \chi_{dDPPC}) + \frac{\Delta H_{dDPPC}}{R}(\chi_{dDPPC}) + (A^L - A^S)(\chi_{dDPPC})(1 - \chi_{dDPPC})^n n^{-n}(n + 1)^{n+1}}{\frac{\Delta H_{DLPC}}{RT_{m-DLPC}}(1 - \chi_{dDPPC}) + \frac{\Delta H_{dDPPC}}{RT_{m-dDPPC}}(\chi_{dDPPC})} \quad (10)$$

Here, $\Delta H_{DLPC}$ and $\Delta H_{dDPPC}$ are the melting enthalpies for DLPC and dDPPC, respectively, and $T_{m-DLPC}$ and $T_{m-dDPPC}$ are their respective melting temperatures. In this equation, $n$ represents a nonlinear exponent of the excess enthalpy and provides the necessary asymmetry between the solidus and liquidus boundaries. The parameters $A^S$ and $A^L$ describe the maximum excess enthalpy for solidus and liquidus, respectively. Equation (10) describes the temperature vs. composition curve for the case when the free energy of the solidus and liquidus is equal. On the phase diagram, this curve lies approximately mid-way between the solidus and liquidus boundaries. We found that values of $n = 2$, $A^L = 112 \text{ cal} \cdot \text{mol}^{-1}$ and $A^S = 172 \text{ cal} \cdot \text{mol}^{-1}$ reasonably describe the boundaries as shown in the Figure 4C, which are not far from the values reported by Van Dijk et al. earlier ($n = 1.65$, $A^L = 75 \text{ cal} \cdot \text{mol}^{-1}$ and $A^S = 250 \text{ cal} \cdot \text{mol}^{-1}$). These values indicate that the maximum excess enthalpy difference between the liquidus and solidus, proportional to the difference $A^L - A^S$, is much smaller in our analysis than that found in the interpretation of Van Dijk et al. The polynomial representing the excess enthalpy ($(\chi_{dDPPC})(1 - \chi_{dDPPC})^n$), on the other hand, only differs minimally between $n = 2$ and $n = 1.65$.

SANS is a technique that can be used to extract structural information in addition to detecting phase separation in vesicles. Anghel et al. [18] studied the dDPPC:DLPC system above and below the miscibility transition using 30 nm in diameter vesicles to capture the characteristic peak from the domains by SANS. As in the present study, their goal was to extract structural information. They approached the analysis using analytical models such as round domains and stripes but were unable to show consistency between their analysis and the boundaries from the phase diagram reported by Van Dijk et al. Thus, using analytical shapes such as circular domains and stripes, which, according to GUV imaging (see References [22,23] and Figure S9.1 in Section S9 of the Supplementary Information), were unlikely, was an important limitation of their analysis. Using a bead vesicle as a real-space model to fit the SANS data from vesicles exhibiting phase separation allowed for a free-form analysis tool that yielded domains showing a phase (liquidus) evolving
into a network of elongated dendrite-like features, as shown in Figure 5A, and which are similar to those found in GUVs (see References [22,23] and Figure S9.1 in Section S9 of the Supplementary Information). Notwithstanding, we also pursued using round domains to analyze the SANS data using a script in Igor™, where we could generate round domain configurations (varying the number of domains but keeping the volume fractions constant). In the high-temperature regime (14.9 °C ≤ *T* ≤ 24 °C), round domains produced larger intensity curves than the data suggested. Lower intensity is attained in the free-form approach by domain structures that “finger” out, as shown in Figure S11.1 in Section S11 of the Supplementary Information. For the lower-temperature regime (9 °C ≤ *T* ≤ 12 °C), the congruency between the data and the calculated scattering for round domains was found to be significantly better. This congruency is partly explained by the ability to form dendrite structures with many small round domains (see Figure S11.1 in Section S11 of the Supplementary Information). In conclusion, given the robustness of the elongated features obtained through MONSA over repeated fits and the images of solidus and liquidus phases of the DPPC:DLPC system in GUVs and by comparing with what can be achieved with round domains, the real-space bead vesicle configurations obtained with MONSA appear physically reasonable.

It is known that this system exhibits a membrane thickness mismatch between the different domains [46]; that is, the solidus domains are expected to be thicker, due to a higher content of dDPPC, than the liquidus domains. In addition, it is known that membrane mismatch can drive domain size [17], where an increase in domain size is associated with an increase in membrane mismatch between domains. This latter result from Heberle et al. suggests that the decreasing characteristic domain size of the liquidus phase is correlated with a decrease in membrane mismatch between solidus and liquidus with decreasing temperature. This is reasonable since, as the temperature lowers, the composition of the solidus changes to include an increasing fraction of DLPC, which would result in the thinning of the solidus phase.

**5. Conclusions**

In closing, we have shown that the use of complementary techniques, which in this case consisted of calorimetry, densitometry and SANS, was necessary to obtain a revised phase diagram for the DPPC:DLPC system. We find that the lower ($\chi_{DPPC} < 0.6$) and higher ($\chi_{DPPC} > 0.7$) solidus boundary branches connect, contrary to previous reports [3,23,28]. Consistency between these separate methods revealed that the solidus phase was composed of gel-phase lipids only at low temperatures, and at higher temperatures, DLPC was found to take intermediate states between its gel- and fluid-phase volumes. This behavior for DLPC was found to be consistent with its unusual gel-to-fluid transition temperature behavior [41]. In addition, we obtained domain structure information using a free-form real-space bead vesicle model and found the domains' shapes to be consistent with a modulated phase [40], whose characteristic size decreases with decreasing temperature. This behavior correlates with a decreasing membrane mismatch between the solidus and liquidus phases due to a concomitant thinning of the solidus phase due to an increase in DLPC content and a thickening of the liquidus phase due to an increase in DPPC content [17].

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/membranes13030323/s1, Figure S1.1: Modification to extrusion metal cylinder; Figure S1.2: New Era Syringe Pump modification; Figure S2.1: SANS Contrast Match series; Figure S3.1: Heat Capacity Traces; Figure S5.1: Densitometry measurements of H2O and D2O; Figure S5.2: Comparison between the molecular volumes obtained in this study and previous studies; Figure S5.3: Errors for the molecular volumes; Figure S6.1: SANS traces for 50 nm vesicles of a 1:1 ratio of hydrogenated DPPC and hydrogenated DLPC in 100% D2O; Figure S6.2: Plot of fit parameters shown in Table S6.1; Figure S7.1: MONSA and SASView vesicle scattering calculations for a monodisperse population of vesicles; Figure S8.1: Volume fraction difference between expected obtained from the MONSA fits for three temperatures: 9C, 11.9C and 14.9C as well as corresponding
chi-square; Figure S8.2: Volume fraction difference between expected obtained from the MONSA fits for three temperatures: 18C, 22C and 24C as well as corresponding chi-square; Figure S8.3: Molecular volumes for the 1:1 dDPPC:DLPC mixture obtained by calorimetry, densitometry and MONSA fits to the SANS data and corresponding phase diagram; Figure S9.1: Fluorescence microscopy images of dDPPC:DLPC 1:1 GUVs at 30C and 24C; Figure S10.1: Correlation peak positions in the I vs Q data; Figure S11.1: Generated scattering curves from round domains and comparison to MONSA fits; Table S2.1: Scattering lengths used in this study; Table S3.1: Latent heats, and solidus and liquidus transition temperatures; Table S5.1: Parameters used to produce molecular volumes from inverted calorimetry traces; Table S5.2: Molecular Volume equations; Table S6.1: Fit parameters from SAS view for SANS traces in Figure S6.1; Table S7.1: ILL D22 SANS instrument parameters; Table S8.1: Volume fraction of the liquidus, $v_f$, as a function of temperature. References [3,22,30,33,34,36,37,41,47–57] are cited in Supplementary Materials.

**Author Contributions:** Conceptualization, U.P.-S. and L.P.; methodology, U.P.-S. and L.P.; software, U.P.-S. and L.P.; validation, U.P.-S. and L.P.; formal analysis, N.K and U.P.-S.; investigation, N.K, U.P.-S. and L.P.; resources, U.P.-S. and L.P.; data curation, U.P.-S. and L.P.; writing—original draft preparation, N.K.; writing—review and editing, N.K., U.P.-S. and L.P.; visualization, N.K. and U.P.-S.; supervision, U.P.-S.; project administration, U.P.-S.; funding acquisition, U.P.-S. All authors have read and agreed to the published version of the manuscript.

**Funding:** This study was funded by NSF, CAREER award DMR-1753238.

**Institutional Review Board Statement:** This article does not contain any studies with human participants or animals performed by any of the authors.

**Data Availability Statement:** The data presented in this study are available in the Supplementary Material.

**Acknowledgments:** The authors thank the Institute Laue Langevin for providing neutron beam time under proposals 9-13-428 and 9-13-510 U.P.-S. acknowledges travel support from the ILL to promote scientific collaboration with L. Porcar, as well as support for density measurements of lipids mixtures studied herein through equipment made available by the Partnership for Soft Condensed Matter (PSCM). U.P.-S. gratefully acknowledges the support from the NSF, CAREER award DMR-1753238.

**Conflicts of Interest:** Authors Krzyzanowski, Porcar and Perez-Salas declare no conflict of interest.

## References

1. Coskun, U.; Simons, K. Cell Membranes: The Lipid Perspective. *Structure* **2011**, *19*, 1543–1548. [CrossRef]

2. Feigenson, G.W. Phase behavior of lipid mixtures. *Nat. Chem. Biol.* **2006**, *2*, 560–563. [CrossRef] [PubMed]

3. Vandijck, P.W.M.; Kaper, A.J.; Oonk, H.A.J.; Degier, J. Miscibility Properties of Binary Phosphatidylcholine Mixtures—Calorimetric Study. *Biochim. Biophys. Acta (BBA) Biomembr.* **1977**, *470*, 58–69. [CrossRef]

4. Lee, A.G. Lipid phase transitions and phase diagrams. II. Mixtures involving lipids. *Biochim. Biophys. Acta (BBA) Rev Biomembr.* **1977**, *472*, 285–344. [CrossRef]

5. Lee, A.G. Lipid phase transitions and phase diagrams. I. Lipid phase transitions. *Biochim. Biophys. Acta (BBA) Rev Biomembr.* **1977**, *472*, 237–281. [CrossRef]

6. Lee, A.G. Calculation of Phase-Diagrams for Non-Ideal Mixtures of Lipids, and a Possible Nonrandom Distribution of Lipids in Lipid Mixtures in Liquid-Crystalline Phase. *Biochim. Biophys. Acta (BBA) Biomembr.* **1978**, *507*, 433–444. [CrossRef]

7. Almeida, P.F. Thermodynamics of lipid interactions in complex bilayers. *Biochim. Biophys. Acta (BBA) Biomembr.* **2009**, *1788*, 72–85. [CrossRef]

8. Simons, K.; Ikonen, E. Functional rafts in cell membranes. *Nature* **1997**, *387*, 569–572. [CrossRef]

9. Simons, K.; Toomre, D. Lipid rafts and signal transduction. *Nat. Rev. Mol. Cell Biol.* **2000**, *1*, 31–39. [CrossRef]

10. Lingwood, D.; Simons, K. Lipid Rafts As a Membrane-Organizing Principle. *Science* **2009**, *327*, 46–50. [CrossRef]

11. Toulmay, A.; Prinz, W.A. Direct imaging reveals stable, micrometer-scale lipid domains that segregate proteins in live cells. *J. Cell Biol.* **2013**, *202*, 35–44. [CrossRef]

12. Rayermann, S.P.; Rayermann, G.E.; Cornell, C.E.; Merz, A.J.; Keller, S.L. Hallmarks of Reversible Separation of Living, Unperturbed Cell Membranes into Two Liquid Phases. *Biophys. J.* **2017**, *113*, 2425–2432. [CrossRef]

13. Levental, I.; Levental, K.R.; Heberle, F.A. Lipid Rafts: Controversies Resolved, Mysteries Remain. *Trends Cell Biol.* **2020**, *30*, 341–353. [CrossRef]

14. Di Pasquale, M.; Nguyen, M.H.L.; Castillo, S.R.; Heberle, F.A.; Marquardt, D. Identifying Membrane Lateral Organization by Contrast-Matched Small Angle Neutron Scattering. *Methods Mol. Biol.* **2012**, *2402*, 163–177.

15. Pan, J.; Heberle, F.; Petruzielo, R.S.; Katsaras, J. Using small-angle neutron scattering to detect nanoscopic lipid domains. *Chem. Phys. Lipids* **2013**, *170–171*, 19–32. [CrossRef]
16. Nickels, J.; Chatterjee, S.; Stanley, C.; Qian, S.; Cheng, X.; Myles, D.A.A.; Standaert, R.F.; Elkins, J.G.; Katsaras, J. The in vivo structure of biological membranes and evidence for lipid domains. *PLoS Biol.* **2017**, *15*, e2002214. [CrossRef]

17. Heberle, F.A.; Petruzielo, R.S.; Pan, J.; Drazba, P.; Kučerka, N.; Standaert, R.F.; Feigenson, G.W.; Katsaras, J. Bilayer Thickness Mismatch Controls Domain Size in Model Membranes. *J. Am. Chem. Soc.* **2013**, *135*, 6853–6859. [CrossRef]

18. Anghel, V.N.P.; Kučerka, N.; Pencer, J.; Katsaras, J. Scattering from laterally heterogeneous vesicles. II. The form factor. *J. Appl. Crystallogr.* **2007**, *40*, 513–525. [CrossRef]

19. Anghel, V.N.P.; Bolmatov, D.; Katsaras, J. Models for randomly distributed nanoscopic domains on spherical vesicles. *Phys. Rev. E* **2018**, *97*, 062405. [CrossRef]

20. Anghel, V.N.P.; Bolmatov, D.; Katsaras, J.; Pencer, J. Domains on a Sphere: Neutron Scattering, Models, and Mathematical Formalism. *Chem. Phys. Lipids* **2019**, *222*, 47–50. [CrossRef]

21. Bolmatov, D.; McClintic, W.T.; Taylor, G.; Stanley, C.B.; Do, C.; Collier, C.P.; Leonenko, Z.; Lavrentovich, M.O.; Katsaras, J. Deciphering Melatonin-Stabilized Phase Separation in Phospholipid Bilayers. *Langmuir* **2019**, *35*, 12236–12245. [CrossRef] [PubMed]

22. Feigenson, G.W.; Buboltz, J.T. Ternary Phase Diagram of Dipalmitoyl-PC/Dilauroyl-PC/Cholesterol: Nanoscopic Domain Formation Driven by Cholesterol. *Biophys. J.* **2001**, *80*, 2775–2788. [CrossRef] [PubMed]

23. Fidorra, M.; Garcia, A.; Ipsen, J.H.; Hartel, S.; Bagatolli, L.A. Lipid domains in giant unilamellar vesicles and their correspondence with equilibrium thermodynamic phases: A quantitative fluorescence microscopy imaging approach. *Biochim. Biophys. Acta (BBA) Biomembr.* **2009**, *1788*, 2142–2149. [CrossRef] [PubMed]

24. Dewhurst, C. GRASP: Graphical Reduction and Analysis SANS Program for Matlab. 2007. Available online: https://www.ill.eu/users/support-labs-infrastructure/software-scientific-tools/grasp (accessed on 25 February 2023).

25. Feigin, L.A.S.; Svergun, D.I. *Structure Analysis by Small-Angle X-Ray and Neutron Scattering*; Springer: Berlin/Heidelberg, Germany, 1987.

26. Petruzielo, R.S.; Heberle, F.A.; Drazba, P.; Katsaras, J.; Feigenson, G.W. Phase behavior and domain size in sphingomyelin-containing lipid bilayers. *Biochim. Biophys. Acta (BBA) Biomembr.* **2013**, *1828*, 1302–1313. [CrossRef] [PubMed]

27. Schmidt, G.; Knoll, W. On the peritectic phase behavior of dimyristoyllecithin/distearoyllecithin mixtures. A densitometric study. *Chem. Phys. Lipids* **1986**, *39*, 329–339. [CrossRef]

28. Heberle, F.A.; Buboltz, J.T.; Stringer, D.; Feigenson, G.W. Fluorescence methods to detect phase boundaries in lipid bilayer mixtures. *Biochim. Biophys. Acta (BBA) Mol. Cell Res.* **2005**, *1746*, 186–192. [CrossRef]

29. Katsaras, J.; Epand, R.F. Absence of chiral domains in mixtures of dipalmitoylphosphatidylcholine molecules of opposite chirality. *Phys. Rev. E* **1997**, *55*, 3751–3753. [CrossRef]

30. Pencer, J.; Mills, T.; Anghel, V.N.P.; Krueger, S.; Epand, R.; Katsaras, J. Detection of submicron-sized raft-like domains in membranes by small-angle neutron scattering. *Eur. Phys. J. E* **2005**, *18*, 447–458. [CrossRef]

31. Nagle, J.F.; Cognet, P.; Dupuy, F.G.; Tristram-Nagle, S. Structure of gel phase DPPC determined by X-ray diffraction. *Chem. Phys. Lipids* **2019**, *218*, 168–177. [CrossRef]

32. Nagle, J.F.; Tristram-Nagle, S. Structure of lipid bilayers. *Biochim. Biophys. Acta (BBA) Rev Biomembr.* **2000**, *1469*, 159–195. [CrossRef]

33. Nagle, J.F.; Venable, R.M.; Maroclo-Kemmerling, E.; Tristram-Nagle, S.; Harper, P.E.; Pastor, R.W. Revisiting Volumes of Lipid Components in Bilayers. *J. Phys. Chem. B* **2019**, *123*, 2697–2709. [CrossRef]

34. Kučerka, N.; Nieh, M.-P.; Katsaras, J. Fluid phase lipid areas and bilayer thicknesses of commonly used phosphatidylcholines as a function of temperature. *Biochim. Biophys. Acta (BBA) Biomembr.* **2011**, *1808*, 2761–2771. [CrossRef]

35. Ebel, H.; Grabitz, P.; Heimburg, T. Enthalpy and Volume Changes in Lipid Membranes. I. The Proportionality of Heat and Volume Changes in the Lipid Melting Transition and Its Implication for the Elastic Constants. *J. Phys. Chem. B* **2001**, *105*, 7353–7360. [CrossRef]

36. Lemmon, E.W. *Thermophysical Properties of Fluid Systems*; National Institute of Standards and Technology: Gaithersburg, MD, USA, 2009.

37. Svergun, D. Restoring Low Resolution Structure of Biological Macromolecules from Solution Scattering Using Simulated Annealing. *Biophys. J.* **1999**, *76*, 2879–2886. [CrossRef]

38. Doucet, J.H.; Alina, G.; Bakker, J.; Bouwman, W.; Butler, P.; Campbell, K.; Gonzales, M.; Heenan, R.; Jackson, A.; Juhas, P.; et al. SasView, version 4.2.2. Released on 20 May 2019. Available online: https://www.sasview.org/ (accessed on 15 February 2023).

39. Chandler, D. *Introduction to Modern Statistical Mechanics*; Oxford University Press: Oxford, UK, 1987.

40. Amazon, J.J.; Feigenson, G.W. Lattice simulations of phase morphology on lipid bilayers: Renormalization, membrane shape, and electrostatic dipole interactions. *Phys. Rev. E* **2014**, *89*, 022702. [CrossRef]

41. Mabrey, S.; Sturtevant, J.M. Investigation of phase transitions of lipids and lipid mixtures by sensitivity differential scanning calorimetry. *Proc. Natl. Acad. Sci. USA* **1976**, *73*, 3862–3866. [CrossRef]

42. Wilkinson, D.A.; Nagle, J.F. Dilatometric study of binary mixtures of phosphatidylcholines. *Biochemistry* **1979**, *18*, 4244–4249. [CrossRef]

43. Gent, M.P.N.; Ho, C. Fluorine-19 nuclear magnetic resonance studies of lipid phase transitions in model and biological membranes. *Biochemistry* **1978**, *17*, 3023–3038. [CrossRef]
44. Losada-Perez, P.; Mertens, N.; de Medio-Vasconcelos, B.; Slenders, E.; Leys, J.; Peeters, M.; van Grinsven, B.; Gruber, J.; Glorieux, C.; Pfeiffer, H.; et al. Phase Transitions of Binary Lipid Mixtures: A Combined Study by Adiabatic Scanning Calorimetry and Quartz Crystal Microbalance with Dissipation Monitoring. *Adv. Cond. Matter. Phys.* **2015**, *2015*, 479318. [CrossRef]

45. Knoll, W.; Ibel, K.; Sackmann, E. Small-angle neutron scattering study of lipid phase diagrams by the contrast variation method. *Biochemistry* **1981**, *20*, 6379–6383. [CrossRef]

46. Lin, W.-C.; Blanchette, C.; Ratto, T.; Longo, M. Lipid Domains in Supported Lipid Bilayer for Atomic Force Microscopy. In *Methods in Membrane Lipids*; Dopico, A.M., Ed.; Humana Press: Totowa, NJ, USA, 2007; pp. 503–513.

47. Sears, V.F. Neutron scattering lengths and cross sections. *Neutron News* **1992**, *3*, 26–37. [CrossRef]

48. Morrow, M.R.; Davis, J.H. Calorimetric and Nuclear-Magnetic-Resonance Study of the Phase-Behavior of Dilauroylphosphatidylcholine Water. *Biochim. Biophys. Acta (BBA) Biomembr.* **1987**, *904*, 61–70. [CrossRef]

49. Zaytseva, Y.V.; Adichtchev, S.V.; Surovtsev, N.V. Raman study of temperature-induced hydrocarbon chain disorder in saturated phosphatidylcholines. *Chem. Phys. Lipids* **2020**, *230*, 104926. [CrossRef] [PubMed]

50. Finegold, L.; Shaw, W.A.; Singer, M.A. Unusual Phase Properties of Dilauryl Phosphatidylcholine (C12pc). *Chem. Phys. Lipids* **1990**, *53*, 177–184. [CrossRef]

51. Hatta, I.; Matuoka, S.; Singer, M.A.; Finegold, L. A new liquid crystalline phase in phosphatidylcholine bilayers as studied by X-ray diffraction. *Chem. Phys. Lipids* **1994**, *69*, 129–136. [CrossRef]

52. Ashkar, R.; Nagao, M.; Butler, P.D.; Woodka, A.C.; Sen, M.K.; Koga, T. Tuning Membrane Thickness Fluctuations in Model Lipid Bilayers. *Biophys. J.* **2015**, *109*, 106–112. [CrossRef]

53. Svergun, D.I.; Nierhaus, K.H. A map of protein-rRNA distribution in the 70 S Escherichia coli ribosome. *J. Biol. Chem.* **2000**, *275*, 14432–14439. [CrossRef]

54. Manalastas-Cantos, K.; Konarev, P.V.; Hajizadeh, N.R.; Kikhney, A.G.; Petoukhov, M.V.; Molodenskiy, D.S.; Panjkovich, A.; Mertens, H.D.T.; Gruzinov, A.; Borges, C.; et al. ATSAS 3.0: Expanded functionality and new tools for small-angle scattering data analysis. *J. Appl. Crystallogr.* **2021**, *54*, 343–355. [CrossRef]

55. Svergun, D.I.; Pethoukov, M.V. MONSA Manual. 2000. Available online: https://www.embl-hamburg.de/biosaxs/manuals/monsa.html#output (accessed on 22 February 2023).

56. Estes, D.J.; Mayer, M. Electroformation of giant liposomes from spin-coated films of lipids. *Colloids Surf. B Biointerfaces* **2005**, *42*, 115–123. [CrossRef]

57. Blosser, M.C.; Starr, J.B.; Turtle, C.W.; Ashcraft, J.; Keller, S.L. Minimal effect of lipid charge on membrane miscibility phase behavior in three ternary systems. *Biophys. J.* **2013**, *104*, 2629–2638. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.