Next Article in Journal
Unravelling Contributions of Astrocytic Connexin 43 to the Functional Activity of Brain Neuron–Glial Networks under Hypoxic State In Vitro
Next Article in Special Issue
Increased Tropism of Extracellular Vesicles Derived from Palmitic Acid-Treated Hepatocytes to Activated Hepatic Stellate Cells
Previous Article in Journal
Interaction Energy between Two Separated Charged Spheres Surrounded Inside and Outside by Electrolyte
Previous Article in Special Issue
New Aspects of Bilayer Lipid Membranes for the Analysis of Ion Channel Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Capturing the Liquid-Crystalline Phase Transformation: Implications for Protein Targeting to Sterol Ester-Rich Lipid Droplets

by
R. Jay Braun
and
Jessica M. J. Swanson
*
Department of Chemistry, University of Utah, Salt Lake City, UT 84112-0850, USA
*
Author to whom correspondence should be addressed.
Membranes 2022, 12(10), 949; https://doi.org/10.3390/membranes12100949
Submission received: 2 September 2022 / Revised: 23 September 2022 / Accepted: 24 September 2022 / Published: 28 September 2022
(This article belongs to the Special Issue Functional Analysis and Applications of Membrane Lipids)

Abstract

:
Lipid droplets are essential organelles that store and traffic neutral lipids. The phospholipid monolayer surrounding their neutral lipid core engages with a highly dynamic proteome that changes according to cellular and metabolic conditions. Recent work has demonstrated that when the abundance of sterol esters increases above a critical concentration, such as under conditions of starvation or high LDL exposure, the lipid droplet core can undergo an amorphous to liquid-crystalline phase transformation. Herein, we study the consequences of this transformation on the physical properties of lipid droplets that are thought to regulate protein association. Using simulations of different sterol-ester concentrations, we have captured the liquid-crystalline phase transformation at the molecular level, highlighting the alignment of sterol esters in alternating orientations to form concentric layers. We demonstrate how ordering in the core permeates into the neutral lipid/phospholipid interface, changing the magnitude and nature of neutral lipid intercalation and inducing ordering in the phospholipid monolayer. Increased phospholipid packing is concomitant with altered surface properties, including smaller area per phospholipid and substantially reduced packing defects. Additionally, the ordering of sterol esters in the core causes less hydration in more ordered regions. We discuss these findings in the context of their expected consequences for preferential protein recruitment to lipid droplets under different metabolic conditions.

Graphical Abstract

1. Introduction

Lipid droplets (LDs) are energy storage organelles that contain, traffic, and buffer fluctuations in the availability of neutral lipids, such as triacylglycerols (TG) and sterol esters (SE). Acting as the hubs of lipid metabolism, they not only store and dispense neutral lipids for energy and membrane building, but also harbor metabolic enzymes that participate in lipolysis and lipogenesis [1,2,3]. Additionally, LDs can interact with transcription factors and proteins involved in the innate immune system, expanding their significance to other facets of the cellular lifecycle [4,5]. Malfunctioning of LD metabolic pathways or alteration of LD proteomes can lead to an onslaught of diseases due to lipid imbalances. This includes, but is not limited to, obesity, non-alcoholic fatty liver disease, lipotoxicity, and cancer [6].
LDs originate in the endoplasmic reticulum (ER) membrane, where primarily TGs accumulate between the bilayer leaflets. At a certain threshold, there are enough neutral lipids to nucleate and form a budding LD that eventually separates from the ER [3]. The newly formed LD with a diameter of 100 nm–100 µm can reattach to the ER through a membrane bridge, diffuse through the cytosol and/or interact with other organelles. The LD lifecycle is largely regulated by the proteins embedded in or associated with the phospholipid (PL) monolayer surrounding the neutral lipid core. This monolayer distinguishes LDs from organelles bound by PL bilayers [2]. LD proteins are grouped into two classes depending on their origin. Class I proteins, recently coined ERTOLD proteins, translocate to LDs from the ER bilayer either during budding or upon reattachment via membrane bridges. They often contain membrane-embedded hairpin-motifs. Class II (CYTOLD) proteins target the LD monolayer from the cytosol typically via amphipathic helices that interact with the LD monolayer [2]. Some of these proteins can bind to either LD monolayers or ER bilayers, as both membranes are thought to share the same PL composition [7], although subtle variations may occur following ER detachment. However, proteins can preferentially target one or the other contingent on the metabolic state of the cell and the lifecycle stage of the LD, such as an expanding LD under growth conditions [2,8].
Identifying how proteins dynamically target LDs to regulate lipid trafficking and metabolism is an emerging area of research. It is somewhat challenging as LDs do not utilize the dedicated protein-targeting machinery used by other organelles, such as those found in peroxisomes, the ER, or mitochondria [9,10]. Various studies have proposed that targeting mechanisms instead rely on the organelle’s unique physical properties [11,12,13,14,15]. Early molecular dynamics (MD) simulations suggested that a substantial amount of interdigitation between the TGs and the PL monolayer led to increased solvent exposure of hydrophobic groups called packing defects [16]. Subsequent all-atom MD simulations spanning 8 µs revealed that certain LD properties require more than 1 µs to reach equilibration [14]. With sufficient sampling, 5–8% of the monolayer became composed of TG molecules that had intercalated into and aligned with the PLs. These surface-oriented TGs (SURF-TGs) not only act as chemically unique membrane components, but they introduce chemically-distinct packing defects with an exposed glycerol group in place of a PL head group. SURF-TGs, along with the substantial level of interdigitation between neutral lipid and PL tails, significantly increase the size and lifetime of observed packing defects in a LD monolayer compared to a bilayer of the same composition.
Class II targeting was linked to physical properties by Prevost and co-workers who conducted MD simulations and in vitro experiments on several amphipathic helices. They found that multiple large hydrophobic groups were required for LD targeting, and that these groups interact with the larger, more persistent packing defects found in LDs monolayers [14]. Following MD work from Kim and co-workers explored growing/budding LDs, which are known to preferentially recruit specific proteins such as CCTα. With an applied surface tension of 15 mN/m per leaflet (representative of a growing LD), the number of SURF-TGs increased to ~18% of the monolayer composition, revealing an interesting interplay between SURF-TG molecules and surface tension. As the monolayer is stretched, more TG molecules intercalate with the PLs, decreasing surface tension. However, as worse surfactants there is a limit to how much SURF-TGs can decrease the surface tension. Thus, for each applied surface tension there is a corresponding average SURF-TG composition. The area per lipid (APL) for the PLs in the LD monolayer also increased by 90.8% due to this surface tension (compared to only 16.4% for a bilayer of the same composition), substantially increasing the number and size of packing defects. It was revealed that the increased packing defects upon expansion, whether in the monolayer or bilayer, enabled preferential targeting of CCTα. Interestingly, the SURF-TG glycerols also provided a preferred protein recruitment pathway via hydrogen bonding with a conserved tryptophan [11].
Class I targeting was also linked to physical properties by Olarte and coworkers who used MD simulations and experimental work to study the ERTOLD protein GPAT4. They found a thermodynamic driving force for targeting by sequence-specific residues that repositioned from unstable regions of the PL bilayer to more stable regions of the LD monolayer, with some forming hydrogen bonds with TG glycerol head groups and/or select water molecules just below the monolayer surface [13]. Although this collection of work provides a general understanding of how LD physical properties enable class I and II protein targeting, these studies have focused on TG-rich LDs. The impact of the neutral lipid core composition remains largely unknown.
Despite carrying numerous types of neutral lipids, LDs are typically categorized as either TG- or SE-dominant, depending on the cell type or the metabolic conditions [7,17,18,19]. In many cases, cells harbor LDs that are TG-rich, as TGs are a primary metabolic substrate [20]. However, some cell types contain SE-dominant LDs; for example, macrophage/foam cells contain a sterol-rich environment due to interaction with low density lipoproteins (LDL). Once esterified, those sterols are stored in the LD core [21]. In adrenal cortex cells, LDs are enriched in SEs which are used as precursors for steroid hormone biosynthesis [22]. Furthermore, recent experimental work demonstrated that certain conditions such as stress or starvation increase the SE ratios in HeLa, Huh7, and yeast cells [17,18,19]. Mahamid and coworkers exposed HeLa cells to stress, through starvation and mitotic arrest, which increased the relative ratio and induced a transition of SEs to a liquid-crystalline smectic phase. A smectic phase is characterized by molecules aligning themselves in an orderly fashion in stacked layers [17,23]. Through microscopy, they observed multiple concentric rings of smectic SEs below the PL monolayer and surrounding an amorphous phase of the remaining TGs in the middle of the LD core. It was suggested that this phase transition and separation of lipids is directly related to the cellular and metabolic conditions and alters LD dynamics. It was further hypothesized that the increase of SEs was likely due to consumption of TG metabolites, and that the concentric smectic rings alter surface properties in a manner that alters occupancy of certain LD-proteins while restricting the accessibility of lipids in the core [17]. A subsequent experimental study demonstrated that acute glucose restriction in yeast cells initiated the same SE liquid-crystalline smectic phase transformation in LDs. Proteomic studies additionally confirmed the LD proteome is altered when the LD is enriched with SEs [19,24]; for example, when the SE smectic phase change is induced, Erg6 and the model class I peptide LiveDrop re-localized to the ER, while TG lipases Tgl4-mNg, Tgl3-mNg, and Tgl1-mNg remained on the LDs. It was suggested that glucose restriction promotes the mobilization of TGs from LDs through lipolysis to provide energy for peroxisomes and the mitochondria, consequently inducing a phase change that remodels the LD proteome [19].
Additionally, in vitro experimental work conducted by Chorlay and Thiam established that LD cores containing different neutral lipid compositions have varying affinities for amphipathic-helix proteins; for example, their model amphipathic helix had a higher affinity for LDs with a TG core, as opposed to a squalene core [12]. Mejhert and coworkers recently revealed various proteins that preferentially target either SE- or TG-rich LDs [4]; for example, the transcription factor Max-like protein X (MLX) seems to only target TG-rich LDs, while Adipose triglyceride lipase (ATGL) can localize on either TG or SE-rich LDs. This is curious, as MLX translocates to the nucleus to control the expression of multiple target genes involved in glucose metabolism. Thus, association of MLX with TG-rich LDs, which blocks its translocation to the nucleus, alters glucose metabolism [4]. This implies that cellular conditions that alter the neutral lipid core composition toward SE-rich LDs can modulate metabolic pathways by altering the LD proteome. This raises the question of how LDs with SE- or TG-rich core compositions can recruit specific proteins, while selectively releasing others.
In this work, we conduct all-atom MD simulations of SE-rich model LDs. Our findings reveal how the neutral lipid core composition can alter the LD surface and core properties. We demonstrate how different SE:TG ratios (50:50 vs. 90:10) impart different physical properties. Our simulations span a timescale of 8 µs and capture, for the first time, the amorphous to smectic liquid crystalline phase transformation in the 90:10 SE:TG LD at physiological temperature. The liquid crystalline SEs display distinct ordered layers in the smectic phase. Simulations of the 50:50 SE to TG ratio system retain an amorphous core, but nonetheless present altered surface properties compared to pure-TG LDs. Core SEs in both systems are less surface-active, intercalating far less into the PL monolayer than TGs. Moreover, TG molecules intercalate less (forming less SURF-TG) with increasing SE concentration. Interdigitation of the neutral lipids with the PL tail region is also altered, decreasing in depth and changing in character. As TG intercalation is replaced with ordered SE intercalation in the 90:10 LD, there is a significant drop in area-per-lipid (APL) and increase in PL tail packing and monolayer order. The reduced APL and number of SURF-TGs results in a monolayer with fewer packing defects exposed to the cytosol. Additionally, the hydration of the LD core decreases compared to that in a TG-rich LD, partially due to increased order. Our findings reveal how an SE-rich core composition affects the dynamic physical properties of LDs. We close with a discussion of the potential implications for how core composition might influence protein targeting and the LD lifecycle.

2. Materials and Methods

2.1. Parameterization of CHYO

We used the CHARMM36 force field for our simulations, due to its history of extensive parametrization and agreement of experimental results for membrane dynamics [25,26]. The neutral lipids triacylglycerol (TG) and cholesterol oleate (CHYO) were used as model lipids for the LD neutral lipid core. CHARMM36 force field parameters were previously created for TG [14]; however, the parameters for CHYO did not exist prior to the time of our simulations, meaning parameterization of the molecule was required. Force fields parameters for both cholesterol and oleic acid had previously been established, so we were faced with the task of developing parameters that would unite the two molecules with an ester linkage (Figure 1a). The initial step was to implement dihedral parameters that are engaged between the two molecule groups as a product of the new ester bond. Quantum mechanical (QM) potential energy surface (PES) scans using the Gaussian 09 software [27] package were performed on these dihedrals using the MP2/cc-pVDZ basis set [28,29,30] as described by previous literature in fitting parameters for CHARMM36 lipid molecules, such as sphingomyelin [31]. Larger basis sets and treatments of electron correlation have previously been demonstrated to have minimal effect on conformational energies [31]. Molecular mechanics (MM) PES scans using the GROMACS version 2019.4 simulation engine [32] were then conducted with the un-parameterized dihedrals, which had their force constants set to 0. These were then subsequently fitted to the QM PES scans using the Monte Carlo Simulated Annealing (MCSA) method, providing force constants for the previously undescribed dihedral parameters [33]. After fitting, the newly parameterized MM scans agreed with the QM scans, except for a slight underestimation of energy on the peaks (Figure 1b). Non-bonded parameters were produced from CHARMM General Forcefield Generator (CGenFF), a tool used in CHARMM force fields to fit partial charges among other parameters [34]. The ester bond and angle types that were needed to connect the cholesterol and oleic acid were already previously defined among the CHARMM36 parameters. To validate the behavior of the newly parameterized cholesteryl oleate, we benchmarked against the only available bulk-phase experimental data: density and phase change properties. At 296 K, CHYO has been experimentally demonstrated to have a density of 960 kg/m3 [35]. At 1 µs, our simulation of bulk isotropic CHYO under NPT conditions reached a value of 940 ± 9 kg/m3. Additionally, pure CHYO, or solutions that contain high ratios of CHYO are also known to undergo a phase transition to a smectic phase at physiological temperatures [36]. Therefore, another objective was to simulate a phase change, which further validate parameterization efforts. The results for phase changes are discussed in-depth above in the Section 3. The topology of CHYO can be found at (https://github.com/jaybraunjr/CHYO-topology).

2.2. System Setup

In this work, two separate LD systems were studied: one with a 90:10 CHYO:TG ratio, and another with a 50:50 CHYO:TG ratio. The ratio of neutral lipids in the LD core evolves throughout the LD lifecycle depending on cell type and metabolic conditions [17,37]. Therefore, we chose ratios to represent equally mixed and extensive SE saturation scenarios. The 90:10 system was also selected because this CHYO:TG ratio was experimentally demonstrated to undergo an isotropic-to-smectic phase transition in a bulk mixture [38]. Like previous studies [14,16], the LD systems were constructed into trilayer models to allow for atomistic simulations to be conducted with computational efficiency. For a trilayer setup, a bilayer is separated by the top and bottom leaflets, and a neutral lipid core is inserted inside. Water is placed on top and bottom which represents the cytosol. Periodic boundary conditions (PBC) allow the monolayer and core to behave as if the systems were continuous. First, for each LD system, a bilayer was created using the CHARMM-GUI membrane builder [39], that contained a heterogeneous population of 88:37:10 amount of the zwitterionic phospholipids: 3-palmitoyl-2-oleoyl-D-glycero-1-phosphatidylcholine (POPC), 2,3-dioleoyl-D-glycero-1- phosphatidylethanolamine (DOPE), and phosphatidylinositol (SAPI), respectively, which is representative of the ER bilayer composition [40]. The next step was to create a neutral lipid core that matches the x–y dimensions of the bilayer, while being thick enough in the z-dimension to capture bulk properties and behaviors. Previous simulations [14] of pure-TG LDs indicated that using a 4 nm thick core led to artefacts such as residual ordering, due to the minimal size in which the TG molecules were allowed to move. Doubling the z dimension to create a system with an 8 nm core has been demonstrated to be sufficient to exhibit bulk properties of the neutral lipids [14].
To create the neutral lipid core, a simulation box with dimensions of approximately 10 × 10 × 4 nm was filled with the neutral lipids CHYO and TG using random placement with the molecular dynamics package Packmol [41]. A total of 242 neutral lipids were placed in the 90:10 ratio system, and 220 in the 50:50 ratio system. The number of lipids in the respective systems was selected to fill the dimensions of 4 nm thick in the z axis, and 10 nm in the x–y dimensions to accommodate for the size of the monolayer. The neutral lipid boxes were then equilibrated for 150 ns at 310 K in NPT and periodic boundary conditions to allow the system to settle. This box was than doubled in the z-dimension to create an 8 nm thick core, and it equilibrated again for 100 ns at 310 K in NPT with PBC applied. Subsequently, the bilayer created from the CHARMM-GUI output was equilibrated for 10 ns. The bilayer leaflets were separated manually to create monolayers in which the equilibrated neutral lipid box was inserted. An extra 1 nm of space in the z-dimensions were allowed between each of the monolayers and neutral lipid box. Neutral lipids that were 2 Å or less from the PLs were removed, and the extra space was reduced with 0.1 ns of NPT simulation. The simulation boxes were solvated with the TIP3 water model [42], and 0.15 M NaCL solution.
The final systems were then simulated for 8 µs using the GROMACS version 2019.4, with a 2 fs timestep. The temperatures were set to 310 K (physiological temperature) using the Nose-Hoover [43,44] thermostat and temperature coupling time constant of 1 ps. The particle mesh Ewald (PME) algorithm [45] was used to calculate the long-range electrostatic interactions with a cutoff of 1.0 nm. Lennard–Jones pair interactions were cutoff at 12 Å with a force-switching function between 8 and 12 Å, and pressure was maintained semi-isotropically using the Parrinello-Rahman barostat [46]. The pressure was set to 1.0 bar, with a compressibility of 4.5 × 10 5 bar−1, and a coupling constant of 5.0 ps. The hydrogen bonds were constrained with the LINCS algorithm [47]. Trajectories were extracted every 100 ps. We analyzed the trajectories using MDAnalysis [48], Visual Molecular Dynamics (VMD) [49], GROMACS tools, and in-house Python scripting.

2.3. Order Parameters

2.3.1. Liquid Crystal Ordering

Liquid crystals display properties of both liquids and solids. Essentially, these are substances that are not as ordered as a solid yet have some degree of alignment. Therefore, determining the order of the molecules in bulk is useful as a method to quantify respective crystalline-like properties. Liquid crystal meso-phases, such as the nematic and smectic phases, are characterized by molecules aligning in a preferred direction, attaining orientational ordering [50,51]. The nematic order parameter, which we will refer to as the orientational order parameter (S), quantifies this degree of ordering and alignment in a material. This is defined as
S = 3 2 cos 2 θ 1 2
where θ is the angle between the long-axis molecular vector of the molecules and the system director. The system director is the common axis in which the molecules collectively tend to point. The brackets indicate that it is an ensemble average of all the CHYO molecules taken over all the timesteps. In an isotropic liquid, the average of the cosine terms is 0, and therefore the order parameter is equal to 0. On the other hand, when S = 1 , it indicates a perfect crystalline system.

2.3.2. Phospholipid Tail Order Parameters

The tail order parameters are calculated with Equation (1) with θ being the angle between the position vector of a carbon atom of an acyl chain bonded to a hydrogen with respect to the system normal (z-axis). A larger value of S indicates increased ordering in the membrane. Again, the brackets indicate an ensemble average [52].

2.4. Classifying Surface Molecules

The neutral lipids were classified into separate categories based on whether they behave as membrane components or lipids within the LD core. This designation was based on the position of certain atoms with respect to the average z-position of the PL tails of the respective monolayer leaflet. The average z-position of the PL tails was calculated by identifying the center-of-mass for the position of all the carbons in the respective tail regions, and averaging their z-values. If all six oxygens of the TG were above the average z-position of the PL tails, then it is considered a SURF-TG. For CHYO, we made this classification (SURF-CHYO) if both ester oxygens were above the average z-position of the PL tails. In some cases, only one to several of the designated oxygens would arise past the average z-position of the PL tails; however, in these cases most the occurrences were short lived (<5 ns), and the neutral lipids would return to the bulk of the core. Neutral lipids that exist within the core are referred to as CORE-TG/CHYO.

2.5. Interdigitation

The degree of interdigitation between the neutral lipids and PLs was calculated in a manner consistent with previous studies [14,16,53]. Interdigitation describes the degree of overlap between neutral lipids and PLs, and it provides us with insight regarding occurrences of membrane-surface properties. This first requires calculating the density profiles of both neutral lipids and PLs with respect to the z-axis, which provides us with the overlap parameter:
ρ o v z = 4 ρ N L z ρ P L z ρ N L z + ρ P L z 2
This overlap parameter ( ρ o v z ) ranges from 0 to 1, with 0 meaning that the density profile is equal to zero with no overlap, and 1 meaning there is full overlap ( ρ N L z = ρ P L z ) . To calculate the degree of interdigitation ( λ o v ) , the overlap parameter is then integrated along the z-axis over the whole box:
λ o v = 0 L ρ o v z d z
where L is the z-dimension of the simulation box, λ o v is then described in units of nanometers, providing us with an area common to the two density profiles. The interdigitation was calculated for CORE-TG/CHYO, SURF-TG/CHYO, and the two groups combined for total interdigitation. While useful in previous work [14,16], λ o v does not effectively differentiate the magnitude of interdigitation for different lipid species in the same system. We therefore also calculated the relative density ( ρ R e l ) for the PL overlap region, i.e., the region in which the PLs and neutral lipids overlap. This allows for the quantification of different neutral lipid species that interdigitate in the monolayer. The equation for relative density is ρ r e l = ρ s p e c i e s ρ t o t a l , where ρ t o t a l is the total NL + PL density in the overlapping region.

2.6. Packing Defect Analysis

We employed a Cartesian-based algorithm to calculate the lipid packing defects, adapted from previous work [16,54,55,56]. For a monolayer leaflet, the lipid atoms whose positions were greater than the threshold z t h r were projected onto a 1 Å spacing two-dimensional (x–y) grid. The threshold z t h r was set to the average z-position of the phosphorous atoms of the leaflet minus 2 nm, which therefore encompasses most atoms of the respective leaflet. If a grid point overlaps with an atom (center of atom and grid point is less than atom’s radius plus half of the diagonal grid, ( 3 / 2   Å ), the z-position of the atom and the atom type are saved in the grid point. If the grid point overlaps with a polar PL head-group, then it is not considered a defect. If the grid point overlaps with PL-acyl groups or the neutral lipids, then it is considered a defect, with reference to whichever atom has the highest z-position. For each of these defect types, the neighboring elementary defects are clustered into one. If the clustered defect contains N elementary defects, it is considered to have a defect size of N  Å 2. The probability of finding a defect with size N was computed and fit to an exponential decay function: P N = c e N / π , where c is the normalization constant and π is the packing defect constant. The packing defect constant represents how slowly the decay function falls off. A higher packing defect constant would generally mean that there is a higher probability of finding a larger defect. The script for this packing defect analysis can be found at (https://github.com/ksy141/SMDAnalysis).

3. Results

3.1. Liquid Crystalline Smectic Phase Transformation

Ordering of LD core: To characterize the influence of the neutral lipid composition, we created LD models with a 50:50 and 90:10 SE:TG ratio. The SE we implemented was cholesteryl oleate (CHYO), one of the more common SEs found within cells [57]. Consistent with previous MD studies, we simulated our 90:10 and 50:50 CHYO:TG systems for 8 µs to ensure adequate equilibration of physical properties [14]. After 3 µs, there was a notable ordering in the 90:10 system that was not present in the 50:50 system. Although the neutral lipids were originally given random orientations in both systems (consistent with the isotropic phase) (Figure S1), those in the 90:10 system aligned to form two distinct layers with clear translational order that is representative of a liquid crystalline smectic phase (Figure 2a-left and Figure S1) (Video S1). We calculated the orientational order parameter ( S ) along with two and three-dimensional radial distribution functions (RDFs) between the center of mass (COM) of the CHYO molecule sterol groups to quantify ordering. An S value of 0 indicates no ordering (isotropic) and a value of 1 indicates perfect ordering (crystalline) (see Section 2). The 90:10 simulation started at a negative value of −0.18, meaning the initial random placement of the CHYO was slightly orthogonal to the system director. Over 3 µs, however, S steadily increased and converged to 0.81 ± 0.02, consistent with an aligned ordering of CHYO molecules along the system director (Figure 2b) [58,59]. This value is consistent with previous experiments and simulations exploring smectic orientational order parameters, which tend to fall in the range of 0.7 0.95 [58,60].
Comparing CHYO RDFs from for the first and fifth microseconds of the simulation also demonstrates ordering as a function of the simulation time (Figure 3). The 3D RDF from the first microsecond has a single peak simply due to co-localization within the LD core (Figure 3a). That from the fifth microsecond has three distinct peaks due to CHYO layering. Similarly, the 2D RDF in the xy plane transitions from no peaks (isotropic ordering) to four distinct peaks (Figure 3b). The smectic-phase CHYO molecules align with their oleate tails alternating orientation in an approximate up/down pattern (Figure S2), as long ago predicted based on crystal structures of sterol esters [36,61]. This is further supported by approximately equal proportions of “up and down” oriented CHYO molecules per smectic layer (Figure S2).
In contrast to the 90:10 results, the 50:50 CHYO:TG system retains an orientational order parameter around zero throughout the simulation (Figure S3). This indicates an amorphous phase with no ordering occurring over the 8 µs simulation. Finally, the spacing between CHYO layers in the 90:10 LD system was calculated based on the density profiles of the CHYO smectic phase. The peaks of the profile correspond to the location of the CHYO oxygens (Figure 3c) due to the overlapping orientation of the tails and side chains (Figure S2). The spacing of the peaks is determined to be 3.451 ± 0.035 nm (Figure 3c), which is consistent with experimental measurements of CHYO in LDs (ranging from 3.4 to 3.6 nm) [17,18,19].

3.2. Physical Properties of SE-Rich LDs

Surface-oriented neutral lipids: SURF-TGs are triacylglycerols that have intercalated into the PL monolayer from the neutral lipid core to behave as chemically unique membrane components (see Section 1 and Figure S4) [11,14]. Structurally, SURF-TG molecules are similar to PLs with their glycerol moieties oriented toward the water of the cytosol and their acyl tails extending towards the neutral lipid core. Kim and Swanson reported SURF-TGs occupied 5–8% of the PL monolayer of a pure-TG LD simulation [14]. Our simulations indicate that after an equilibration period of 4 µs, ~2.5% of the 50:50 system monolayer and less than 1% of the 90:10 system monolayer are composed of SURF-TGs (Figure 4). CHYO does not appear to be as surface active as TG. CHYO molecules can enter the monolayer for short periods of time but are readily pushed back to the core. They are also relatively few in number, with an average of only 1.3 and 0.21 CHYO molecules (equivalent to 0.96% and 0.15%) existing at the monolayer surface for the 50:50 the 90:10 systems, respectively. Thus, the presence of SE in the LD core reduces the number of SURF-TG proportional to the decreased concentration of TG, and does not intercalate itself to form a surface-oriented SE.

3.3. Interdigitation

The overlap of neutral lipids with PL tails can also alter surface properties. To quantify this for SE-rich LDs, we calculated both the overlap distance (λOV), a traditionally used but approximate quantification of degree of interdigitation, and the relative densities ( ρ r e l ), which more accurately reflect the magnitude of interdigitation (see Section 2). TG molecules that are found explicitly within the core will herein be called CORE-TG to differentiate them from SURF-TG. Since so few CHYO molecules intercalate into the PL monolayer, we do not distinguish CORE and SURF CHYO.
Starting with the 50:50 LD, the density profiles (Figure 5a) show that TG contributes more density than CHYO in the middle of the LD (~41% CHYO: 59% TG). This is simply a consequence of TG being a larger molecule (885.5 g/mol versus 651.1 g/mol or 42.4%:57.6% of the neutral lipid mass for CHYO:TG). Moving to the region of PL overlap; however, the density profiles from intercalating CHYO and CORE-TG are almost overlapping (Figure 5a) and the total relative density for CHYO increases (~44%:56% CHYO:CORE-TG). Including the PL density this is ~19% CHYO: 24% TG: 57% PL (Figure S5). This suggests that CHYO molecules intercalate with lipid tails more easily than the CORE-TG molecules, which is perhaps surprising given that TG has three colocalized oleoyl tails to CHYO’s one. The overlap distance, (λOV Table S1), which is the integrated overlap parameter ρ o v z over z (Figure 5a bottom), shows equivalent interdigitation distances for CHYO and CORE-TG. Additionally, we find that the overlap parameter profiles for both CORE-TG and SURF-TG are very similar to what was previously reported for pure-TG LDs [14], suggesting the presence of SEs in an amorphous LD core does not significantly change the nature of TG intercalation into the PL monolayer.
For the 90:10 LD, the overlap distance time-dependence (Figure S6) shows that λOV of both CHYO and CORE-TG decrease in the 90:10 system as the core undergoes the liquid-crystalline phase transition. Thus, ordering of the liquid smectic CHYO layers influences interdigitation. Once ordered, the density in the middle of the LD (Figure 5b), which would be 87%:13% for CHYO:TG based on relative masses alone, is now slightly increased in CHYO (~88%:12%). Consistently, the relative density intercalated into the PL tails is slightly decreased in CHYO (85%:15%), which is equivalently ~36.5% CHYO: 6.5% TG: 57% PL when including PL density (Figure S5). Notably, however, the total magnitude of intercalation is not decreased compared to the 50:50 LD. In both systems neutral lipid intercalation makes up ~42–43% of the density in the overlap region. Thus, the liquid smectic ordering retains neutral lipid intercalation and pushes some CORE-TG into the PL overlap region and out of the ordered smectic rings. However, most of the TG molecules are expected to be sequestered into the amorphous LD core under the smectic rings as resolved in microscopy [17,18]. Interestingly, the overlap parameter also shows that the intercalating CORE-TG in the 90:10 system is shifted slightly down into the LD core and away from the PL-aqueous interface.
Lastly, focusing on SURF-TG molecules, they intercalate closer to the PL headgroups with both density and overlap-parameter profiles (Figure 5) that are consistent with previous findings for pure TG-LDs [14]. The reduced density magnitudes are expected from the decreased frequency in the 50:50 and 90:10 systems described above. Thus, the phase transition to liquid smectic seems to decrease the number of SURF-TGs, but not to change their alignment in the PL monolayer.
Collectively, for the 50:50 LD we find that: (1) CHYO interdigitates slightly more than CORE-TG but with similar depth, and (2) the presence of 50% CHYO decreases the magnitude of both SURF-TG and CORE-TG intercalation, but does not change their depth profiles when compared to pure-TG LDs. In contrast, for the 90:10 LD we can say that ordering of the smectic phase: (1) retains the magnitude of intercalation density, (2) decreases the depth of CHYO intercalation, (3) significantly reduces the number of SURF-TG molecules, and (4) shifts some of the CORE-TG into the PL monolayer (or amorphous core), albeit with a slightly shifted depth profile. The potential biological significance of these changes will be discussed below (see Section 4).

3.4. Phospholipid Properties

To characterize the packing and degree of order in the PL monolayer, we calculated the area per lipid (APL) and tail order parameters of the PLs. APL is calculated by dividing surface area of a membrane by the number of components. For bilayers, the exact relationship between APL and lipid order can be nuanced [62], but APL generally decreases with increased lipid packing density, saturated tails and smaller lipid headgroups. The increased packing density of PLs was also previously demonstrated to negatively correlate with the binding affinity of amphipathic helices to artificial LDs [12], potentially due to altered packing defects (further discussed below). The APL for the 90:10 system converges to a reduced value of 0.619 nm2 after the core gains order (Figure 6), while the 50:50 system converges to a larger value of 0.67 nm2. For reference, the previously reported pure-TG LD had an APL of 0.694 nm2, which is significantly larger than that of the ER bilayer (0.628 nm2) due to the intercalation of SURF-TG [14]. There is a positive correlation between the APL and depth of interdigitation (λOV) in the 90:10 system, as indicated by the Pearson correlation coefficient of r = 0.803 (Figure 6b). As CHYO gains order, its intercalation depth into the PL interface decreases (Figure S6) and the monolayer packs more efficiently, thereby decreasing the APL (Figure 6a). In contrast, for the 50:50 CHYO:TG system, which does not have any ordering within the core, the reduced APL compared to a pure-TG LD is largely due to the reduced amount of SURF-TG in the monolayer (2.5% vs. 5–8%) [14]. It was previously shown that APL in the pure-TG LD is strongly correlated with the number of SURF-TG. With three oleoyl tails, intercalated SURF-TG molecules take up a substantial amount of room. For the 50:50 system, there remains a correlation between APL and the number of SURF-TG, but it is less significant than it was in a pure-TG LD. This is likely due to the interplay between reduced SURF-TG and the altered intercalation described above.
Tail order parameters (SCD) quantify the degree of order in the phospholipid tails. For POPC and DOPE, the 90:10 CHYO:TG system has significantly increased SCD values compared to the 50:50 CHYO:TG system (Figure 7) and the ER bilayer (Figure S7). This increased order in the 90:10 system is expected given the significant decrease in APL combined with an ordered core interface; the PL tails have less area to move in as they pack more tightly together. Interestingly, the 50:50 system, which also has a smaller APL due to fewer SURF-TG molecules, has a slight increase in fluidity in the PL tail region compared with a pure-TG LD. This is likely due to replacement of approximately half of the CORE-TG intercalation (3 oleoyl tails) with CORE-CHYO intercalation (1 oleoyl tail). Collectively, the APL and order parameters indicate that the 90:10 system has a significantly more ordered and tightly packed PL monolayer than a pure-TG LD, while the 50:50 monolayer has a more moderate reduction in the APL due to loss of SURF-TG and slightly less ordered PL tails due to roughly half of the CORE-TG intercalation being replaced by CHYO intercalation.

3.5. Packing Defects

Packing defects play a crucial role in recruiting amphipathic helices (class II proteins) to both bilayer and monolayer membranes [2,11,15,56,63]. Bilayers have PL-acyl defects, which are defects that expose the hydrophobic acyl tails of the PLs to the aqueous cytosol. LD monolayers additionally exhibit defects caused by the neutral lipid core. Therefore, we began by distinguishing packing defects based on their molecular composition; PL-acyl, TG-acyl, and CHYO-acyl are defects in which the respective acyl tails are exposed to the cytosol, while TG-glyc are exposed TG-glycerol groups, and CHYO-side are defects caused by either the sterol rings or sidechain of the cholesterol group.
The 50:50 CHYO:TG system has significantly more packing defects as easily visualized when looking down on the monolayers from a bird’s eye view (Figure 8a,b). Quantifying this difference with the probability of observing a packing defect as a function of defect size (Figure 8c) verifies that the 50:50 CHYO:TG system has larger defects. The dominant defect-type for all systems are the PL-acyl defects. To quantify the different types of defects, we calculated packing defect constants ( π ) (Figure 9), which describe the decay rate of the probability of finding a defect with size N fit to a function P N = c e N / π (see Section 2). A higher packing defect constant indicates that there is a higher probability of finding larger defects. In the 50:50 CHYO:TG system, the TG-glycerol and TG-acyl defects have slightly higher defect constants than CHYO (Figure 9c), which is sensible given the presence of 2.5% SURF-TG and nearly identical intercalation of CHYO and TG. Surprisingly, the defects from TG and CHYO are nearly identical in the 90:10 LD, despite the fact that CHYO is present in 9-fold excess and makes up 85% of the intercalated neutral lipid mass (Figure 9c). This indicates that the few remaining SURF-TG (<1%) and intercalated TG molecules contribute much more significantly to the packing defects than CHYO, potentially due to the broader 3-tail structure of TG compared to streamlined 1-tail structure of CHYO.
For a broader interpretation, we also calculated the overall probability of any type of defect occurring (Figure 9a). The 90:10 system has a defect constant of 13.89 ± 0.21 Å2 while the 50:50 system has a defect constant of 26.44 ± 0.25 Å2. To further compare large packing defects, we then calculated the probability of defects that occur over the size of 50 Å2. The 50:50 CHYO:TG system is 3.66 times more likely to exhibit these large defects than the 90:10 system (Figure 9b). To summarize, the presence of CHYO up to 50% decreases the magnitude of packing defects compared to a pure-TG LD, but not substantially since intercalating CHYO still contribute to packing defects. However, upon transformation to the smectic phase with 90% CHYO, packing defects decrease substantially, with the small relative number of remaining TG molecules contributing more to defects than CHYO. In both cases, PL-acyl defects are the dominant defect type, but even these are reduced with increasing CHYO concentrations.

3.6. Hydration of the LD Core

Our simulations revealed some water molecules penetrate into the neutral lipid core, consistent with previous work on pure-TG LDs [14]. These results suggest that a higher ratio of CHYO is inversely associated to hydration (Figure 10), with the 90:10 LD and 50:50 LD systems converging to water densities of 0.0025 g cm 3 and 0.004 g cm 3   , respectively. This is reduced from the previously reported 0.008 g cm 3 for pure-TG LDs [14]. It is important to note that these magnitudes of hydration are likely overestimated since the standard CHARMM36 parameters used in our model for TG are not ideal for bulk hydration properties and interfacial (water/lipid) surface tension [64,65]. The central limitation in these, and all parameters developed to date, is lack of polarizability, which is essential to capture proper neutral lipid interactions in both nonpolar core and polar interfacial environments; for example, bulk TG is hydrophobic, but at a water interface the glycerol moieties should have larger partial charges in response to proximal water. The standard CHARMM36 parameters are on the polar side and thus slightly over-hydrate bulk TG, which was experimentally measured at 1.8 × 10−3 g/mL [66]. However, these parameters correctly capture the APL for the monolayer, which is approximately 15% larger than that of the ER bilayer [12,14]. This is, in large part, due to the presence of SURF-TGs, which behave as worse surfactants than PLs consistent with the small surface tension for LD monolayers [12,67,68]. Updates to the additive forcefields of TG that correct bulk hydration and interfacial water/lipid surface tension do not correctly represent the APL of a LD monolayer, likely due to increased hydrophobicity that retains TG in the LD core [65,69].
The surface properties and APL are, in part, what make LDs unique organelles, as membrane properties are responsible for recruiting proteins. Since our primary focus is on the LD surface properties, the standard CHARMM36 parameters were selected as most appropriate for this study. Additionally, although the hydration profiles reported herein are slightly overestimated, we anticipate they correctly capture some degree of hydration in the LD core that is enhanced just below the PL monolayer and that decreases with increasing SE/CHYO concentrations. Upon formation of the liquid-crystalline smectic phase, any retained core water would localize with the oxygens, pulling back from the PL-neutral lipid interface. These predictions should be compared to experimental analysis of LD core hydration, which would be a very helpful addition to the field, as well as future studies with polarizable force fields.

4. Discussion

Recent studies on TG-rich LDs have revealed unique physical properties that are at least partially responsible for protein recruitment. However, there is growing recognition that SE-rich LDs are different and that these differences likely have important implications for metabolic cellular processes and health [17,18,19,21,22]. Notably, it has been shown that when LDs are transitioned to SE-rich cores there is a shift in their proteome as they exclude some proteins while retaining association with others [4,19,70]. Our objective is to understand how this differential protein recruitment occurs. Using simulations spanning 8 µs, we have herein calculated the physical properties for LDs containing two ratios of SEs to TGs. Our findings provide clues into the mechanisms that would differentiate protein targeting to LDs dependent on their core composition.
Most notably, we observed the anticipated phase change in a 90:10 CHYO:TG LD with the bulk CHYO gaining order over the span of ~3 µs. The resulting orientational order parameter (S) of 0.81 and lattice spacing between layers of 3.45 nm is indicative of a highly ordered smectic phase [71,72]. This spacing is also consistent with the 3.4–3.6 nm SE smectic spacing reported from cryogenic electron microscopy [17,18,19]. The phase change in our study occurred at physiological temperature, further supporting the recent experimental results suggesting SE-induced phase changes do occur in biological systems [17,18,19]. The liquid-crystalline smectic phase within LD cores is distinctive from an amorphous phase of TG-rich LDs, and many of the unique surface properties discussed next stem from this ordering. It should be noted, however, that different species of SEs undergo phase-transitions at different temperatures and/or ratios, depending on the saturation and length of the fatty acid that is attached through the cholesterol ester linkage. In fact, some SEs do not undergo an isotropic-smectic phase transition at all [36]. The effect different types of SEs impose on LDs could pose an interesting question in future work.
It is first important to note that the PL monolayer surrounding a LD of any composition is very different from one leaflet in a PL bilayer. This is partially due to the significant interdigitation of neutral lipids (NLs) into the PL tails complicating the interface. It is also contributed to by surface active TGs (SURF-TGs), which fully intercalate and align with the PLs and act as unique membrane components. We find notable changes in both intercalation and surface-active neutral lipids with increasing SE concentrations.
First, the abundance of surface-active neutral lipids decreases—moderately for the 50:50 LD and significantly the 90:10 LD (Figure 4). CHYO does not seem to be surface active like TG, so there is no SURF-CHYO equivalent of SURF-TG. This is consistent with the decreased solubility of CHYO in PL membranes compared to TG [73,74,75]. Additionally, the number of SURF-TG are reduced from 5–8% in a pure-TG LD to 2.5% in the 50:50 CHYO:TG LD and less than 1% in the 90:10 LD. This will not only impact proteins that directly interact with SURF-TG molecules (discussed with packing defects below), but may also have implications for growth and shrinkage. SURF-TGs act as dynamic membrane components that can modulate surface tension in expanding/budding LDs. Thus, SE-rich LDs may be more resistant to expansion. The CHYO LDs in our simulations are likely comparable to LDs during stressed or starved conditions [17,19]. Therefore, shrinkage instead of expansion would likely occur, as TGs would be actively consumed as metabolites. The decreased surface tension would, in turn, decrease the number of SURF-TGs intercalated into the LD monolayer, consistent with our observed decrease in SURF-TGs with increasing SE concentration.
Focusing on intercalation, there is a decrease in the depth of intercalation from CORE-TG and CHYO molecules into the PL tails with increasing SE concentration (Figure S6). Intercalation from CORE-TG is reduced from a λOV of 1.83 nm for a pure-TG LD [14] to 1.657 nm for the 50:50 systems, and 1.188 for the 90:10 LD (Table S1). As expected, some of this CORE-TG intercalation is replaced by CHYO intercalation. In fact, the density profiles (Figure 5) show that the single CHYO tails intercalate slightly more readily than the three colocalized TG tails with nearly identical intercalation depth in the 50:50 LD (Figure 5a bottom, Table S1). This replacement of TG with CHYO in the PL tail region is a contributing factor to the shift in monolayer packing (APL), order and defects discussed next. For the 90:10 LD, CHYO makes up slightly less than 90% of the intercalation. This is a result of a small number of CORE-TG molecules being pushed into the PL tail region and out of the ordered smectic rings. Similar to CORE-TG, the phase transformation also decreases the depth of CHYO intercalation (Figure S6, Table S1). However, the total magnitude of intercalation density remains the same as that observed in the 50:50 LD. Thus, the more pronounced changes in physical properties for the 90:10 LD system discussed below are a consequence of the altered depth and shifted composition (mostly now CHYO tails) of neutral lipid intercalation, not due to the exclusion of tail intercalation.
Turning to APL, as the 90:10 LD core gains order, the APL for the PL monolayer also decreases (Figure 6a). A decrease in APL is generally associated with increased PL packing, but the exact causes can be nuanced. In this case, multiple factors contribute. The first is the decrease in SURF-TG which were previously shown to strongly correlate with the APL in a pure-TG LD [14]. The drop in APL from 0.694 nm2 for a pure-TG LD to 0.67 nm2 for the 50:50 LD to 0.619 nm2 for the 90:10 LD is consistent with their relative abundance of SURF-TG (5–8%, 2.5% and <1%, respectively). However, reduced depth of neutral lipid interdigitation (λOV Table S1) also contributes, as supported by the correlation between these two variables (Figure 6b). Lastly, the shift in composition to increased CHYO intercalation combined with reduced conformational flexibility of the smectic phase are expected to be large contributing factors to the decreased APL and attendant increase in PL packing. In the smectic phase, the CHYO molecules are aligned with their tails and side chains oriented in a parallel fashion (Figure S2), which leads to a more uniform interface with the PL tails. This interface of the interdigitating smectic-CHYO molecules induces increased ordering in the PL tails as indicated by the increased tail order parameters (Figure 7).
With increasing order and decreasing APL, one would anticipate a change in packing defects. This is highly relevant since experiments and MD studies have demonstrated that both APL and packing defects are correlated with the binding affinities of class II amphipathic helix proteins, such as CCT α and the amphipathic helix from model protein caveolin 1 [11,12]. Defects enable protein recruitment through exposure of PL-acyl tails and neutral lipids that provide ideal interfaces for class II hydrophobic residues. The small drop in APL for the 50:50 LD, and more significant drop for the 90:10 LD system are consistent with the moderate and significant reduction in packing defects in the 50:50 and 90:10 systems, respectively (Figure 8 and Figure 9). Not only are there less defects, but they are also much smaller, with the probability of finding a defect over 50 Å2 reduced 3-fold in the 90:10 CHYO:TG system when compared to TG-rich systems (Figure 9b).
The defect distributions also change. The 50:50 CHYO:TG and pure-TG LDs actually have a similar magnitude of packing defects (26.44 ± 0.25 and 28.62 ± 0.45 Å2, respectively), with some of the TG defects of course being replaced by CHYO defects in the 50:50 system (Figure 9c) [14]. However, TG molecules still contribute more to the packing defects than CHYO molecules in the 50:50 system, likely due to the remaining 2.5% SURF-TG. In the 90:10 smectic LD, the decrease in PL and TG defects is more significant with CHYO contributing far less than 90% of the neutral lipid defects, indicating again that it is less effective in creating defects. These changes in the types of defects present will likely have implications for which proteins are attracted to the LD monolayer; for example, the preferential interactions between glycerol defects and a tryptophan in CCTα, a protein segment that is recruited during LD expansion, could be minimal for a SE-rich LD [11]. Moreover, the presence of CHYO defects may enable entirely different protein interactions yet to be discovered. Proteins that associate to LDs with high SE ratios likely have an association mechanism that circumvents the need to interact with large membrane packing defects, whether it be specific sequences of charged amino acids that interact instead with PL headgroups, or membrane embedded motifs that burry and anchor into the PLs and core. Understanding these association mechanisms for model proteins and peptides will be a valuable focus for future experimental and computational research.
Hydration of the core has previously been proposed to help stabilize certain proteins, such as the class I ERTOLD model peptide LiveDrop [13]. Specific residues of this membrane-spanning peptide that are emersed in the LD core have charged and polar residues, seemingly unfit for a hydrophobic neutral lipid environment. Coordination with water molecules that are also hydrogen bonded to TG glycerol oxygens could provide the requisite polar interactions to enable such conformations. We find that the degree of hydration for the 50:50 and 90:10 CHYO:TG LDs is decreased compared to that of a pure-TG LD [14]. This makes sense given that TG has a total of six oxygens in its glycerol group, which create a polar assemblage for water interactions. Sterol-esters like CHYO on the other hand, have two oxygens, meaning fewer contact points for water. Additionally, the more densely packed monolayer of the 90:10 CHYO:TG LD may make it more difficult for waters to enter past the membrane, which could introduce kinetic barriers for polar protein residues to drag water in. Once inside the core, water molecules will have less mobility (i.e., stabilizing entropic freedom) due to the ordering and density of the smectic CHYOs. This is apparent in the hydration profile for the 90:10 system where the water is predominantly co-localized with the ordered ester oxygens of CHYO (Figure 11 left). This ordering may also have implications for the stabilization of class I proteins since the waters are pulled further down into the core, making them potentially inaccessible for transmembrane residues close to the PL/core interface.
To that point, transmembrane motifs are thought to coordinate and make hydrogen bonds with the TG-glycerol oxygens as well [12]. With CHYO-oleate tails and sidechains existing as the primary chemical component directly below the PL tails for the 90:10 system, potential polar interactions would be significantly reduced. This may contribute to the destabilization of certain proteins as the smectic SE-hydrocarbon tails and chains fill in below the PLs during the phase transition (Figure S2b) [19].

5. Conclusions

In summary, our results provide insight into the surface, interfacial, and core physical properties of SE-rich LDs, which vary greatly from those of TG-rich LDs. Namely, with a high ratio of SEs, the LD core undergoes a phase change to a liquid-crystalline smectic phase. This, in turn, induces ordering in the PL monolayer, reducing the depth that neutral lipids interdigitate and number of SURF-TG molecules, increasing lipid packing and order parameters, decreasing the APL and packing defects, and structuring fewer water molecules in the LD core. CHYO is fundamentally a different molecule from TG and accordingly influences the LD PL monolayer in different ways. The CHYO tails interdigitate just as, if not more readily as TG tails, but they contribute less to packing defects. At lower concentrations, they induce less order in the PL tail region than pure TG does (compare 50:50 to pure-TG in Figure S7), potentially due to CYO’s single oleate tail. However, at sufficiently high concentration the alignment of CHYO into smectic rings translates ordering into the monolayer, significantly reducing the packing defects that play a pivotal role in class II CYTOLD protein association. Collectively, the changes induced by even 50% SE concentration (e.g., changing the percentage and chemical composition of defects) are expected to be sufficient for influencing protein association. Even subtle shifts in the kon and koff rates of competing proteins may be enough to discriminate association under conditions of competition and crowding. However, the significant changes induced by the amorphous to liquid-crystalline smectic phase transformation under high SE concentrations will have a profound impact on the LD proteome. Understanding the detailed mechanisms of protein association to LDs and how they are impacted by neutral lipid composition will be important topics for future studies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/membranes12100949/s1, Figure S1: (a) 90:10 CHYO:TG system, and (b), the 50:50 CHYO:TG system during the 1st µs; Figure S2: The orientational order parameter (S) for the 50:50 CHYO:TG LD averages 0, indicating no ordering has occurred; Figure S3: Representation of SURF-TGs in the PL monolayer. The SURF-TGs behave as membrane components; Figure S4: (a) Interdigitation as a function of time for CORE-TG, and (b) CORE-CHYO. (c) Interdigitation depth as a function of time for total-TG, and (d) Average interdigitation after initial equilibration of 3 µs to allow for CHYO ordering to occur. (d) Interdigitation total-TG and (e) total-CHYO. (f) Average CORE interdigitation depth after initial equilibration of 3 µs; Figure S5: The interdigitation of the neutral lipids and phospholipids. (a) ρ r e l for the 90:10 system and (b) ρ r e l for the 50:50 system; Figure S6: PL tail order parameters for (a) DOPE and (b) POPC; Figure S7: Calculation of the orientation of CHYO molecules per smectic layer (a), and a screenshot of CHYO in one of the layers displaying the alternating orientation (b); Table S1: The average magnitude of interdigitation of neutral lipids for multiple LD species; Video S1: CHYO ordering in LD core over 8 µs.

Author Contributions

Conceptualization, R.J.B. and J.M.J.S.; methodology, R.J.B.; validation, R.J.B. and J.M.J.S.; formal analysis, R.J.B.; investigation, R.J.B.; writing—original draft preparation, R.J.B.; writing—review and editing, J.M.J.S.; supervision, J.M.J.S.; funding acquisition, J.M.J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Institute of General Medicine of the National Institutes of Health under award number R35GM143117 and computational support from the Extreme Science and Engineering Discovery Environment supported by the National Science Foundation (Grant No. ACI-1548562) under allocation MCB200018 as well as the Center for High Performance Computing at the University of Utah.

Data Availability Statement

All parameters required for reproducing the simulation data are contained herein. Specific input files, trajectories and analysis are available upon request.

Acknowledgments

The authors thank Siyoung Kim, Rich Pastor and Jeffrey Klauda for helpful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, C.; Liu, P. The New Face of the Lipid Droplet: Lipid Droplet Proteins. Proteomics 2019, 19, 1700223. [Google Scholar] [CrossRef] [PubMed]
  2. Olarte, M.J.; Swanson, J.M.J.; Walther, T.C.; Farese, R.V. The CYTOLD and ERTOLD pathways for lipid droplet–protein targeting. Trends Biochem. Sci. 2021, 47, 39–51. [Google Scholar] [CrossRef] [PubMed]
  3. Walther, T.C.; Chung, J.; Farese, R.V. Lipid Droplet Biogenesis. Annu. Rev. Cell Dev. Biol. 2017, 33, 491–510. [Google Scholar] [CrossRef]
  4. Mejhert, N.; Kuruvilla, L.; Gabriel, K.R.; Elliott, S.D.; Guie, M.-A.; Wang, H.; Lai, Z.W.; Lane, E.A.; Christiano, R.; Danial, N.N.; et al. Partitioning of MLX-Family Transcription Factors to Lipid Droplets Regulates Metabolic Gene Expression. Mol. Cell 2020, 77, 1251–1264.e9. [Google Scholar] [CrossRef]
  5. Bosch, M.; Sánchez-Álvarez, M.; Fajardo, A.; Kapetanovic, R.; Steiner, B.; Dutra, F.; Moreira, L.; López, J.A.; Campo, R.; Marí, M.; et al. Mammalian lipid droplets are innate immune hubs integrating cell metabolism and host defense. Science 2020, 370, eaay8085. [Google Scholar] [CrossRef] [PubMed]
  6. Onal, G.; Kutlu, O.; Gozuacik, D.; Emre, S.D. Lipid Droplets in Health and Disease. Lipids Health Dis. 2017, 16, 128. [Google Scholar] [CrossRef]
  7. Bartz, R.; Li, W.-H.; Venables, B.; Zehmer, J.K.; Roth, M.R.; Welti, R.; Anderson, R.G.W.; Liu, P.; Chapman, K.D. Lipidomics reveals that adiposomes store ether lipids and mediate phospholipid traffic. J. Lipid Res. 2007, 48, 837–847. [Google Scholar] [CrossRef]
  8. Meyers, A.; Weiskittel, T.M.; Dalhaimer, P. Lipid Droplets: Formation to Breakdown. Lipids 2017, 52, 465–475. [Google Scholar] [CrossRef]
  9. McGilvray, P.T.; Anghel, S.A.; Sundaram, A.; Zhong, F.; Trnka, M.J.; Fuller, J.R.; Hu, H.; Burlingame, A.L.; Keenan, R.J. An ER translocon for multi-pass membrane protein biogenesis. Elife 2020, 9, e56889. [Google Scholar] [CrossRef]
  10. Kim, P.K.; Hettema, E.H. Multiple pathways for protein transport to peroxisomes. J. Mol. Biol. 2015, 427, 1176–1190. [Google Scholar] [CrossRef] [Green Version]
  11. Kim, S.; Oh, M.I.; Swanson, J.M.J. Stressed Lipid Droplets: How Neutral Lipids Relieve Surface Tension and Membrane Expansion Drives Protein Association. J. Phys. Chem. B 2021, 125, 5572–5586. [Google Scholar] [CrossRef] [PubMed]
  12. Chorlay, A.; Thiam, A.R. Neutral lipids regulate amphipathic helix affinity for model lipid droplets. J. Cell Biol. 2020, 219, e201907099. [Google Scholar] [CrossRef]
  13. Olarte, M.J.; Kim, S.; Sharp, M.E.; Swanson, J.M.J.; Farese, R.V.; Walther, T.C. Determinants of Endoplasmic Reticulum-to-Lipid Droplet Protein Targeting. Dev. Cell 2020, 54, 471–487.e7. [Google Scholar] [CrossRef] [PubMed]
  14. Kim, S.; Swanson, J.M.J. The Surface and Hydration Properties of Lipid Droplets. Biophys. J. 2020, 119, 1958–1969. [Google Scholar] [CrossRef]
  15. Prévost, C.; Sharp, M.E.; Kory, N.; Lin, Q.; Voth, G.A.; Farese, R.V.; Walther, T.C. Mechanism and Determinants of Amphipathic Helix-Containing Protein Targeting to Lipid Droplets. Dev. Cell 2018, 44, 73–86.e4. [Google Scholar] [CrossRef] [PubMed]
  16. Bacle, A.; Gautier, R.; Jackson, C.L.; Fuchs, P.F.J.; Vanni, S. Interdigitation between Triglycerides and Lipids Modulates Surface Properties of Lipid Droplets. Biophys. J. 2017, 112, 1417–1430. [Google Scholar] [CrossRef]
  17. Mahamid, J.; Tegunov, D.; Maiser, A.; Arnold, J.; Leonhardt, H.; Plitzko, J.M.; Baumeister, W. Liquid-crystalline phase transitions in lipid droplets are related to cellular states and specific organelle association. Proc. Natl. Acad. Sci. USA 2019, 116, 16866–16871. [Google Scholar] [CrossRef]
  18. Czabany, T.; Wagner, A.; Zweytick, D.; Lohner, K.; Leitner, E.; Ingolic, E.; Daum, G. Structural and biochemical properties of lipid particles from the yeast Saccharomyces cerevisiae. J. Biol. Chem. 2008, 283, 17065–17074. [Google Scholar] [CrossRef]
  19. Rogers, S.; Gui, L.; Kovalenko, A.; Reetz, E.; Nicastro, D.; Henne, W.M. Liquid-crystalline lipid phase transitions in lipid droplets selectively remodel the LD proteome. BioRxiv 2021. [Google Scholar] [CrossRef]
  20. Watt, M.J.; Spriet, L.L.; Watt, M.J. Triacylglycerol lipases and metabolic control: Implications for health and disease. Am. J. Physiol.-Endocrinol. Metab. 2010, 299, 162–168. [Google Scholar] [CrossRef] [Green Version]
  21. Ayyappan, J.P.; Paul, A.; Goo, Y.H. Lipid droplet-associated proteins in atherosclerosis (Review). Mol. Med. Rep. 2016, 13, 4527–4534. [Google Scholar] [CrossRef] [PubMed]
  22. Yu, J.; Zhang, L.; Li, Y.; Zhu, X.; Xu, S.; Zhou, X.-M.; Wang, H.; Zhang, H.; Liang, B.; Liu, P. The Adrenal Lipid Droplet is a New Site for Steroid Hormone Metabolism. Proteomics 2018, 18, e1800136. [Google Scholar] [CrossRef] [PubMed]
  23. Saupel, A. Recent Results in the Field of Liquid Crystals. Angew. Chem. Int. Ed. Engl. 1968, 7, 97–112. [Google Scholar] [CrossRef]
  24. Khor, V.K.; Ahrends, R.; Lin, Y.; Shen, W.-J.; Adams, C.M.; Roseman, A.N.; Cortez, Y.; Teruel, M.N.; Azhar, S.; Kraemer, F.B. The proteome of cholesteryl-ester-enriched versus triacylglycerol-enriched lipid droplets. PLoS ONE 2014, 9, e105047. [Google Scholar] [CrossRef] [PubMed]
  25. MacKerell, A.D., Jr.; Bashford, D.; Bellott, M.L.; Dunbrack, R.L., Jr.; Evanseck, J.D.; Field, M.J.; Fischer, S.; Gao, J.; Guo, H.; Ha, S.; et al. All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins. J. Phys. Chem. B 1998, 102, 3586–3616. [Google Scholar] [CrossRef] [PubMed]
  26. Klauda, J.B.; Venable, R.M.; Freites, J.A.; O’Connor, J.W.; Tobias, D.J.; Mondragon-Ramirez, C.; Vorobyov, I.; MacKerell, A.D., Jr.; Pastor, R.W. Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid Types. J. Phys. Chem. B 2010, 114, 7830–7843. [Google Scholar] [CrossRef]
  27. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 09, revision A.02; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  28. Dunning, T.H. Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen. J. Chem. Phys. 1989, 90, 1007–1023. [Google Scholar] [CrossRef]
  29. Head-Gordon, M.; Pople, J.A.; Frisch, M.J. MP2 energy evaluation by direct methods. Chem. Phys. Lett. 1988, 153, 503–506. [Google Scholar] [CrossRef]
  30. Moller, C.; Plesset, M.S. Note on an Approximation Treatment for Many-Electron Systems. Phys. Rev. 1934, 46, 618–622. [Google Scholar] [CrossRef]
  31. Venable, R.M.; Sodt, A.J.; Rogaski, B.; Rui, H.; Hatcher, E.; MacKerell, A.D.; Pastor, R.W.; Klauda, J.B. CHARMM all-atom additive force field for sphingomyelin: Elucidation of hydrogen bonding and of positive curvature. Biophys. J. 2014, 107, 134–145. [Google Scholar] [CrossRef] [Green Version]
  32. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. Gromacs: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1–2, 19–25. [Google Scholar] [CrossRef]
  33. Guvench, O.; MacKerell, A.D. Automated conformational energy fitting for force-field development. J. Mol. Model. 2008, 14, 667–679. [Google Scholar] [CrossRef] [PubMed]
  34. Vanommeslaeghe, K.; Raman, E.P.; MacKerell, A.D. Automation of the CHARMM General Force Field (CGenFF) II: Assignment of Bonded Parameters and Partial Atomic Charges. J. Chem. Inf. Model. 2012, 52, 3155–3168. [Google Scholar] [CrossRef] [PubMed]
  35. Miller, K.W.; Small, D.M. Triolein-cholesteryl oleate-cholesterol-lecithin emulsions: Structural models of triglyceride-rich lipoproteins. Biochemistry 1983, 22, 443–451. [Google Scholar] [CrossRef]
  36. Small, D.M. Handbook of Lipid Research; Plenum Press: New York, NY, USA, 1986; Volume 4, pp. 149–180. [Google Scholar]
  37. Fujimoto, T.; Parton, R.G. Not just fat: The structure and function of the lipid droplet. Cold Spring Harb. Perspect. Biol. 2011, 3, a004838. [Google Scholar] [CrossRef]
  38. Croll, D.H.; Small, D.M.; Hamilton, J.A. Molecular motions and thermotropic phase behavior of cholesteryl esters with triolein. Biochemistry 1985, 24, 7971–7980. [Google Scholar] [CrossRef]
  39. Jo, S.; Kim, T.; Iyer, V.G.; Im, W. CHARMM-GUI: A web-based graphical user interface for CHARMM. J. Comput. Chem. 2008, 29, 1859–1865. [Google Scholar] [CrossRef]
  40. Yang, Y.; Lee, M.; Fairn, G.D. Phospholipid subcellular localization and dynamics. J. Biol. Chem. 2018, 293, 6230–6240. [Google Scholar] [CrossRef]
  41. Martinez, L.; Andrade, R.; Birgin, E.G.; Martínez, J.M. PACKMOL: A package for building initial configurations for molecular dynamics simulations. J. Comput. Chem. 2009, 30, 2157–2164. [Google Scholar] [CrossRef]
  42. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.W.; Klein, M.L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  43. Nosé, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 1984, 81, 511–519. [Google Scholar] [CrossRef]
  44. Hoover, W.G. Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A 1985, 31, 1695–1697. [Google Scholar] [CrossRef] [PubMed]
  45. Essmann, U.; Perera, L.; Berkowitz, M.L.; Darden, T.; Lee, H.; Pedersen, L.G. A smooth particle mesh Ewald method. J. Chem. Phys. 1995, 103, 8577–8593. [Google Scholar] [CrossRef]
  46. Parrinello, M.; Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. J. Appl. Phys. 1981, 52, 7182–7190. [Google Scholar] [CrossRef]
  47. Hess, B. P-LINCS: A parallel linear constraint solver for molecular simulation. J. Chem. Theory Comput. 2008, 4, 116–122. [Google Scholar] [CrossRef] [PubMed]
  48. Gowers, R.J.; Linke, M.; Barnoud, J.; Reddy, T.J.E.; Melo, M.N.; Seyler, S.L.; Domanski, J.; Dotson, D.L.; Buchoux, S.; Kenney, I.M.; et al. MDAnalysis: A Python Package for the Rapid Analysis of Molecular Dynamics Simulations. In Proceedings of the Python in Science Conference, Austin, TX, USA, 11–17 July 2016; pp. 102–109. [Google Scholar] [CrossRef]
  49. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  50. Mcmillan, W.L. Simple Molecular Model for the SmecticAPhase of Liquid Crystals. Phys. Rev. A 1971, 4, 1238–1246. [Google Scholar] [CrossRef]
  51. Singh, G.; Agra-Kooijman, D.M.; Fisch, M.R.; Vengatesan, M.R.; Song, J.-K.; Kumar, S. Orientational Order in the Nematic and Heliconical Nematic Liquid Crystals. arXiv 2015, arXiv:1510.08362. [Google Scholar]
  52. Schindler, H.; Seelig, J. Order parameters in relation to thermodynamic properties of a phospholipid bilayer. Statistical mechanical interpretation. Biochemistry 1975, 14, 2283–2287. [Google Scholar] [CrossRef]
  53. Das, C.; Noro, M.G.; Olmsted, P.D. Simulation studies of stratum corneum lipid mixtures. Biophys. J. 2009, 97, 1941–1951. [Google Scholar] [CrossRef] [Green Version]
  54. Gautier, R.; Bacle, A.; Tiberti, M.L.; Fuchs, P.F.; Vanni, S.; Antonny, B. PackMem: A Versatile Tool to Compute and Visualize Interfacial Packing Defects in Lipid Bilayers. Biophys. J. 2018, 115, 436–444. [Google Scholar] [CrossRef] [PubMed]
  55. Cui, H.; Lyman, E.; Voth, G.A. Mechanism of membrane curvature sensing by amphipathic helix containing proteins. Biophys. J. 2011, 100, 1271–1279. [Google Scholar] [CrossRef]
  56. Vanni, S.; Vamparys, L.; Gautier, R.; Drin, G.; Etchebest, C.; Fuchs, P.F.; Antonny, B. Amphipathic lipid packing sensor motifs: Probing bilayer defects with hydrophobic residues. Biophys. J. 2013, 104, 575–584. [Google Scholar] [CrossRef] [PubMed]
  57. Degirolamo, C.; Shelness, G.S.; Rudel, L.L. LDL cholesteryl oleate as a predictor for atherosclerosis: Evidence from human and animal studies on dietary fat. J. Lipid Res. 2009, 50, S434–S439. [Google Scholar] [CrossRef] [PubMed]
  58. Poll, K.; Sims, M.T. An insight into de Vries behaviour of smectic liquid crystals from atomistic molecular dynamics simulations. J. Mater. Chem. C 2020, 8, 13040–13052. [Google Scholar] [CrossRef]
  59. Das, P.; Basu, S.; Sinha, R.K.; Das, U. Dipole-induced chiral smectic-C phase in a eutectic mixture of cholesterol esters. Chem. Phys. Lett. 2005, 410, 417–422. [Google Scholar] [CrossRef]
  60. Haege, C.; Jagiella, S.; Giesselmann, F. Molecular Electron Density Distribution and X-Ray Diffraction Patterns of Smectic A Liquid Crystals—A Simulation Study. ChemPhysChem 2019, 20, 2466–2472. [Google Scholar] [CrossRef]
  61. Burks, C.; Engelman, D.M. Cholesteryl myristate conformation in liquid crystalline mesophases determined by neutron scattering. Proc. Natl. Acad. Sci. USA 1981, 78, 6863–6867. [Google Scholar] [CrossRef]
  62. Venable, R.M.; Brown, F.L.H.; Pastor, R.W. Mechanical properties of lipid bilayers from molecular dynamics simulation. Chem. Phys. Lipids 2015, 192, 60–74. [Google Scholar] [CrossRef]
  63. Nuscher, B.; Kamp, F.; Mehnert, T.; Odoy, S.; Haass, C.; Kahle, P.J.; Beyer, K. α-synuclein has a high affinity for packing defects in a bilayer membrane: A thermodynamics study. J. Biol. Chem. 2004, 279, 21966–21975. [Google Scholar] [CrossRef] [Green Version]
  64. Kim, S.; Voth, G.A. Physical Characterization of Triolein and Implications for Its Role in Lipid Droplet Biogenesis. J. Phys. Chem. B 2021, 125, 6874–6888. [Google Scholar] [CrossRef] [PubMed]
  65. Campomanes, P.; Prabhu, J.; Zoni, V.; Vanni, S. Recharging your fats: CHARMM36 parameters for neutral lipids triacylglycerol and diacylglycerol. Biophys. Rep. 2021, 1, 100034. [Google Scholar] [CrossRef] [PubMed]
  66. Platford, R.F. The octanol-water partitioning of some hydrophobic and hydrophilic compounds. Chemosphere 1983, 12, 1107–1111. [Google Scholar] [CrossRef]
  67. Miller, K.W.; Small, D.M. Surface-to-Core and Interparticle Equilibrium Distributions of Triglyceride-rich Lipoprotein Lipids. J. Biol. Chem. 1983, 258, 13772–13784. [Google Scholar] [CrossRef]
  68. Chorlay, A.; Forêt, L.; Thiam, A.R. Origin of gradients in lipid density and surface tension between connected lipid droplet and bilayer. Biophys. J. 2021, 120, 5491–5503. [Google Scholar] [CrossRef] [PubMed]
  69. Kim, S.; Swanson, J.M.J.; Voth, G.A. Computational Studies of Lipid Droplets. J. Phys. Chem. B 2022, 126, 2145–2154. [Google Scholar] [CrossRef]
  70. Hsieh, K.; Lee, Y.K.; Londos, C.; Raaka, B.M.; Dalen, K.T.; Kimmel, A.R. Perilipin family members preferentially sequester to either triacylglycerol-specific or cholesteryl-esterspecific intracellular lipid storage droplets. J. Cell Sci. 2012, 125, 4067–4076. [Google Scholar] [CrossRef]
  71. Leadbetter, A.J.; Norris, E.K. Distribution functions in three liquid crystals from X-ray diffraction measurements. Mol. Phys. 1979, 38, 669–686. [Google Scholar] [CrossRef]
  72. Engelman, D.M.; Humlan, G.M. Molecular Organization of the Cholesteryl Ester Droplets in the Fatty Streaks of Human Aorta. J. Clin. Investig. 1976, 58, 997–1007. [Google Scholar] [CrossRef]
  73. Hamilton, J.A.; Miller, K.W.; Small, D.M. Solubilization of triolein and cholesteryl oleate in egg phosphatidylcholine vesicles. J. Biol. Chem. 1983, 258, 12821–12826. [Google Scholar] [CrossRef]
  74. Spooner, P.J.R.; Small, D.M. Effect of Free Cholesterol on Incorporation of Triolein in Phospholipid Bilayers. Biochemistry 1986, 26, 5820–5825. [Google Scholar] [CrossRef] [PubMed]
  75. Chaban, V.V.; Khandelia, H. Distribution of neutral lipids in the lipid droplet core. J. Phys. Chem. B 2014, 118, 11145–11151. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Parameterization of CHYO (a) required fitting of two dihedrals at the ester linkage. The dihedral potential energy scans (b) for QM and fitted MM are in agreement, except the MM tends to underestimate energy at the peaks.
Figure 1. Parameterization of CHYO (a) required fitting of two dihedrals at the ester linkage. The dihedral potential energy scans (b) for QM and fitted MM are in agreement, except the MM tends to underestimate energy at the peaks.
Membranes 12 00949 g001
Figure 2. (a)-left 90:10 and (a)-right 50:50 ratio CHYO:TG LD systems after 3 μs. The neutral lipids are multicolored to allow for identification of order and disorder (b) Nematic order parameter of CHYO in the 90:10 ratio system, converging to a value of 0.81 ± 0.02.
Figure 2. (a)-left 90:10 and (a)-right 50:50 ratio CHYO:TG LD systems after 3 μs. The neutral lipids are multicolored to allow for identification of order and disorder (b) Nematic order parameter of CHYO in the 90:10 ratio system, converging to a value of 0.81 ± 0.02.
Membranes 12 00949 g002
Figure 3. Three and two-dimensional RDFs (a,b) of the 90:10 ratio system. The reference is the CHYO center of mass. (c) The spacing between smectic layers is 3.451 ± 0.035 nm, comparable to experimental results.
Figure 3. Three and two-dimensional RDFs (a,b) of the 90:10 ratio system. The reference is the CHYO center of mass. (c) The spacing between smectic layers is 3.451 ± 0.035 nm, comparable to experimental results.
Membranes 12 00949 g003
Figure 4. (a) SURF-TGs are not as abundant in SE-rich LDs. Around 2.5% of the 50:50 CHYO:TG ratio system monolayer is composed of SURF-TGs, and less than 1% for the 90:10 ratio system. (b) The average number (‘count’) of SURF-TGs throughout the trajectory after equilibration, with standard deviation for error bars.
Figure 4. (a) SURF-TGs are not as abundant in SE-rich LDs. Around 2.5% of the 50:50 CHYO:TG ratio system monolayer is composed of SURF-TGs, and less than 1% for the 90:10 ratio system. (b) The average number (‘count’) of SURF-TGs throughout the trajectory after equilibration, with standard deviation for error bars.
Membranes 12 00949 g004
Figure 5. Density profiles (top) and overlap parameters (bottom) for the 50:50 CHYO:TG system (a), and 90:10 CHYO:TG systems (b).
Figure 5. Density profiles (top) and overlap parameters (bottom) for the 50:50 CHYO:TG system (a), and 90:10 CHYO:TG systems (b).
Membranes 12 00949 g005
Figure 6. (a) The APL for the 90:10 CHYO:TG ratio LD significantly drops as the system’s core undergoes a phase transition. (b) The APL positively correlates with interdigitation depth, λOV.
Figure 6. (a) The APL for the 90:10 CHYO:TG ratio LD significantly drops as the system’s core undergoes a phase transition. (b) The APL positively correlates with interdigitation depth, λOV.
Membranes 12 00949 g006
Figure 7. Tail order parameters for (a) DOPE and (b) POPC. In both cases, the 90:10 CHYO:TG ratio LD system is more ordered, due to increased packing concomitate with decreased APL.
Figure 7. Tail order parameters for (a) DOPE and (b) POPC. In both cases, the 90:10 CHYO:TG ratio LD system is more ordered, due to increased packing concomitate with decreased APL.
Membranes 12 00949 g007
Figure 8. Birds-eye-view of PL monolayer of (a) 90:10 and (b) 50:50 CHYO:TG systems. Light and dark blue are the PL head and acyl groups, while green and yellow are TG glycerol and acyl groups. Red and purple indicate PL-acyl and TG-glycerol defects that are larger than 15 Å2, respectively. (c) The probability of finding defects of a given sizes fit to exponential curves highlighting the increased propensity of large defects in the 50:50 LD.
Figure 8. Birds-eye-view of PL monolayer of (a) 90:10 and (b) 50:50 CHYO:TG systems. Light and dark blue are the PL head and acyl groups, while green and yellow are TG glycerol and acyl groups. Red and purple indicate PL-acyl and TG-glycerol defects that are larger than 15 Å2, respectively. (c) The probability of finding defects of a given sizes fit to exponential curves highlighting the increased propensity of large defects in the 50:50 LD.
Membranes 12 00949 g008
Figure 9. (a) The packing defect constants for the 90:10, 50:50, and a pure TG LDs show that increasing CHYO concentration leads to less defects. (b) The probability of finding a defect over 50 Å 2 shows the same trend. (c) The respective types of defects demonstrate that PL-acyl are the dominant defects for all 3 systems, while 10% TG induce defects almost as much as the 90% CHYO in the 90:10 ratio LD.
Figure 9. (a) The packing defect constants for the 90:10, 50:50, and a pure TG LDs show that increasing CHYO concentration leads to less defects. (b) The probability of finding a defect over 50 Å 2 shows the same trend. (c) The respective types of defects demonstrate that PL-acyl are the dominant defects for all 3 systems, while 10% TG induce defects almost as much as the 90% CHYO in the 90:10 ratio LD.
Membranes 12 00949 g009
Figure 10. The 90:10 ratio system has the least amount of water in the core, while the pure-TG LD has the most.
Figure 10. The 90:10 ratio system has the least amount of water in the core, while the pure-TG LD has the most.
Membranes 12 00949 g010
Figure 11. Hydration profiles of the LD systems. The cartoon above the plots shows water (dark blue) phospholipids (cyan), and neutral lipid oxygens (red). The plots display the water density as a function of the z-axis of the systems. The 90:10 ratio system (left) has ordering of the CHYO oxygens, leading to ordering of the hydration profile. The 50:50 ratio system (right) has a more dispersed hydration profile.
Figure 11. Hydration profiles of the LD systems. The cartoon above the plots shows water (dark blue) phospholipids (cyan), and neutral lipid oxygens (red). The plots display the water density as a function of the z-axis of the systems. The 90:10 ratio system (left) has ordering of the CHYO oxygens, leading to ordering of the hydration profile. The 50:50 ratio system (right) has a more dispersed hydration profile.
Membranes 12 00949 g011
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Braun, R.J.; Swanson, J.M.J. Capturing the Liquid-Crystalline Phase Transformation: Implications for Protein Targeting to Sterol Ester-Rich Lipid Droplets. Membranes 2022, 12, 949. https://doi.org/10.3390/membranes12100949

AMA Style

Braun RJ, Swanson JMJ. Capturing the Liquid-Crystalline Phase Transformation: Implications for Protein Targeting to Sterol Ester-Rich Lipid Droplets. Membranes. 2022; 12(10):949. https://doi.org/10.3390/membranes12100949

Chicago/Turabian Style

Braun, R. Jay, and Jessica M. J. Swanson. 2022. "Capturing the Liquid-Crystalline Phase Transformation: Implications for Protein Targeting to Sterol Ester-Rich Lipid Droplets" Membranes 12, no. 10: 949. https://doi.org/10.3390/membranes12100949

APA Style

Braun, R. J., & Swanson, J. M. J. (2022). Capturing the Liquid-Crystalline Phase Transformation: Implications for Protein Targeting to Sterol Ester-Rich Lipid Droplets. Membranes, 12(10), 949. https://doi.org/10.3390/membranes12100949

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop