skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on August 26, 2026

Title: Progress in Lipid Droplet Simulations: How Polarizability Shapes Triacylglycerol Behavior in Bulk and Interfacial Environments
Triacylglycerols (TG) are the primary neutral lipids in lipid droplets (LDs), organelles responsible for lipid storage, metabolism, and signaling. Molecular dynamics (MD) simulations have provided valuable insight into LD structure, but fixed-charge force fields struggle to capture TG behavior across both hydrophobic cores and polar interfaces. Here, we develop and evaluate a polarizable TG model using the Drude2023 lipid force field and benchmark its performance against experimental measurements of bulk density, TG−water interfacial tension, core hydration, and monolayer expansion. The Drude model accurately reproduces the experimental properties and captures key monolayer features such as surface-oriented TGs (SURF-TGs) and chemically distinct membrane packing defects. Compared to fixed-charge models such as C36-standard and C36-cutoff, the Drude polarizable model is the only force field able to capture the dual nature of TG at polar−nonpolar interfaces like the LD monolayer and more homogeneous hydrophobic environments, like the LD core. However, C36-standard is consistent with the Drude results for the LD monolayer, while C36-cutoff is consistent with the decreased hydration in the LD core. Even with large applied surface tensions, C36-cutoff does not produce Drude-like LD monolayer properties. These results highlight the importance of dynamic polarizability and establish Drude2023 as a more reliable framework for simulating TG in heterogeneous systems like LDs.  more » « less
Award ID(s):
2341008
PAR ID:
10650560
Author(s) / Creator(s):
;
Publisher / Repository:
Journal of Chemical Theory and Computation
Date Published:
Journal Name:
Journal of Chemical Theory and Computation
Volume:
21
Issue:
16
ISSN:
1549-9618
Page Range / eLocation ID:
8235 to 8248
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The activation of transcription factor Max-Like Protein x (MLX) is modulated by competition between active dimerization and inactive association with cytosolic lipid droplets (LDs). However, LD association has been shown to depend on the neutral lipid composition. This work explores the mechanism by which MLX specifically targets LDs rich in triacylglycerol (TG) over those with abundant sterol esters (SE). We compare the association ensembles for a potential minimal targeting sequence, an amphipathic helix-loop-helix hairpin, and the full dimerization and cytoplasmic localization domain (DCD), finding the latter requires larger packing defects and quantifiably alters LD membrane properties. Surprisingly, direct interactions with TG neutral lipids are not observed for either sequence. Instead, targeting to SE-rich LDs is blocked for both sequences by insufficient packing defects. We additionally explore the full mechanism of hairpin association, aiming to understand sequence-specific features that enable strong membrane association. We find that there are multiple association pathways, but that each involves a catch, dive, snorkeling, and embedding phase. The combination of multiple catch and dive residues placed on opposing ends of amphipathic helices lengthens the catch phase, greatly enhancing association in a manner that resembles kinetic selection. Once bound, locking interhelical interactions block dissociation. Collectively, our findings suggest that in addition to relative binding affinities, both kinetics and altered surface properties due to protein association could influence competition within the LD proteome. SignificanceThe transcription factor Max-Like Protein x (MLX plays) a central role in metabolic regulation by responding to nutrient status and, simultaneously, neutral lipid composition. This work reveals how MLX selectively targets triacylglycerol-rich lipid droplets (LDs) through sequence-specific interactions with packing defects. We show that LD surface modulates MLX binding and that MLX in turn alters monolayer properties, highlighting a dynamic interplay between protein association and membrane properties. These findings provide new insight into how protein localization and function may be regulated at LD surfaces, with implications for nutrient sensing and, more broadly, transcriptional control relevant to metabolic and disease states. 
    more » « less
  2. Abstract Lipid droplets (LDs) are dynamic organelles that contain an oil core mainly composed of triglycerides (TAG) that is surrounded by a phospholipid monolayer and LD-associated proteins called perilipins (PLINs). During LD biogenesis, perilipin 3 (PLIN3) is recruited to nascent LDs as they emerge from the endoplasmic reticulum. Here, we analyze how lipid composition affects PLIN3 recruitment to membrane bilayers and LDs, and the structural changes that occur upon membrane binding. We find that the TAG precursors phosphatidic acid and diacylglycerol (DAG) recruit PLIN3 to membrane bilayers and define an expanded Perilipin-ADRP-Tip47 (PAT) domain that preferentially binds DAG-enriched membranes. Membrane binding induces a disorder to order transition of alpha helices within the PAT domain and 11-mer repeats, with intramolecular distance measurements consistent with the expanded PAT domain adopting a folded but dynamic structure upon membrane binding. In cells, PLIN3 is recruited to DAG-enriched ER membranes, and this requires both the PAT domain and 11-mer repeats. This provides molecular details of PLIN3 recruitment to nascent LDs and identifies a function of the PAT domain of PLIN3 in DAG binding. 
    more » « less
  3. Abstract Lipid droplets (LDs) are organelles that are necessary for eukaryotic and prokaryotic metabolism and energy storage. They have a unique structure consisting of a spherical phospholipid monolayer encasing neutral lipids such as triacylglycerol (TAG). LDs have garnered increased interest for their implications in disease and for drug delivery applications. Consequently, there is an increased need for tools to study their structure, composition, and dynamics in biological contexts. In this work, we utilize CHARMM‐GUIMembrane Builderto simulate and analyze LDs with and without a plant LD protein, oleosin. The results show thatMembrane Buildercan generate biologically relevant all‐atom LD systems with relatively short equilibration times using a new TAG library having optimized headgroup parameters. TAG molecules originally inserted into a lipid bilayer aggregate in the membrane center, forming a TAG‐only core flanked by two monolayers. The TAG‐only core thickness stably grows with increasing TAG mole fraction. A 70 % TAG system has a core that is thick enough to house oleosin without its interactions with the distal leaflet or disruption of its secondary structure. We hope thatMembrane Buildercan aid in the future study of LD systems, including their structure and dynamics with and without proteins. 
    more » « less
  4. Electrostatic interactions are fundamental to biomolecular structure, stability, and function. While these interactions are traditionally modeled using fixed-charge force fields, such approaches are not transferable among di↵erent molecular environments. Polarizable force fields, such as DRUDE, address this limitation by explicitly incorpo- rating polarization e↵ect. However, their performance does not uniformly surpass that of nonpolarizable force fields, since multiple factors such as bonded terms, dihedral correction maps, and solvent screening also modulate biomolecular dynamics. In this work, we study the Im7 protein to evaluate the structural and dynamic behaviors of non-polarizable (CHARMM36m) and polarizable (DRUDE2019) force fields relative to NMR experiments. Our simulations show that DRUDE better stabilizes ↵-helices than CHARMM36m, including shorter ones that contain helix-breaking residues. However, both force fields underestimate loop dynamics, particularly in the loop I region, mainly due to restricted dihedral angle sampling. Moreover, salt bridge analysis reveals that DRUDE and CHARMM36m preferentially stabilize di↵erent salt bridges, driven by ionic interactions, charge screening by the environment, and neighboring residue flex- ibility Additionally, the latest DRUDE2019 variant, featuring updated NBFIX and NBTHOLE parameters for ion-protein interactions, demonstrated improved accuracy in modeling Na+-protein interactions. These findings are further supported by simu- lations of CBD1, a protein with a -sheet and flexible loops, which exhibited similar trends of stable structured regions and restricted loop dynamics across both force fields. These findings highlight the need to balance bonded and non-bonded interactions along with dihedral correction maps while incorporating polarization e↵ects to improve the accuracy of force fields to model protein structure and dynamics. 
    more » « less
  5. null (Ed.)
    Lipid droplets (LDs) are ubiquitously expressed organelles; the only intracellular organelles that contain a lipid monolayer rather than a bilayer. Proteins localize and bind to this monolayer as they do to intracellular lipid bilayers. The mechanism by which cytosolic LD binding proteins recognize, and bind, to this lipid interface remains poorly understood. Amphipathic α-helix bundles form a common motif that is shared between cytosolic LD binding proteins (e.g., perilipins 2, 3, and 5) and apolipoproteins, such as apoE and apoLp-III, found on lipoprotein particles. Here, we use pendant drop tensiometry to expand our previous work on the C-terminal α-helix bundle of perilipin 3 and the full-length protein. We measure the recruitment and insertion of perilipin 3 at mixed lipid monolayers at an aqueous-phospholipid-oil interface. We find that, compared to its C-terminus alone, the full-length perilipin 3 has a higher affinity for both a neat oil/aqueous interface and a phosphatidylcholine (PC) coated oil/aqueous interface. Both the full-length protein and the C-terminus show significantly more insertion into a fully unsaturated PC monolayer, contrary to our previous results at the air-aqueous interface. Additionally, the C-terminus shows a preference for lipid monolayers containing phosphatidylethanolamine (PE), whereas the full-length protein does not. These results strongly support a model whereby both the N-terminal 11-mer repeat region and C-terminal amphipathic α-helix bundle domains of perilipin 3 have distinct lipid binding, and potentially biological roles. 
    more » « less