Kinetic Evidence Suggests Spinodal Phase Separation in Stratum

Apr 4, 2014 - barrier to skin permeability, has been extensively studied, only limited data are extant concerning the kinetic mechanism for the format...
0 downloads 0 Views 2MB Size
Article pubs.acs.org/JPCB

Kinetic Evidence Suggests Spinodal Phase Separation in Stratum Corneum Models by IR Spectroscopy Richard Mendelsohn,*,† Ibrahim Selevany,† David J. Moore,† M. Catherine Mack Correa,‡ Guangru Mao,‡ Russel M. Walters,‡ and Carol R. Flach† †

Department of Chemistry, Newark College, Rutgers University, Newark, New Jersey 07102, United States Johnson & Johnson Consumer Companies, Inc., 199 Grandview Road, Skillman, New Jersey 08558, United States



ABSTRACT: Although lipid structure in models for the stratum corneum (SC), the main barrier to skin permeability, has been extensively studied, only limited data are extant concerning the kinetic mechanism for the formation of domains, lamellar phases, and lipid packing motifs. Such information would be of substantial interest in the characterization of the effects of disease states which disrupt the barrier. Kinetic IR spectroscopy measurements probed the temporal sequence of molecular events producing ordered structures in a three-component SC model of equimolar ceramide[NS] (cer[NS]), perdeuterated stearic acid-d35 (SA-d35), and cholesterol. Samples, heated above Tm, were quenched to 31 °C, and then spectra were recorded at ∼15 min intervals for a total of 20−150 h. IR provides unique molecular structure information about headgroup H-bonding, lipid packing, and lipid chain order. The following sequence for phase separation was observed: (1) Formation of ceramide amide H-bonds from disordered forms to ordered structures (0.5−4 h); (2) appearance of ordered ceramide chains with some orthorhombically packed structures (0.5−8 h); and (3) phase separation of large orthorhombic domains of SA-d35 (4−10 h). A spinodal decomposition mechanism, defined by continuous composition changes during the phase separation, suggests a qualitative description for these events.



INTRODUCTION The primary permeability barrier in human skin resides within the ∼15 μm thick outermost layer, the stratum corneum (SC), which is composed of anucleated corneocytes embedded in a continuous lamellar lipid matrix. Lipid structural organization in the SC has been widely studied with physical methods (e.g., see refs 1−10), and the results have begun to provide a basis for understanding barrier function with the potential for elucidating how skin lipid organization is altered in human disease states. Beginning approximately 20 years ago, smallangle X-ray scattering has been used6,9,11 to identify the existence of two lamellar phases with periodicities of ∼6 and 13 nm, respectively termed the short (SPP) and long (LPP) periodicity phases. In addition to the lamellar organization, the lateral organization of SC lipids has been studied. IR spectroscopy12−17 and X-ray diffraction18−20 measurements on SC and SC lipid models reveal the presence of a substantial fraction of orthorhombic chain packing. Fatty acids are considered to be important for formation of this phase. Thermotropic studies on SC also reveal the occurrence of an orthorhombic → hexagonal packing transition near the physiological temperature of skin. This transition is followed at higher, non physiological temperatures by a hexagonal → disorder transition in the lipid chains.5,12,21,22 Yet, structural studies alone are insufficient to provide a complete description of skin organization/function. Thus, the SC is a dynamic structure, actively responding to the environment, and turns over in an approximately two week © 2014 American Chemical Society

period. Simply stated, SC lipids end up highly ordered in a healthy skin barrier. They probably did not start that way earlier in the temporal sequence, and they can be perturbed from this state. In spite of the need for kinetic studies of SC assembly, investigations of the mechanisms of the development of lamellar phases, lipid packing motifs, and domains are substantially less well elaborated than SC structural studies. Yet such information would be of interest in a variety of areas within biophysics in general, and within skin biophysics in particular. For instance, functions controlled by biomembranes such as lipid or protein transport, intracellular signaling, and regulation of integral membrane proteins have been suggested to involve ordered domains termed “lipid rafts” (e.g., refs 23 and 24). In skin biophysics, the occurrence of ordered lipid domains is a central feature in descriptions of SC lipid organization such as Forslind’s domain mosaic model25 which depicts the barrier lipids as segregated into ordered structures surrounded by grain boundaries formed from lipids in a more fluid (liquid crystalline) state that may provide control of hydration and transport levels. Elucidation at the molecular level of the physical mechanisms and kinetics of lipid domain formation would also help define potential pathways for permeation of hydrophobic species as well as provide clues as to why such a variety of ceramides and additional chemical species are present in the skin barrier. Understanding the Received: January 28, 2014 Revised: April 2, 2014 Published: April 4, 2014 4378

dx.doi.org/10.1021/jp501003c | J. Phys. Chem. B 2014, 118, 4378−4387

The Journal of Physical Chemistry B

Article

formation and dissipation kinetics of domains formed from particular ceramide and fatty acid species would aid in defining biological processes such as epidermal desquamation, which includes degradation of layered structures in the intercellular spaces. In addition, Oren et al. showed26 that antimicrobial peptides important in the skin’s ability to repair itself are found within the SC extracellular matrix. The role of enzymes and the dynamics of lipids at the SG−SC interface in controlling in the formation of the barrier has been reviewed by Elias.27 His group also28 investigated the colocalization of phospholipases with polar lipids, suggested to induce fusion of lamellar body discs as a prelude to formation of the SC. The presence of variable lipid domains may also provide one method by which enzyme activity is controlled at different depths within the SC.29,30 Finally, in a medical context, abnormal lamellar structures and disrupted lateral order have been associated with disrupted barrier function in patients with psoriasis, ichthyosis, or atopic dermatitis.31−34 Spinodal decomposition,35 extensively applied to the study of alloy36−38 and polymer phase separation,39−44 may provide a general framework for understanding the formation of complex layered structures, but has not, to our knowledge, been considered in skin biophysics. The closest reported studies with molecules related to those in the SC are vibrational (IR and Raman) spectroscopy,45−49 and small-angle neutron and X-ray scattering50 investigations of binary n-alkane mixtures quenched from melted states. In particular, Snyder and co-workers46−49 tracked the kinetics of demixing by utilizing a variety of IR spectral features to evaluate chain conformational ordering, packing motifs, and domain formation. The dependence of demixing rates on temperature and isotopic composition of alkanes in the C28−C36 length range have been probed. In the SC, at least 11 subclasses of ceramides have been characterized and coexist with relatively large levels of cholesterol, free fatty acids, and relatively minor amounts of cholesterol sulfate, glucosylceramides, and cholesterol esters. The current work is directed toward developing IR spectroscopy for studying lateral phase separation and the early stages of ordered lamellar structure formation in skin lipid models. The complexity of the SC lipid structure will necessitate examination of a large number of skin lipid mixtures over a wide range of compositions and temperatures to evaluate the possible relevance of spinodal mechanisms. For these initial studies, we have selected two previously studied ternary equimolar mixtures, namely ceramide [NS]/ stearic acid/cholesterol and ceramide[AS]/stearic acid/cholesterol (see Figure 1 for structures) for which currently available physical measurements on the same or similar systems51,52 at least begin to define the existence regime for various phases. We are well aware that most studies of phase separation in molecular systems utilize relatively simple binary mixtures whose phase diagrams are or can be readily established. Without such data for the current ternary systems, we must at these early stages forego a quantitative approach, but are guided instead by the criterion of Cahn,53 an initial developer of the theory, who noted explicitly that “To establish that the composition of the phases change continuously during phase separation would be a direct proof” for a spinodal mechanism to be present. The power of IR spectroscopy for these structurally oriented kinetic measurements is clear from the data presented below. The approach directly tracks temporal rearrangements of Hbonds and formation of separate orthorhombic domains from

Figure 1. Chemical structures of skin lipid constituents used in the current work and H-bonded dimer structure with a trans geometry for SA-d35.

the stearic acid and ceramide constituents. The kinetic sequence of events in the early stages of phase separation in each of the aforementioned mixtures is thereby determined. It is also demonstrated that the addition of a single −OH group in going from [NS] → [AS] ceramide introduces substantial differences in the observed kinetics, a finding of consequence for understanding the possible roles of the different ceramide species. We have previously reported preliminary measurements for similar systems54 that focused primarily on SA domain formation. The current work elaborates upon the analysis and introduces a variety of additional IR spectral parameters that provide mechanistic information about the phase separation events.



EXPERIMENTAL SECTION Cholesterol (Sigma Chemical Co., St. Louis, MO) was of stated purity greater than 99%. Chain perdeuterated stearic acid (SAd35), isotopic purity of 99%, was obtained from CDN Isotopes (Pointe Claire, Quebec, Canada). Porcine brain ceramide [NS] (nonhydroxy fatty acid sphingosine ceramide) was purchased from Avanti Polar Lipids, Inc. (Alabaster, AL) and was >99% pure. Ceramide [AS] (α-hydroxy sphingosine ceramide) from bovine brain was purchased from MATREYA, LLC (Pleasant Gap, PA) and was of stated purity greater than 98%. The chain length compositions (in %) of the major lipid components (>7%) of the two ceramides are as follows: Ceramide [NS] C18:0 (11), C22:0 (10), C24:0 (24), C24:1 (31); Ceramide [AS] C18:0 (24), C22:0 (8), C24:0 (35), C24:1 (17). All substances were used without further purification. Equimolar mixtures of the desired ceramide, cholesterol, and SA-d35 were prepared by codissolving the lipids in a chloroform/methanol mixture (9:1 v:v), removing bulk liquid under a gentle stream of N2, and placing the samples under high vacuum for 3 h. The lipids were subsequently hydrated in excess citrate buffer (pH 5.5) by repeated cycles of heating in a water bath (80 °C), cooling and mixing by vortex action. The sample was then sandwiched between ZnSe windows and 4379

dx.doi.org/10.1021/jp501003c | J. Phys. Chem. B 2014, 118, 4378−4387

The Journal of Physical Chemistry B

Article

placed in a thermostated IR transmission cell (Harrick Scientific, Ossining, NY). The sample was warmed for one additional hour at temperatures greater that the major order− disorder phase transition ∼(65 °C) and quenched within 5 min to 31 °C (most experiments) or 22 °C as desired. When the required temperature was achieved, spectral acquisition with a Thermo-Nicolet 6700 spectrometer (Thermo Electron Corporation, Madison, WI), equipped with a Mercury−Cadmium− Telluride (MCT) detector, was initiated. Data acquisition continued for the desired period with a spectrum generally being acquired every 15 min for the first 12 h and every 15 or 30 min for the duration of the experiment. Absorbance spectra were generated from 128 coadded, Fourier transformed, interferograms collected at ∼2 cm−1 resolution. Occasional experiments were conducted in which spectra were acquired every 7.5 min. The spectrometer and sample compartments were purged with N2 gas for the duration of the experiment to remove residual water vapor. The cer[NS]-containing ternary mixture was examined in nine separate experiments while the ceramide [AS]-containing sample was examined in three separate experiments. Representative data are shown in the figures. Spectra were processed with GRAMS/32, Thermo Fischer Scientific (baseline correction, spectral subtraction, Fourier deconvolution), ISys 3.1, Malvern Instruments (peak intensity), or software developed at the National Research Council of Canada (peak position). Parameters used for particular operations are noted in the figure captions.

Figure 2. (A) Temporal evolution of the IR CO stretching region in a cer[NS]/cholesterol/SA-d35 equimolar mixture quenched from 65 to 31 °C as described in the text. Top: overlaid spectra acquired every 15 min for the first 3 h immediately following the quench highlighting the evolution of the sharp bands in the ceramide Amide I region. Bottom: overlaid spectra acquired every hour ranging from 3 to 20 h after the quench highlighting the delayed downshift in the stearic acid CO frequency. The arrow indicates increasing time within each spectral set. (B) IR intensities as a function of time for the SA-d35 (1694 cm−1) band (red) and for the cer[NS] Amide I (1644 cm−1) band (black).



RESULTS A. Phase Separation Events: Headgroup Reorganization, Chain Ordering, and Orthorhombic Phase Formation in Ceramide [NS]/Cholesterol/Stearic Acid-d35. i. Kinetics of Headgroup H-Bond Rearrangements. Figure 2A depicts the temporal evolution of IR spectra post quenching in the spectral region 1510−1760 cm−1. As indicated in the figure, this spectral region encompasses the Amide I (CO stretch) and Amide II (mixed C−N stretch and N−H in-plane bend, ∼1550 cm−1) vibrations of the cer[NS] amide group, as well as the CO stretching mode of SA-d35. The first 3 h of spectral data acquired at 15 min intervals are presented in the top and those acquired from 3 to 20 hours in 1 h intervals are offset and displayed in the bottom of Figure 2A. Significant changes with different time dependencies were noted in each region following the quench. The Amide I contour appears initially as a relatively broad band (halfwidth >50 cm−1) centered at ∼1640 cm−1. Such features in H-bonded systems are common in condensed phase spectra (e.g., thermally denatured proteins) and reflect contours inhomogeneously broadened by a wide distribution of H-bond lengths and geometries to the Amide CO groups immediately following the thermal quenching. As the system evolves temporally, remarkable changes take place from the initial contour. The broad spectral feature is resolved into four relatively sharp features with progressively increasing intensities assigned to three Amide I components at ∼1644 and 1625 cm−1 along with a shoulder at ∼1611 cm−1 and a weak feature of unknown structural origin at 1594 cm−1. A shift in the Amide II contour (∼1550 cm−1) is also evident. The Amide I contour (Figure 2A), at post quenching times longer than a few hours, thus depicts a well-ordered structure, with a much reduced range of H-bond geometries for each of the three identified Amide I sub-bands. The dominant H-bonded Amide I structure at long times presumably arises from the 1644 cm−1

band, the most intense in this spectral region, assuming the extinction coefficients for each of the aforementioned features to be similar. The CO stretching contour of the SA-d35 initially reveals a main peak near 1710 cm−1 and a shoulder at ∼1734 cm−1 (Figure 2A, top). This contour evolves into a broad single band centered at ∼1694 cm−1 (Figure 2A, bottom). The time evolution of the 1644 and 1694 cm−1 bands, plotted in Figure 2B, demonstrate quite different kinetics. The ordering process for the ceramide Amide CO, which commences close to the quenching time, precedes the stearic acid headgroup structure changes which become significant only at ∼4 h post quenching and diminish in rate after ∼7 h. We note that in repeat experiments a range of delay times (from 4 h, arises from H-bonds in the well-known dimer structure shown in Figure 1. The dimer most likely exhibits the trans geometry also shown in the figure. This structure is common for fatty acid crystals with chain lengths greater than 10 carbons. ii. Kinetics of Ceramide Chain Ordering. In addition to the temporal evolution of the polar headgroup H-bonded structure, insight into the chain behavior of the cer[NS] during the post quench period is revealed through the time-dependence of the asymmetric methylene stretching frequency (υasym(CH2)). The ordering of the chains as deduced from changes in this spectral parameter is shown in Figure 4. The process starts within a few moments of quenching and continues for ∼7−8 h. Several possibilities exist for the origin of the chain ordering. First, the quenching process induces chain ordering from a conformationally disordered state, the remnants of which may be observed, as follows: the CH2 asymmetric stretching frequencies at the earliest times appear at ∼2919.7 cm−1, a frequency which indicates that the chains are already significantly ordered, as the frequency range for this mode typically lies between ∼2916−2917 cm−1 (gel phase) and at ∼2925 cm−1 (liquid crystal, disordered phase). The initial value of ∼2919.7 represents chains that are only slightly disordered. As chain ordering proceeds, the frequency is reduced from

Figure 5. A cartoon depicting the early stages of lamellar phase formation (red, ceramide molecules, two chains; black, stearic acid). Stage 1: Disordered structures during and immediately following quenching. Stage 2: Formation of bilayer structures, some segregation of the SA-d35 with (interlayer) H-bond dimeric structures, and poorly packed cer[NS] with some disordered chain ends. Stage 3: more extensive orthorhombic domains of SA-d35, better packed and chainlength matched ceramide chains (also with some orthorhombic structure). Whether Stage 4 (well separated lamellae) occurs is unknown at present. Cholesterol and water molecules are omitted from this representation. 4381

dx.doi.org/10.1021/jp501003c | J. Phys. Chem. B 2014, 118, 4378−4387

The Journal of Physical Chemistry B

Article

structure, some segregation of the SA-d35 with (interlayer) Hbond dimeric structures (discussed above), and poorly packed cer[NS] with some disordered chain ends. At short times post quench (Stage 2), the longer chains of the cer[NS] (depicted in red) extend beyond both the shorter ceramide chains and the SA-d35 chains (depicted in black), or the cholesterol structure (not included in the cartoon), and are presumably not optimally packed. Thus, the protruding ends of the longer chains initially bend toward the lamellar surface. The subsequent chain interaction (Figure 5, Stage 3) is suggested to result from coalescence between the saturated chains on different molecules and results in a reduced fraction of chains with disordered end regions, possibly coupled with interdigitation. Whether Stage 4 (well-separated lamellae) occurs is unknown at present and will require kinetic diffraction studies for detection. A final comment about the cartoon in Figure 5 is that cholesterol (mostly silent in the current IR analysis) and water molecules have been omitted. It is noted that examination of the water O−H stretching modes in our bulk phase samples compared with samples prepared by an alternative method (“spray technique”) reveal the clear presence of bulk phase water in all our current preparations. Although the precise location of water cannot be ascertained in the current experiments, it seems reasonable to assume that the majority of the water resides in the headgroup (interbilayer) region. Experiments quantifying the amount of water in the samples will be presented at a later date. The observed ceramide chain ordering (Figure 4) is accompanied by formation of orthorhombic structures. For this purpose, it is common (for a review, see ref 58) to monitor packing-induced changes in the scissoring mode contour at ∼1467 cm−1. However, the presence of overlapped cholesterol vibrations in this spectral region renders the contour too complex to unambiguously resolve. As an alternative approach, the CH2 rocking mode contour is shown in Figure 6A. The main CH2 rocking component at ∼721 cm−1 develops a low intensity shoulder at ∼726 cm−1 beginning at times shortly after quenching. In alkanes and phospholipids the rocking mode splittings are indicative of the presence of orthorhombic perpendicular packing. The same interpretation is suggested here for ceramides. The kinetics of the process is tracked in Figure 6B as shoulder intensity at 726 cm−1. An initial fairly early temporal evolution with a half-life of ∼2 h is followed by a slower process that continues for the duration of the experiment (∼20 h). Thus, Stage 3 (Figure 5) depicts better packed and chain-length matched ceramide chains (some with orthorhombic structure) as well as the formation of extensive large orthorhombic domains of SA-d35 as discussed below. iii. Kinetics of Stearic Acid-d35 Phase Separation. Stage 3 in the kinetic sequence is accompanied by temporal evolution of the splitting of the SA-d35 CD2 scissoring mode contour. Typical experimental data are shown in Figure 7A. The origins of the three components of this contour have been thoroughly analyzed by Snyder and co-workers48,49 and applied to lipid systems.59−61 The central peak near 1088 cm−1 which appears as a single band at short times following the quench, diminishes in intensity at long times, broadens, and develops shoulders and finally two new peaks at 1086 and 1092 cm−1. The 1088 cm−1 feature arises from populations of SA-d35 molecules with either of the following properties:

Figure 6. (A) Temporal evolution of the CH2 rocking region for the sample discussed in the caption to Figure 2. The shoulder that arises temporally at 726 cm−1 is assigned to orthorhombically packed cer[NS] phases. Spectra were baseline corrected between 707 and 732 cm−1. (B) The intensity variation of the 726 cm−1 band as a function of time following the quench (red). Also shown is the intensity variation of the 1085 cm−1 peak (black, measured from baseline corrected (1075−1097 cm−1) deconvolved spectra with γ = 2.2 and 10% Fourier smoothing), the low frequency component of the CD2 contour arising from SA-d35 packed in an orthorhombic phase. The lower X-axis scale is logarithmic.

(a) Molecules that possess either disordered chains or conformationally ordered, hexagonally packed chains. In each instance, these do not undergo the vibrational interchain interactions characteristic of orthorhombically packed chains. (b) Molecules with chains that are orthorhombically packed but are physically isolated from other orthorhombic chains with the same isotopic composition. The frequency difference between the (isolated) CH2 and CD2 scissoring modes (∼380 cm−1) precludes strong interchain vibrational coupling between them. Thus, interactions between isotopically different chains (H, D) in the current case do not contribute to the splitting of the scissoring contour. The observed splitting therefore results only from interactions between perdeuterated SA chains in orthorhombic phases. Beginning at 3−4 h following the quench, the scissoring doublet arising from perdeuterated C18 chains packed in an orthorhombic motif appears with progressively increasing intensity at ∼1086 and 1092 cm−1 (Figure 7B). Prior to revealing discrete new spectral features, the halfwidth of the original scissoring contour increases. Snyder et al.46,48 have analyzed these types of spectral changes during lateral phase separation in mixtures of alkanes of unequal chain length in which one component is deuterated. Two different simulation models have been utilized. Initially,48 it was suggested that the observed splitting between the high and low frequency 4382

dx.doi.org/10.1021/jp501003c | J. Phys. Chem. B 2014, 118, 4378−4387

The Journal of Physical Chemistry B

Article

indicating that a similar level of hexagonal or isolated, orthorhombic SA-d35 molecules remain. B. Effects of Quenching Temperature and Ceramide Structure on Phase Separation Kinetics. i. Effect of Temperature Changes. As shown in Figure 8, a reduction in

Figure 7. (A) Temporal evolution of the splitting of the SA-d35 CD2 scissoring mode contour for the sample discussed in the caption to Figure 2. At short times following the quench, a single band arising from ordered, hexagonally packed chains appears as a single feature at 1088 cm−1, while the progressive formation of an orthorhombically packed phase is characterized by a doublet with components at 1086 and 1092 cm−1. (B) Temporal intensity variation of the 1086 cm−1 peak, the low frequency component of the CD2 scissoring contour arising from SA-d35 packed in an orthorhombic phase (blue). Also plotted is the temporal intensity variation of the 1088 cm−1 peak (red), the central component of the CD2 contour arising from SA-d35 packed in a hexagonal phase or in an orthorhombic phase physically isolated from other orthorhombic chains with the same isotopic composition. Intensities are measured from deconvolved spectra (see Figure 6B caption for details).

Figure 8. (A) Effect of quenching temperature on orthorhombic phase formation kinetics of SA-d35, as tracked by the disappearance of the hexagonal marker band at 1088 cm−1 at 22C (red) and at 31C (blue).The lower X-axis scale is logarithmic. Intensities are measured from deconvolved spectra (see Figure 6B caption for details). (B) Effect of quenching temperature on evolution kinetics of H-bond rearrangement of the cer[NS] Amide I peak at 1644 cm−1 at 22C (red) and at 31C (blue). The lower X-axis scale is logarithmic. Intensity measurements were made after the spectra were baseline corrected (1489−1771 cm−1).

components of the contour represents the average number of chains in domains ranging from the smallest domains detectable (∼2−3 chains) to domains of ∼100 molecules in size, which produced the maximal splitting. More recently,46 the band contours have been modeled as linear combinations of deuterated and proteated species in varying mixtures of known concentrations, assuming that the overall IR scissoring contour is described by a domain size distribution function. The splitting (preceded by broadening) observed for SA-d35 in the current experiments increases from the minimum detectable to a maximum of ∼6 cm−1 as noted. This value corresponds closely to that observed52 for pure SA-d35 in an orthorhombic packed state and so suggests, according to Snyder’s first interpretation,48 the presence of large domains with more than 100 chains and a maximum of a few non-SA-d35 chains, since the presence of the latter (i.e., either cholesterol or ceramide) would reduce the splitting from the maximum. Once the maximal splitting is attained following the quench, the 1086 and 1092 cm−1 doublet increases in intensity, without additional frequency shifts. Thus, at times greater than about 4 h, a relatively pure orthorhombic SA-d35 phase has separated laterally, i.e. in an intrabilayer (as deduced from the scissoring band frequencies) and, with the onset of formation of lamellar structures, in an interbilayer manner, as deduced from the frequencies of the SA CO stretching modes. We note that, after ∼30 h, 20% of the central peak intensity is retained

the quenching temperature from 31 to 22 °C significantly slows the kinetics of stearic acid orthorhombic phase formation. The data are presented with a logarithmic time axis for viewing convenience. The rate of orthorhombic phase formation (Figure 8A) is tracked by the disappearance of the hexagonal/disordered marker band at 1088 cm−1. Lowering the quenching temperature to 22 °C slows the ordering kinetics by about an order of magnitude. The cer[NS] headgroup rearrangement from an initially disordered form to an Amide I structure with fixed H-bond geometries and a strong feature at 1644 cm−1 is also slowed to a similar degree (Figure 8B). ii. Effects of Ceramide[NS] → Ceramide[AS] Substitution on Phase Separation Kinetics. The replacement of cer[NS] with cer[AS] significantly changes the kinetics of various steps in the phase separation sequence. Thus, in Figure 9, addition of a single −OH group to the polar region of the ceramide is seen to significantly speed up the changes in the stearic acid Hbonding kinetics. The intensity of the H-bonded SA CO peak at ∼1695 cm−1 in the cer[AS]/SA-d35/cholesterol mixture begins to form immediately upon quenching, and the initial rapid phase of the process is completed after ∼2 h. In contrast, as noted previously, the intensity of the equivalent band in the cer[NS]-containing ternary system depicts a clear initial 4 h lag time followed by a more cooperative transition that spans the time interval ∼4−7 h. 4383

dx.doi.org/10.1021/jp501003c | J. Phys. Chem. B 2014, 118, 4378−4387

The Journal of Physical Chemistry B

Article

2), namely the spectral response to alterations in the Hbonding state of the polar headgroups. Significant differences have been noted in Amide I H-bonding patterns between the ternary systems studied here, compared with IR spectra of pure ceramides.65,66 The early steps in ordered lipid lamellae formation observed herein for the first time involve organization of the cer[NS] polar headgroup leading to the formation of three ordered Hbonds with similar kinetics, each characterized by a different Amide I frequency. Without isotopic labeling of the SA acid function along with molecular modeling, it is not possible to definitively assign a particular Amide I mode to a particular Hbonded structure. The observed conformational ordering and orthorhombic phase formation within the cer[NS] chains occurs on a time scale similar to changes in the ceramide polar headgroup region (compare Figures 2B, 4, and 6B). In contrast to the relatively rapid formation of specific Hbonds involving the ceramide, formation of the (predominantly trans) H-bonded dimer in SA (suggested from the frequency shift in the acid carboxyl stretch) commences after a 3−4 h time lag (Figure 2B). SA chain ordering follows a similar time scale. The SA dimer formation along with chain conformational ordering presumably requires a physical separation between the hydrophobic regions in the two molecules so their chains are effectively in separate monolayers and therefore provides clear evidence of the lamellar structure in this system. This process is accompanied temporally by phase separation of the majority of the SA chains into essentially pure orthorhombic domains. B. A Spinodal Mechanism of Phase Separation? Spinodal decomposition provides a framework to be considered as a possible mechanism for the formation of lamellar phases and domains within lamellae in SC models. Homogeneous solutions, such as the one phase regions of solutions possessing a miscibility gap, may develop fluctuations in chemical composition upon supercooling into the spinodal region of the phase diagram. These fluctuations are small at first but grow until domains of equilibrium composition develop. This differs from the alternate phase separation mechanism of nucleation and growth in which the phase separated material has the equilibrium composition required from the phase diagram. The molecules closest structurally to those studied here for which spinodal decomposition has been suggested are binary alkane mixtures differing in chain length and for which one component has been perdeuterated. Gilbert et al.50 have examined binary paraffin mixtures of C30H(D)62 with C36D(H)74 using time-resolved small-angle neutron scattering, as well as binary mixtures of CnH2n+2/C36D74 for 28 ≤ n ≤ 31. A single mixing process resulting in alternating lamellae was suggested. The phase separation time increased as chain length increased and decreased with increasing chain length mismatch. The results were consistent with interplay between C36 conformational defects and screw motion in the individual alkane chains. In the current skin models the formation of ordered H-bonds in ceramide and SA headgroups define additional types of steps in the lateral reorganization. A limitation in the IR studies of alkanes is the inability of the method, due to the lateral nature of the interchain interactions, to reveal the presence of stacked lamellae of a single species at long times. In the current experiments, headgroup CO frequency shifts provide some structural information with regard to lamellar organization. The current study is complicated by several factors compared to alkane mixtures and is probably best viewed as a

Figure 9. Effect of the addition of a single −OH group to the polar region of the ceramide component on the SA-d35 H-bonding kinetics, as tracked by formation of the 1695 cm−1 band in ceramide[AS] (red), compared to cer[NS] (black). Intensity measurements were made after the spectra were baseline corrected (1489−1771 cm−1).



DISCUSSION A. Sequence of Kinetic Events: A Possible Mechanism for the Onset of Barrier Structure. Lipid organization in human SC and in relevant models for the lipid barrier has been widely studied. In a recent review,10 the authors note the importance of characterizing both the lamellar structure (SPP and LPP phases) as well as the lipid phase behavior parallel to the plane of the lamellae. As is evident from the current experiments, IR spectroscopy is most profitably used to track local component concentrations and structures that characterize the time evolution of microscopic compositional heterogeneity, e.g. domain formation. As noted in the Introduction, structurally oriented studies of the kinetic mechanism of skin barrier formation are sparse. Yet, kinetics plays an essential role in barrier function. Hints as to the formation of metastable phases and the duration of reequilibration kinetics have appeared in the literature. Lafleur62 noted that metastable phases form in an equimolar mixture of bovine brain Type IlI ceramide (aka human type 2, NFA, NS)/ cholesterol/palmitic acid and that the time needed to reform the equilibrium state was “several weeks” at 5 °C and “rapidly” at 40 °C. More recently from this laboratory, Moore et al.52 tracked the kinetics of orthorhombic domain formation in SC lipid models and catalogued kinetic regimes with a power law of the form P = ctα, where P is some (measured) property related to domain size and “c” is a constant. In an application to human SC,63 the partial reformation of the lipid orthorhombic phase was tracked in isolated human SC following heating of the tissue to ∼55 °C and quenching to either 25 or 30 °C. Gorcea et al.64 have most recently used methods similar to those in the current study to qualitatively track lipid domain formation in normal and ceramide deficient SC lipid models. A very small fraction of hexagonally packed (or more disordered) fatty acid was observed in the ceramide deficient models compared to the normal SC model. The biological importance of kinetic behavior has been considered in a recent review by Menon et al.1 in which the dynamic nature of the SC in response to barrier perturbation was demonstrated. These kinetic processes, addressed to date mainly by indirect approaches, appear sufficiently important to justify greater effort toward a more quantitative understanding with more direct physically oriented approaches. IR spectroscopy is particularly effective for the study of chain systems since spectral alterations induced by structural changes are sensitive to both lipid packing and conformational order. An additional element of IR sensitivity is illustrated here (Figure 4384

dx.doi.org/10.1021/jp501003c | J. Phys. Chem. B 2014, 118, 4378−4387

The Journal of Physical Chemistry B

Article

three-component SC model of equimolar cer[NS], SA-d35, and cholesterol. IR provides unique molecular structure information about headgroup H-bonding, lipid packing, and chain order. The following temporal sequence for phase separation/lamellar structure formation was observed: 1) reorganization of ceramide amide H-bonds from disordered forms to ordered structures; 2) appearance of ordered ceramide chains with some orthorhombically packed structures; and 3) phase separation of large orthorhombically packed domains of SAd35 with H-bonded dimeric headgroups. A spinodal decomposition mechanism, defined by continuous composition changes during the phase separation, is suggested as a qualitative description for these events.

demonstration of the feasibility of the approach. The complicating factors are, first, the use of a ternary system with the ceramide component possessing a heterogeneous distribution of alkyl chains. This heterogeneous ceramide species was of course selected to provide a commonly used SC barrier model. Most spinodal experiments on alkanes and polymers involve binary mixtures. The inherent structural complexity of skin models will require extensive experimental work (DSC, IR, etc.) to generate a complete phase diagram from which to define the spinodal regime. Second, the sensitivity of IR to the size of orthorhombic domains that form is limited to 100 chains or so when the splitting of the methylene scissoring contour is used as the spectroscopic probe.48 At that point, the full splitting of the orthorhombic components of the methylene scissoring contour is observed. Thus, detection of the formation of large domains of alternating lamellae lies outside the sensitivity range of the IR technique and must be addressed by scattering measurements. According to Cahn’s stated criteria53 for direct proof of a spinodal mechanism, the compositions of the phases change continuously during phase separation. In several figures6−9 it is evident that the SA phase separates continuously into relatively pure orthorhombically packed domains at ∼3−8 h following the quench. In the early stages of this process, the peak separation between the components of the contour increases to its maximum difference of ∼6 cm−1. Since the splitting is known be sensitive to the number of deuterated chains adjacent to each other, the composition is obviously changing during this stage of the phase separation. These observations are consistent with the Cahn criterion for a spinodal process. A final issue of potential interest to skin biology relevant to the current experiments is the question of why so many ceramide species have evolved and seem to be essential for the correct formation of skin structure. The data in Figure 9 show that addition of a single −OH to the cer[NS] structure has dramatic consequences for the kinetics of barrier formation and likely for the final barrier structure. Several possible causes must be considered. The presence of an additional −OH functional group in Cer[AS] evidently alters the propensity and numbers of H-bonds that form not only to the ceramide headgroup but also to the fatty acid headgroup. This in turn may alter the overall level of SC hydration. Finally, slight differences in chain length distribution (see the Experimental Section) of the acid chain between the two naturally derived ceramides may also influence the kinetics. Further studies are required to delineate the relative importance of each of these factors. However, it seems reasonable even at this early stage to tentatively suggest that each ceramide species will have unique kinetic behavior both for the formation of lamellar structures and for the phase separation processes including the formation of orthorhombic phases. Given that, during the two week SC turnover time between the secretion of lamellar body lipids at the SG-SC border and exfoliation, a proficient barrier is established in healthy skin, the current work begins to suggest models for the complex kinetics and structural rearrangements responsible for competent barrier function.



AUTHOR INFORMATION

Corresponding Author

*Telephone: 973-353-5613. Fax: 973-353-1264. E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS I.S. received student support from Rutgers University. J & J Consumer Companies, Inc. generously provided support for this study. The underpinnings for spectral interpretations were illuminated by discussions with the late R. G. Snyder.



REFERENCES

(1) Menon, G. K.; Cleary, G. W.; Lane, M. E. The Structure and Function of the Stratum Corneum. Int. J. Pharm. 2012, 435, 3−8. (2) McIntosh, T. J.; Stewart, T. E.; Downing, D. T. X-ray Diffraction Analysis of Isolated Skin Lipids: Reconstitution of Intercellular Lipid Domains. Biochemistry 1996, 35, 3649−3653. (3) Mojumdar, E. H.; Groen, D.; Gooris, G. S.; Barlow, D. J.; Lawrence, M. J.; Deme, B.; Bouwstra, J. A. Localization of Cholesterol and Fatty Acids in a Model Lipid Membrane: A Neutron Diffraction Approach. Biophys. J. 2013, 105, 911−918. (4) Bouwstra, J. A.; Honeywell-Nguyen, P. L. Skin Structure and Mode of Action of Vesicles. Adv. Drug Delivery Rev. 2002, 54 (Supp 1), S41−55. (5) Gay, C. L.; Guy, R. H.; Golden, G. M.; Mak, V. H. W.; Francoeur, M. L. Characterization of Low-Temperature (i.e. < 65 °C) Lipid Transitions in Human Stratum Corneum. J. Invest. Dermatol. 1994, 103, 233−239. (6) Groen, D.; Gooris, G. S.; Bouwstra, J. A. New Insights into Stratum Corneum Lipid Organization by X-ray Diffraction Analysis. Biophys. J. 2009, 97, 2242−2249. (7) Hill, J., Jr.; Wertz, P. W. Molecular Models of the Interlamellar Lipid Lamellae from Epidermal Stratum Corneum. Biochim. Biophys. Acta 2003, 1616, 1121−126. (8) Mao, G.; VanWyck, D.; Xiao, X.; Mack Correa, M. C.; Gunn, E.; Flach, C. R.; Mendelsohn, R.; Walters, R. M. Oleic Acid Disorders Stratum Corneum Lipids in Langmuir Monolayers. Langmuir. 2013, 29, 4857−4865. (9) Groen, D.; Gooris, G. S.; Barlow, D. J.; Lawrence, M. J.; van Mechelen, J. B.; Deme, B.; Bouwstra, J. A. Disposition of Ceramide in Model Lipid Membranes Determined by Neutron Diffraction. Biophys. J. 2011, 100, 1481−1489. (10) Bouwstra, J. A.; Gooris, G. S. The Lipid Organization in Human Stratum Corneum and Model Systems. Open Dermatology Journal 2010, 4, 10−13. (11) White, S. H.; Mirejovsky, D.; King, G. I. Structure of Lamellar Lipid Domains and Corneocyte Envelopes of Murine Stratum Corneum. An X-ray Diffraction Study. Biochemistry 1988, 27, 3725− 3732.



CONCLUSIONS The current experiments provide a framework for investigating the sequence of dynamic events involved in the formation and dissipation of the SC barrier in normal and pathological states. Kinetic IR spectroscopy measurements probed the temporal sequence of molecular events producing ordered structures in a 4385

dx.doi.org/10.1021/jp501003c | J. Phys. Chem. B 2014, 118, 4378−4387

The Journal of Physical Chemistry B

Article

(12) Damien, F.; Boncheva, M. The Extent of Orthorhombic Lipid Phases in the Stratum Corneum Determines the Barrier Efficiency of Human Skin In Vivo. J. Invest. Dermatol. 2010, 130, 611−614. (13) Ongpipattanakul, B.; Francoeur, M. L.; Potts, R. O. Polymorphism in Stratum Corneum Lipids. Biochim. Biophys. Acta 1994, 1190, 115−122. (14) Moore, D. J.; Rerek, M. E.; Mendelsohn, R. FTIR Spectroscopy Studies of the Conformational Order and Phase Behavior of Ceramides. J. Phys. Chem. B 1997, 101, 8933−8940. (15) Moore, D. J.; Rerek, M. E.; Mendelsohn, R. Lipid Domains and Orthorhombic Phases in Model Stratum Corneum: Evidence from Fourier Transform Infrared Spectroscopy Studies. Biochem. Biophys. Res. Commun. 1997, 231, 797−801. (16) Janssens, M.; Gooris, G. S.; Bouwstra, J. A. Infrared Spectroscopy Studies of Mixtures Prepared with Synthetic Ceramides Varying in Head Group Architecture: Coexistence of Liquid and Crystalline Phases. Biochim. Biophys. Acta 2009, 1788, 732−742. (17) Krill, S. L.; Knutson, K.; Higuchi, W. I. The Stratum Corneum Lipid Thermotropic Phase Behavior. Biochim. Biophys. Acta 1992, 1112, 281−286. (18) Bouwstra, J. A.; Gooris, G. A.; Weerheim, A.; Kempeneer, J.; Ponec, M. Characterization of Stratum Corneum Structure in Reconstructed Epidermis by X-Ray Diffraction. J. Lipid Res. 1995, 36, 496−504. (19) Chen, Y.-L.; Wiedmann, T. S. Human Stratum Corneum Lipids have a Distorted Orthorhombic Packing at the Surface of Cohesive Failure. J. Invest. Derm. 1996, 107, 15−19. (20) Bouwstra, J. A.; Gooris, G. S.; Dubbelnaar, F. E. R.; Ponec, M. Phase Behavior of Lipid Mixtures Based on Human Ceramides: Coexistence of Crystalline and Liquid Phases. J. Lipid Res. 2001, 42, 1759−1770. (21) Golden, G. M.; Guzek, D. B.; Harris, R. R.; McKie, J. E.; Potts, R. O. Lipid Thermotropic Transitions in Human Stratum Corneum. J. Invest. Derm. 1986, 86, 255−259. (22) Mendelsohn, R.; Flach, C. R.; Moore, D. J. Determination of Molecular Conformation and Permeation in Skin via IR Spectroscopy, Microscopy and Imaging. Biochim. Biophys. Acta 2006, 1758, 923−933. (23) Brown, D. A.; London, E. Structure and Function of Sphingolipid- and Cholesterol-Rich Membrane Rafts. J. Biol. Chem. 2001, 275, 17221−17224. (24) Subczynski, W. K.; Kusumi, A. Dynamics of Raft Molecules in the Cell and Artificial Membranes: Approaches by Pulse EPR Spin Labeling and Single Molecule Optical Microscopy. Biochim. Biophys. Acta 2003, 1610, 231−243. (25) Forslind, B. A Domain Mosaic Model of the Skin Barrier. Acta Derm. Venereol. 2004, 74, 1−6. (26) Oren, A.; Ganz, T.; Liu, L.; Meerloo, T. In Human Epidermis, Beta-Defensin 2 is Packaged in Lamellar Bodies. Exp. Mol. Pathol. 2003, 74, 180−182. (27) Elias, P. M. Epidermal Barrier Function: Intercellular Lamellar Lipid Structures, Origin, Composition and Metabolism. J. Controlled Release 1991, 15, 199−208. (28) Elias, P. M.; Menon, G. K.; Grayson, S.; Brown, B. E. Membrane Structural Alterations in Murine Stratum Corneum: Relationship to the Localization of Polar Lipids and Phospholipases. J. Invest. Derm. 1988, 91, 3−10. (29) Sondell, B.; Thornell, L. E.; Stigbrand, T.; Egelrud, T. Immunolocalization of Stratum Corneum Chymotryptic Enzyme in Human Skin and Oral Epithelium with Monoclonal Antibodies: Evidence of a Proteinase Specifically Expressed in Keratinizing Squamous Epithelia. J. Histochem. Cytochem. 1994, 42, 459−65. (30) Vičanová, J.; Boelsma, E.; Mommaas, A. M.; Kempenaar, J. A.; Forslind, B.; Pallon, J.; Egelrud, T.; Koerte, H. K.; Ponec, M. Normalization of Epidermal Calcium Distribution Profile in Reconstructed Human Epidermis Is Related to Improvement of Terminal Differentiation and Stratum Corneum Barrier Formation. J. Invest. Derm. 1998, 111, 97−106. (31) Wohlrab, J.; Vollmann, A.; Wartewig, S.; Marsch, W. C.; Neubert, R. Noninvasive Characterization of Human Stratum

Corneum of Undiseased Skin of Patients with Atopic Dermatitis and Psoriasis as Studied by Fourier Transform Raman Spectroscopy. Biopolymers 2001, 62, 141−146. (32) Lavrijsen, A. P. M.; Bouwstra, J. A.; Gooris, G. S.; Weerheim, A.; Boddé, H. E.; Maria, P. Reduced Skin Barrier Function Parallels Abnormal Stratum Corneum Lipid Organization in Patients with Lamellar Ichthyosis. J. Invest. Derm. 1995, 105, 619−624. (33) Menon, G.; Ghadially, R. Morphology of Lipid Alterations in the Epidermis: A Review. Microsc. Res. Technol. 1997, 37, 180−192. (34) Ghadially, R.; Reed, J. T.; Elias, P. M. Stratum Corneum Structure and Function Correlates with Phenotype in Psoriasis. J. Invest. Derm. 1996, 107, 558−564. (35) Langer, J. S. Theory of Spinodal Decomposition in Alloys. Annals of Physics 1971, 65, 53−86. (36) Ditchek, B.; Schwartz, L. H. Applications of Spinodal Alloys. Annu. Rev. Mater. Sci. 1979, 9, 219−253. (37) Kubota, K.; Kuwhara, N. Spinodal Decomposition in a Binary Mixture. Phys. Rev. Lett. 1992, 68, 197−200. (38) Snyder, H. L.; Meakin, P. Phase Separation Dynamics: Comparison of Experimental Results. J. Chem. Phys. 1983, 79, 5588−5594. (39) Karim, A.; Slawecki, T. M.; Kumar, S. K.; Douglas, J. F.; Satija, S. K.; Han, C. C.; Russell, T. P.; Liu, Y.; Overney, R.; Sokolov, J.; Rafailovich, M. H. Phase-Separation Induced Surface Patterns in Thin Polymer Blend Films. Macromolecules 1998, 31, 857−862. (40) Bates, F. S. Polymer-Polymer Phase Behavior. Science 1991, 251, 898−905. (41) He, M.; Liu, Y.; Feng, Y.; Jiang, M.; Han, C. C. Spinodal Decomposition in a Hydrogen-Bonded Polymer Blend. Macromolecules 1991, 24, 464−473. (42) Higgins, J. S.; Lipson, J. E. G.; White, R. P. A Simple Approach to Polymer Mixture Miscibility. Philos. Trans. R. Soc. London, Ser. A 2010, 368, 1009−1025. (43) Fan, J.; Han, T.; Haataja, M. Hydrodynamic Effects on Spinodal Decomposition Kinetics in Planar Lipid Bilayer Members. J. Chem. Phys. 2010, 133, 23510−23519. (44) Le, T. D.; Olson, U.; Wennerstrom, H.; Uhrmeister, P.; Rathke, B.; Strey, R. Binodal and Spinodal Curves of an L3 (Sponge) Phase. Phys. Chem. Chem. Phys. 2001, 3, 4346−4354. (45) White, J. W.; Wei, Z. P.; Epperson, J. E.; Snyder, R. G. The Mechanism of Microphase Separation in Paraffin Hydrocarbon Mixtures. Mol. Phys. 1997, 1, 1039−1058. (46) Snyder, R. G.; Clavell-Grunbaum, D.; Strauss, H. L. Spinodal Phase Separation of Unstable Solid-State Binary n-Alkane Mixtures. J. Phys. Chem. B 2007, 111, 13957−13966. (47) Snyder, R. G.; Conti, G.; Strauss, H. L.; Dorset, D. L. Thermally-Induced Mixing in Partially Microphase Segregated Binary n-Alkane Crystals. J. Chem. Phys. 1993, 97, 7342−7350. (48) Snyder, R. G.; Goh, M. C.; Srivatsovoy, V. J. P.; Strauss, H. L.; Dorset, D. L. Measurement of the Growth Kinetics of Microdomains in Binary n-Alkane Solid Solutions by Infrared Spectroscopy. J. Phys. Chem. 1992, 96, 10008−10019. (49) Snyder, R. G.; Goh, M. C.; Srivatsovoy, V. J. P.; Cates, D. A.; Strauss, H. L.; White, J. L.; Dorset, D. L. Hydrogen/Deuterium Isotope Effects on Microphase Separation in Unstable Crystalline Mixtures of Binary n-Alkanes. J. Chem. Phys. 1994, 98, 674−684. (50) Gilbert, E. P.; Reynolds, P. A.; Thiyagarajan, P.; Wozniak, D. G.; White, J. W. Microphase Separation Kinetics in n-Alkane Mixtures. Phys. Chem. Chem. Phys. 1999, 1, 2715−2724. (51) Chen, X.; Kwak, S.; Lafleur, M.; Bloom, M.; Kitson, N.; Thewalt, J. Fatty Acids Influence “Solid” Phase Formation in Models of Stratum Corneum Intercellular Membranes. Langmuir 2007, 23, 5548−5556. (52) Garidel, P.; Folting, B.; Schaller, I.; Kerth, A. The Microstructure of the Stratum Corneum Lipid Barrier: Mid-Infrared Spectroscopic Studies of Hydrated Ceramide: Palmitic Acid: Cholesterol Model Systems. Biophys. Chem. 2010, 150, 144−156. (53) Cahn, J. W. Phase Separation by Spinodal Decomposition in Isotropic Systems. J. Chem. Phys. 1965, 42, 93−99. 4386

dx.doi.org/10.1021/jp501003c | J. Phys. Chem. B 2014, 118, 4378−4387

The Journal of Physical Chemistry B

Article

(54) Moore, D. J.; Snyder, R.; Rerek, M. E.; Mendelsohn, R. Kinetics of Membrane Raft Formation: Fatty Acid Domains in Stratum Corneum Lipid Models. J. Phys. Chem. B 2006, 110, 2378−2386. (55) Umemura, J. Change of Molecular Configuration in Crystalline Decanoic Acid as Studied by Infrared Spectra and Normal Coordinate Analysis. J. Chem. Phys. 1978, 68, 42−48. (56) Hayashi, S.; Umemura, J. Infrared Spectroscopic Evidence for the Coexistence of Two Molecular Configurations in Crystalline Fatty Acids. J. Chem. Phys. 1975, 63, 1732−1740. (57) Maroncelli, M.; Strauss, H. L.; Snyder, R. G. Structure of the nAlkane Binary Solid n-C19H40/n-C21H44 by Infrared Spectroscopy and Calorimetry. J. Phys. Chem. 1985, 89, 5260−5267. (58) Mendelsohn, R.; Moore, D. J. Vibrational Spectroscopic Studies of Lipid Domains in Biomembranes and Model Systems. Chem. Phys. Lipids 1998, 96, 141−157. (59) Moore, D. J.; Rerek, M. E.; Mendelsohn, R. Lipid Domains and Orthorhombic Phases in Model Stratum Corneum: Evidence from Fourier Transform Infrared Spectroscopy Studies. Biochem. Biophys. Res. Commun. 1997, 231, 797−801. (60) Chen, H.-C.; Mendelsohn, R.; Rerek, M. E.; Moore, D. J. Fourier Transform Infrared Spectroscopy and Differential Scanning Calorimetry Studies of Fatty Acid Homogeneous Ceramide 2. Biochim. Biophys. Acta 2000, 1468, 293−303. (61) Shimomura, M.; Song, K.; Rabolt, J. F. Spectroscopic Determination of Diffusion in Langmuir−Blodgett Films. 1. Interlayer Diffusion of Cadmium Arachidate and Its Deuterated Analog Through a Polymer Interface. Langmuir 1992, 8, 887−893. (62) Lafleur, M. Phase Behaviour of Model Stratum Corneum Lipid Mixtures: An Infrared Spectroscopy Investigation. Can. J. Chem. 1998, 76, 1501−1511. (63) Pensack, R. D.; Michniak, B.; Moore, D. J.; Mendelsohn, R. Infrared Kinetic/Structural Studies of Barrier Reformation in Intact Stratum Corneum Following Thermal Perturbation. Appl. Spectrosc. 2006, 60, 1399−1404. (64) Gorcea, M.; Hadgraft, J.; Moore, D. J.; Lane, M. E. Fourier Transform Infrared Spectroscopy Studies of Lipid Domain Formation in Normal and Ceramide Deficient Stratum Corneum Lipid Models. Int. J. Pharm. 2012, 435, 63−68. (65) Moore, D. J.; Rerek, M. E.; Mendelsohn, R. FTIR Spectroscopy Studies of the Conformational Order and Phase Behavior of Ceramides. J. Phys. Chem. B 1997, 101, 8933−8940. (66) Rerek, M. E.; Chen, H.-C.; Markovic, B.; VanWyck, D.; Garidel, P.; Mendelsohn, R.; Moore, D. J. Phytosphingosine and Sphingosine Ceramide Headgroup Hydrogen Bonding: Structural Insights through Thermotropic Hydrogen/Deuterium Exchange. J. Phys. Chem. B 2001, 105, 9355−9362.

4387

dx.doi.org/10.1021/jp501003c | J. Phys. Chem. B 2014, 118, 4378−4387