Mathematical Methods for Images and SurfacesView this Special Issue
Characterization of the Lateral Distribution of Fluorescent Lipid in Binary-Constituent Lipid Monolayers by Principal Component Analysis
Lipid lateral organization in binary-constituent monolayers consisting of fluorescent and nonfluorescent lipids has been investigated by acquiring multiple emission spectra during measurement of each force-area isotherm. The emission spectra reflect BODIPY-labeled lipid surface concentration and lateral mixing with different nonfluorescent lipid species. Using principal component analysis (PCA) each spectrum could be approximated as the linear combination of only two principal vectors. One point on a plane could be associated with each spectrum, where the coordinates of the point are the coefficients of the linear combination. Points belonging to the same lipid constituents and experimental conditions form a curve on the plane, where each point belongs to a different mole fraction. The location and shape of the curve reflects the lateral organization of the fluorescent lipid mixed with a specific nonfluorescent lipid. The method provides massive data compression that preserves and emphasizes key information pertaining to lipid distribution in different lipid monolayer phases. Collectively, the capacity of PCA for handling large spectral data sets, the nanoscale resolution afforded by the fluorescence signal, and the inherent versatility of monolayers for characterization of lipid lateral interactions enable significantly enhanced resolution of lipid lateral organizational changes induced by different lipid compositions.
Biomembranes contain a variety of different lipids with a broad array of physicochemical properties. Such lipid variety has stimulated ideas regarding the existence of nonrandom lateral organizational states, including the so-called raft microdomains [1–7]. Model membranes, that is, bilayer vesicles and monolayer films, have proven useful for dissection of the lipid lateral distributional tendencies because simple combinations are possible and adjustments to lipid composition are relatively straightforward. With lipid monolayer systems, imaging of the monolayer surfaces to obtain direct insights into lipid lateral organizational states at macroscopic levels, that is, micron resolution, can be achieved using epifluorescence and/or Brewster angle microscopy [7–10]. These approaches have provided fundamental insights into the role(s) played by line tension and other parameters in stabilizing lipid macrodomains. However, their resolution capabilities are constrained to the micron range, limiting their effectiveness for detection of microdomains.
Recently, we showed that nanoscale changes in lipid packing and lateral organization can be detected in mixed monolayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1-palmitoyl-2-[7-(4, 4-difluoro-1,3,5,7-tetramethyl-4-bora-3a,4a-diaza-s-indacene-8-yl) heptanoyl]-sn-glycero-3-phosphocholine (Me4BODIPY-PC) by direct monitoring of changes in fluorescence emission from multiple spectra acquired simultaneously with the force-area isotherms [11, 12]. Me4BODIPY-PC mimics the behavior of fluid-phase PCs with unsaturated acyl chains and does not perturb lipid packing when the fluorophore concentration is kept low (e.g., ≤1 mole%). At higher mole fractions (e.g., 10 or 20 mole%), the probe begins to exert its own influence on the system rather than serving as a completely “impartial” reporter, a feature observed for virtually all probe molecules. Spectral characterization of Me4BODIPY-PC incorporated into lipid monolayers and bilayers indicates retention of excitation and emission transitions associated with monomeric forms of other BODIPY derivatives (i.e., absorption and emission wavelength maxima () near 498 nm and 506–515 nm, resp.). Me4BODIPY-PC also displays substantial broadening of its emission spectrum as a function of increasing fluorophore concentration due to elevated fluorescence near 570 nm relative to the 510–520 nm region. Although the photophysical mechanism accounting for the emission broadening may differ from the excited state dimer fluorescence of dimethyl-BODIPY derivatives (see Discussion and ), the response serves as an effective indicator of Me4BODIPY-PC local concentration and lateral distribution within lipid monolayers and bilayers. We recently showed that alterations in Me4BODIPY-PC packing density in POPC monolayers, achieved by lateral compression at constant Me4BODIPY-PC mole fraction or by variation of Me4BODIPY-PC mole fraction, were detectable in the emission spectra even though corresponding force-area isotherms indicated no deviation from ideal mixing behavior . The findings suggested that direct monitoring of fluorescence spectral changes could provide a highly sensitive tool for studying lipid lateral interactions in monolayers. However, to extract the information from the multiple emission spectra, that is, as many as 600 spectra per isotherm, an effective means is needed to analyze the data.
When the input data is extremely large and/or contains redundancy, then transformation into a reduced representation provides an effective means for preservation of essential information. Several data reduction methods are available [13–17] such as principal component analysis (PCA), semidefinite embedding, partial least square, and multidimensional scaling. PCA was designed by Karl Pearson  to compress and identify unknown trends in a multidimensional data set. Data dimensionality is reduced by performing a covariance analysis. Herein, we show that PCA provides an effective way to analyze fluorescence spectra of binary constituent lipid monolayers. By using this method, each spectrum can be represented by one point on a plane. We find that points, belonging to the same lipid constituents and experimental conditions, form a curve on the plane where each point belongs to a different mole fraction. The location and shape of the curve are characteristics to the lateral organization of the respective monolayer. Use of phosphatidylcholine (PC) or sphingomyelin (SM), containing either saturated-monounsaturated or disaturated hydrocarbon chains, as the nonfluorescent lipid in the binary constituent monolayers enables evaluation of the role played by matrix lipid phase state and composition in controlling the lateral mixing of the fluorescent lipid by PCA.
2. Materials and Methods
1-Palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1,2-dipalmitoylphosphatidylcholine (DPPC) were purchased from Avanti Polar Lipids (Alabaster AL). Palmitoyl sphingomyelin (PSM) and oleoyl sphingomyelin (OSM) were produced by reacylation of lyso-SM with the desired fatty acyl residue and purified as described previously [18, 19]. 1-Palmitoyl-2-[7-(4,4-difluoro-1,3,5,7-tetramethyl-4-bora-3a, 4a-diaza-s-indacene-8-yl)heptanoyl]-sn- glycero-3-phosphocholine (Me4BODIPY-PC) was synthesized and purified as described earlier .
2.1. Fluorescence Spectra of Binary-Constituent Lipid Monolayers
The Langmuir film balance platform was outfitted for acquisition of fluorescence emission intensity as a function of wavelength while simultaneously measuring surface pressure (π) and dipole potential (V) as a function of lipid cross-sectional molecular area (A) of the lipid monolayer [11, 12]. Briefly, BODIPY lipid films were excited at a incident angle using 488 nm unpolarized light from an argon-ion laser (Model 2122-45L, JDS Uniphase, San Jose, CA) equipped with a model-3 light-intensity controller and a fiber optic coupler (Model HPUC-23-488-S-3, FAC-2BL; Oz Optics, Nepean, ON, Canada). Fluorescence emission was collected perpendicular to the interface at a distance of ~1 cm using a fiber optic spectrometer (Model PC2000-ISA, Ocean Optics, Dunedin, FL) equipped with an L2 lens and 200-m slit. A 500-nm-long pass filter (500EFLP, Omega Optical, Brattleboro, VT) was mounted between the emission collimator and the detector to reduce scattered excitation light. Fluorescence emission spectral intensities were collected each second. Although monolayer compression was continuous during the spectral data acquisition cycle, the fractional change in lipid concentration during each acquisition cycle was ≤0.0073. Control emission spectra were unaffected by gas phase, that is, air or argon, or by 0.01% sodium azide in the subphase buffer.
Subphase buffer was maintained at C via a thermostated, circulating water bath and was produced using water purified by reverse osmosis, activated charcoal adsorption, mixed-bed deionization, then passed through a Milli-Q UV Plus System (Millipore Corp., Bedford, MA), and filtered through a 0.22 m Millipak 40 membrane. Subphase buffer contained 10 mM potassium phosphate (pH 6.6), 100 mM NaCl, and 0.01% NaN3 and was kept stored under argon, cleaned by passage through a seven-stage series filtration set-up consisting of an Alltech activated charcoal gas purifier, a LabClean filter, and a series of Balston disposable filters consisting of two adsorption (carbon) and three filter units (93% and 99.99% efficiency at 0.1 m). The film balance was housed in an isolated laboratory supplied with clean air by a Bioclean Air Filtration system equipped with charcoal and HEPA filters and was kept under humidified argon in a separate enclosure. Other features contributing to isotherm reproducibility include automated lipid spreading via a modified HPLC autoinjector, automated surface cleaning by multiple barrier sweeps between runs, and highly accurate, reproducible setting of the subphase level by an automated aspirator. Glassware was acid cleaned, rinsed with purified water, and then with hexane/ethanol (95: 5) before use.
Lipid monolayers were formed by automated spreading (51.67 L aliquots) of mixtures made from stock solutions dissolved in toluene/ethanol (5:6) or hexane/isopropanol/water (70: 30: 2.5). After spreading on the subphase surface and a delay period of 4 min, lipid films were compressed at a rate of ≤4 /molecule/min. Surface pressure and area calibration of the film balance were performed as detailed previously . Phospholipid concentrations, including BPC, were determined by a modified Bartlett assay . Solvent purity was verified by dipole potential measurements using a 210Po ionizing electrode .
2.2. Principal Component Analysis of the Fluorescence Spectra
Here we provide a short description of PCA. The measured fluorescence spectra can be represented by S, an nxm matrix, where m is the number of measured spectra. Each spectrum is measured at the same n wavelengths. Sij is the intensity of fluorescence emission measured at the ith wavelength of the jth spectrum. The covariance matrix C, a symmetric mxm matrix, can be created from the S matrix as follows:
where is the average intensity of the jth spectrum. There are m eigenvalues and m respective eigenvectors for the covariance matrix. The eigenvalues of C are real positive numbers and are sorted in descending order in the eigenvalue vector, . From the eigenvectors, one can create an mxm matrix, V where each column represents an eigenvector. In V, the eigenvectors are sorted as in , that is, the eigenvector the eigenvector in the first column has the largest eigenvalue; in the second column, the next largest eigenvalue, and so forth . By multiplying the S and V matrices and transposing the product, one obtains the mxn final data matrix, P = (SV)T. We refer to the rows of the final data matrix as principal vectors. The kth element of the ith principal vector is
In the case of our spectra, the first two eigenvalues in are much larger than the rest of the eigenvalues, as shown in Table 2. Thus each fluorescent spectrum can be approximated by the linear combination of the first two principal vectors, that is
or if the intensity of each spectra at the first wavelength is zero, then
Fluorescence spectra of lipid monolayers containing Me4BODIPY-PC as one of the two lipid constituents were acquired while simultaneously measuring the surface pressure versus average molecular area isotherm (Figure 1, upper panel). In every case, Me4BODIPY-PC was one lipid constituent (Figure 1, lower panel), while the type and amount of the second lipid constituent (POPC, DPPC, OSM, and PSM) varied. This experimental strategy provided insights into the effects of lipid phase state and composition on Me4BODIPY-PC lateral mixing. Figure 2 serves as an example, showing the fluorescence spectra of Me4BODIPY-PC/POPC (Figures 2(a) and 2(c)) and Me4BODIPY-PC/DPPC (Figures 2(b) and 2(d)) monolayers obtained at two representative surface pressures of 5 and 30 mN/m. It is noteworthy that surface pressures in the 30–35 mN/m region mimic conditions found in biomembranes [22, 23]. For each panel in Figure 2, the four spectra correspond to the following Me4BODIPY-PC mole fractions: 0.01 (red), 0.1 (blue), 0.2 (green), 1 (black).
Interpretation of the Me4BODIPY-PC spectral responses in Figure 2 is facilitated by consideration of the phase behavior of the lipids comprising the binary-constituent lipid monolayers. POPC and Me4BODIPY-PC both are known to display fluid phase behavior at all surface pressures below monolayer collapse (e.g., ~45 mN/m), while DPPC displays fluid phase behavior only at low surface pressures (e.g., 8 mN/m ) and undergoes a two-dimensional phase transition resulting in ordered lipid hydrocarbon chains at high surface pressures (e.g., [11, 18, 24]). At 5 mN/m (Figures 2(c) and 2(d)), both lipid constituents display fluid-phase behavior and the emission spectra are similar in intensity and shape at equivalent Me4BODIPY-PC mole fractions. However, at 30 mN/m (Figures 2(a) and 2(b)), only Me4BODIPY-PC and POPC remain fluid; whereas, DPPC is gel-like and chain ordered [11, 18, 24]. Under the latter condition, dramatic changes are observed in Me4BODIPY-PC/DPPC spectral intensity and shape, that is, quenching and broadening, compared to equivalent Me4BODIPY-PC mole fractions in the Me4BODIPY-PC/POPC mixed monolayers.
With BODIPY fluorophores, dilute noninteracting monomers exhibit a narrow emission peak centered in the 506–515 nm range after excitation [11, 12, 25–27]. At sufficiently high BODIPY concentrations, additional absorption and emission peaks are observed, often reflecting dimer emission and/or homo resonance energy transfer (Förster distance = 57 ). Johansson and colleagues have shown that emission of ground-state dimethyl-BODIPY dimers, denoted (J-dimer) occurs when the BODIPY rings are oriented in planes with their transition dipoles aligned at . dimers absorb near 570 nm and emit near 630 nm . Energy transfer to the ground-state dimers from excited-state monomers is typically responsible for the emission peak observed near 630 nm. A second type of ground-state dimer, denoted (H-dimer) and characterized by sandwich-like stacking of the BODIPY rings, results in parallel alignment of the transition dipoles and absorption near 477 nm, but produces no fluorescence emission after excitation [25, 26]. In the case of Me4-BODIPY-PC, we find no evidence of the 620–630 nm peak characteristic of dimethyl-BODIPY dimer fluorescence but rather the presence of an emission shoulder (~570 nm) at high Me4BODIPY-PC surface concentrations. Substantial spectral broadening occurs because of increased fluorescence near 570 nm relative to the 510–520 nm region. Thus, detection of surface concentration changes and lateral heterogeneity in the mixing of Me4BODIPY-PC in lipid monolayers and bilayers becomes evident.
In addition to formation of the emission shoulder, intensity decreases occur in the ~510 nm emission peaks at 20 mole% compared to 10 mole% Me4BODIPY-PC (Figures 2(a), 2(c), and 2(d)). This response could reflect nonfluorescent relaxation occurring for excited monomers by a mechanism involving Förster energy transfer between excited state monomer and ground state dimers aligned differently than dimers [25, 26]. In contrast to the above cases, the spectra in Figure 2(b) are broad at every Me4BODIPY-PC mole fraction, that is, fluorophore dimers form even at very low probe concentration, suggesting highly nonideal mixing and possible lateral phase separation of Me4BODIPY-PC in the DPPC monolayers at 30 mN/m (Figure 2(b)). Under such conditions, the specific area of Me4BODIPY-PC (30 mN/m; 80 /molecule)  is considerably larger than the DPPC specific area (46.2 /moL) [18, 24]. Me4BODIPY-PC molecules are unable to fit into the tight lattice of DPPC molecules. Breaking this lattice structure at many points is not energetically feasible. Thus, even at low concentration, Me4BODIPY-PC molecules are expected to segregate, forming clusters in the monolayer. In the clusters, the local probe concentration is high, where dimers (or excimers) may form [25–28], resulting in spectral broadening (Figure 2(b); red curve). The preceding discussion of spectral and force-area data facilitates interpretation of the results of the principal component analyses that follow.
In general, analogous monolayer phase behavior is exhibited by 18:1 sphingomyelin (SM) as by POPC in that 18:1-SM remains fluid at all surface pressures. However, at 30 mN/m, 16:0 SM exists as a mixture of coexisting gel-like chained-ordered and fluid chain-disordered phases, in contrast to DPPC which exists as a gel-like, chain-ordered phase (e.g, [18, 24]). It is noteworthy that the spectra of Me4BODIPY-PC/18:1-SM and Me4BODIPY-PC/16:0-SM (not shown) are broadened at surface pressures of 5 and 30 mN/m as well as at average molecular cross-sectional areas of 64 /molecule and 75 /molecule with increasing probe concentration, reflecting subtle differences in lipid structure.
Table 1 shows that a total of 52 fluorescence spectra were simultaneously analyzed by PCA using the spectra of 13 binary mixtures taken at two lateral pressures and two specific areas. We determined the eigenvalues and eigenvectors of the respective covariance matrix, C, by reducing the symmetric matrix to tridiagonal form [16, 17]. The first six elements of the sorted eigenvalue vector, are shown in Table 2.
Since the first two eigenvalues are much larger than the rest of the 52 eigenvalues, one is able to characterize each fluorescence spectrum by only the first two principal vectors (first two rows of the final data matrix, P). This is a massive dimensionality reduction from 52 to 2. In Figure 3, the components of the first two principal vectors are plotted against the respective wavelengths. The first principal vector (red curve) shows strong similarity to spectra measured at low Me4BODIPY-PC mole fractions with the exception of the spectrum (red curve) shown in Figure 2(b). The peak location (510 nm) of the first principal vector and its half width (40 nm) are similar to that of spectra taken at 1 mole% Me4BODIPY-PC concentration in Figures 2(a), 2(c), and 2(d). As was pointed out above, these spectra originate from fluorescence emission of the excited monomer probe. The second principal vector (Figure 2, blue curve) resembles higher-order fluorescence, that is, dimer or excimer, because it peaks at 572 nm where the spectral shoulder appears at high Me4BODIPY-PC concentrations. Note that each spectrum taken at high Me4BODIPY-PC concentration is expected to result in both monomer and dimer (or excimer) fluorescence because dimer formation from monomers is a reversible interaction [25–28].
In order to demonstrate the validity of the PCA approach, the measured intensities of the 52 fluorescence spectra were plotted against the calculated intensities using (4) (Figure 4). A straight line could be fitted to the points with high confidence. The slope and the intercept of the fitted line are and , respectively, while the correlation coefficient of the fit is 0.9997. Consequently, the spectra of BODIPY-labeled, two-constituent lipid monolayers can be approximated by the linear combinations of the first two principal vectors given by (4).
In the case of the jth spectrum, the coefficients of the linear combinations are the jth elements of the first two eigenvectors, and . There is a strong positive correlation between the value of and the peak height of the jth spectrum, while correlates with the half-width of the spectrum. In Figure 5(a), the elements of the first eigenvector are plotted against the respective peak height. The correlation coefficient between and the peak height is 0.96. In Figure 5(b), the elements of the second eigenvector are plotted against the half-width of the respective spectrum. The correlation coefficient between and the half-width is 0.91. Note that in Figure 5(b), the point that deviates most from the general trend, with coordinates (261,0.5), belongs to the spectrum of pure Me4BODIPY-PC monolayer at 30 mN/m.
According to (4), the jth elements of the first two eigenvectors, and , fully specify the jth spectrum, while the first two principal vectors are characteristics of all the measured spectra. Thus, geometrically each spectrum can be represented by one point in a plane with coordinate axes: and . For example, in the case of the jth spectrum, the coordinates of the respective point are and . Figures 6(a), 6(b), 6(c), and 6(d) show the planes with points representing different spectra, measured at surface pressures of 5 mN/m and 30 mN/m, and at average molecular cross-sectional areas of 75 /molecule and 64 /molecule, respectively. In each panel, points belonging to monolayers containing the same lipid constituents are connected by the same color line: Me4BODIPY-PC/POPC (black); Me4BODIPY-PC/DPPC (green); Me4BODIPY-PC/18:1-SM (blue); Me4BODIPY-PC/16:0-SM (red). The zig-zag nature of the line connecting the points illustrates how the characteristic point representing the spectrum travels as the mole fraction of Me4BODIPY-PC increases. When the probe’s mole fraction approaches zero, the characteristic point of each line is expected to approach the origin (, ). This is the case because, with decreasing Me4BODIPY-PC concentration, the fluorescence intensity at any wavelength approaches zero, , and thus in (4), both and should approach zero. (Note that the principal vectors are independent of the Me4BODIPY-PC concentration.) Thus in spite of the fact that there is no spectrum at zero fluorophore concentration, each line can be expected to start from the origin.
It was pointed out above that the and coordinates of the characteristic point correlate with the half-width and height of the spectrum, respectively, and these spectral properties depend on dimer (or excimer) formation by the fluorophore. Since the fluorescence spectrum depends on the surface concentration, that is, mole fraction, of the fluorescent lipid in the respective monolayer, one can expect the location and shape of the lines in the plane to reflect changes in Me4BODIPY-PC surface concentration and lateral organization. It is noteworthy that two distinctly different line shapes are observed in Figure 6.(1)The initial slope of the line is negative, starting from zero fluorophore content, if the monolayer is in fluid or fluid/gel mixed phase. This situation occurs for all lipids shown in Figure 6, except for the green curve (DPPC) in Figure 6(b) (30 mN/m). At low-mole fraction, the Me4BODIPY-PC is expected to mix ideally with the fluid phase lipids and to be monomerically dispersed. With increasing fluorophore concentration, the number of monomers increases (i.e., the peak height and increases), and the half-width of the spectrum slightly decreases (i.e., decreases). Between 10 and 20 mole% Me4BODIPY-PC, dimers are expected to substantially increase in the monolayer, resulting in an emission shoulder at 570 nm (i.e., the half-width and increases). On the other hand, dimer formation makes possible nonfluorescent energy transfer between excited state monomer and ground state dimer (i.e., the peak height and decreases).(2)The initial slope of the line is positive, starting from zero fluorophore content, if the monolayer is in gel phase (Figure 6(b), green curve). Even at low-mole fraction, Me4BODIPY-PC mixes nonideally with gel phase DPPC, and forms aggregates. Thus increasing monomer and dimer fluorescence can be detected simultaneously from low-fluorophore concentration, that is, and increase simultaneously with increasing peak height and half-width. The response suggests that dimer abundance can be affected by the lipid lateral organizational state.
Under fluid-phase monolayer conditions, curves are highly similar for same-type lipids (e.g., POPC and DPPC; OSM and PSM). Nonetheless, detectable differences exist between the lipid types, consistent with subtle differences in the lateral mixing of Me4BODIPY-PC in each lipid type. The differences in curves generated from spectra acquired at 30 mN/m imply that both lipid type and phase state affect Me4BODIPY-PC lateral distribution under lipid packing conditions approximating biomembranes. Currently underway are comprehensive studies involving: more lipid mixing ratios, different BODIPY-labeled lipids, and other nonfluorescent lipid components.
Collectively, our data demonstrate the usefulness of PCA for deciphering fluorescence spectra obtained from lipid monolayers with respect to the lateral distribution of lipid fluorophore. The exceptional promise shown by PCA for handling large spectral data sets, the nanoscale resolution afforded by the fluorescence signal, and the inherent versatility of lipid monolayers for characterization of lipid lateral interactions combine to provide a new approach expected to facilitate investigation of lipid lateral organizational changes induced by differing lipid compositions and by proteins.
The authors are grateful for the continuing support of NIH/NIGMS 45928, NIH/NHLBI 49180, and the Hormel Foundation. This work was also supported by Russian Foundation for Basic Research Grant 09-04-00313 and by contract NIH/NIAID HHSN266200500021C. They thank Maureen M. Momsen for technical help with the monolayer experiments.
T. E. Thompson and T. W. Tillack, “Organization of glycosphingolipids in bilayers and plasma membranes of mammalian cells,” Annual Review of Biophysics and Biophysical Chemistry, vol. 14, pp. 361–386, 1985.View at: Google Scholar
R. E. Brown, “Sphingolipid organization in biomembranes: what physical studies of model membranes reveal,” Journal of Cell Science, vol. 111, no. 1, pp. 1–9, 1998.View at: Google Scholar
H. M. McConnell, “Structures and transitions in lipid monolayers at the air-water interface,” Annual Review of Physical Chemistry, vol. 42, no. 1, pp. 171–195, 1991.View at: Google Scholar
H. Möhwald, “Phospholipid monolayers,” in Phospholipids Handbook, G. Cevc, Ed., chapter 16, pp. 579–601, Marcel Dekker, New York, NY, USA, 1993.View at: Google Scholar
M. Dahim, N. K. Mizuno, X.-M. Li, W. E. Momsen, M. M. Momsen, and H. L. Brockman, “Physical and photophysical characterization of a BODIPY phosphatidylcholine as a membrane probe,” Biophysical Journal, vol. 83, no. 3, pp. 1511–1524, 2002.View at: Google Scholar
I. Guyon, S. Gunn, M. Nikravesh, and L. A. Zadeh, Eds., Feature Extraction: Foundations and Applications, vol. 207 of Studies in Fuzziness and Soft Computing, Springer, Berlin, Germany, 2006.
G. H. Golub and C. F. van Loan, Matrix Computations, The Johns Hopkins University Press, 3rd edition, 1996.
W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling, “Eigensystems,” in Numerical Recipes, chapter 11, 1989.View at: Google Scholar
K. Pearson, “On lines and planes of closest fit to systems of points in space,” Philosophical Magazine, vol. 2, pp. 559–572, 1901.View at: Google Scholar
L. I. Smith, “A tutorial on principal component analysis,” 2002, http://www.cs.otago.ac.nz/cosc453/student_tutorials/principal_components.pdf.View at: Google Scholar
J. M. Smaby, V. S. Kulkarni, M. Momsen, and R. E. Brown, “The interfacial elastic packing interactions of galactosylceramides, sphingomyelins, and phosphatidylcholines,” Biophysical Journal, vol. 70, no. 2, pp. 868–877, 1996.View at: Google Scholar
X.-M. Li, J. M. Smaby, M. M. Momsen, H. L. Brockman, and R. E. Brown, “Sphingomyelin interfacial behavior: the impact of changing acyl chain composition,” Biophysical Journal, vol. 78, no. 4, pp. 1921–1931, 2000.View at: Google Scholar
R. E. Brown and H. L. Brockman, “Using monomolecular films to characterize lipid lateral interactions,” Methods in Molecular Biology, vol. 398, pp. 41–58, 2007.View at: Google Scholar
R. A. Demel, W. S. M. Geurts Van Kessel, R. F. A. Zwaal, B. Roelofsen, and L. L. M. Van Deenen, “Relation between various phospholipase actions on human red cell membranes and the interfacial phospholipid pressure in monolayers,” Biochimica et Biophysica Acta, vol. 406, no. 1, pp. 97–107, 1975.View at: Google Scholar
F. Bergström, I. Mikhalyov, P. Hägglöf, R. Wortmann, T. Ny, and L. B.-Å. Johansson, “Dimers of dipyrrometheneboron difluoride (BODIPY) with light spectroscopic applications in chemistry and biology,” Journal of the American Chemical Society, vol. 124, no. 2, pp. 196–204, 2002.View at: Publisher Site | Google Scholar
J. B. Birks, “Excimer fluorescence. II. Lifetime studies of pyrene solutions,” Proceedings of the Royal Society A, vol. 275, pp. 575–588, 1963.View at: Google Scholar