Toward Chemically Resolved Computer Simulations of Dynamics and Remodeling of Biological Membranes

Variable Non-Gaussian Transport of Nanoplastic on Supported Lipid Bilayers in Saline Conditions

Nanoplastic-lipid interaction is vital to understanding the nanoscale mechanism of plastic adsorption and aggregation on a lipid membrane surface. However, a single-particle mechanistic picture of the nanoplastic transport process on a lipid surface remains unclear. Here, we report a salt-dependent non-Gaussian transport mechanism of polystyrene particles on a supported 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) lipid bilayer surface. Particle stickiness on the POPC surface increases with salt concentration, where the particles stay longer at the surface and diffuse to shorter distances. Additionally, a non-Gaussian diffusion state dominates the transport process at high salt concentrations. Our current study provides insight into the transport mechanism of polystyrene (PS) particles on supported lipid membranes, which is essential to understanding fundamental questions regarding the adsorption mechanisms of nanoplastics on lipid surfaces.

Exploring the Molecular Dynamics of a Lipid-A Vesicle at the Atom Level: Morphology and Permeation Mechanism

Lipid-A was previously shown to spontaneously aggregate into a vesicle via the hybrid particle field approach. We assess the validity of the proposed vesiculation mechanism by simulating the resulting lipid-A vesicle at the atom level. The spatial confinement imposed by the vesicle geometry on the conformation and packing of lipid-A induces significant heterogeneity of physical properties in the inner and outer leaflets. It also induces tighter molecular packing and lower acyl chain order compared to the lamellar arrangement. Around 5% of water molecules passively permeates the vesicle membrane inward and outward. The permeation is facilitated by interactions with water molecules that are transported across the membrane by a network of electrostatic interactions with the hydrogen bond donors/acceptors in the N-acetylglucosamine ring and upper region of the acyl chains of lipid-A. The permeation process takes place at low rates but still at higher frequencies than observed for the lamellar arrangement of lipid-A. These findings not only substantiate the proposed lipid-A vesiculation mechanism but also reveal the complex structural dynamics of an important nonlamellar arrangement of lipid-A.

Combining particle and field-theoretic polymer models with multi-representation simulations

Particle-based and field-theoretic simulations are both widely used methods to predict the properties of polymeric materials. In general, the advantages of each method are complementary. Field-theoretic simulations are preferred for polymers with high molecular weights and can provide direct access to chemical potentials and free energies, which makes them the method-of-choice for calculating phase diagrams. The trade-off is that field-theoretic simulations sacrifice the molecular details present in particle-based simulations, such as the configurations of individual molecules and their dynamics. In this work, we describe a new approach to conduct "multi-representation" simulations that efficiently map between particle-based and field-theoretic simulations. Our approach involves the construction of formally equivalent particle-based and field-based models, which are then simulated subject to the constraint that their spatial density profiles are equal. This constraint provides the ability to directly link particle-based and field-based simulations and enables calculations that can switch between one representation to the other. By switching between particle/field representations during a simulation, we demonstrate that our approach can leverage many of the advantages of each representation while avoiding their respective limitations. Although our method is illustrated in the context of complex sphere phases in linear diblock copolymers, we anticipate that it will be useful whenever free energies, rapid equilibration, molecular configurations, and dynamic information are all simultaneously desired.

On the equivalence of the hybrid particle-field and Gaussian core models

Hybrid particle-field molecular dynamics is a molecular simulation strategy, wherein particles couple to a density field instead of through ordinary pair potentials. Traditionally considered a mean-field theory, a momentum and energy-conserving hybrid particle-field formalism has recently been introduced, which was demonstrated to approach the Gaussian Core model potential in the grid-converged limit. Here, we expand on and generalize the correspondence between the Hamiltonian hybrid particle-field method and particle-particle pair potentials. Using the spectral procedure suggested by Bore and Cascella, we establish compatibility to any local soft pair potential in the limit of infinitesimal grid spacing. Furthermore, we document how the mean-field regime often observed in hybrid particle-field simulations is due to the systems under consideration, and not an inherent property of the model. Considering the Gaussian filter form, in particular, we demonstrate the ability of the Hamiltonian hybrid particle-field model to recover all structural and dynamical properties of the Gaussian Core model, including solid phases, a first-order phase transition, and anomalous transport properties. We quantify the impact of the grid spacing on the correspondence, as well as the effect of the particle-field filtering length scale on the emergent particle-particle correlations.

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.

Slip-Spring Hybrid Particle-Field Molecular Dynamics for Coarse-Graining Branched Polymer Melts: Polystyrene Melts as an Example

The topology of chains significantly modifies the dynamical properties of polymer melts. Here, we extend a recently developed efficient simulation method, namely the slip-spring hybrid particle-field (SS-hPF) model, to study the structural and dynamical properties of branched polymer melts over large spatial-temporal scales. In the coarse-grained SS-hPF simulation of polymers, the bonded potentials are derived by iterative Boltzmann inversion from the underlying fine-grained model. The nonbonded potentials are computed from a density functional field instead of pairwise interactions used in standard molecular dynamics simulations, which increases the computational efficiency by a factor of 10-20. The entangled dynamics is lost due to the soft-core nature of density functional field interactions. It is recovered by a multichain slip-spring model that is rigorously parametrized from existing experimental or simulation data. To quantitatively predict the relaxation and diffusion of branched polymers, which are dominated by arm retraction rather than chain reptation, the slip-spring algorithm is augmented to improve the polymer dynamics near the branch point. Multiple dynamical observables, e.g., diffusion coefficients, arm relaxations, and tube survival probabilities, are characterized in an example coarse-grained model of symmetric and asymmetric star-shaped polystyrene melts. Consistent dynamical behaviors are identified and compared with theoretical predictions. With a single rescaling factor, the prediction of diffusion coefficients agrees well with the available experimental measurements. In this work, an efficient approach is provided to build chemistry-specific coarse-grained models for predicting the dynamics of branched polymers.

Investigating the structural properties of hydrophobic solvent-rich lipid bilayers

In vitro reconstitutions of lipid membranes have proven to be an indispensable tool to rationalize their molecular complexity and to understand their role in countless cellular processes. However, amongst the various techniques used to reconstitute lipid bilayers in vitro, several approaches are not solvent-free, but rather contain residual hydrophobic solvents in between the two bilayer leaflets, generally as a consequence of the procedure used to generate the bilayer. To what extent the presence of these hydrophobic solvents modifies bilayer properties with respect to native, solvent-free, conditions remains an open question that has important implications for the appropriate interpretation of numerous experimental observations. Here, we thorouhgly characterize hydrophobic solvent-rich lipid bilayers using atomistic molecular dynamics simulations. Our data indicate that while the presence of hydrophobic solvents at high concentrations, such as hexadecane, has a significant effect on membrane thickness, their effects on surface properties, membrane order and lateral stress are quite moderate. Our results corroborate the validity of in vitro approaches as model systems for the investigations of biological membranes but raise a few cautionary aspects that must be considered when investigating specific membrane properties.

Biomolecular modeling thrives in the age of technology

The biomolecular modeling field has flourished since its early days in the 1970s due to the rapid adaptation and tailoring of state-of-the-art technology. The resulting dramatic increase in size and timespan of biomolecular simulations has outpaced Moore's law. Here, we discuss the role of knowledge-based versus physics-based methods and hardware versus software advances in propelling the field forward. This rapid adaptation and outreach suggests a bright future for modeling, where theory, experimentation and simulation define three pillars needed to address future scientific and biomedical challenges.

Estimating the accuracy of the MARTINI model towards the investigation of peripheral protein–membrane interactions

In this article, we investigate the ability of the MARTINI CG force field, specifically the 3 open-beta version, to reproduce known experimental observations regarding the membrane binding behavior of 12 peripheral membrane proteins and peptides.

Out of Sight, Out of Mind: The Effect of the Equilibration Protocol on the Structural Ensembles of Charged Glycolipid Bilayers

Molecular dynamics (MD) simulations represent an essential tool in the toolbox of modern chemistry, enabling the prediction of experimental observables for a variety of chemical systems and processes and majorly impacting the study of biological membranes. However, the chemical diversity of complex lipids beyond phospholipids brings new challenges to well-established protocols used in MD simulations of soft matter and requires continuous assessment to ensure simulation reproducibility and minimize unphysical behavior. Lipopolysaccharides (LPS) are highly charged glycolipids whose aggregation in a lamellar arrangement requires the binding of numerous cations to oppositely charged groups deep inside the membrane. The delicate balance between the fully hydrated carbohydrate region and the smaller hydrophobic core makes LPS membranes very sensitive to the choice of equilibration protocol. In this work, we show that the protocol successfully used to equilibrate phospholipid bilayers when applied to complex lipopolysaccharide membranes occasionally leads to a small expansion of the simulation box very early in the equilibration phase. Although the use of a barostat algorithm controls the system dimension and particle distances according to the target pressure, fluctuation in the fleeting pressure occasionally enables a few water molecules to trickle into the hydrophobic region of the membrane, with spurious solvent buildup. We show that this effect stems from the initial steps of NPT equilibration, where initial pressure can be fairly high. This can be solved with the use of a stepwise-thermalization NVT/NPT protocol, as demonstrated for atomistic MD simulations of LPS/DPPE and lipid-A membranes in the presence of different salts using an extension of the GROMOS forcefield within the GROMACS software. This equilibration protocol should be standard procedure for the generation of consistent structural ensembles of charged glycolipids starting from atomic coordinates not previously pre-equilibrated. Although different ways to deal with this issue can be envisioned, we investigated one alternative that could be readily available in major MD engines with general users in mind.

Supramolecular Packing Drives Morphological Transitions of Charged Surfactant Micelles

The shape and size of self-assembled structures upon local organization of their molecular building blocks are hard to predict in the presence of long-range interactions. Combining small-angle X-ray/neutron scattering data, theoretical modelling, and computer simulations, sodium dodecyl sulfate (SDS), over a broad range of concentrations and ionic strengths, was investigated. Computer simulations indicate that micellar shape changes are associated with different binding of the counterions. By employing a toy model based on point charges on a surface, and comparing it to experiments and simulations, it is demonstrated that the observed morphological changes are caused by symmetry breaking of the irreducible building blocks, with the formation of transient surfactant dimers mediated by the counterions that promote the stabilization of cylindrical instead of spherical micelles. The present model is of general applicability and can be extended to all systems controlled by the presence of mobile charges.

Can Polarity-Inverted Surfactants Self-Assemble in Nonpolar Solvents?

We investigate the self-assembly process of a surfactant with inverted polarity in water and cyclohexane using both all-atom and coarse-grained hybrid particle-field molecular dynamics simulations. Unlike conventional surfactants, the molecule under study, proposed in a recent experiment, is formed by a rigid and compact hydrophobic adamantane moiety, and a long and floppy triethylene glycol tail. In water, we report the formation of stable inverted micelles with the adamantane heads grouping together into a hydrophobic core and the tails forming hydrogen bonds with water. By contrast, microsecond simulations do not provide evidence of stable micelle formation in cyclohexane. Validating the computational results by comparison with experimental diffusion constant and small-angle X-ray scattering intensity, we show that at laboratory thermodynamic conditions the mixture resides in the supercritical region of the phase diagram, where aggregated and free surfactant states coexist in solution. Our simulations also provide indications as to how to escape this region to produce thermodynamically stable micellar aggregates.

Large-scale simulation of biomembranes incorporating realistic kinetics into coarse-grained models

Biomembranes are two-dimensional assemblies of phospholipids that are only a few nanometres thick, but form micrometre-sized structures vital to cellular function. Explicit molecular modelling of biologically relevant membrane systems is computationally expensive due to the large number of solvent particles and slow membrane kinetics. Coarse-grained solvent-free membrane models offer efficient sampling but sacrifice realistic kinetics, thereby limiting the ability to predict pathways and mechanisms of membrane processes. Here, we present a framework for integrating coarse-grained membrane models with continuum-based hydrodynamics. This framework facilitates efficient simulation of large biomembrane systems with large timesteps, while achieving realistic equilibrium and non-equilibrium kinetics. It helps to bridge between the nanometer/nanosecond spatiotemporal resolutions of coarse-grained models and biologically relevant time- and lengthscales. As a demonstration, we investigate fluctuations of red blood cells, with varying cytoplasmic viscosities, in 150-milliseconds-long trajectories, and compare kinetic properties against single-cell experimental observations.

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.

Insight into the formation and permeability of ionic liquid unilamellar vesicles by molecular dynamics simulation

Unilamellar vesicles in solution could open up new horizons for reaction and material delivery, but the formation mechanism especially for the permeability of the small molecule through the vesicle membrane is still unknown. In this study, the formation and permeability of the unilamellar vesicles formed by the ionic liquid 1-dodecyl-3-methylimidazolium salicylate ([C12mim][Sal]) have been investigated by molecular dynamics simulation. Starting from a random distribution of ionic liquids, the entire process of vesicle formation could be observed on a nanosecond time scale, during which planar and cup-like structures are formed at the intermediate stage. Energy analysis reveals that the electrostatic interactions between cations and anions play a dominant role in forming and stabilizing the vesicle. Radial density distribution functions indicate that the final stable vesicle is a spherical bilayer structure. Besides, it was found that the structure of vesicles is maintained with the increase of temperature, while the water molecules in the vesicles could be completely exchanged quickly. These results suggest that vesicles may be beneficial for the enrichment or release of molecules.

Aggregation of Lipid A Variants: A Hybrid Particle-Field Model

Lipid A is one of the three components of bacterial lipopolysaccharides constituting the outer membrane of Gram-negative bacteria, and is recognized to have an important biological role in the inflammatory response of mammalians. Its biological activity is modulated by the number of acyl-chains that are present in the lipid and by the dielectric medium, i.e., the type of counter-ions, through electrostatic interactions. In this paper, we report on a coarse-grained model of chemical variants of Lipid A based on the hybrid particle-field/molecular dynamics approach (hPF-MD). In particular, we investigate the stability of Lipid A bilayers for two different hexa- and tetra-acylated structures. Comparing particle density profiles along bilayer cross-sections, we find good agreement between the hPF-MD model and reference all-atom simulation for both chemical variants of Lipid A. hPF-MD models of constituted bilayers composed by hexa-acylated Lipid A in water are stable within the simulation time. We further validate our model by verifying that the phase behavior of Lipid A/counterion/water mixtures is correctly reproduced. In particular, hPF-MD simulations predict the correct self-assembly of different lamellar and micellar phases from an initially random distribution of Lipid A molecules with counterions in water. Finally, it is possible to observe the spontaneous formation and stability of Lipid A vesicles by fusion of micellar aggregates.

Nanoparticles at Biomimetic Interfaces: Combined Experimental and Simulation Study on Charged Gold Nanoparticles/Lipid Bilayer Interfaces

The poor understanding of the interaction of nanomaterials with biologically relevant interfaces is recognized as one of the major issues currently limiting the development of nanomedicine. The central purpose of this study is to compare experimental (confocal microscopy, fluorescence correlation spectroscopy, X-ray reflectivity) and computational (molecular dynamics simulations) results to thoroughly describe the interaction of cationic gold nanoparticles (AuNPs) with mixed zwitterionic/anionic lipid membranes. The adhesion of AuNPs to the lipid membrane is investigated on different length scales from a structural and dynamical point of view; with this approach, a series of complex phenomena, spanning from lipid extraction, localized membrane disruption, lateral phase separation, and slaved diffusion, are characterized and interpreted from a molecular level to macroscopic observations.

Molecular structure and multi-body potential of mean force in silica-polystyrene nanocomposites

We perform a systematic application of the hybrid particle-field molecular dynamics technique [Milano, et al., J. Chem. Phys., 2009, 130, 214106] to study interfacial properties and potential of mean force (PMF) for separating nanoparticles (NPs) in a melt. Specifically, we consider Silica NPs bare or grafted with Polystyrene chains, aiming to shed light on the interactions among free and grafted chains affecting the dispersion of NPs in the nanocomposite. The proposed hybrid models show good performances in catching the local structure of the chains, and in particular their density profiles, documenting the existence of the "wet-brush-to-dry-brush" transition. By using these models, the PMF between pairs of ungrafted and grafted NPs in Polystyrene matrix are calculated. Moreover, we estimate the three-particle contribution to the total PMF and its role in regulating the phase separation on the nanometer scale. In particular, the multi-particle contribution to the PMF is able to give an explanation of the complex experimental morphologies observed at low grafting densities. More in general, we propose this approach and the models utilized here for a molecular understanding of specific systems and the impact of the chemical nature of the systems on the composite final properties.

Molecular dynamics simulations of lipid nanodiscs

A lipid nanodisc is a discoidal lipid bilayer stabilized by proteins, peptides, or polymers on its edge. Nanodiscs have two important connections to structural biology. The first is associated with high-density lipoprotein (HDL), a particle with a variety of functionalities including lipid transport. Nascent HDL (nHDL) is a nanodisc stabilized by Apolipoprotein A-I (APOA1). Determining the structure of APOA1 and its mimetic peptides in nanodiscs is crucial to understanding pathologies related to HDL maturation and designing effective therapies. Secondly, nanodiscs offer non-detergent membrane-mimicking environments and greatly facilitate structural studies of membrane proteins. Although seemingly similar, natural and synthetic nanodiscs are different in that nHDL is heterogeneous in size, due to APOA1 elasticity, and gradually matures to become spherical. Synthetic nanodiscs, in contrast, should be homogenous, stable, and size-tunable. This report reviews previous molecular dynamics (MD) simulation studies of nanodiscs and illustrates convergence and accuracy issues using results from new multi-microsecond atomistic MD simulations. These new simulations reveal that APOA1 helices take 10-20 μs to rearrange on the nanodisc, while peptides take 2 μs to migrate from the disc surfaces to the edge. These systems can also become kinetically trapped depending on the initial conditions. For example, APOA1 was trapped in a biologically irrelevant conformation for the duration of a 10 μs trajectory; the peptides were similarly trapped for 5 μs. It therefore remains essential to validate MD simulations of these systems with experiments due to convergence and accuracy issues. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Intramolecular structural parameters are key modulators of the gel-liquid transition in coarse grained simulations of DPPC and DOPC lipid bilayers

The capability of coarse-grained models based on the MARTINI mapping to reproduce the gel-liquid phase transition in saturated and unsaturated model lipids was investigated. We found that the model is able to reproduce a lower critical temperature for 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) with respect to 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC). Nonetheless, the appearance of a gel phase for DOPC is strictly dependent on the intramolecular parameters chosen to model its molecular structure. In particular, we show that the bending angle at the coarse-grained bead corresponding to the unsaturated carbon-carbon bond acts as an order parameter determining the temperature of the phase transition. Structural analysis of the molecular dynamics simulations runs evidences that in the gel phase, the packing of the lipophilic tails of DOPC assume a different conformation than in the liquid phase. In the latter phase, the DOPC geometry resembles that of the relaxed free molecule. DPPC:DOPC mixtures show a single phase transition temperature, indicating that the observation of a phase separation between the two lipids requires the simulation of systems with sizes much larger than the ones used here.

Biomembrane solubilization mechanism by Triton X-100: a computational study of the three stage model

The solubilization mechanism of lipid membranes in the presence of Triton X-100 (TX-100) is investigated at molecular resolution using hybrid particle field–self consistence field simulations.