Seipin transmembrane segments critically function in triglyceride nucleation and lipid droplet budding from the membrane
Lipid droplets (LDs) are organelles formed in the endoplasmic reticulum (ER) to store triacylglycerol (TG) and sterol esters. The ER protein seipin is key for LD biogenesis. Seipin forms a cage-like structure, with each seipin monomer containing a conserved hydrophobic helix and two transmembrane (TM) segments. How the different parts of seipin function in TG nucleation and LD budding is poorly understood. Here, we utilized molecular dynamics simulations of human seipin, along with cell-based experiments, to study seipin’s functions in protein–lipid interactions, lipid diffusion, and LD maturation. An all-atom simulation indicates that seipin TM segment residues and hydrophobic helices residues located in the phospholipid tail region of the bilayer attract TG. Simulating larger, growing LDs with coarse-grained models, we find that the seipin TM segments form a constricted neck structure to facilitate conversion of a flat oil lens into a budding LD. Using cell experiments and simulations, we also show that conserved, positively charged residues at the end of seipin’s TM segments affect LD maturation. We propose a model in which seipin TM segments critically function in TG nucleation and LD growth.
eLife, 2022, 11, e75808.
Backmapping with Mapping and Isomeric Information
I present a powerful and flexible backmapping tool named Multiscale Simulation Tool (mstool) that converts a coarse-grained (CG) system into all-atom (AA) resolution and only requires AA to CG mapping and isomeric information (cis/trans/chiral). The backmapping procedure includes two simple steps: a) the AA atoms are randomly placed near the corresponding CG beads according to the provided mapping scheme. b) the energy minimization is performed with two modifications in the forcefield. First, the nonbonded interactions are replaced with a cosine function to ensure numerical stability. Second, additional dihedrals are imposed to maintain the molecule’s isomeric properties. To test the tool’s simplicity and robustness, I have demonstrated a backmapping of the four-bead CG lipid model into the corresponding AA resolution (resolution increased by a factor of 34) without using an intermediate resolution. The tool is freely available at github.com/ksy141/mstool.
J. Phys. Chem. B., submitted. (invited article)
Computational Studies of Lipid Droplets
Lipid droplets (LDs) are intracellular organelles whose primary function is energy storage. Known to emerge from the endoplasmic reticulum (ER) bilayer, LDs have a unique structure with a core consisting of neutral lipids, triacylglycerol (TG) or sterol esters (SE), surrounded by a phospholipid (PL) monolayer and decorated by proteins that come and go throughout their complex lifecycle. In this Feature Article, we review recent developments in computational studies of LDs, a rapidly growing area of research. We highlight how molecular dynamics (MD) simulations have provided valuable molecular-level insight into LD targeting and LD biogenesis. Additionally, we review the physical properties of TG from different force fields compared with experimental data. Possible future directions and challenges are discussed.
J. Phys. Chem. B., 2022, 126, 11, 2145-2154. (Feature Article)
Key Factors Governing Initial Stages of Lipid Droplet Formation
Lipid droplets (LDs) are neutral lipid storage organelles surrounded by a phospholipid (PL) monolayer. LD biogenesis from the endoplasmic reticulum is driven by phase separation of neutral lipids, overcoming surface tension and membrane deformation. However, the core biophysics of the initial steps of LD formation remains relatively poorly understood. Here, we use a tunable, phenomenological coarse–grained model to study triacylglycerol (TG) nucleation in a bilayer membrane. We show that PL rigidity has a strong influence on TG lensing and membrane remodeling: when membrane rigidity increases, TG clusters remain more planar with high anisotropy but a minor degree of phase nucleation. This finding is confirmed by advanced sampling simulations that calculate nucleation free energy as a function of the degree of nucleation and anisotropy. We also show that asymmetric tension, controlled by the number of PL molecules on each membrane leaflet, determines the budding direction. A TG lens buds in the direction of the monolayer containing excess PL molecules to allow for better PL coverage of TG, consistent with the reported experiments. Finally, two governing mechanisms of the LD growth, Ostwald ripening and merging, are observed. Taken together, this study characterizes the interplay between two thermodynamic quantities during the initial LD phases, the TG bulk free energy and membrane remodeling energy.
J. Phys. Chem. B., 2022, 126, 2, 453-462.
Key Factors Governing Initial Stages of Lipid Droplet Formation
Lipid droplets (LDs) are neutral lipid-storing organelles surrounded by a phospholipid (PL) monolayer. At present, how LDs are formed in the endoplasmic reticulum (ER) bilayer is poorly understood. In this study, we present a revised all-atom (AA) triolein (TG) model, the main constituent of the LD core, and characterize its properties in a bilayer membrane to demonstrate the implications of its behavior in LD biogenesis. In bilayer simulations, TG resides at the surface, adopting PL-like conformations (denoted in this work as SURF-TG). Free energy sampling simulation results estimate the barrier for TG relocating from the bilayer surface to the bilayer center to be ∼2 kcal/mol in the absence of an oil lens. SURF-TG is able to modulate membrane properties by increasing PL ordering, decreasing bending modulus, and creating local negative curvature. The other neutral lipid, dioleoyl-glycerol (DAG), also reduces the membrane bending modulus and populates negative curvature regions. A phenomenological coarse-grained (CG) model is also developed to observe larger-scale SURF-TG-mediated membrane deformation. CG simulations confirm that TG nucleates between the bilayer leaflets at a critical concentration when SURF-TG is evenly distributed. However, when one monolayer contains more SURF-TG, the membrane bends toward the other leaflet, followed by TG nucleation if a concentration is higher than the critical threshold. The central conclusion of this study is that SURF-TG is a negative curvature inducer, as well as a membrane modulator. To this end, a model is proposed in which the accumulation of SURF-TG in the luminal leaflet bends the ER bilayer toward the cytosolic side, followed by TG nucleation.
J. Phys. Chem. B., 2021, 125, 25, 6874-6888.