Mesoscale simulation of biomembranes with FreeDTS

Confinement Induces Morphological and Topological Transitions in Multivesicles

The study of self-assembly in confined spaces has gained significant attention among amphiphilic superstructures and colloidal design. The additional complexity introduced by interactions between contents and their containers, along with the effects of shape and lipid mixing, makes multivesicular bodies an interesting subject of study. Despite its promising applications in biomedicine, such as drug delivery and biomimetic materials, much remains unexplored. Here we investigate the effects of confinement on vesicles with varying lipid tail lengths. We first analyze the morphological changes of single spherical vesicles undergoing dehydration, which leads to a prolate-to-oblate transition. Our findings reveal that reductions in water content induce changes of shape while minimally affecting the surface area needed to maintain the hydration layer of lipid phosphate groups. Additionally, using extensive coarse-grained molecular dynamics simulations, we explore how vesicles confined within other vesicles evolve through topological changes into unexpected structures, mainly influenced by the lipid hydrocarbon lengths. Our results highlight the interplay between confinement, curvature-induced lipid sorting, and lipid-mixing entropy, leading to exquisitely self-assembled superstructures.

Delineating the shape of COat Protein complex-II coated membrane bud

Curvature-generating proteins that direct membrane trafficking assemble on the surface of lipid bilayers to bud transport intermediates, which move protein and lipid cargoes from one cellular compartment to another. However, it remains unclear what controls the overall shape of the membrane bud once curvature induction has begun. In vitro experiments showed that excessive concentrations of the COPII protein Sar1 promoted the formation of membrane tubules from synthetic vesicles, while COPII-coated transport intermediates in cells are generally more spherical or lobed in shape. To understand the origin of these morphological differences, we employ atomistic, coarse-grained (CG), and continuum mesoscopic simulations of membranes in the presence of multiple curvature-generating proteins. We first characterize the membrane-bending ability of amphipathic peptides derived from the amino terminus of Sar1, as a function of interpeptide angle and concentration using an atomistic bicelle simulation protocol. Then, we employ CG simulations to reveal that Sec23 and Sec24 control the relative spacing between Sar1 protomers and form the inner-coat unit through an attachment with Sar1. Finally, using dynamical triangulated surface simulations based on the Helfrich Hamiltonian, we demonstrate that the uniform distribution of spacer molecules among curvature-generating proteins is crucial to the spherical budding of the membrane. Overall, our analyses suggest a new role for Sec23, Sec24, and cargo proteins in COPII-mediated membrane budding process in which they act as spacers to preserve a dispersed arrangement of Sar1 protomers and help determine the overall shape of the membrane bud.

There and back again: bridging meso- and nano-scales to understand lipid vesicle patterning

We describe a complete methodology to bridge the scales between nanoscale molecular dynamics and (micrometer) mesoscale Monte Carlo simulations in lipid membranes and vesicles undergoing phase separation, in which curving molecular species are furthermore embedded. To go from the molecular to the mesoscale, we notably appeal to physical renormalization arguments enabling us to rigorously infer the mesoscale interaction parameters from its molecular counterpart. We also explain how to deal with the physical timescales at stake at the mesoscale. Simulating the as-obtained mesoscale system enables us to equilibrate the long wavelengths of the vesicles of interest, up to the vesicle size. Conversely, we then backmap from the meso- to the nano-scale, which enables us to equilibrate in turn the short wavelengths down to the molecular length-scales. By applying our approach to the specific situation of patterning a vesicle membrane, we show that macroscopic membranes can thus be equilibrated at all length-scales in achievable computational time offering an original strategy to address the fundamental challenge of timescale in simulations of large bio-membrane systems.

Diffusion Analyses along Mean and Gaussian-Curved Membranes with CurD

Curved cellular membranes are both abundant and functionally relevant. While novel tomography approaches reveal the structural details of curved membranes, their dynamics pose an experimental challenge. Curvature especially affects the diffusion of lipids and macromolecules, yet neither experiments nor continuum models distinguish geometric effects from those caused by curvature-induced changes in membrane properties. Molecular simulations could excel here, yet despite community interest toward curved membranes, tools for their analysis are still lacking. Here, we satisfy this demand by introducing CurD, our novel and openly available implementation of the Vertex-oriented Triangle Propagation algorithm to the study of lipid diffusion along membranes with mean and/or Gaussian curvature. This approach, aided by our highly optimized implementation, computes geodetic distances significantly faster than conventional implementations of path-finding algorithms. Our tool, applied to coarse-grained simulations, allows for the first time the analysis of curvature effects on diffusion at size scales relevant to physiological processes such as endocytosis. Our analyses with different membrane geometries reveal that Gaussian curvature plays a surprisingly small role on lipid motion, whereas mean curvature; i.e., the packing of lipid headgroups largely dictates their mobility.