Structure of fully hydrated fluid phase lipid bilayers with monounsaturated chains

Structure of the gel phase of diC22:1PC lipid bilayers determined by x-ray diffraction

High-resolution x-ray data are reported for the ordered phases of long-chain di-monounsaturated C22:1 phosphocholine lipid bilayers. Similar to PC lipids that have saturated chains, diC22:1PC has a subgel phase and a gel phase, but dissimilarly, we find no ripple phase. Our quantitative focus is on the structure of the gel phase. We have recorded 17 lamellar orders, indicating a very well-ordered structure. Fitting to a model provides the phases of the orders. The Fourier construction of the electron density profile has two well-defined headgroup peaks and a very sharp and deep methyl trough. The wide-angle scattering exhibits two Bragg rods that provide the area per molecule. They have an intensity pattern quite different than that of lipids with saturated chains. Models of chain packing indicate that ground state chain configurations are tilted primarily toward next nearest neighbors with an angle that is also consistent with the modeling of the electron density profile. Wide-angle modeling also indicates broken mirror symmetry between the monolayers. Our wide-angle results and our electron density profile together leads to the hypothesis that the sn-1 and sn-2 chains have equivalent penetration depths in contrast to the gel phase structure of lipids with saturated hydrocarbon chains.

Exhaustive in silico design and screening of novel antipsychotic compounds with improved pharmacodynamics and blood-brain barrier permeation properties

Antipsychotic drugs or neuroleptics are widely used in the treatment of psychosis as a manifestation of schizophrenia and bipolar disorder. However, their effectiveness largely depends on the blood-brain barrier (BBB) permeation (pharmacokinetics) and drug-receptor pharmacodynamics. Therefore, in this study, we developed and implemented the in silico pipeline to design novel compounds (n = 260) as leads using the standard drug scaffolds with improved PK/PD properties from the standard scaffolds. As a result, the best candidates (n = 3) were evaluated in molecular docking to interact with serotonin and dopamine receptors. Finally, haloperidol (HAL) derivative (1-(4-fluorophenyl)-4-(4-hydroxy-4-{4-[(2-phenyl-1,3-thiazol-4-yl)methyl]phenyl}piperidin-1-yl)butan-1-one) was identified as a "magic shotgun" lead compound with better affinity to the 5-HT2A, 5-HT1D, D2, D3, and 5-HT1B receptors than the control molecule. Additionally, this hit substance was predicted to possess similar BBB permeation properties and much lower toxicological profiles in comparison to HAL. Overall, the proposed rational drug design platform for novel antipsychotic drugs based on the BBB permeation and receptor binding might be an invaluable asset for a medicinal chemist or translational pharmacologist.Communicated by Ramaswamy H. Sarma.

Effect of Cholesterol on the Structure and Composition of Glyco-DIBMA Lipid Particles

Different amphiphilic co-polymers have been introduced to produce polymer-lipid particles with nanodisc structure composed of an inner lipid bilayer and polymer chains self-assembled as an outer belt. These particles can be used to stabilize membrane proteins in solution and enable their characterization by means of biophysical methods, including small-angle X-ray scattering (SAXS). Some of these co-polymers have also been used to directly extract membrane proteins together with their associated lipids from native membranes. Styrene/maleic acid and diisobutylene/maleic acid are among the most commonly used co-polymers for producing polymer-lipid particles, named SMALPs and DIBMALPs, respectively. Recently, a new co-polymer, named Glyco-DIBMA, was produced by partial amidation of DIBMA with the amino sugar N-methyl-d-glucosamine. Polymer-lipid particles produced with Glyco-DIBMA, named Glyco-DIBMALPs, exhibit improved structural properties and stability compared to those of SMALPs and DIBMALPs while retaining the capability of directly extracting membrane proteins from native membranes. Here, we characterize the structure and lipid composition of Glyco-DIBMALPs produced with either 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) or 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). Glyco-DIBMALPs were also prepared with mixtures of either POPC or DMPC and cholesterol at different mole fractions. We estimated the lipid content in the Glyco-DIBMALPs and determined the particle structure and morphology by SAXS. We show that the Glyco-DIBMALPs are nanodisc-like particles whose size and shape depend on the polymer/lipid ratio. This is relevant for designing nanodisc particles with a tunable diameter according to the size of the membrane protein to be incorporated. We also report that the addition of >20 mol % cholesterol strongly perturbed the formation of Glyco-DIBMALPs. Altogether, we describe a detailed characterization of the Glyco-DIBMALPs, which provides relevant inputs for future application of these particles in the biophysical investigation of membrane proteins.

Confinement effect on hydrolysis in small lipid vesicles

In living organisms most chemical reactions take place within the confines of lipid-membrane bound compartments, while confinement within the bounds of a lipid membrane is thought to be a key step in abiogenesis. In previous work we demonstrated that confinement in the aqueous cavity of a lipid vesicle affords protection against hydrolysis, a phenomenon that we term here confinement effect (C _e) and that we attributed to the interaction with the lipid membrane. Here, we show that both the size and the shape of the cavity of the vesicle modulate the C _e. We link this observation to the packing of the lipid following changes in membrane curvature, and formulate a mathematical model that relates the C _e to the radius of a spherical vesicle and the packing parameter of the lipids. These results suggest that the shape of the compartment where a molecule is located plays a major role in controlling the chemical reactivity of non-enzymatic reactions. Moreover, the mathematical treatment we propose offers a useful tool for the design of vesicles with predictable reaction rates of the confined molecules, e.g., drug delivery vesicles with confined prodrugs. The results also show that a crude form of signal transduction, devoid of complex biological machinery, can be achieved by any external stimuli that drastically changes the structure of the membrane, like the osmotic shocks used in the present work.

A skin pharmacokinetics study of permeation enhancers: The root cause of dynamic enhancement effect on in vivo drug permeation

Skin pharmacokinetics (SPK) of permeation enhancers can answer the question of why enhancement effects different at the kinetic level. Herein, SPK of permeation enhancers were classified into two categories, namely, lateral elimination (elimination to surrounding stratum corneum (SC)) and longitudinal elimination (elimination to deep epidermal (EP)). They were evaluated with a specific parameter for permeation enhancers, diffusion ratio (DR_SC-EP), according to results of tissue-distribution test, molecular dynamic (MD) simulation, and confocal laser scanning microscopy (CLSM). The linear relationship between k_e-enahcer and Δ C_max-drug (R² = 0.92), MRT_enhancer and Δ T_max-drug (R² = 0.97), AUC_t-enhancer and Δ AUC_t-drug (R² = 0.90) suggesting that SPK of permeation enhancers precisely controlled dynamic process of drug permeation in vivo. The molecular mechanisms of the dynamic effect of SPK process on drug transdermal behaviors were characterized by modulated-temperature differential scanning calorimetry (MTDSC), dielectric spectroscopy, small-angle X-ray scattering (SAXS), solid-state NMR. Permeation enhancers with high molecular weight (M.W.) and high polar surface area (P.S.A.) had good compatibility and strong interaction strength with SC, leading their lateral-elimination behavior, causing their low DR_SC-EP and resulting in low k_e-enhancer, long MRT_enhancer, and large AUC_t-enhancer. Consequently, skin barrier can be rapidly opened fast and to a great extent. In summary, compared with SPK of permeation enhancers with longitudinal elimination, SPK of permeation enhancers with lateral elimination can enable more sustainable and greater drug permeation. The information about SPK of permeation enhancers offered a criterion to estimate its permeation-enhancement effect on the drug and its subsequent application in transdermal formulations.

Incorporation and localisation of alkanes in a protomembrane model by neutron diffraction

Protomembranes at the origin of life were likely composed of short-chain lipids, readily available on the early Earth. Membranes formed by such lipids are less stable and more permeable under extreme conditions, so a novel membrane architecture was suggested to validate the accuracy of this assumption. The model membrane includes the presence of a layer of alkanes in the mid-plane of the protomembrane in between the two monolayer leaflets and lying perpendicular to the lipid acyl chains. Here, we investigated such a possibility experimentally for membranes formed by the short-chain phospholipid 1,2-didecanoyl-sn-glycero-3-phophocholine, including or not the alkanes eicosane, squalane or triacontane by means of neutron membrane diffraction and contrast variation. We found strong indications for incorporation of two of the three alkanes in the membrane mid-plane through the determination of neutron scattering length density profiles with hydrogenated vs deuterated alkanes and membrane swelling at various relative humidities indicating a slightly increased bilayer thickness when the alkanes are incorporated into the bilayers. The selectivity of the incorporation points out the role of the length of the n-alkanes with respect to the capacity of the membrane to incorporate them.

Coarse-Grained Molecular Dynamics of pH-Sensitive Lipids

pH-sensitive lipids represent a class of lipids that can be protonated and destabilized in acidic environments, as they become positively charged in response to low-pH conditions. They can be incorporated into lipidic nanoparticles such as liposomes, which are able to change their properties and allow specific drug delivery at the acidic conditions encountered in some pathological microenvironments. In this work, we used coarse-grained molecular-dynamic simulations to study the stability of neutral and charged lipid bilayers containing POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and various kinds of ISUCA ((F)2-(imidazol-1-yl)succinic acid)-derived lipids, which can act as pH-sensitive molecules. In order to explore such systems, we used a MARTINI-derived forcefield, previously parameterized using all-atom simulation results. We calculated the average area per lipid, the second-rank order parameter and the lipid diffusion coefficient of both lipid bilayers made of pure components and mixtures of lipids in different proportions, under neutral or acidic conditions. The results show that the use of ISUCA-derived lipids disturbs the lipid bilayer structure, with the effect being particularly marked under acidic conditions. Although more-in depth studies on these systems must be carried out, these initial results are encouraging and the lipids designed in this research could be a good basis for developing new pH-sensitive liposomes.

The Effect of Cholesterol in SOPC Lipid Bilayers at Low Temperatures

We study the behavior of lipid bilayers composed of SOPC (1-stearoyl-2-oleoyl-sn-glycero-3-phosphocholine) with different concentrations of cholesterol, ranging from 10 mol% to 50 mol% at 273 K. To this end, we carry out extensive atomistic molecular dynamic simulations with the aid of the Slipid force field aiming at computing basic bilayer parameters, as well as thermodynamic properties and structural characteristics. The obtained results are compared to available relevant experimental data and the outcome of atomistic simulations performed on bilayers composed of analogous phospholipids. Our results show a good quantitative, as well as qualitative, agreement with the main trends associated with the concentration increase in cholesterol. Moreover, it comes out that a change in the behavior of the bilayer is brought about at a concentration of about 30 mol% cholesterol. At this very concentration, some of the bilayer properties are found to exhibit a saturation and a significant long-range ordering of the lipid molecules in the membrane shows up.

How Does the Study MD of pH-Dependent Exposure of Nanoparticles Affect Cellular Uptake of Anticancer Drugs?

The lack of knowledge about the uptake of NPs by biological cells poses a significant problem for drug delivery. For this reason, designing an appropriate model is the main challenge for modelers. To address this problem, molecular modeling studies that can describe the mechanism of cellular uptake of drug-loaded nanoparticles have been conducted in recent decades. In this context, we developed three different models for the amphipathic nature of drug-loaded nanoparticles (MTX-SS-γ-PGA), whose cellular uptake mechanism was predicted by molecular dynamics studies. Many factors affect nanoparticle uptake, including nanoparticle physicochemical properties, protein-particle interactions, and subsequent agglomeration, diffusion, and sedimentation. Therefore, the scientific community needs to understand how these factors can be controlled and the NP uptake of nanoparticles. Based on these considerations, in this study, we investigated for the first time the effects of the selected physicochemical properties of the anticancer drug methotrexate (MTX) grafted with hydrophilic-γ-polyglutamic acid (MTX-SS-γ-PGA) on its cellular uptake at different pH values. To answer this question, we developed three theoretical models describing drug-loaded nanoparticles (MTX-SS-γ-PGA) at three different pH values, such as (1) pH 7.0 (the so-called neutral pH model), (2) pH 6.4 (the so-called tumor pH model), and (3) pH 2.0 (the so-called stomach pH model). Exceptionally, the electron density profile shows that the tumor model interacts more strongly with the head groups of the lipid bilayer than the other models due to charge fluctuations. Hydrogen bonding and RDF analyses provide information about the solution of the NPs with water and their interaction with the lipid bilayer. Finally, dipole moment and HOMO-LUMO analysis showed the free energy of the solution in the water phase and chemical reactivity, which are particularly useful for determining the cellular uptake of the NPs. The proposed study provides fundamental insights into molecular dynamics (MD) that will allow researchers to determine the influence of pH, structure, charge, and energetics of NPs on the cellular uptake of anticancer drugs. We believe that our current study will be useful in developing a new model for drug delivery to cancer cells with a much more efficient and less time-consuming model.

Production of Lipid Constructs by Design via Three-Dimensional Nanoprinting

Atomic force microscopy (AFM) in conjunction with microfluidic delivery was utilized to produce three-dimensional (3D) lipid structures following a custom design. While AFM is well-known for its spatial precision in imaging and 2D nanolithography, the development of AFM-based nanotechnology into 3D nanoprinting requires overcoming the technical challenges of controlling material delivery and interlayer registry. This work demonstrates the concept of 3D nanoprinting of amphiphilic molecules such as 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). Various formulations of POPC solutions were tested to achieve point, line, and layer-by-layer material delivery. The produced structures include nanometer-thick disks, long linear spherical caps, stacking grids, and organizational chiral architectures. The POPC molecules formed stacking bilayers in these constructions, as revealed by high-resolution structural characterizations. The 3D printing reached nanometer spatial precision over a range of 0.5 mm. The outcomes reveal the promising potential of our designed technology and methodology in the production of 3D structures from nanometer to continuum, opening opportunities in biomaterial sciences and engineering, such as in the production of 3D nanodevices, chiral nanosensors, and scaffolds for tissue engineering and regeneration.

Determining the Bending Rigidity of Free-Standing Planar Phospholipid Bilayers

We describe a method to determine membrane bending rigidity from capacitance measurements on large area, free-standing, planar, biomembranes. The bending rigidity of lipid membranes is an important biological mechanical property that is commonly optically measured in vesicles, but difficult to quantify in a planar, unsupported system. To accomplish this, we simultaneously image and apply an electric potential to free-standing, millimeter area, planar lipid bilayers composed of DOPC and DOPG phospholipids to measure the membrane Young's (elasticity) modulus. The bilayer is then modeled as two adjacent thin elastic films to calculate bending rigidity from the electromechanical response of the membrane to the applied field. Using DOPC, we show that bending rigidities determined by this approach are in good agreement with the existing work using neutron spin echo on vesicles, atomic force spectroscopy on supported lipid bilayers, and micropipette aspiration of giant unilamellar vesicles. We study the effect of asymmetric calcium concentration on symmetric DOPC and DOPG membranes and quantify the resulting changes in bending rigidity. This platform offers the ability to create planar bilayers of controlled lipid composition and aqueous ionic environment, with the ability to asymmetrically alter both. We aim to leverage this high degree of compositional and environmental control, along with the capacity to measure physical properties, in the study of various biological processes in the future.

Structure of POPC Lipid Bilayers in OPLS3e Force Field

It is crucial for molecular dynamics simulations of biomembranes that the force field parameters give a realistic model of the membrane behavior. In this study, we examined the OPLS3e force field for the carbon-hydrogen order parameters S_CH of POPC (1-palmitoyl-2-oleoylphosphatidylcholine) lipid bilayers at varying hydration conditions and ion concentrations. The results show that OPLS3e behaves similarly to the CHARMM36 force field and relatively accurately follows the experimentally measured S_CH for the lipid headgroup, the glycerol backbone, and the acyl tails. Thus, OPLS3e is a good choice for POPC bilayer simulations under many biologically relevant conditions. The exception are systems with an abundancy of ions, as similarly to most other force fields OPLS3e strongly overestimates the membrane-binding of cations, especially Ca²⁺. This leads to undesirable positive charge of the membrane surface and drastically lowers the concentration of Ca²⁺ in the surrounding solvent, which might cause issues in systems sensitive to correct charge distribution profiles across the membrane.

Stepwise remodeling and subcompartment formation in individual vesicles by three ESCRT-III proteins

The endosomal sorting complex required for transport (ESCRT) is a multi-protein machinery involved in several membrane remodeling processes. Different approaches have been used to resolve how ESCRT proteins scission membranes. However, the underlying mechanisms generating membrane deformations are still a matter of debate. Here, giant unilamellar vesicles, microfluidic technology, and micropipette aspiration are combined to continuously follow the ESCRT-III-mediated membrane remodeling on the single-vesicle level for the first time. With this approach, we identify different mechanisms by which a minimal set of three ESCRT-III proteins from Entamoeba histolytica reshape the membrane. These proteins modulate the membrane stiffness and spontaneous curvature to regulate bud size and generate intraluminal vesicles even in the absence of ATP. We demonstrate that the bud stability depends on the protein concentration and membrane tension. The approaches introduced here should open the road to diverse applications in synthetic biology for establishing artificial cells with several membrane compartments.

Extensive Molecular Dynamics Simulations Disclosed the Stability of mPGES-1 Enzyme and the Structural Role of Glutathione (GSH) Cofactor

A deep in silico investigation of various microsomal prostaglandin E₂ synthase-1 (mPGES-1) protein systems is here reported using molecular dynamics (MD) simulations. Firstly, eight different proteins models (Models A-H) were built, starting from the active enzyme trimer system (Model A), namely that bound to three glutathione (GSH) cofactor molecules, and then gradually removing the GSHs (Models B-H), simulating each of them for 100 ns in explicit solvent. The analysis of the obtained data disclosed the structural role of GSH in the chemical architecture of mPGES-1 enzyme, thus suggesting the unlikely displacement of this cofactor, in accordance with experimentally determined protein structures co-complexed with small molecule inhibitors. Afterwards, Model A was submitted to microsecond-scale molecular dynamics simulation (total simulation time=10 μs), in order to shed light about the dynamical behaviour of this enzyme at atomic level and to obtain further structural features and protein function information. We confirmed the structural stability of the enzyme machinery, observing a conformational rigidity of the protein, with a backbone RMSD of ∼3 Å along the simulation time, and highlighting the strong active contribution of GSH molecules due to their active role in packing the protein chains through a tight binding at monomer interfaces. Furthermore, the focused analysis on R73 residue disclosed its role in solvent exchange events, probably excluding its function as route for GSH to enter towards the endoplasmic reticulum membrane, in line with the recently reported function of cap domain residues F44-D66 as gatekeeper for GSH entrance into catalytic site.

Zwitterionic Dendrimersomes: A Closer Xenobiotic Mimic of Cell Membranes

Building functional mimics of cell membranes is an important task toward the development of synthetic cells. So far, lipid and amphiphilic block copolymers are the most widely used amphiphiles with the bilayers by the former lacking stability while membranes by the latter are typically characterized by very slow dynamics. Herein, a new type of Janus dendrimer containing a zwitterionic phosphocholine hydrophilic headgroup (JD^PC ) and a 3,5-substituted dihydrobenzoate-based hydrophobic dendron is introduced. JD^PC self-assembles in water into zwitterionic dendrimersomes (z-DSs) that faithfully recapitulate the cell membrane in thickness, flexibility, and fluidity, while being resilient to harsh conditions and displaying faster pore closing dynamics in the event of membrane rupture. This enables the fabrication of hybrid DSs with components of natural membranes, including pore-forming peptides, structure-directing lipids, and glycans to create raft-like domains or onion vesicles. Moreover, z-DSs can be used to create active synthetic cells with life-like features that mimic vesicle fusion and motility as well as environmental sensing. Despite their fully synthetic nature, z-DSs are minimal cell mimics that can integrate and interact with living matter with the programmability to imitate life-like features and beyond.

Molecular recognition of a membrane-anchored HIV-1 pan-neutralizing epitope

Antibodies against the carboxy-terminal section of the membrane-proximal external region (C-MPER) of the HIV-1 envelope glycoprotein (Env) are considered as nearly pan-neutralizing. Development of vaccines capable of producing analogous broadly neutralizing antibodies requires deep understanding of the mechanism that underlies C-MPER recognition in membranes. Here, we use the archetypic 10E8 antibody and a variety of biophysical techniques including single-molecule approaches to study the molecular recognition of C-MPER in membrane mimetics. In contrast to the assumption that an interfacial MPER helix embodies the entire C-MPER epitope recognized by 10E8, our data indicate that transmembrane domain (TMD) residues contribute to binding affinity and specificity. Moreover, anchoring to membrane the helical C-MPER epitope through the TMD augments antibody binding affinity and relieves the effects exerted by the interfacial MPER helix on the mechanical stability of the lipid bilayer. These observations support that addition of TMD residues may result in more efficient and stable anti-MPER vaccines.

Polystyrene-modified carbon nanotubes: Promising carriers in targeted drug delivery

To design drug-delivery agents for therapeutic and diagnostic applications, understanding the mechanisms by which covalently functionalized carbon nanotubes penetrate and interact with cell membranes is of great importance. Here, we report all-atom molecular dynamics results from polystyrene and carboxyl-terminated polystyrene-modified carbon nanotubes and show their translocation behavior across a model lipid bilayer together with their potential to deliver a molecule of the drug ibuprofen into the cell. Our results indicate that functionalized carbon nanotubes are internalized by the membrane in hundreds of nanoseconds and that drug loading increases the internalization speed further. Both loaded and unloaded tubes cross the closest leaflet of the bilayer by nonendocytic pathways, and for the times studied, the drug molecule remains trapped inside the pristine tube while remaining attached at the end of polystyrene-modified tube. On the other hand, carboxyl-terminated polystyrene functionalization allows the drug to be completely released into the lower leaflet of the bilayer without imposing damage to the membrane. This study shows that polystyrene functionalization is a promising alternative and facilitates drug delivery as a benchmark case.

Phospholipid headgroups govern area per lipid and emergent elastic properties of bilayers

Phospholipid bilayers are liquid-crystalline materials whose intermolecular interactions at mesoscopic length scales have key roles in the emergence of membrane physical properties. Here we investigated the combined effects of phospholipid polar headgroups and acyl chains on biophysical functions of membranes with solid-state 2H NMR spectroscopy. We compared the structural and dynamic properties of phosphatidylethanolamine and phosphatidylcholine with perdeuterated acyl chains in the solid-ordered (so) and liquid-disordered (ld) phases. Our analysis of spectral lineshapes of 1,2-diperdeuteriopalmitoyl-sn-glycero-3-phosphoethanolamine (DPPE-d62) and 1,2-diperdeuteriopalmitoyl-sn-glycero-3-phosphocholine (DPPC-d62) in the so (gel) phase indicated an all-trans rotating chain structure for both lipids. Greater segmental order parameters (SCD) were observed in the ld (liquid-crystalline) phase for DPPE-d62 than for DPPC-d62 membranes, while their mixtures had intermediate values irrespective of the deuterated lipid type. Our results suggest the SCD profiles of the acyl chains are governed by methylation of the headgroups and are averaged over the entire system. Variations in the acyl chain molecular dynamics were further investigated by spin-lattice (R1Z) and quadrupolar-order relaxation (R1Q) measurements. The two acyl-perdeuterated lipids showed distinct differences in relaxation behavior as a function of the order parameter. The R1Z rates had a square-law dependence on SCD, implying collective mesoscopic dynamics, with a higher bending rigidity for DPPE-d62 than for DPPC-d62 lipids. Remodeling of lipid average and dynamic properties by methylation of the headgroups thus provides a mechanism to control the actions of peptides and proteins in biomembranes.

Small-Angle Neutron Scattering for Studying Lipid Bilayer Membranes

Small-angle neutron scattering (SANS) is a powerful tool for studying biological membranes and model lipid bilayer membranes. The length scales probed by SANS, being from 1 nm to over 100 nm, are well-matched to the relevant length scales of the bilayer, particularly when it is in the form of a vesicle. However, it is the ability of SANS to differentiate between isotopes of hydrogen as well as the availability of deuterium labeled lipids that truly enable SANS to reveal details of membranes that are not accessible with the use of other techniques, such as small-angle X-ray scattering. In this work, an overview of the use of SANS for studying unilamellar lipid bilayer vesicles is presented. The technique is briefly presented, and the power of selective deuteration and contrast variation methods is discussed. Approaches to modeling SANS data from unilamellar lipid bilayer vesicles are presented. Finally, recent examples are discussed. While the emphasis is on studies of unilamellar vesicles, examples of the use of SANS to study intact cells are also presented.

Transmembrane proteins tetraspanin 4 and CD9 sense membrane curvature

Multiple membrane-shaping and remodeling processes are associated with tetraspanin proteins by yet unknown mechanisms. Tetraspanins constitute a family of proteins with four transmembrane domains present in every cell type. Prominent examples are tetraspanin4 and CD9, which are required for the fundamental cellular processes of migrasome formation and fertilization, respectively. These proteins are enriched in curved membrane structures, such as cellular retraction fibers and oocyte microvilli. The factors driving this enrichment are, however, unknown. Here, we revealed that tetraspanin4 and CD9 are curvature sensors with a preference for positive membrane curvature. To this end, we used a biomimetic system emulating membranes of cell retraction fibers and oocyte microvilli by membrane tubes pulled out of giant plasma membrane vesicles with controllable membrane tension and curvature. We developed a simple thermodynamic model for the partitioning of curvature sensors between flat and tubular membranes, which allowed us to estimate the individual intrinsic curvatures of the two proteins. Overall, our findings illuminate the process of migrasome formation and oocyte microvilli shaping and provide insight into the role of tetraspanin proteins in membrane remodeling processes.

Structural and mechanical properties of the red blood cell's cytoplasmic membrane seen through the lens of biophysics

Red blood cells (RBCs) are the most abundant cell type in the human body and critical suppliers of oxygen. The cells are characterized by a simple structure with no internal organelles. Their two-layered outer shell is composed of a cytoplasmic membrane (RBC cm ) tethered to a spectrin cytoskeleton allowing the cell to be both flexible yet resistant against shear stress. These mechanical properties are intrinsically linked to the molecular composition and organization of their shell. The cytoplasmic membrane is expected to dominate the elastic behavior on small, nanometer length scales, which are most relevant for cellular processes that take place between the fibrils of the cytoskeleton. Several pathologies have been linked to structural and compositional changes within the RBC cm and the cell's mechanical properties. We review current findings in terms of RBC lipidomics, lipid organization and elastic properties with a focus on biophysical techniques, such as X-ray and neutron scattering, and Molecular Dynamics simulations, and their biological relevance. In our current understanding, the RBC cm 's structure is patchy, with nanometer sized liquid ordered and disordered lipid, and peptide domains. At the same time, it is surprisingly soft, with bending rigidities κ of 2-4 kBT. This is in strong contrast to the current belief that a high concentration of cholesterol results in stiff membranes. This extreme softness is likely the result of an interaction between polyunsaturated lipids and cholesterol, which may also occur in other biological membranes. There is strong evidence in the literature that there is no length scale dependence of κ of whole RBCs.

Membranolytic Mechanism of Amphiphilic Antimicrobial β-Stranded [KL]n Peptides

Amphipathic peptides can act as antibiotics due to membrane permeabilization. KL peptides with the repetitive sequence [Lys-Leu]_n-NH₂ form amphipathic β-strands in the presence of lipid bilayers. As they are known to kill bacteria in a peculiar length-dependent manner, we suggest here several different functional models, all of which seem plausible, including a carpet mechanism, a β-barrel pore, a toroidal wormhole, and a β-helix. To resolve their genuine mechanism, the activity of KL peptides with lengths from 6–26 amino acids (plus some inverted LK analogues) was systematically tested against bacteria and erythrocytes. Vesicle leakage assays served to correlate bilayer thickness and peptide length and to examine the role of membrane curvature and putative pore diameter. KL peptides with 10–12 amino acids showed the best therapeutic potential, i.e., high antimicrobial activity and low hemolytic side effects. Mechanistically, this particular window of an optimum β-strand length around 4 nm (11 amino acids × 3.7 Å) would match the typical thickness of a lipid bilayer, implying the formation of a transmembrane pore. Solid-state ¹⁵N- and ¹⁹F-NMR structure analysis, however, showed that the KL backbone lies flat on the membrane surface under all conditions. We can thus refute any of the pore models and conclude that the KL peptides rather disrupt membranes by a carpet mechanism. The intriguing length-dependent optimum in activity can be fully explained by two counteracting effects, i.e., membrane binding versus amyloid formation. Very short KL peptides are inactive, because they are unable to bind to the lipid bilayer as flexible β-strands, whereas very long peptides are inactive due to vigorous pre-aggregation into β-sheets in solution.

Lipid/water interface of galactolipid bilayers in different lyotropic liquid-crystalline phases

In this study, carried out using computational methods, the organisation of the lipid/water interface of bilayers composed of galactolipids with both α-linolenoyl acyl chains is analysed and compared in three different lyotropic liquid-crystalline phases. These systems include the monogalactosyldiglyceride (MGDG) and digalactosyldiglyceride (DGDG) bilayers in the lamellar phase, the MGDG double bilayer during stalk phase formation and the inverse hexagonal MGDG phase. For each system, lipid-water and direct and water-mediated lipid-lipid interactions between the lipids of one bilayer leaflet and those of two apposing leaflets at the onset of new phase (stalk) formation, are identified. A network of interactions between DGDG molecules and its topological properties are derived and compared to those for the MGDG bilayer.

The bending rigidity of the red blood cell cytoplasmic membrane

An important mechanical property of cells is the membrane bending modulus, κ. In the case of red blood cells (RBCs) there is a composite membrane consisting of a cytoplasmic membrane and an underlying spectrin network. Literature values of κ are puzzling, as they are reported over a wide range, from 5 kBT to 230 kBT. To disentangle the contribution of the cytoplasmic membrane from the spectrin network, we investigated the bending of red blood cell cytoplasmic membranes (RBCcm) in the absence of spectrin and adenosine triphosphate (ATP). We used a combination of X-ray diffuse scattering (XDS), neutron spin-echo (NSE) spectrometry and Molecular Dynamics (MD) simulations. Our results indicate values of κ of order 4 kBT to 6 kBT, relatively small compared to literature values for most single component lipid bilayers. We suggest two ways this relative softness might confer biological advantage.

Molecular Dynamics Simulation of Transport Mechanism of Graphene Quantum Dots Through Different Cell Membranes

Exploring the mechanisms underlying the permeation of graphene quantum dots (GQDs) through different cell membranes is key for the practical application of GQDs in medicine. Here, the permeation process of GQDs through different lipid membranes was evaluated using molecular dynamics (MD) simulations. Our results showed that GQDs can easily permeate into 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE) lipid membranes with low phospholipid molecule densities but cannot permeate into 1-palmitoyl-2-oleoyl phosphatidylethanolamine (POPE) lipid membranes with high phospholipid densities. Free energy calculation showed that a high-energy barrier exists on the surface of the POPE lipid membrane, which prevents GQDs from entering the cell membrane interior. Further analysis of the POPE membrane structure showed that sparsely arranged phospholipid molecules of the low-density lipid membrane facilitated the entry of GQDs into the interior of the membrane, compared to compactly arranged molecules in the high-density lipid membrane. Our simulation study provides new insights into the transmembrane transport of GQDs.

Lipid chain-driven interaction of a lipidated Src-family kinase Lyn with the bilayer membrane

N-Myristoylation is a process of ubiquitous protein modification, which promotes the interaction of lipidated proteins on cell surfaces, in conjunction with reversible S-palmitoylation. We report the cooperative lipid-lipid interaction of two acyl chains of proteins, which increases the protein-membrane interaction and facilitates selective targeting of membranes containing anionic lipids. Lyn is a member of the Src family kinases distributed on the membrane surface by N-myristoyl and neighbouring S-palmitoyl chain anchors at the unique N-terminus domain. We prepared N-terminal short segments of lipidated Lyn to investigate the behaviour of each acyl chain in the lipid composition-dependent membrane interaction by solid-state nuclear magnetic resonance (NMR) analysis. Solid-state ³¹P-NMR studies revealed that S-palmitoylation of N-myristoylated Lyn peptides increased the interaction between peptides and phospholipid head groups, particularly with the anionic phosphatidylserine-containing bilayers. The solid-state ²H-NMR of Lyn peptides with a perdeutero N-myristoyl chain indicated an increase (0.6-0.8 Å) in the extent of the N-myristoyl chain in the presence of nearby S-palmitoyl chains, probably through the interaction via the acyl chains. The cooperative hydrocarbon chain interaction of the two acyl chains of Lyn increased membrane binding by extending the hydrocarbon chains deeper into the membrane interior, thereby promoting the peptide-membrane surface interaction between the cationic peptide side chains and the anionic lipid head groups. This lipid-driven mechanism by S-palmitoylation promotes the partition of the lipidated proteins to the cytoplasmic surface of the cell membranes and may be involved in recruiting Lyn at the signalling domains rich in anionic lipids.

Interaction of Phthalates with Lipid Bilayer Membranes

Phthalates are esters of phthalic acid, widely used as additives in the manufacture of plastics. They are not covalently linked to polymer chains and can easily leach out, disperse in the environment, and get into contact with living organisms. Several short chain phthalates are classified as endocrine disruptors or hormonal active agents, and have also been reported to promote various kinds of cancer. However, the biological effects of longer chain analogues are less well known. Moreover, little is known on the permeation of phthalates and their metabolites through biological membranes and on their effects on the physical properties of membranes. Here we explore the interaction of a group of phthalates and their main metabolites with model biological membranes. We focus on three industrially relevant phthalates, with acyl chains of different sizes, and their monoester metabolites. We use molecular dynamics simulations to predict the distribution in model membranes, as well as permeabilities and effects on the structural, dynamic, and elastic properties of the membranes. We find that alterations of membrane properties are significant and only weakly affected by the size of acyl chains, suggesting that modifications of molecular size may not be sufficient to reduce the impact of this class of molecules on the environment and health.

Phosphatidylserine and Phosphatidylethanolamine Asymmetry Have a Negligible Effect on the Global Structure, Dynamics, and Interactions of the KRAS Lipid Anchor

The intrinsically disordered C-terminus of the prominent oncogenic protein KRAS-4B (KRAS) selectively interacts and clusters with phosphatidylserine (PS) lipids in the plasma membrane (PM). This 11-residue segment, called tK, contains a polybasic domain (PBD) of six contiguous lysine residues and a farnesylated cysteine. Previous molecular dynamics (MD) simulation studies of tK in phosphatidylcholine (PC)/PS bilayers have suggested that backbone conformational dynamics modulate tK-PS interactions. These simulations have been conducted in symmetric membranes whereas the PM is compositionally asymmetric, with the inner leaflet, where KRAS is localized, being enriched with PS and phosphatidylethanolamine (PE) lipids. To examine if bilayer asymmetry affects tK conformational dynamics and interaction with lipids, we conducted two 10 μs long MD simulations of tK bound to a PC/PS and a PC/PS/PE bilayer in which the PS and PE lipids are distributed in one leaflet. We found that, first, these compositional asymmetries caused differences in acyl chain dynamics between leaflets, but the equilibrium structural and dynamic properties of the two asymmetric bilayers are similar; second, in both systems tK is highly dynamic and samples at least two distinct conformational states; third, PS-tK hydrogen-bonding interactions vary with peptide backbone conformations, and lysine side chains in the PBD predominantly interact with the serine oxygens of PS. These results are in good agreement with previous observations of tK in symmetric membranes. The effects of POPS asymmetry or the presence of POPE on tK are limited to modulating the relative contribution of individual side chains to interactions with lipids and redistributing conformational substates. Additional observations include the larger flexibility of tK in the current simulations, which we attribute to the longer duration of the simulations and the use of the CHARMM36m force field, which more accurately models intrinsically disordered peptides such as tK.

Impact of Surface Polarity on Lipid Assembly under Spatial Confinement

Molecular dynamics (MD) simulations in the MARTINI model are used to study the assembly of 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) molecules under spatial confinement, such as during solvent evaporation from ultrasmall (femtoliter quantity) droplets. The impact of surface polarity on molecular assembly is discussed in detail. To the best of our knowledge, this work represents the first of its kind. Our results reveal that solvent evaporation gives rise to the formation of well-defined stacks of lipid bilayers in a smectic alignment. These smectic mesophases form on both polar and nonpolar surfaces but with a notable distinction. On polar surfaces, the director of the stack is oriented perpendicular to the support surface. By contrast, the stacks orient at an angle on the nonpolar surfaces. The packing of head groups on surfaces and lipid molecular mobility exhibits significant differences as surface polarity changes. The role of glycerol in the assembly and stability is also revealed. The insights revealed from the simulation have a significant impact on additive manufacturing, biomaterials, model membranes, and engineering protocells. For example, POPC assemblies via evaporation of ultrasmall droplets were produced and characterized. The trends compare well with the bilayer stack models. The surface polarity influences the local morphology and structures at the interfaces, which could be rationalized via the molecule-surface interactions observed from simulations.

Development of OPLS-AA/M Parameters for Simulations of G Protein-Coupled Receptors and Other Membrane Proteins

G protein-coupled receptors (GPCRs) and other membrane proteins are valuable drug targets, and their dynamic nature makes them attractive systems for study with molecular dynamics (MD) simulations and free energy approaches. Here, we report the development, implementation, and validation of OPLS-AA/M force field parameters to enable simulations of these systems. These efforts include the introduction of post-translational modifications including lipidations and phosphorylation. We also modify previously reported parameters for lipids to be more consistent with the OPLS-AA force field standard and extend their coverage. These new parameters are validated on a variety of test systems, with the results compared to high-level quantum mechanics calculations, experimental data, and simulations with other force fields. The results demonstrate that the new parameters reliably reproduce the behavior of membrane protein systems.

Effect of encapsulated protein on the dynamics of lipid sponge phase: a neutron spin echo and molecular dynamics simulation study

Lipid membranes are highly mobile systems with hierarchical, time and length scale dependent, collective motions including thickness fluctuations, undulations, and topological membrane changes, which play an important role in membrane interactions. In this work we have characterised the effect of encapsulating two industrially important enzymes, β-galactosidase and aspartic protease, in lipid sponge phase nanoparticles on the dynamics of the lipid membrane using neutron spin echo (NSE) spectroscopy and molecular dynamics (MD) simulations. From NSE, reduced membrane dynamics were observed upon enzyme encapsulation, which were dependent on the enzyme concentration and type. By fitting the intermediate scattering functions (ISFs) with a modified Zilman and Granek model including nanoparticle diffusion, an increase in membrane bending rigidity was observed, with a larger effect for β-galactosidase than aspartic protease at the same concentration. MD simulations for the system with and without aspartic protease showed that the lipids relax more slowly in the system with protein due to the replacement of the lipid carbonyl-water hydrogen bonds with lipid-protein hydrogen bonds. This indicates that the most likely cause of the increase in membrane rigidity observed in the NSE measurements was dehydration of the lipid head groups. The dynamics of the protein itself were also studied, which showed a stable secondary structure of protein over the simulation, indicating no unfolding events occurred.

Probing the Link between Pancratistatin and Mitochondrial Apoptosis through Changes in the Membrane Dynamics on the Nanoscale

Pancratistatin (PST) is a natural antiviral alkaloid that has demonstrated specificity toward cancerous cells and explicitly targets the mitochondria. PST initiates apoptosis while leaving healthy, noncancerous cells unscathed. However, the manner by which PST induces apoptosis remains elusive and impedes the advancement of PST as a natural anticancer therapeutic agent. Herein, we use neutron spin-echo (NSE) spectroscopy, molecular dynamics (MD) simulations, and supporting small angle scattering techniques to study PST's effect on membrane dynamics using biologically representative model membranes. Our data suggests that PST stiffens the inner mitochondrial membrane (IMM) by being preferentially associated with cardiolipin, which would lead to the relocation and release of cytochrome c. Second, PST has an ordering effect on the lipids and disrupts their distribution within the IMM, which would interfere with the maintenance and functionality of the active forms of proteins in the electron transport chain. These previously unreported findings implicate PST's effect on mitochondrial apoptosis.

Physicochemical Properties Altered by the Tail Group of Lipid Membranes Influence Huntingtin Aggregation and Lipid Binding

Huntington's disease is a neurodegenerative disorder caused by an expanded polyglutamine (polyQ) domain within the huntingtin protein (htt) that initiates toxic protein aggregation. Htt directly interacts with membranes, influencing aggregation and spurring membrane abnormalities. These interactions are facilitated by the 17 N-terminal residues (Nt17) that form an amphipathic α-helix implicated in both lipid binding and aggregation. Here, the impact of unsaturation in phospholipid tails on htt-lipid interaction and htt aggregation was determined. There was no correlation between the degree of htt-lipid complexation and the degree of htt aggregation in the presence of each lipid system, indicating that lipid systems with different properties uniquely alter the membrane-mediated aggregation mechanisms. Also, the association between Nt17 and membrane surfaces is determined by complementarity between hydrophobic residues and membrane defects and how easily the peptide can partition into the bilayer. Our results provide critical insights into how membrane physical properties influence downstream htt aggregation.

Distinct lipid membrane-mediated pathways of Tau assembly revealed by single-molecule analysis

The conversion of intrinsically disordered Tau to highly ordered amyloid aggregates is associated with a wide range of neurodegenerative diseases termed tauopathies. The presence of lipid bilayer membranes is a critical factor that accelerates the abnormal aggregation of Tau protein. However, the lipid membrane-induced conformational changes of Tau and the mechanism for the accelerated fibrillation remain elusive. In this study, single-molecule Förster resonance energy transfer (smFRET) and fluorescence correlation spectroscopy (FCS) were applied to detect the conformational changes and intermolecular interactions of full-length Tau in the presence of different concentrations of 1,2-dimyristoyl-sn-glycero-3-phosphatidylserine (DMPS) vesicles. The results show that the conformation of Tau becomes expanded with opening of the N-terminal and C-terminal domains of Tau upon binding to DMPS. At low DMPS concentrations, Tau forms oligomers with a partially extended conformation which facilitates the amyloid fibrillization process. At high DMPS concentrations, Tau monomer binds to lipid membranes in a fully expanded conformation at low density thus inhibiting intermolecular aggregation. Our study reveals the underlying mechanisms by which lipid membranes influence amyloid formation of Tau, providing a foundation for further understanding of the pathogenesis and physiology of the interplay between Tau protein and lipid membranes.

A vesicle microrheometer for high-throughput viscosity measurements of lipid and polymer membranes

Viscosity is a key property of cell membranes that controls mobility of embedded proteins and membrane remodeling. Measuring it is challenging because existing approaches involve complex experimental designs and/or models, and the applicability of some methods is limited to specific systems and membrane compositions. As a result there is scarcity of systematic data, and the reported values for membrane viscosity vary by orders of magnitude for the same system. Here, we show how viscosity of membranes can be easily obtained from the transient deformation of giant unilamellar vesicles. The approach enables a noninvasive, probe-independent, and high-throughput measurement of the viscosity of membranes made of lipids or polymers with a wide range of compositions and phase state. Using this novel method, we have collected a significant amount of data that provides insights into the relation between membrane viscosity, composition, and structure.

Visualizing Lipid Dynamics Role in Infrared Neural Stimulation using Stimulated Raman Scattering

Infrared neural stimulation, or INS, uses pulsed infrared light to yield label-free neural stimulation with broad experimental and translational utility. Despite its robust demonstration, INS's mechanistic and biophysical underpinnings have been the subject of debate for more than a decade. The role of lipid membrane thermodynamics appears to play an important role in how fast IR-mediated heating nonspecifically drives action potential generation. Direct observation of lipid membrane dynamics during INS remains to be shown in a live neural model system. We used hyperspectral stimulated Raman scattering (hsSRS) microscopy to study biochemical signatures of high-speed vibrational dynamics underlying INS in a live neural cell culture model. Findings suggest that lipid bilayer structural changes are occurring during INS in vitro in NG108-15 neuroglioma cells. Lipid-specific signatures of cell SRS spectra varied with stimulation energy and radiant exposure. Spectroscopic observations agree with high-speed ratiometric fluorescence imaging of a conventional lipophilic membrane structure reporter, di-4-ANNEPS. Overall, the presented findings support the hypothesis that INS causes changes in the lipid membrane of neural cells by changing lipid membrane packing order. Furthermore, this work highlights the potential of hsSRS as a method to study biophysical and biochemical dynamics safely in live cells.

Taking the Monte-Carlo gamble: How not to buckle under the pressure!

Consistent buckling distortions of a large membrane patch (200 × 200 Å) are observed during molecular dynamics (MD) simulations using the Monte-Carlo (MC) barostat in combination with a hard Lennard-Jones (LJ) cutoff. The buckling behavior is independent of both the simulation engine and the force field but requires the MC barostat-hard LJ cutoff combination. Similar simulations of a smaller patch (90 × 90 Å) do not show buckling, but do show a small, systematic reduction in the surface area accompanied by ~1 Å thickening suggestive of compression. We show that a mismatch in the way potentials and forces are handled in the dynamical equations versus the MC barostat results in a compressive load on the membrane. Moreover, a straightforward application of elasticity theory reveals that a minimal compression of the linear dimensions of the membrane, inversely proportional to the edge length, is required for buckling, explaining this differential behavior. We recommend always using LJ force or potential-switching when the MC barostat is employed to avoid undesirable membrane deformations.

Deciphering Ethanol-Driven Swelling, Rupturing, Aggregation, and Fusion of Lipid Vesicles Using Coarse-Grained Molecular Dynamics Simulations

Traditionally, liquid ethanol is known to enhance the permeability of lipid membranes and causes vesicle aggregation and fusion. However, how the amphiphilic ethanol molecules perturb the lipid vesicles to facilitate their aggregation or fusion has not been addressed at any level of molecular simulations. Herein, not only have we developed a coarse-grained (CG) model for liquid ethanol, its aqueous mixture, and hydrated lipid membranes for molecular dynamics (MD) simulations, but also utilized it to delineate the aggregation and fusion of lipid vesicles using CG-MD simulations with multimillion particles. We have systematically parametrized the force-field for pure ethanol and its interactions with hydrated POPC and POPE model lipid membranes. In this process, we have successfully reproduced the bulk ethanol structure and concentration-dependent density of aqueous ethanol. To quantify the interaction of ethanol with lipid membranes, we have reproduced the transfer free energy of the ethanol molecule across the hydrated bilayers, and the concentration-dependent distribution of ethanol molecules across the lipid bilayers. After having acceptable force-field parameters for ethanol-membrane interactions, we have checked the effect of ethanol toward the vesicles comprising POPC lipids. We observe a rapid increase in the size of the POPC lipid vesicles with increasing amounts of ethanol up to 30 mol %. We unambiguously observe swelling and decrease in the thickness of the POPC vesicles with increasing amounts of ethanol up to 30 mol %, beyond which the vesicles begin to lose their integrity and rupture at higher mol % of ethanol. The fusion study of two vesicles demonstrates that fused vesicles can be obtained from 20 to 30 mol % of ethanol provided that they are brought closer than a critical distance at a particular mol %. The multivesicle simulations show that along with the increase in the sizes of vesicles the propensity of vesicle aggregation increases as the mol % of ethanol increases.

Comparative Molecular Dynamics Simulation Studies of Realistic Eukaryotic, Prokaryotic, and Archaeal Membranes

We present a comparative all-atom molecular dynamics simulation study of 18 biomembrane systems with lipid compositions corresponding to eukaryotic, bacterial, and archaebacterial membranes together with three single-component lipid bilayers. A total of 105 lipid types used in this study include diverse sterols and glycerol-based lipids with acyl chains of various lengths, unsaturation degrees, and branched or cyclic moieties. Our comparative analysis provides deeper insight into the influences of sterols and lipid unsaturation on the structural and mechanical properties of these biomembranes, including water permeation into the membrane hydrocarbon core. For sterol-containing membranes, sterol fraction is correlated with the membrane thickness, the area compressibility modulus, and lipid order but anticorrelated with the area per lipid and sterol tilt angles. Similarly, for all 18 biomembranes, lipid order is correlated with the membrane thickness and area compressibility modulus. Sterols and lipid unsaturation produce opposite effects on membrane thickness, but only sterols influence water permeation into the membrane. All membrane systems are accessible for public use in CHARMM-GUI Archive. They can be used as templates to expedite future modeling of realistic cell membranes with transmembrane and peripheral membrane proteins to study their structure, dynamics, molecular interactions, and function in a nativelike membrane environment.

Can Urea and Trimethylamine-N-oxide Prevent the Pressure-Induced Phase Transition of Lipid Membrane?

Organisms dwelling in ocean trenches are exposed to the high hydrostatic pressure of ocean water. Increasing pressure can alter the membrane packing density and fluidity and trigger the fluid-to-gel phase transition. To combat environmental stress, the organisms synthesize small polar solutes, which are known as osmolytes. Urea and trimethylamine-N-oxide (TMAO) are two such solutes found in deep-sea creatures. While TMAO stabilizes protein, urea induces protein denaturation. These solutes strongly influence the packing density and membrane fluidity of the lipid bilayer at different conditions. But can these solutes affect the pressure-induced phase transition of the lipid membrane? In the present work, we have studied the effect of these two solutes on pressure-induced fluid-to-gel phase transition based on the all-atom molecular dynamics (MD) simulation approach. A high-pressure-stimulated fluid-to-gel phase transition of the membrane is seen at 800 bar, which is consistent with previous experiments. We have also observed that in the low-pressure region (1-400 bar), urea slightly increases the membrane fluidity where TMAO decreases the same. However, the phase transition pressure remains almost unchanged on the addition of urea while TMAO shifts the phase transition toward a lower pressure. We have found that the hydrogen (H)-bond interaction between lipid and urea plays an important role in preserving the fluidity of the membrane in the low-pressure zone. However, at a higher pressure, both water and urea are excluded from the membrane surface. TMAO is also excluded from the interfacial region of the membrane at all pressures. Exclusion from the membrane surface further triggers the phase transition of the lipid membrane from the fluid to gel phase at a high pressure.

Lipid Membrane Binding and Cell Protection Efficacy of Poly(1,2-butylene oxide)-b-poly(ethylene oxide) Copolymers

Poloxamers consisting of poly(ethylene oxide) (PEO) and poly(propylene oxide) segments can protect cell membranes against various forms of stress. We investigated the role of the hydrophobic block chemistry on polymer/membrane binding and cell membrane protection by comparing a series of poly(butylene oxide)-b-PEO (PBO-b-PEO) copolymers to poloxamer analogues, using a combination of pulsed-field-gradient (PFG) NMR experiments and a lactate dehydrogenase (LDH) cell assay. We found that the more hydrophobic PBO-b-PEO copolymers bound more significantly to model liposomes composed of 1-palmitol-2-oleoyl-glycero-3-phosphocholine (POPC) compared to poly(propylene oxide) (PPO)/PEO copolymers. However, both classes of polymers performed similarly when compared by an LDH assay. These results present an important comparison between polymers with similar structures but with different binding affinities. They also provide mechanistic insight as enhanced polymer/lipid membrane binding did not directly translate to increased cell protection in the LDH assay, and therefore, additional factors need to be considered when trying to achieve greater membrane protection efficacy.

Lipid21: Complex Lipid Membrane Simulations with AMBER

We extend the modular AMBER lipid force field to include anionic lipids, polyunsaturated fatty acid (PUFA) lipids, and sphingomyelin, allowing the simulation of realistic cell membrane lipid compositions, including raft-like domains. Head group torsion parameters are revised, resulting in improved agreement with NMR order parameters, and hydrocarbon chain parameters are updated, providing a better match with phase transition temperature. Extensive validation runs (0.9 μs per lipid type) show good agreement with experimental measurements. Furthermore, the simulation of raft-like bilayers demonstrates the perturbing effect of increasing PUFA concentrations on cholesterol molecules. The force field derivation is consistent with the AMBER philosophy, meaning it can be easily mixed with protein, small molecule, nucleic acid, and carbohydrate force fields.

Specificity of Loxosceles α clade phospholipase D enzymes for choline-containing lipids: Role of a conserved aromatic cage

Spider venom GDPD-like phospholipases D (SicTox) have been identified to be one of the major toxins in recluse spider venom. They are divided into two major clades: the α clade and the β clade. Most α clade toxins present high activity against lipids with choline head groups such as sphingomyelin, while activities in β clade toxins vary and include preference for substrates containing ethanolamine headgroups (Sicarius terrosus, St_βIB1). A structural comparison of available structures of phospholipases D (PLDs) reveals a conserved aromatic cage in the α clade. To test the potential influence of the aromatic cage on membrane-lipid specificity we performed molecular dynamics (MD) simulations of the binding of several PLDs onto lipid bilayers containing choline headgroups; two SicTox from the α clade, Loxosceles intermedia αIA1 (Li_αIA) and Loxosceles laeta αIII1 (Ll_αIII1), and one from the β clade, St_βIB1. The simulation results reveal that the aromatic cage captures a choline-headgroup and suggest that the cage plays a major role in lipid specificity. We also simulated an engineered St_βIB1, where we introduced the aromatic cage, and this led to binding with choline-containing lipids. Moreover, a multiple sequence alignment revealed the conservation of the aromatic cage among the α clade PLDs. Here, we confirmed that the i-face of α and β clade PLDs is involved in their binding to choline and ethanolamine-containing bilayers, respectively. Furthermore, our results suggest a major role in choline lipid recognition of the aromatic cage of the α clade PLDs. The MD simulation results are supported by in vitro liposome binding assay experiments.

Envenomation following bites from recluse spiders (Loxosceles) causes loxoscelism, a necrotic tissue breakdown in mammals, and leads to skin degeneration and systemic reactions in the worst case. Recluse spiders belong to the Sicariidae family which also includes six-eyed sand spiders in the genera Sicarius and Hexopthalma. While sicariid spiders are found natively on all continents except Australia, treatments of loxoscelism are typically antibody based and available in some regions of the Americas. Sphingomyelinase D/phospholipase D enzymes are one of the major toxins in venom of sicariid spiders, and have been divided in two clades called α and β. The activity of α and β clades toxins differs; most α clade toxins present high activity against lipids with choline headgroups (-N (CH₃)₃⁺) such as sphingomyelin, while activities in β clade toxins vary and include preference for substrates containing ethanolamine headgroups (-NH₃⁺). When comparing the structures of two α clade toxins and one β clade toxin, we noticed the presence in the α clade toxins only of a cage consisting of three aromatic amino acids. In this work we used numerical molecular simulations to probe the role of this cage in the preference of α clade toxins for choline head groups over ethanolamine head groups.

Effects of cholesterol on the size distribution and bending modulus of lipid vesicles

The influence of cholesterol fraction in the membranes of giant unilamellar vesicles (GUVs) on their size distributions and bending moduli has been investigated. The membranes of GUVs were synthesized by a mixture of two elements: electrically neutral lipid 1, 2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and cholesterol and also a mixture of three elements: electrically charged lipid 1,2-dioleoyl-sn-glycero-3-phospho-(1′-rac-glycerol) (DOPG), DOPC and cholesterol. The size distributions of GUVs have been presented by a set of histograms. The classical lognormal distribution is well fitted to the histograms, from where the average size of vesicle is obtained. The increase of cholesterol content in the membranes of GUVs increases the average size of vesicles in the population. Using the framework of Helmholtz free energy of the system, the theory developed by us is extended to explain the experimental results. The theory determines the influence of cholesterol on the bending modulus of membranes from the fitting of the proper histograms. The increase of cholesterol in GUVs increases both the average size of vesicles in population and the bending modulus of membranes.

Reconstitution of Membrane-associated Components of a G-protein Signaling Pathway on Membrane-coated Nanoparticles (Lipobeads)

G-protein coupled signaling pathways are organized into multi-protein complexes called signalosomes that are located within and on cellular membranes. We describe the use of silica nanoparticles coated with a unilamellar phospholipid bilayer (lipobeads) to reconstitute the activated photoreceptor G-protein α-subunit (Gtα*) with its cognate effector (phosphodiesterase-6; PDE6) for biochemical and structural studies of the activation mechanism regulating this GPCR signaling pathway. Lipobeads are prepared by resuspending dried-down phospholipid mixtures with monodisperse 70 nm silica particles, followed by extrusion through a 100 nm membrane filter. This uniform and supported liposomal preparation is easily sedimented, permitting the separation of soluble from membrane-associated proteins. Upon loading lipobeads with Gtα* and PDE6, we find that activation of PDE6 catalysis by Gtα* occurs much more efficiently than in the absence of membranes. Chemical cross-linking of membrane-confined proteins allows detection of changes in protein-protein interactions, resulting from G-protein activation of PDE6. The advantages of using lipobeads over partially purified membrane preparations or traditional liposomal preparations are generally applicable to the study of other membrane-confined signal transduction pathways.

Helical structure motifs made searchable for functional peptide design

The systematic design of functional peptides has technological and therapeutic applications. However, there is a need for pattern-based search engines that help locate desired functional motifs in primary sequences regardless of their evolutionary conservation. Existing databases such as The Protein Secondary Structure database (PSS) no longer serves the community, while the Dictionary of Protein Secondary Structure (DSSP) annotates the secondary structures when tertiary structures of proteins are provided. Here, we extract 1.7 million helices from the PDB and compile them into a database (Therapeutic Peptide Design database; TP-DB) that allows queries of compounded patterns to facilitate the identification of sequence motifs of helical structures. We show how TP-DB helps us identify a known purification-tag-specific antibody that can be repurposed into a diagnostic kit for Helicobacter pylori. We also show how the database can be used to design a new antimicrobial peptide that shows better Candida albicans clearance and lower hemolysis than its template homologs. Finally, we demonstrate how TP-DB can suggest point mutations in helical peptide blockers to prevent a targeted tumorigenic protein-protein interaction. TP-DB is made available at http://dyn.life.nthu.edu.tw/design/ .

In Silico-Based Experiments on Mechanistic Interactions between Several Intestinal Permeation Enhancers with a Lipid Bilayer Model

Oral administration of drugs is generally considered convenient and patient-friendly. However, oral administration of biological drugs exhibits low oral bioavailability (BA) due to enzymatic degradation and low intestinal absorption. A possible approach to circumvent the low BA of oral peptide drugs is to coformulate the drugs with permeation enhancers (PEs). PEs have been studied since the 1960s and are molecules that enhance the absorption of hydrophilic molecules with low permeability over the gastrointestinal epithelium. In this study, we investigated the impact of six PEs on the structural properties of a model membrane using molecular dynamics (MD) simulations. The PEs included were the sodium salts of the medium chain fatty acids laurate, caprate, and caprylate and the caprylate derivative SNAC─all with a negative charge─and neutral caprate and neutral sucrose monolaurate. Our results indicated that the PEs, once incorporated into the membrane, could induce membrane leakiness in a concentration-dependent manner. Our simulations suggest that a PE concentration of at least 70-100 mM is needed to strongly affect transcellular permeability. The increased aggregation propensity seen for neutral PEs might provide a molecular-level mechanism for the membrane disruptions seen at higher concentrations in vivo. The ability for neutral PEs to flip-flop across the lipid bilayer is also suggestive of possible intracellular modes of action other than increasing membrane fluidity. Taken together, our results indicate that MD simulations are useful for gaining insights relevant to the design of oral dosage forms based around permeability enhancer molecules.