Computationally efficient prediction of area per lipid

Exploration of lipid bilayer mechanical properties using molecular dynamics simulation

Important biological structures known for their exceptional mechanical qualities, lipid bilayers are essential to many cellular functions. Fluidity, elasticity, permeability, stiffness, tensile strength, compressibility, shear viscosity, line tension, and curvature elasticity are some of the fundamental characteristics affecting their behavior. The purpose of this review is to examine these characteristics in more detail by molecular dynamics simulation, elucidating their importance and the elements that lead to their appearance in lipid bilayers. Comprehending the mechanical characteristics of lipid bilayers is critical for creating medications, drug delivery systems, and biomaterials that interact with biological membranes because it allows one to understand how these materials respond to different stresses and deformations. The influence of mechanical characteristics on important lipid bilayer properties is examined in this review. The mechanical properties of lipid bilayers were clarified through the use of molecular dynamics simulation analysis techniques, including bilayer thickness, stress-strain analysis, lipid bilayer area compressibility, membrane bending rigidity, and time- or ensemble-averaged the area per lipid evaluation. We explain the significance of molecular dynamics simulation analysis methods, providing important new information about the stability and dynamic behavior of the bilayer. In the end, we hope to use molecular dynamics simulation to provide a comprehensive understanding of the mechanical properties and behavior of lipid bilayers, laying the groundwork for further studies and applications. Taken together, careful investigation of these mechanical aspects deepens our understanding of the adaptive capacities and functional roles of lipid bilayers in biological environments.

Influence of phospholipid head and tail molecular structures on cell membrane mechanical response under tension

Biological cell membranes are primarily comprised of a diverse lipid bilayer with multiple phospholipid (lipid) types, each of which is comprised of a hydrophilic headgroup and two hydrophobic hydrocarbon tails. The lipid type determines the molecular structure of head and tail groups, which can affect membrane mechanics at nanoscale and subsequently cell viability under mechanical loading. Hence, using molecular dynamics simulations, the current study investigated seven membrane phospholipids and the effect of their structural differences on physical deformation, mechanoporation damage, and mechanical failure of the membranes under tension. The inspected phospholipids showed similar yield stresses and strains, as well as pore evolution and damage, but significantly different failure strains. In general, failure occurred at a lower strain for lipids with a larger equilibrium area per lipid. The obtained results suggest that larger headgroup structure, greater degree of unsaturation, and tail-length asymmetry influenced the phospholipids' ability to pack against each other, increased the fluidity and equilibrium area per lipid of the membrane, and resulted in lower failure strain. Overall, this study provides insights on how different phospholipid structures affect membrane physical responses at the molecular level and serves as a reference for future studies of more complex membrane systems with intricate biophysical properties.

Towards development of new antimalarial compounds through in silico and in vitro assays

Malaria is one of most widespread infectious disease in world. The antimalarial therapy presents a series of limitations, such as toxicity and the emergence of resistance, which makes the search for new drugs urgent. Thus, it becomes necessary to explore essential and exclusive therapeutic targets of the parasite to achieve selective inhibition. Enoyl-ACP reductase is an enzyme of the type II fatty acid biosynthetic pathway and is responsible for the rate-limiting step in the fatty acid elongation cycle. In this work, we use hierarchical virtual screening and drug repositioning strategies to prioritize compounds for phenotypic assays and molecular dynamics studies. The molecules were tested against chloroquine-resistant W2 strain of Plasmodium falciparum (EC50 between 330.05 and 13.92 µM). Nitrofurantoin was the best antimalarial activity at low micromolar range (EC50 = 13.92 µM). However, a hit compound against malaria must have a biological activity value below 1 µM. A large number of molecules present problems with permeability in biological membranes and reaching an effective concentration in their target's microenvironment. Nitrofurantoin derivatives with inclusions of groups which confer increased lipid solubility (methyl groups, halogens and substituted and unsubstituted aromatic rings) have been proposed. These derivatives were pulled through the lipid bilayer in molecular dynamics simulations. Molecules 14, 18 and 21 presented lower free energy values than nitrofurantoin when crossing the lipid bilayer.

Controlling spatial distribution of functional lipids in a supported lipid bilayer prepared from vesicles

Conjugating biomolecules, such as antibodies, to bioconjugate moieties on lipid surfaces is a powerful tool for engineering the surface of diverse biomaterials, including cells and nanoparticles. We developed supported lipid bilayers (SLBs) presenting well-defined spatial distributions of functional moieties as models for precisely engineered functional biomolecular-lipid surfaces. We used quartz crystal microbalance with dissipation (QCM-D) and atomic force microscopy (AFM) to determine how vesicles containing a mixture of 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC) and 1,2-distearoyl-sn-glycero-3-phosphoethanolamine-N-[azido(polyethylene glycol)-2000] (DSPE-PEG-N3) form SLBs as a function of the lipid phase transition temperature (Tm). Above the DPPC Tm, DPPC/DSPE-PEG-N3 vesicles form SLBs with functional azide moieties on SiO2 substrates via vesicle fusion. Below this Tm, DPPC/DSPE-PEG-N3 vesicles attach to SiO2 intact. Intact DPPC/DSPE-PEG-N3 vesicles on the SiO2 surfaces fuse and rupture to form SLBs when temperature is brought above the DPPC Tm. AFM studies show uniform and complete DPPC/DSPE-PEG-N3 SLB coverage of SiO2 surfaces for different DSPE-PEG-N3 concentrations. As the DSPE-PEG-N3 concentration increases from 0.01 to 6 mol%, the intermolecular spacing of DSPE-PEG-N3 in the SLBs decreases from 4.6 to 1.0 nm. The PEG moiety undergoes a mushroom to brush transition as DSPE-PEG-N3 concentration varies from 0.1 to 2.0 mol%. Via copper-free click reaction, IgG was conjugated to SLB surfaces with 4.6 nm or 1.3 nm inter-DSPE-PEG-N3 spacing. QCM-D and AFM data show; 1) uniform and complete IgG layers of similar mass and thickness on the two types of SLB; 2) a higher-viscosity/less rigid IgG layer on the SLB with 4.6 nm inter-DSPE-PEG-N3 spacing. Our studies provide a blueprint for SLBs modeling spatial control of functional macromolecules on lipid surfaces, including surfaces of lipid nanoparticles and cells.

How living cells are affected during the cold atmospheric pressure plasma treatment

When the electric discharge process is limited by high voltage electrodes shielding, the ionization measure would be controlled to less than one percent and the temperature to less than 37 °C even at atmospheric pressure, so-called cold atmospheric pressure plasma (CAP). CAP has been found to have profound medical applications in association with its reactive oxygen and nitrogen species (ROS/RNS). In this way that during plasma exposure, the subjected medium (e.g. cell cytoplasmic membrane in plasma therapy) interacts with ROS/RNS. Accordingly, a precise study of the mentioned interactions and their consequences on the cells' behavior changes, is necessary. The results lead to the reduction of possible risks and provide the opportunity of optimizing the efficacy of CAP before the development of CAP applications in the field of plasma medicine. In this report molecular dynamic (MD) simulation is used to investigate the mentioned interactions and a proper and compatible comparison with the experimental results is presented. Based on this, the effects of H2O2, NO and O2 on the living cell's membrane are investigated in biological conditions. Our results show that: i) The hydration of phospholipid polar heads would be enhanced associated with the H2O2 presence. ii) A new definition of the surface area assigned to each phospholipid (APL), more reliable and compatible with the physical expectations, is introduced. iii) The long-term behavior of NO and O2 is their penetration into the lipid bilayer and sometimes passing through the membrane into the cell. The latter would be an indication of internal cells' pathways activation leading to modification of cells' function.

The effects of adsorbed benzo(a)pyrene on dynamic behavior of polystyrene nanoplastics through phospholipid membrane: A molecular simulation study

Nanoplastics (NPs) are mainly generated from the decomposition of plastic waste and industrial production, which have attracted much attention due to the potential risk for humans. The ability of NPs to penetrate different biological barriers has been proved, but the understanding of molecular details is very limited, especially for organic pollutant-NP combinations. Here, we investigated the uptake process of polystyrene NPs (PSNPs) combined with benzo(a)pyrene (BAP) molecules by dipalmitoylphosphatidylcholine (DPPC) bilayers by molecular dynamics (MD) simulations. The results showed that the PSNPs can adsorb and accumulate BAP molecules in water phase and then carried BAP molecules to enter DPPC bilayers. At the same time, the adsorbed BAP promoted the penetration of PSNPs into DPPC bilayers effectively by hydrophobic effect. The process of BAP-PSNP combinations penetrating into DPPC bilayers can be summarized into four steps including adhesion on the DPPC bilayer surface, uptake by the DPPC bilayer, BAP molecules detached from the PSNPs, and the PSNPs depolymerized in the bilayer interior. Furthermore, the amount of adsorbed BAP on PSNPs affected the properties of DPPC bilayers directly, especially the fluidity of DPPC bilayers that determine the physiologic function. Obviously, the combined effect of PSNPs and BAP enhanced the cytotoxicity. This work not only presented a vivid transmembrane process of BAP-PSNP combinations and revealed the nature of the effects of adsorbed benzo(a)pyrene on the dynamic behavior of polystyrene nanoplastics through phospholipid membrane, but also provide some necessary information of the potential damage for organic pollutant-nanoplastic combinations on human health at a molecular level.

Study of trioleoylglycerol two-layer and adiposome cross-section mimicking four-layer systems through atomic-level simulations

Adiposomes are artificially prepared lipid droplet (LD)-mimetic structures, which, unlike LDs, do not harbor proteins. The dynamics of interaction between triacylglycerols (TAGs), drug molecule, and phospholipids in adiposomes is currently not well-established. Trioleoylglycerol (TOG) molecule was divided into three parts: two oleoyl tails and one 2-monooleoylglycerol (MOG). Forcefield parameters for two oleoyl tails were adopted from the AMBER18 repository while that of the MOG forcefield was taken from the literature. Charge correction was performed on the MOG forcefield before its utilization. After charge correction, the resulting TOG molecule had zero charge. TOG bilayer (2L) and tetralayer (4L) systems were prepared and simulated. TOG bilayer (2L) systems-modeled from two different initial conformations, the TOG3 conformation and the TOG2:1 conformation-showed that TOG2:1 conformation was more prevailing irrespective of the starting conformation and was subsequently used in further simulations. The hydrated TOG 2L system showed TOG-water solution solubility of 0.051 mol L-1 which is near experimental values. This validated the correct parameterization of the TOG molecule. The simulations of 4L systems showed stable membrane behaviors toward the end of simulations. It was also observed that in the 4L system, the TOG molecules showed the formation of micelles with the drug molecule. Almost six TOGs remained continuously in contact with the drug molecule throughout the simulation. The availability of charge-corrected TOG parameterization is expected to equip future studies with a framework for molecular dynamics simulations of adiposomes and/or LDs at the atomic level.

Biophysical Impact of Lipid A Modification Caused by Mobile Colistin Resistance Gene on Bacterial Outer Membranes

Expression of mobile colistin resistance gene mcr-1 results in the addition of phosphoethanolamine (pEtN) to the lipid A headgroup in the bacterial outer membrane (OM) of Gram-negative bacteria, increasing the resistance to the last-line polymyxins. However, the potential biological consequences of such modification remain unclear. Using coarse-grained molecular simulations with quantitative lipidomics models, we discovered pEtN modification of the lipid A headgroup caused substantial changes to the morphology and physicochemical properties of the OM. Single-lipid level structural and energetic analyses revealed that this modification resulted in lipid A-pEtN adopting an abnormally twisted and slanted conformation with a closer packing state because of strengthened inter-lipid attraction. The consequent accumulation of lipid A-pEtN produced a negative curvature of the OM and altered the membrane's tension, fluidity, and rigidity. Our results provide a key mechanistic connection between mcr-1 expression and biophysical changes in the bacterial OM.

Influence of Cholesterol on the Orientation of the Farnesylated GTP-Bound KRas-4B Binding with Anionic Model Membranes

The Ras family of proteins is tethered to the inner leaflet of the cell membranes which plays an essential role in signal transduction pathways that promote cellular proliferation, survival, growth, and differentiation. KRas-4B, the most mutated Ras isoform in different cancers, has been under extensive study for more than two decades. Here we have focused our interest on the influence of cholesterol on the orientations that KRas-4B adopts with respect to the plane of the anionic model membranes. How cholesterol in the bilayer might modulate preferences for specific orientation states is far from clear. Herein, after analyzing data from in total 4000 ns-long molecular dynamics (MD) simulations for four KRas-4B systems, properties such as the area per lipid and thickness of the membrane as well as selected radial distribution functions, penetration of different moieties of KRas-4B, and internal conformational fluctuations of flexible moieties in KRas-4B have been calculated. It has been shown that high cholesterol content in the plasma membrane (PM) favors one orientation state (OS1), exposing the effector-binding loop for signal transduction in the cell from the atomic level. We confirm that high cholesterol in the PM helps KRas-4B mutant stay in its constitutively active state, which suggests that high cholesterol intake can increase mortality and may promote cancer progression for cancer patients. We propose that during the treatment of KRas-4B-related cancers, reducing the cholesterol level in the PM and sustaining cancer progression by controlling the plasma cholesterol intake might be taken into account in anti-cancer therapies.

Importance of Hydrophilic Groups on Modulating the Structural, Mechanical, and Interfacial Properties of Bilayers: A Comparative Molecular Dynamics Study of Phosphatidylcholine and Ion Pair Amphiphile Membranes

An ion pair amphiphile (IPA), a molecular complex composed of two oppositely charged amphiphiles, is a phospholipid mimic which differs from a phospholipid only in the hydrophilic compositions. Here, we utilized molecular dynamics (MD) simulations to compare the bilayer systems composed of phosphatidylcholines (PC) and alkyltrimethylammonium-alkylsulfate IPAs with various alkyl chain lengths. The membrane properties for both liquid-disordered (Ld) and gel (S) phase bilayers were examined via running simulations above and below the main transition temperatures. The electrostatic attraction between the IPA hydrophilic groups leads to a more ordered molecular packing within both S and Ld phase IPA membranes, as revealed by the molecular area, deuterium order parameter, and gauche conformation analyses. Furthermore, IPA bilayers possess a higher area compressibility modulus, molecular tilt modulus, and effective bending rigidity than PC systems. The variation of hydrophilic groups of IPA also leads to fewer hydrogen bonds on the membrane surface and smaller electrostatic potentials for the biomimetic bilayer. The non-covalently linked head groups of IPA further decouple alkyl tilting and surface water retention. The combined results reveal the importance of hydrophilic groups of amphiphiles on modulating the membrane properties, which also provides insights for designs of biomimetic membranes.

Effect of 5-trans Isomer of Arachidonic Acid on Model Liposomal Membranes Studied by a Combined Simulation and Experimental Approach

Unsaturated fatty acids are found in humans predominantly in the cis configuration. Fatty acids in the trans configuration are primarily the result of human processing (trans fats), but can also be formed endogenously by radical stress. The cis-trans isomerization of fatty acids by free radicals could be connected to several pathologies. Trans fats have been linked to an increased risk of coronary artery disease; however, the reasons for the resulting pathogenesis remain unclear. Here, we investigate the effect of a mono-trans isomer of arachidonic acid (C20:4-5trans, 8cis, 11cis, 14cis) produced by free radicals in physiological concentration on a model erythrocyte membrane using a combined experimental and theoretical approach. Molecular Dynamics (MD) simulations of two model lipid bilayers containing arachidonic acid and its 5-trans isomer in 3 mol% were carried out for this purpose. The 5-trans isomer formation in the phospholipids was catalyzed by HOCH₂CH₂S· radicals, generated from the corresponding thiol by γ-irradiation, in multilamellar vesicles of SAPC. Large unilamellar vesicles were made by the extrusion method (LUVET) as a biomimetic model for cis-trans isomerization. Atomic Force Microscopy and Dynamic Light Scattering were used to measure the average size, morphology, and the z-potential of the liposomes. Both results from MD simulations and experiments are in agreement and indicate that the two model membranes display different physicochemical properties in that the bilayers containing the trans fatty acids were more ordered and more rigid than those containing solely the cis arachidonic acid. Correspondingly, the average size of the liposomes containing trans isomers was smaller than the ones without.

Computer modelling studies of the bilayer/water interface

This review summarises high resolution studies on the interface of lamellar lipid bilayers composed of the most typical lipid molecules which constitute the lipid matrix of biomembranes. The presented results were obtained predominantly by computer modelling methods. Whenever possible, the results were compared with experimental results obtained for similar systems. The first and main section of the review is concerned with the bilayer-water interface and is divided into four subsections. The first describes the simplest case, where the interface consists only of lipid head groups and water molecules and focuses on interactions between the lipid heads and water molecules; the second describes the interface containing also mono- and divalent ions and concentrates on lipid-ion interactions; the third describes direct inter-lipid interactions. These three subsections are followed by a discussion on the network of direct and indirect inter-lipid interactions at the bilayer interface. The second section summarises recent computer simulation studies on the interactions of antibacterial membrane active compounds with various models of the bacterial outer membrane. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Force Field Development for Lipid Membrane Simulations

With the rapid development of computer power and wide availability of modelling software computer simulations of realistic models of lipid membranes, including their interactions with various molecular species, polypeptides and membrane proteins have become feasible for many research groups. The crucial issue of the reliability of such simulations is the quality of the force field, and many efforts, especially in the latest several years, have been devoted to parametrization and optimization of the force fields for biomembrane modelling. In this review, we give account of the recent development in this area, covering different classes of force fields, principles of the force field parametrization, comparison of the force fields, and their experimental validation. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.