Membrane stress profiles from self-consistent field theory

The diffuse interface description of fluid lipid membranes captures key features of the hemifusion pathway and lateral stress profile

Topological transitions of lipid membranes are ubiquitous in key biological processes for cell life, like neurotransmission, fertilization, morphogenesis, and viral infections. Despite this, they are not well understood due to their multiscale nature, which limits the use of molecular models and calls for a mesoscopic approach such as the celebrated Canham-Helfrich one. Unfortunately, such a model cannot handle topological transitions, hiding the crucial involved forces and the appearance of the experimentally observed hemifused intermediates. In this work, we describe the membrane as a diffuse interface preserving the Canham-Helfrich elasticity. We show that pivotal features of the hemifusion pathway are captured by this mesoscopic approach, e.g. a (meta)stable hemifusion state and the fusogenic behavior of negative monolayer spontaneous curvatures. The membrane lateral stress profile is calculated as a function of the elastic rigidities, yielding a coarse-grained version of molecular models findings. Insights into the fusogenic mechanism are reported and discussed.

Reaction-driven assembly: controlling changes in membrane topology by reaction cycles

Chemical reaction cycles are prototypical examples how to drive systems out of equilibrium and introduce novel, life-like properties into soft-matter systems. We report simulations of amphiphilic molecules in aqueous solution. The molecule's head group is permanently hydrophilic, whereas the reaction cycle switches the molecule's tail from hydrophilic (precursor) to hydrophobic (amphiphile) and vice versa. The reaction cycle leads to an arrest in coalescence and results in uniform vesicle sizes that can be controlled by the reaction rate. Using a continuum description and particle-based simulation, we study the scaling of the vesicle size with the reaction rate. The chemically active vesicles are inflated by precursor, imparting tension onto the membrane and, for specific parameters, stabilize pores.

HylleraasMD: A Domain Decomposition-Based Hybrid Particle-Field Software for Multiscale Simulations of Soft Matter

We present HylleraasMD (HyMD), a comprehensive implementation of the recently proposed Hamiltonian formulation of hybrid particle-field molecular dynamics. The methodology is based on a tunable, grid-independent length-scale of coarse graining, obtained by filtering particle densities in reciprocal space. This enables systematic convergence of energies and forces by grid refinement, also eliminating nonphysical force aliasing. Separating the time integration of fast modes associated with internal molecular motion from slow modes associated with their density fields, we enable the first time-reversible, energy-conserving hybrid particle-field simulations. HyMD comprises the optional use of explicit electrostatics, which, in this formalism, corresponds to the long-range potential in particle-mesh Ewald. We demonstrate the ability of HyMD to perform simulations in the microcanonical and canonical ensembles with a series of test cases, comprising lipid bilayers and vesicles, surfactant micelles, and polypeptide chains, comparing our results to established literature. An on-the-fly increase of the characteristic coarse-grain length significantly speeds up dynamics, accelerating self-diffusion and leading to expedited aggregation. Exploiting this acceleration, we find that the time scales involved in the self-assembly of polymeric structures can lie in the tens to hundreds of picoseconds instead of the multimicrosecond regime observed with comparable coarse-grained models.

Dynamic Processes and Mechanical Properties of Lipid-Nanoparticle Mixtures

In this study, we investigate the dynamic processes and mechanical properties of lipid nanoparticle mixtures in a melt via dissipation particle dynamic simulation. By investigating the distribution of nanoparticles in lamellar and hexagonal lipid matrices in equilibrium state and dynamic processes, we observe that the morphology of such composites depends not only on the geometric features of the lipid matrix but also on the concentration of nanoparticles. The dynamic processes are also demonstrated by calculating the average radius of gyration, which indicates the isotropic conformation of lipid molecules in the x-y plane and that the lipid chains are stretched in the z direction with the addition of nanoparticles. Meanwhile, we predict the mechanical properties of lipid-nanoparticle mixtures in lamellar structures by analyzing the interfacial tensions. Results show that the interfacial tension decreased with the increase in nanoparticle concentration. These results provide molecular-level information for the rational and a priori design of new lipid nanocomposites with ad hoc tailored properties.

Soft Matter under Pressure: Pushing Particle–Field Molecular Dynamics to the Isobaric Ensemble

Hamiltonian hybrid particle-field molecular dynamics is a computationally efficient method to study large soft matter systems. In this work, we extend this approach to constant-pressure (NPT) simulations. We reformulate the calculation of internal pressure from the density field by taking into account the intrinsic spread of the particles in space, which naturally leads to a direct anisotropy in the pressure tensor. The anisotropic contribution is crucial for reliably describing the physics of systems under pressure, as demonstrated by a series of tests on analytical and monatomic model systems as well as realistic water/lipid biphasic systems. Using Bayesian optimization, we parametrize the field interactions of phospholipids to reproduce the structural properties of their lamellar phases, including area per lipid, and local density profiles. The resulting model excels in providing pressure profiles in qualitative agreement with all-atom modeling, and surface tension and area compressibility in quantitative agreement with experimental values, indicating the correct description of long-wavelength undulations in large membranes. Finally, we demonstrate that the model is capable of reproducing the formation of lipid droplets inside a lipid bilayer.

Nonequilibrium Processes in Polymer Membrane Formation: Theory and Experiment

Porous polymer and copolymer membranes are useful for ultrafiltration of functional macromolecules, colloids, and water purification. In particular, block copolymer membranes offer a bottom-up approach to form isoporous membranes. To optimize permeability, selectivity, longevity, and cost, and to rationally design fabrication processes, direct insights into the spatiotemporal structure evolution are necessary. Because of a multitude of nonequilibrium processes in polymer membrane formation, theoretical predictions via continuum models and particle simulations remain a challenge. We compiled experimental observations and theoretical approaches for homo- and block copolymer membranes prepared by nonsolvent-induced phase separation and highlight the interplay of multiple nonequilibrium processes─evaporation, solvent-nonsolvent exchange, diffusion, hydrodynamic flow, viscoelasticity, macro- and microphase separation, and dynamic arrest─that dictates the complex structure of the membrane on different scales.

Asymmetric Lipid Membranes under Shear Flows: A Dissipative Particle Dynamics Study

We investigate the phase behavior of the asymmetric lipid membranes under shear flows, using the dissipative particle dynamics simulation. Two cases, the weak and strong shear flows, are considered for the asymmetric lipid microstructures. Three typical asymmetric structures, the membranes, tubes, and vesicle, are included in the phase diagrams, where the effect of two different types of lipid chain length on the formation of asymmetric membranes is evaluated. The dynamic processes are demonstrated for the asymmetric membranes by calculating the average radius of gyration and shape factor. The result indicates that different shear flows will affect the shape of the second type of lipid molecules; the shape of the first type of lipid molecules is more stable than that of the second type of lipid molecules. The mechanical properties are investigated for the asymmetric membranes by analyzing the interface tension. The results reveal an absolute pressure at the junctions of different types of particles under the weak shear flow; the other positions are almost in a state of no pressure; there is almost no pressure inside the asymmetric lipid membrane structure under the strong shear flow. The findings will help us to understand the potential applications of asymmetric lipid microstructures in the biological and medical fields.

Virtual bending method to calculate bending rigidity, saddle-splay modulus, and spontaneous curvature of thin fluid membranes

A method to calculate the bending rigidity κ, saddle-splay modulus κ[over ¯], and spontaneous curvature C_{0} of a fluid membrane is proposed. Virtual work for the bending deformations into cylindrical and spherical shapes is calculated for a flat membrane. This method does not require a force decomposition, unlike the existing stress-profile method. The first derivative of the deformation gives κC_{0} and is a discrete form of the first moment of the stress profile. The second derivatives give κ and κ[over ¯] and include the variance terms of the first derivatives, which are not accounted for in the stress-profile method. This method is examined for a solvent-free meshless membrane model and a dissipative-particle-dynamics two-bead amphiphilic molecular model. It is concluded that κ and κ[over ¯] of a thin membrane can be accurately calculated, whereas for a thick membrane or one with an explicit solvent, a further extension to include the volume-fluctuation effects is required for an accurate estimation. The amplitude of the volume-fluctuation effects can be evaluated using the parameter dependence in the present method.

Regulation Mechanism of Bubbling Deformation and Fracture Toughness of the Membrane by Asymmetric Phospholipids: A Model System Study

Dynamic regulation of the deformation modulus and fracture toughness of a membrane is critical to organelles and cells for matching their conflicting needs of resilient and fractured behaviors. These properties implement the protection of the function in the normal condition and the fission function in the endocytosis condition of a membrane. Naturally, a membrane contains phospholipids that have different hydrophilic and hydrophobic group length. The diffusion and aggregation of the phospholipids with asymmetry of the hydrophilic-hydrophobic ratio on the membrane play a key role in regulating the mechanical behaviors passively to the external force. In present work, the effects of the asymmetry of phospholipids on the bubbling deformation and fracture toughness of the membrane to external stretching are investigated in a model system. A disk-shaped micelle formed from the blend of symmetric and asymmetric diblock copolymers in a selective solvent is considered as the membrane sheet. Its mechanically responsive behaviors are investigated by self-consistent field theory. By analyzing the evolution of different components during the stretching process, the mechanism of formation of the bubbling structure is revealed. Moreover, the fracture toughness depending on the asymmetry of the phospholipids is determined quantitatively.

Predicting asymmetric phospholipid microstructures in solutions

Asymmetric phospholipid microstructures, such as asymmetric phospholipid membranes, have potential applications in biological and medicinal processes. Here, we used the dissipative particle dynamics simulation method to predict the asymmetric phospholipid microstructures in aqueous solutions. The asymmetric phospholipid membranes, tubes and vesicles are determined and characterized by the chain density distributions and order parameters. The phase diagrams are constructed to evaluate the effects of the chain length on the asymmetric structure formations at equilibrium states, while the average radius of gyration and shape factors are calculated to analyze the asymmetric structure formations in the non-equilibrium processes. Meanwhile, we predicted the mechanical properties of the asymmetric membranes by analyzing the spatial distributions of the interface tensions and osmotic pressures in solutions.

Hybrid particle-field molecular dynamics under constant pressure

Hybrid particle-field methods are computationally efficient approaches for modeling soft matter systems. So far, applications of these methodologies have been limited to constant volume conditions. Here, we reformulate particle-field interactions to represent systems coupled to constant external pressure. First, we show that the commonly used particle-field energy functional can be modified to model and parameterize the isotropic contributions to the pressure tensor without interfering with the microscopic forces on the particles. Second, we employ a square gradient particle-field interaction term to model non-isotropic contributions to the pressure tensor, such as in surface tension phenomena. This formulation is implemented within the hybrid particle-field molecular dynamics approach and is tested on a series of model systems. Simulations of a homogeneous water box demonstrate that it is possible to parameterize the equation of state to reproduce any target density for a given external pressure. Moreover, the same parameterization is transferable to systems of similar coarse-grained mapping resolution. Finally, we evaluate the feasibility of the proposed approach on coarse-grained models of phospholipids, finding that the term between water and the lipid hydrocarbon tails is alone sufficient to reproduce the experimental area per lipid in constant-pressure simulations and to produce a qualitatively correct lateral pressure profile.

Generalization of the swelling method to measure the intrinsic curvature of lipids

Via computer simulation of a coarse-grained model of two-component lipid bilayers, we compare two methods of measuring the intrinsic curvatures of the constituting monolayers. The first one is a generalization of the swelling method that, in addition to the assumption that the spontaneous curvature linearly depends on the composition of the lipid mixture, incorporates contributions from its elastic energy. The second method measures the effective curvature-composition coupling between the apposing leaflets of bilayer structures (planar bilayers or cylindrical tethers) to extract the spontaneous curvature. Our findings demonstrate that both methods yield consistent results. However, we highlight that the two-leaflet structure inherent to the latter method has the advantage of allowing measurements for mixed lipid systems up to their critical point of demixing as well as in the regime of high concentration (of either species).