Revisiting Lipid Nanoparticle Composition and Structure: A Critical Take on Simulation Approaches

Marius F.W. Trollmann; Paolo Rossetti; Rainer A. B\"ockmann

REVIEW 3 major objections 3 minor 7 references

Revisiting Lipid Nanoparticle Composition and Structure: A Critical Take on Simulation Approaches

T0 review · 3 major / 3 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A simulation ensemble that fixes lipid monolayer area artificially blocks the pH-dependent phase transition in lipid nanoparticles, so the criticized study's composition conclusions are not trustworthy.

desk verdict Worth taking seriously as a comment, but it stops short of proving that the NPzT ensemble is what suppresses the phase transition. read the letter →

arxiv 2412.02731 v1 pith:WUVWUP7S submitted 2024-12-03 cond-mat.soft physics.bio-ph

Marius F.W. Trollmann , Paolo Rossetti , Rainer A. Böckmann This is my paper · ORCID

classification cond-mat.softphysics.bio-ph

keywords lipidnanoparticlespH-dependentphasetransitionmoleculardynamicssimulationNPzTensembleelectrostaticenergydirectcoexistencemRNAdeliveryartifacts

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This comment challenges a recent molecular dynamics study of lipid nanoparticles (LNPs) that reported minimal structural differences between neutral and acidic pH in direct coexistence simulations. The authors argue that the study's NPzT ensemble, which fixes the monolayer area and couples pressure only in the z-direction, artificially prevents the pH-dependent phase transition in which protonated aminolipids form a bilayer. They support this by estimating that the reported charge density of 0.4–0.5 $e_0$/nm$^3$ would give an electrostatic energy of about $10^{11}$–$10^{13}$ kJ/mol for realistic LNP sizes, a physically improbable value. If the critique is right, the earlier conclusions about how LNP composition affects structure and RNA release remain the better supported picture.

What carries the argument

The central object is the NPzT ensemble (constant number of particles, pressure coupling only along the z-axis, fixed lateral area). It restrains the simulated monolayer from changing its in-plane area, which is exactly the degree of freedom that must change when protonated aminolipids drive the system from a hydrophobic-core/monolayer arrangement to a lipid bilayer at low pH. By removing that degree of freedom, the ensemble creates an energetic barrier to the phase transition, which the paper identifies as the source of the criticized study's divergent results.

What would settle it

Run the same direct coexistence system from the criticized study under isotropic pressure coupling (allowing lateral area fluctuations) at low pH. If a lipid bilayer phase does not form, the NPzT ensemble is not the cause and the central claim fails.

Watch

Extended reading notes

Core claim

The paper's central claim is that the NPzT ensemble, which fixes the monolayer area and couples only the z-direction to a barostat, artificially blocks the pH-dependent phase transition in lipid nanoparticle (LNP) direct coexistence simulations. The criticized study reports charge densities of about 0.4–0.5 $e_0$/nm$^3$ at low pH; for 100–200 nm LNPs this implies an electrostatic energy of roughly $10^{11}$ to $10^{13}$ kJ/mol (with $\epsilon_r \approx 10$), which the authors call improbable and equivalent to 20–140 times the energy of TNT. They argue this unphysical charge buildup stems from the ensemble preventing the lateral area contraction that accompanies protonation, so the simulated structures are not an unbiased representation of LNP core–shell architecture.

Load-bearing premise

The critique assumes that the NPzT pressure coupling, rather than differences in system size, lipid composition, or force field parameters, is what prevents the pH-driven phase transition in the criticized study; no control simulation with isotropic coupling is run to verify this.

Editorial extensions

If this is right

If the NPzT ensemble artificially blocks the phase transition, then the criticized study's finding of minimal pH-dependent structural changes is not a reliable description of LNP behavior.
The reported charge densities (~0.4–0.5 $e_0$/nm$^3$) imply electrostatic energies that are physically implausible for 100–200 nm particles, suggesting the simulated structures are not representative of real LNPs.
Differences between LNP formulations, such as the appearance of DSPC-water micellar structures, may be artifacts of the fixed-area constraint rather than composition-dependent biology.
The energetic barrier created by the ensemble may artificially suppress RNA escape under neutral pH, so conclusions about RNA release mechanisms need re-examination.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

A control simulation under isotropic pressure coupling with the identical system would directly test the ensemble-bias hypothesis; the comment does not run such a test, so the causal role of the NPzT ensemble remains asserted rather than demonstrated.
The electrostatic energy estimate assumes a uniform charge distribution and a dielectric constant of about 10; using a different dielectric model could change the magnitude, so the estimate is an order-of-magnitude argument rather than a precise calculation.
Similar area-constrained or semi-isotropic ensembles are common in LNP simulations, so the critique may generalize beyond the single target study to other direct coexistence approaches.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Quick take: this is a comment, not a new simulation study. It raises a legitimate concern about fixed-area pressure coupling (NPzT) in Garaizar et al.'s LNP direct-coexistence simulations, but the causal claim is asserted, not demonstrated.

What is genuinely new here is the specific argument that the NPzT ensemble biases the direct-coexistence setup. The earlier literature, including the authors' own work, sees a pH-dependent phase transition; Garaizar et al. do not. Pointing at the ensemble as the likely cause is a fresh and plausible criticism. The back-of-the-envelope electrostatic energy estimate (10^11–10^13 kJ/mol, equivalent to tens or hundreds of times TNT) is a useful red flag. Even with the crude assumptions, it makes the protonated-charge-in-core picture look physically suspicious. The writing is appropriately cautious: 'may stem from' and 'likely hinders' are honest hedges, not overclaims.

The soft spots are in proportion to the format. There is no control simulation. The authors could have rerun the same system under isotropic pressure coupling to show the phase transition is recovered; without that, system size, composition, force field, or timescale differences remain alternative explanations. The energy estimate depends on an assumed dielectric constant (epsilon_r ~ 10) and a uniform charge distribution; it's an order-of-magnitude argument, not a quantitative proof. The comment also leans heavily on the authors' own prior simulations (refs 1-4) as the expected behavior, which creates a mild circularity: the discrepancy is interpreted through the lens of their own results. That is not disqualifying—the area-constraint mechanism is a plausible independent physical argument—but it does mean the comment is best read as a hypothesis with supporting circumstantial evidence.

For a published comment, this is solid enough. It identifies a concrete methodological concern in a prominent paper and gives a physically reasonable mechanism. A serious referee would probably ask for a direct test, but the argument deserves to be in the literature so the community can see the exchange. I would bring it to our reading group and would be inclined to accept it for peer review, ideally with a request for a matched control simulation or at least a more explicit statement that the causal role of NPzT is a hypothesis.

Referee Report

3 major / 3 minor

Summary. This manuscript is a comment on Garaizar et al. (PNAS 2024), which used direct-coexistence simulations to study lipid nanoparticle (LNP) structure at different pH. The comment argues that Garaizar et al.'s NPzT ensemble, which fixes the lateral monolayer area and couples pressure only in the z-direction, artificially suppresses a pH-dependent transition from a core-monolayer to a bilayer phase. It supports this by noting that Garaizar et al. observe protonated aminolipids uniformly distributed in the core with a charge density of approximately 0.4-0.5 e0/nm3, and estimates the electrostatic energy of such a distribution as 10^11-10^13 kJ/mol, which it calls improbable. The comment concludes that the Garaizar et al. simulations 'may not provide conclusive insights' into how LNP composition affects activity.

Significance. If the criticism is correct, the manuscript identifies a nontrivial ensemble artifact that could invalidate the main conclusions of a high-profile PNAS simulation study and would have implications for the design of future LNP simulations. The manuscript is short, clearly written, and its central claim is falsifiable. Its main strengths are the concrete identification of a plausible mechanism (fixed-area NPzT coupling) and an order-of-magnitude energetic argument that the reported charge distribution is unrealistic. However, the causal link between the NPzT ensemble and the missing phase transition is not demonstrated; the comment relies on comparing with the authors' own earlier simulations that used different system setups, and it does not provide a control simulation or a quantitative isolation of the ensemble effect. The manuscript is therefore a useful critical hypothesis but not yet a fully supported demonstration.

major comments (3)

[Full text, third paragraph] The electrostatic-energy estimate of 10^11-10^13 kJ/mol assumes a uniform charge distribution, a specific dielectric constant (epsilon_r ≈ 10), and no explicit counterion screening. In a direct-coexistence simulation with explicit water and ions, the effective electrostatic penalty could be far lower than the estimate for an isolated charged sphere. Because this estimate is the primary quantitative evidence that the Garaizar et al. configuration is 'improbable,' the authors should either include counterions in the calculation or clearly state the estimate as a zero-order bound, not a definitive proof.
[Full text, fourth paragraph] The central claim that the NPzT ensemble 'artificially block[s]' the pH-dependent transition is not directly tested. The comment does not run a matched control simulation under isotropic pressure coupling (or NPxyT) with the same system, force field, and composition as Garaizar et al. Differences in system size, lipid composition, simulation time, or force-field parameters could equally explain the discrepancy with refs [1-4]. The hedged language ('may stem from,' 'likely hinders') is honest, but the conclusion that Garaizar et al.'s results are inconclusive rests on this untested causal assumption. I recommend adding either a control simulation or an explicit statement that the causal role of NPzT is a hypothesis requiring further work.
[Full text, references [1-4] and paragraphs 1-2] The expected behavior—the pH-dependent phase transition—is taken from the authors' own prior simulations and from other works using different setups (e.g., different box geometries, force fields, and simulation lengths). The manuscript does not establish that those systems are interchangeable with the Garaizar et al. system beyond the common direct-coexistence geometry. This creates a circularity risk: the disagreement with refs [1-4] is interpreted as an error in Garaizar et al. without proving that the comparison systems are controlled. A specific observable, such as the area per lipid or a lipid-order parameter, should be compared quantitatively across the models to justify the comparison.

minor comments (3)

[Full text, third paragraph] The derivation of the energy range (10^11-10^13 kJ/mol) is not shown; a formula or a one-line calculation would help readers reproduce the estimate and understand the scaling with LNP size.
[Full text, fourth paragraph] The phrase 'systemic bias' should be 'systematic bias,' and 'high-perfomance' in the author affiliation should be 'high-performance.' Also, the figure caption contains 'from, together' which should read 'form, together.'
[Figure 1 caption] The notation 'N PzT' appears with an extra space; for consistency with the text 'NPzT' is preferable.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the NPzT critique rests on an independent electrostatic-energy estimate and ensemble mechanics, with only a minor supporting self-citation.

full rationale

This is a comment critiquing another group's simulation protocol, not a derivation with fitted parameters or a prediction. The expected pH-dependent phase transition is introduced as background and supported by refs [1]–[4]; although ref [1] is the authors' own prior work, refs [2]–[4] are independent, so the self-citation is not the sole load-bearing support. The central causal claim—that fixing the monolayer area and coupling only the z-direction (NPzT) suppresses the core-to-bilayer reorganization—is argued from the mechanics of the ensemble and from an order-of-magnitude electrostatic-energy estimate (10^11–10^13 kJ/mol) computed from the reported charge density, not from reusing the target result. The absence of a matched control simulation under isotropic pressure coupling is a genuine limitation in proving causation, but it is a correctness/experimental-design concern, not circularity. No equation or fitted value is renamed as a prediction, and no uniqueness theorem or ansatz is imported via self-citation. Hence no step reduces to its own input by construction.

Assumptions & free parameters 1 free parameters · 2 assumptions · 0 invented entities

The critique rests on two assumptions about simulation ensembles and on a single hand-picked dielectric constant for the energy estimate. No new entities are introduced, and no free parameters are fitted to data beyond the dielectric constant, which is assumed.

free parameters (1)

dielectric constant epsilon_r = 10
Assumed for the electrostatic energy estimate of the LNP core charge distribution. The resulting energy range spans three orders of magnitude, and the dielectric constant is chosen by hand without a clear justification.

assumptions (2)

domain assumption The NPzT ensemble fixes the lateral area of the monolayer, preventing the natural area change associated with protonation.
Central to the critique; asserted in the paragraph describing the pressure coupling scheme. This assumption is not directly verified in the comment.
domain assumption Unbiased simulations (e.g., with isotropic pressure coupling) show a pH-dependent phase transition from a core-monolayer to a bilayer structure.
The comment takes the phase-transition picture from earlier simulations (refs [1-4]) as the correct reference behavior, without presenting new data or addressing potential differences in model parameters.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Revisiting Lipid Nanoparticle Composition and Structure: A Critical Take on Simulation Approaches." pith.science (2026). https://pith.science/paper/WUVWUP7S

@misc{pith2026241202731,
  author       = {Pith},
  title        = {Pith review of: Revisiting Lipid Nanoparticle Composition and Structure: A Critical Take on Simulation Approaches},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WUVWUP7S}},
  note         = {Machine review of arXiv:2412.02731}
}

read the original abstract

Comment to article published in Proc. Natl. Acad. Sci. U. S. A.: Garaizar, A. et al. 'Toward understanding lipid reorganization in RNA lipid nanoparticles in acidic environments.' Proc. Natl. Acad. Sci. U. S. A. 121, e2404555121 (2024)

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

7 extracted references · 7 canonical work pages

[1]

Trollmann, M. F. W. & B¨ ockmann, R. A. mRNA lipid nanoparticle phase transition. Biophys. J. 121, 3927–3939 (2022). 2

work page 2022
[2]

Philipp, J. et al. pH-dependent structural transitions in cationic ionizable lipid mesophases are critical for lipid nanoparticle function. Proc. Natl. Acad. Sci. U. S. A. 120, e2310491120 (2023)

work page 2023
[3]

Atomistic insights into organization of RNA-loaded lipid nanoparticles

Palonc´ yov´ a, M.et al. Atomistic insights into organization of RNA-loaded lipid nanoparticles. J. Phys. Chem. B127, 1158–1166 (2023)

work page 2023
[4]

Ramezanpour, M. et al. Ionizable amino lipid interactions with POPC: implica- tions for lipid nanoparticle function. Nanoscale 11, 14141–14146 (2019)

work page 2019
[5]

Yanez Arteta, M. et al. Successful reprogramming of cellular protein production through mRNA delivered by functionalized lipid nanoparticles. Proc. Natl. Acad. Sci. U. S. A.115, E3351–E3360 (2018)

work page 2018
[6]

Kjølbye, L. R. et al. Martini 3 building blocks for lipid nanoparticle design. ChemRxiv (2024)

work page 2024
[7]

Garaizar, A. et al. Toward understanding lipid reorganization in RNA lipid nanoparticles in acidic environments. Proc. Natl. Acad. Sci. U. S. A. 121, e2404555121 (2024). 3

work page 2024

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.

[1] [1]

Trollmann, M. F. W. & B¨ ockmann, R. A. mRNA lipid nanoparticle phase transition. Biophys. J. 121, 3927–3939 (2022). 2

work page 2022

[2] [2]

Philipp, J. et al. pH-dependent structural transitions in cationic ionizable lipid mesophases are critical for lipid nanoparticle function. Proc. Natl. Acad. Sci. U. S. A. 120, e2310491120 (2023)

work page 2023

[3] [3]

Atomistic insights into organization of RNA-loaded lipid nanoparticles

Palonc´ yov´ a, M.et al. Atomistic insights into organization of RNA-loaded lipid nanoparticles. J. Phys. Chem. B127, 1158–1166 (2023)

work page 2023

[4] [4]

Ramezanpour, M. et al. Ionizable amino lipid interactions with POPC: implica- tions for lipid nanoparticle function. Nanoscale 11, 14141–14146 (2019)

work page 2019

[5] [5]

Yanez Arteta, M. et al. Successful reprogramming of cellular protein production through mRNA delivered by functionalized lipid nanoparticles. Proc. Natl. Acad. Sci. U. S. A.115, E3351–E3360 (2018)

work page 2018

[6] [6]

Kjølbye, L. R. et al. Martini 3 building blocks for lipid nanoparticle design. ChemRxiv (2024)

work page 2024

[7] [7]

Garaizar, A. et al. Toward understanding lipid reorganization in RNA lipid nanoparticles in acidic environments. Proc. Natl. Acad. Sci. U. S. A. 121, e2404555121 (2024). 3

work page 2024