Carbohydrate-derived amphiphilic macromolecules: a biophysical structural characterization and analysis of binding behaviors to model membranes.

The design and synthesis of enhanced membrane-intercalating biomaterials for drug delivery or vascular membrane targeting is currently challenged by the lack of screening and prediction tools. The present work demonstrates the generation of a Quantitative Structural Activity Relationship model (QSAR) to make a priori predictions. Amphiphilic macromolecules (AMs) "stealth lipids" built on aldaric and uronic acids frameworks attached to poly(ethylene glycol) (PEG) polymer tails were developed to form self-assembling micelles. In the present study, a defined set of novel AM structures were investigated in terms of their binding to lipid membrane bilayers using Quartz Crystal Microbalance with Dissipation (QCM-D) experiments coupled with computational coarse-grained molecular dynamics (CG MD) and all-atom MD (AA MD) simulations. The CG MD simulations capture the insertion dynamics of the AM lipophilic backbones into the lipid bilayer with the PEGylated tail directed into bulk water. QCM-D measurements with Voigt viscoelastic model analysis enabled the quantitation of the mass gain and rate of interaction between the AM and the lipid bilayer surface. Thus, this study yielded insights about variations in the functional activity of AM materials with minute compositional or stereochemical differences based on membrane binding, which has translational potential for transplanting these materials in vivo. More broadly, it demonstrates an integrated computational-experimental approach, which can offer a promising strategy for the in silico design and screening of therapeutic candidate materials

Construction of a synaptic membrane for in silico endocannabinoid investigations

The mechanism of a chemical reaction cannot be assumed prior to identifying the solvent, thus to accurately model membrane-bound cannabinoid receptors, details of the lipid bilayer must be considered as well. Traditionally thought of as a permeable barrier, the membrane also serves as a medium for hydrophobic ligand delivery; the Reggio group has published on the entry of sn-2-arachidonoylglycerol into the cannabinoid CB2 receptor through the trans-membrane region. Currently receptor and ligand trajectories are calculated in a single-phospholipid bilayer, however biological membranes are comprised of many lipid species. Herein nanosecond timescale molecular dynamics studies of various membrane compositions are reported, culminating to a final model mimicking a synaptic plasma membrane in phospholipid composition. Structural changes to phospholipids in increasingly heterogeneous compositions are discussed, along with interactions between endocannabinoid N-arachidonoylglycine (NAGly) and membrane components. The final model provides a significant contrast to a pure phosphatidylcholine bilayer, in addition to proposing new investigations of ligand diffusion, phospholipid clustering, and similar hydrogen bonding studies

The ELBA Force Field for Coarse-Grain Modeling of Lipid Membranes

A new coarse-grain model for molecular dynamics simulation of lipid membranes is presented. Following a simple and conventional approach, lipid molecules are modeled by spherical sites, each representing a group of several atoms. In contrast to common coarse-grain methods, two original (interdependent) features are here adopted. First, the main electrostatics are modeled explicitly by charges and dipoles, which interact realistically through a relative dielectric constant of unity (). Second, water molecules are represented individually through a new parametrization of the simple Stockmayer potential for polar fluids; each water molecule is therefore described by a single spherical site embedded with a point dipole. The force field is shown to accurately reproduce the main physical properties of single-species phospholipid bilayers comprising dioleoylphosphatidylcholine (DOPC) and dioleoylphosphatidylethanolamine (DOPE) in the liquid crystal phase, as well as distearoylphosphatidylcholine (DSPC) in the liquid crystal and gel phases. Insights are presented into fundamental properties and phenomena that can be difficult or impossible to study with alternative computational or experimental methods. For example, we investigate the internal pressure distribution, dipole potential, lipid diffusion, and spontaneous self-assembly. Simulations lasting up to 1.5 microseconds were conducted for systems of different sizes (128, 512 and 1058 lipids); this also allowed us to identify size-dependent artifacts that are expected to affect membrane simulations in general. Future extensions and applications are discussed, particularly in relation to the methodology's inherent multiscale capabilities

The order-disorder transition in model lipid bilayers is a first-order hexatic to liquid phase transition

We characterize the order-disorder transition in a model lipid bilayer using molecular dynamics simulations. We find that the ordered phase is hexatic. In particular, in-plane structures possess a finite concentration of 5-7 disclination pairs that diffuse throughout the plane of the bilayer, and further, in-plane structures exhibit long-range orientational order and short-range translational order. In contrast, the disordered phase is liquid. The transition between the two phases is first order. Specifically, it exhibits hysteresis, and coexistence exhibits an interface with capillary scaling. The location of the interface and its spatial fluctuations are analyzed with a spatial field constructed from a rotational-invariant for local 6-fold orientational order. As a result of finite interfacial tension, there necessarily exist associated forces of assembly between membrane-bound solutes that pre-melt the ordered phase.Comment: Addressed the comments from colleagues, corrected typos, clarified text, updated references. The new draft also contains new results relating to the hexatic phas

Determination of biomembrane bending moduli in fully atomistic simulations.

The bilayer bending modulus (Kc) is one of the most important physical constants characterizing lipid membranes, but precisely measuring it is a challenge, both experimentally and computationally. Experimental measurements on chemically identical bilayers often differ depending upon the techniques employed, and robust simulation results have previously been limited to coarse-grained models (at varying levels of resolution). This Communication demonstrates the extraction of Kc from fully atomistic molecular dynamics simulations for three different single-component lipid bilayers (DPPC, DOPC, and DOPE). The results agree quantitatively with experiments that measure thermal shape fluctuations in giant unilamellar vesicles. Lipid tilt, twist, and compression moduli are also reported

Accurate Binding of Sodium and Calcium to a POPC Bilayer by Effective Inclusion of Electronic Polarization

Binding affinities and stoichiometries of Na+ and Ca2+ ions to phospholipid bilayers are of paramount significance in the properties and functionality of cellular membranes. Current estimates of binding affinities and stoichiometries of cations are, however, inconsistent due to limitations in the available experimental and computational methods. In this work, we improve the description of the binding details of Na+ and Ca2+ ions to a 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayer by implicitly including electronic polarization as a mean field correction, known as the electronic continuum correction (ECC). This is applied by scaling the partial charges of a selected state-of-the-art POPC lipid model for molecular dynamics simulations. Our improved ECC-POPC model reproduces not only the experimentally measured structural parameters for the ion-free membrane, but also the response of lipid headgroup to a strongly bound cationic amphiphile, as well as the binding affinities of Na+ and Ca2+ ions. With our new model, we observe on the one side negligible binding of Na+ ions to POPC bilayer, while on the other side stronger interactions of Ca2+ primarily with phosphate oxygens, which is in agreement with the previous interpretations of the experimental spectroscopic data. The present model results in Ca2+ ions forming complexes with one to three POPC molecules with almost equal probabilities, suggesting more complex binding stoichiometries than those from simple models used to interpret the NMR data previously. The results of this work pave the way to quantitative molecular simulations with realistic electrostatic interactions of complex biochemical systems at cellular membranes.Peer reviewe