Interaction of C60 fullerenes with asymmetric and curved lipid membranes: a molecular dynamics study

Hydrophobic Matching Dictates over the Linear Rule of Mixtures in Binary Lipid Membranes

Biological membranes feature heterogeneous mixtures of lipids with different head and tail characteristics. Their biophysical properties are dictated by the intimate interaction among different constituent lipids. Previous studies suggest that the membrane area-per-lipid (APL) deviates from the linear rule of mixtures (LRM) for binary lipid membranes, but the underlying mechanism remains elusive. Our molecular dynamics (MD) simulations of binary lipid membranes consisting of lipids with different tail characteristics reveal a competitive mechanism whereby lipids tend to deform each other to minimize the hydrophobic mismatch between their tails. Depending on the relative tail lengths and saturation levels, this may result in an either positive or negative deviation of APL from the LRM. As lipid packing plays an essential role in membrane fusion and peptide-membrane binding, our findings may help guide the selection of lipids for the effective rational design of nanoliposomes and membrane-targeting peptides.

Plasma membrane lipid bilayer is druggable: Selective delivery of gemcitabine-squalene nano-medicine to cancer cells

Up to now the lipid bilayers were rarely considered as targets in cancer therapy despite pronounced differences in lipid composition between plasma membranes of benign and malignant cells. In this study we demonstrate that the lipid bilayer of the plasma membrane is druggable and suitable for facilitating selective delivery of amphiphilic gemcitabine-squalene nanomedicines to cancer cells. Data from radioactive assays, fluorescent membrane probes and molecular dynamics simulations provide evidence of selective accumulation of gemcitabine-squalene in the plasma membranes with disrupted lipid asymmetry and its subsequent preferential uptake by malignant cells. This causes pronounced cytotoxicity on cancer cells in comparison to their benign counterparts originating from the same tissue.

Current Trends and Changes in Use of Membrane Molecular Dynamics Simulations within Academia and the Pharmaceutical Industry

There has been an almost exponential increase in the use of molecular dynamics simulations in basic research and industry over the last 5 years, with almost a doubling in the number of publications each year. Many of these are focused on neurological membranes, and biological membranes in general, applied to the medical industry. A smaller portion have utilized membrane simulations to answer more basic questions related to the function of specific proteins, chemicals or biological processes. This review covers some newer studies, alongside studies from the last two decades, to determine changes in the field. Some of these are basic, while others are more profound, such as multi-component embedded membrane machinery. It is clear that many facets of the discipline remain the same, while the focus on and uses of the technology are broadening in scope and utilization as a general research tool. Analysis of recent literature provides an overview of the current methodologies, covers some of the recent trends or advances and tries to make predictions of the overall path membrane molecular dynamics will follow in the coming years. In general, the overview presented is geared towards the general scientific community, who may wish to introduce the use of these methodologies in light of these changes, making molecular dynamic simulations more feasible for general scientific or medical research.

Amphiphilic anti-SARS-CoV-2 drug remdesivir incorporates into the lipid bilayer and nerve terminal membranes influencing excitatory and inhibitory neurotransmission

Remdesivir is a novel antiviral drug, which is active against the SARS-CoV-2 virus. Remdesivir is known to accumulate in the brain but it is not clear whether it influences the neurotransmission. Here we report diverse and pronounced effects of remdesivir on transportation and release of excitatory and inhibitory neurotransmitters in rat cortex nerve terminals (synaptosomes) in vitro. Direct incorporation of remdesivir molecules into the cellular membranes was shown by FTIR spectroscopy, planar phospholipid bilayer membranes and computational techniques. Remdesivir decreases depolarization-induced exocytotic release of L-[¹⁴C] glutamate and [³H] GABA, and also [³H] GABA uptake and extracellular level in synaptosomes in a dose-dependent manner. Fluorimetric studies confirmed remdesivir-induced impairment of exocytosis in nerve terminals and revealed a decrease in synaptic vesicle acidification. Our data suggest that remdesivir dosing during antiviral therapy should be precisely controlled to prevent possible neuromodulatory action at the presynaptic level. Further studies of neurotropic and membranotropic effects of remdesivir are necessary.

The interplay between surface-functionalized gold nanoparticles and negatively charged lipid vesicles

The comprehensive understanding of the interactions between gold nanoparticles (AuNPs) and phospholipid vesicles has important implications in various biomedical applications; however, this is not yet well understood. Here, coarse-grained molecular dynamics (CGMD) simulations were performed to study the interactions between functionalized AuNPs and negatively charged lipid vesicles, and the effects of the surface chemistry and surface charge density (SCD) of AuNPs were analyzed. It is revealed that AuNPs with different surface ligands adhere to the membrane surface (anionic AuNPs) or get into the vesicle bilayer (hydrophobic and cationic AuNPs). Due to the loose arrangement of lipid molecules, AuNPs penetrate curved vesicle membranes more easily than planar lipid bilayers. Cationic AuNPs present three different interaction modes with the vesicle, namely insertion, partial penetration and complete penetration, which are decided by the SCD difference. Both hydrophobic interaction and electrostatic interaction play crucial roles in the interplay between cationic AuNPs and lipid vesicles. For the cationic AuNP with a low SCD, it gets into the lipid bilayer without membrane damage through the hydrophobic interaction, and it is finally stabilized in the hydrophobic interior of the vesicle membrane in a thermodynamically stable "snorkeling" configuration. For the cationic AuNP with a high SCD, it crosses the vesicle membrane and gets into the vesicle core through a membrane pore induced by strong electrostatic interaction. In this process, the membrane structure is destroyed. These findings provide a molecular-level understanding of the interplay between AuNPs and lipid vesicles, which may further expand the application of functional AuNPs in modern biomedicine.

Curvature increases permeability of the plasma membrane for ions, water and the anti-cancer drugs cisplatin and gemcitabine

In this work the permeability of a model asymmetric plasma membrane, for ions, water and the anti-cancer drugs cisplatin and gemcitabine is studied by means of all-atom molecular dynamics simulations. It is shown for the first time that permeability of the highly curved membrane increases from one to three orders of magnitude upon membrane bending depending on the compound and the sign of curvature. Our results suggest that the membrane curvature could be an important factor of drug translocation through the membrane.

The asymmetry of plasma membranes and their cholesterol content influence the uptake of cisplatin

The composition of the plasma membrane of malignant cells is thought to influence the cellular uptake of cisplatin and to take part in developing resistance to this widespread anti-cancer drug. In this work we study the permeation of cisplatin through the model membranes of normal and cancer cells using molecular dynamics simulations. A special attention is paid to lipid asymmetry and cholesterol content of the membranes. The loss of lipid asymmetry, which is common for cancer cells, leads to a decrease in their permeability to cisplatin by one order of magnitude in comparison to the membranes of normal cells. The change in the cholesterol molar ratio from 0% to 33% also decreases the permeability of the membrane by approximately one order of magnitude. The permeability of pure DOPC membrane is 5-6 orders of magnitude higher than one of the membranes with realistic lipid composition, which makes it as an inadequate model for the studies of drug permeability.

Dependence of fullerene aggregation on lipid saturation due to a balance between entropy and enthalpy

It is well-known that fullerenes aggregate inside lipid membranes and that increasing the concentration may lead to (lethal) membrane rupture. It is not known, however, how aggregation and rupture depend on the lipid type, what physical mechanisms control this behavior and what experimental signatures detect such changes in membranes. In this paper, we attempt to answer these questions with molecular simulations, and we show that aggregation and membrane damage depend critically on the degree of saturation of the lipid acyl chains: unsaturated bonds, or "kinks", impose a subtle but crucial compartmentalization of the bilayer into core and surface regions leading to three distinct fullerene density maxima. In contrast, when the membrane has only fully saturated lipids, fullerenes prefer to be located close to the surface under the head groups until the concentration becomes too large and the fullerenes begin clustering. No clustering is observed in membranes with unsaturated lipids. The presence of "kinks" reverses the free energy balance; although the overall free energy profiles are similar, entropy is the dominant component in unsaturated bilayers whereas enthalpy controls the fully saturated ones. Fully saturated systems show two unique signatures: 1) membrane thickness behaves non-monotonously while the area per lipid increases monotonously. We propose this as a potential reason for the observations of low fullerene concentrations being effective against bacteria. 2) The fullerene-fullerene radial distribution function (RDF) shows splitting of the second peak indicating the emergence short-range order and the importance of the second-nearest neighbor interactions. Similar second peak splitting has been reported in metal glasses.

Effect of C60 on the phase transition behavior of a lipid bilayer under high pressure

Interactions between fullerenes and cells and effects on the main transition of lipid bilayers have attracted much attention in biophysics in recent years. By employing coarse-grained molecular dynamics simulations, we obtained the temperature–pressure phase diagrams of a dipalmitoylphosphatidylcholine bilayer, which exhibits a gel phase and a fluid phase, with variation of the C₆₀versus lipid ratios. The simulation results show that the critical area per lipid at the fluid–gel main phase transition boundary increases with the increasing ratios of C₆₀. A critical area per lipid of 0.594 ± 0.002 nm² is obtained when the ratio of C₆₀ reaches 6.4% while that of the pure bilayer case is 0.572 ± 0.002 nm². The main transition temperature, T_m, remains almost unchanged with the addition of C₆₀ below a ratio of 4.7%, while a 2 K decrease of T_m is observed at a ratio of 6.4% under various pressures. Consequently, the presence of C₆₀ in the bilayer, with the ratio of C₆₀ less than 4.7%, will not influence the main transition behavior of the bilayer even under pressure as high as 1500 bar. The radial distribution function analyses suggest that the presence of C₆₀ produces no impact on the radial distribution of the lipids in the bilayers. The lateral density profiles show that the incorporation of C₆₀ with relatively high ratios stabilizes the thickness of the bilayer.

By employing coarse-grained molecular dynamics simulations, we obtained the temperature–pressure phase diagrams of a dipalmitoylphosphatidylcholine bilayer, which exhibits a gel phase and a fluid phase, with variation of the C₆₀versus lipid ratios.

The influence of curvature on the properties of the plasma membrane. Insights from atomistic molecular dynamics simulations

In this work we conduct a systematic analysis of the influence of curvature on various properties of a realistic model of mammalian plasma membrane with asymmetric lipid content of monolayers and a realistic concentration of cholesterol. In order to do this we developed new technique for simulating membranes with the global membrane curvature restricted to any desirable value while keeping free lateral diffusion of the lipids and without introducing artifacts or perturbing the membrane structure. We show that the thickness of the hydrophobic core of the concave monolayer decreases by approximately 1.3 Å in comparison to that of the flat membrane, while the thickness of the convex monolayer does not change. The order parameter of the lipid tails decreases significantly in the certain layers of the curved membrane. The area per lipid increases in the convex monolayer and decreases in the concave monolayer. The cholesterol inclination angle decreases when the curvature of a particular monolayer changes from concave to convex. The amount of cholesterol in the minor fraction located between the membrane leaflets is zero in the membrane with positive curvature and increases to 1.7% in the flat membrane and to 2.5% in the membrane with negative curvature.

Interaction between charged nanoparticles and vesicles: coarse-grained molecular dynamics simulations

Interactions between charged nanoparticles and curved zwitterionic lipid vesicles with different surface charge densities (ρ): insertion, repulsion, adsorption, and penetration.