Acyl chain order parameter profiles in phospholipid bilayers: computation from molecular dynamics simulations and comparison with 2H NMR experiments

Manufacturing process of liposomal Formation: A coarse-grained molecular dynamics simulation

A method of producing liposomes has been previously developed using a continuous manufacturing technology that involves a co-axial turbulent jet in co-flow. In this study, coarse-grained molecular dynamics (CG-MD) simulations were used to gain a deeper understanding of how the self-assembly process of liposomes is affected by the material attributes (such as the concentration of ethanol) and the process parameters (such as temperature), while also providing detailed information on a nano-scale molecular level. Specifically, the CG-MD simulations yield a comprehensive internal view of the structure and formation mechanisms of liposomes containing DPPC, DPPG, and cholesterol molecules. The importance of this work is that structural details on the molecular level are proposed, and such detail is not possible to obtain through experimental studies alone. The assessment of structural properties, including the area per lipid, diffusion coefficient, and order parameters, indicated that a thicker bilayer was observed at higher ethanol concentrations, while a thinner bilayer was present at higher temperatures. These conditions led to more water penetrating the interior of the bilayer and an unstable structure, as indicated by a larger contact area between lipids and water, and a higher coefficient of lipid lateral diffusion. However, stable liposomes were found through these evaluations at lower ethanol concentrations and/or lower process temperatures. Furthermore, the CG-MD model was further compared and validated with experimental and computational data including liposomal bilayer thickness and area per lipid measurements.

Interfacial Glucose to Regulate Hydrated Lipid Bilayer Properties: Influence of Concentrations

A series of atomistic molecular dynamics (MD) simulations were carried out with a hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayer with the variation of glucose concentrations from 0 to 30 wt % in the presence of 0.3 M NaCl. The study suggested that although the thickness of the lipid bilayer dropped significantly with the increase in glucose concentration, it expanded laterally at high glucose levels due to the intercalation of glucose between the headgroups of adjacent lipids. We adopted the surface assessment via the grid evaluation method to compute the deviation of the bilayer's key structural features for the different amounts of glucose present. This suggested that the accumulation of glucose molecules near the headgroups influences the local lipid bilayer undulation and crimping of the lipid tails. We find that the area compressibility modulus increases with the glucose level, causing enhanced bilayer rigidity arising from the slow lateral diffusion of lipids. The restricted lipid motion at high glucose concentrations controls the sustainability of the curved bilayer surface. Calculations revealed that certain orientations of CO → of interfacial glucose with the PN → of lipid headgroups are preferred, which helps the glucose to form direct hydrogen bonds (HBs) with the lipid headgroups. Such lipid-glucose (LG) HBs relax slowly at low glucose concentrations and exhibit a higher lifetime, whereas fast structural relaxation of LG HBs with a shorter lifetime was noticed at a higher glucose level. In contrast, lipid-water (LW) HBs exhibited a higher lifetime at a higher glucose level, which gradually decreased with the glucose level lowering. The study interprets that the glucose concentration-driven LW and LG interactions are mutually inclusive. Our detailed analysis will exemplify small saccharide concentration-driven membrane stabilizing efficiency, which is, in general, helpful for drug delivery study.

Antimicrobial peptide interactions with bacterial cell membranes

Antimicrobial peptides (AMPs) are potential alternatives for common antibiotics because of their greater activity and efficiency against a broad range of viruses, bacteria, fungi, and parasites. In this project, two antimicrobial peptides including magainin 2 and protegrin 1 with α-helix and β-sheet secondary structures were selected to investigate their interactions with different lipid bilayers such as 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoserine (POPS), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), POPC/POPG (7:3), POPC/POPS (7:3), POPG/POPE(1:3), and POPG/POPE(3:1). The obtained structures of the AMPs illustrated that protegrin 1 cannot maintain its secondary structure in the solution phase in contrast to magainin 2. The head groups of the lipid units play a key role in the stability of the lipid bilayers. The head parts of the lipid membranes by increasing the internal H-bond contribute to membrane compactness. The POPG and POPS units inside the POPC/POPG and POPC/POPS membranes increase the order of the POPC units. The cationic residues of the AMPs form remarkable electrostatic interactions with the negatively charged membrane surfaces, which play a key role in the stabilization process of the peptide secondary structures. The Arg residues of protegrin 1 and the Gly1, Lys4, Lys10, Lys11, Lys14, and Glu19 of the magainin 2 have the most important roles in the complexation process. The values of Gibbs binding energies (ΔG) indicate that the complexation process between AMPs and different bacterial membranes is favorable from the thermodynamic viewpoint and AMPs could form stable complexes with the lipid bilayers. As a result of ΔG values, protegrin 1 forms a more stable complex with POPG/POPE(3:1), while the α-helix has more affinity to the POPG/POPE(1:3) bacterial membranes. Therefore, it can be considered that β-sheet and α-helix AMPs are more effective against gram-positive and gram-negative bacteria, respectively. The results of this study can provide useful details about the antimicrobial peptide interactions with the bacterial cell, which can be employed for designing new antimicrobial materials with greater efficiency.Communicated by Ramaswamy H. Sarma.

Molecular dynamics simulation of apolipoprotein E3 lipid nanodiscs

Nanodiscs are binary discoidal complexes of a phospholipid bilayer circumscribed by belt-like helical scaffold proteins. Using coarse-grained and all-atom molecular dynamics simulations, we explore the stability, size, and structure of nanodiscs formed between the N-terminal domain of apolipoprotein E3 (apoE3-NT) and variable number of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) molecules. We study both parallel and antiparallel double-belt configurations, consisting of four proteins per nanodisc. Our simulations predict nanodiscs containing between 240 and 420 DMPC molecules to be stable. The antiparallel configurations exhibit an average of 1.6 times more amino acid interactions between protein chains and 2 times more ionic contacts, compared to the parallel configuration. With one exception, DMPC order parameters are consistently larger in the antiparallel configuration than in the parallel one. In most cases, the root mean square deviation of the positions of the protein backbone atoms is smaller in the antiparallel configuration. We further report nanodisc size, thickness, radius of gyration, and solvent accessible surface area. Combining all investigated parameters, we hypothesize the antiparallel protein configuration leading to more stable and more rigid nanodiscs than the parallel one.

[Atomic Force Microscopy to Measure the Mechanical Property of Nanosized Lipid Vesicles and Its Applications]

Nanoparticles, including liposomes and lipid nanoparticles, have garnered global attention due to their potential applications in pharmaceuticals, vaccines, and gene therapies. These particles enable targeted delivery of new drug modalities such as highly active small molecules and nucleic acids. However, for widespread use of nanoparticle-based formulations, it is crucial to comprehensively analyze their characteristics to ensure both efficacy and safety, as well as enable consistent production. In this context, this review focuses on our research using atomic force microscopy (AFM) to study liposomes and lipid nanoparticles. Our work significantly contributes to the capability of AFM to measure various types of liposomes in an aqueous medium, providing valuable insights into the mechanical properties of these nanoparticles. We discuss the applications of this AFM technique in assessing the quality of nanoparticle-based pharmaceuticals and developing membrane-active peptides.

Structure and Elasticity of Mitochondrial Membranes: A Molecular Dynamics Simulation Study

Mitochondria are known as the powerhouse of the cell because they produce energy in the form of adenosine triphosphate. They also have other crucial functions such as regulating apoptosis, calcium homeostasis, and reactive oxygen species production. To perform these diverse functions, mitochondria adopt specific structures and frequently undergo dynamic shape changes, indicating that their mechanical properties play an essential role in their functions. To gain a detailed understanding at the molecular level of the structure and mechanical properties of mitochondria, we carry out atomistic molecular dynamics simulations for three inner mitochondrial membranes and three outer mitochondrial membrane models. These models take into account variations in cardiolipin and cholesterol concentrations as well as the symmetry/asymmetry between the two leaflets. Our simulations allow us to calculate various structural quantities and the bending, twisting, and tilting elastic moduli of the membrane models. Our results indicate that the structures of the inner and outer mitochondrial membranes are quite similar and do not depend much on the variation in lipid compositions. However, the bending modulus of the membranes increases with increasing concentrations of cardiolipin or cholesterol but decreases with a membrane asymmetry. Notably, we found that the dipole potential of the membrane increases with an increasing cardiolipin concentration. Finally, possible roles of cardiolipin in regulating ion and proton currents and maintaining the cristate are discussed in some details.

Anomalous lateral diffusion of lipids during the fluid/gel phase transition of a lipid membrane

A lipid membrane undergoes a phase transition from fluid to gel phase upon changing external thermodynamic conditions, such as decreasing temperature and increasing pressure. Extremophilic organisms face the challenge of preventing this deleterious phase transition. The main focus of their adaptive strategy is to facilitate effective temperature sensing through sensor proteins, relying on the drastic changes in packing density and membrane fluidity during the phase transition. Although the changes in packing density parameters due to the fluid/gel phase transition are studied in detail, the impact on membrane fluidity is less explored in the literature. Understanding the lateral diffusive dynamics of lipids in response to temperature, particularly during the fluid/gel phase transition, is albeit crucial. Here we have simulated the phase transition of a single component lipid membrane composed of dipalmitoylphosphatidylcholine (DPPC) lipids using a coarse-grained (CG) model and studied the changes of the structural and dynamical properties. It is observed that near the phase transition point, both fluid and gel phase domains coexist together. The dynamics remains highly non-Gaussian for a long time even when the mean square displacement reaches the Fickian regime at a much earlier time. This Fickian yet non-Gaussian diffusion (FnGD) is a characteristic of a highly heterogeneous system, previously observed for the lateral diffusion of lipids in raft mimetic membranes having liquid-ordered and liquid-disordered phases co-existing together. We have analyzed the molecular trajectories and calculated the jump-diffusion of the lipids, stemming from sudden jump translations, using a translational jump-diffusion (TJD) approach. An overwhelming contribution of the jump-diffusion of the lipids is observed suggesting anomalous diffusion of lipids during fluid/gel phase transition of the membrane. These results are important in unravelling the intricate nature of lipid diffusion during the phase transition of the membrane and open up a new possibility of investigating the most significant change of membrane properties during phase transition, which can be effectively sensed by proteins.

Enhanced Uptake of Anticancer C6-Pep Dimer in a Model Membrane Caused by Differential pKa in Acidic Microenvironment

Acidic tumor microenvironment (TME) presents a challenge for the action of antitumor drugs by acting as an additional barrier for the passive crossing of the cell membrane by chemotherapic agents playing a critical role in the proliferation of tumor cells. Anticancer lipopeptide C6-Pep dimer containing the leucine zipper motif shows an increased uptake into the model tumor membrane in TME, and application of external heat might lead to the uncoiling of the zipper, which could result in cell lysis. This work investigated the cause of this increased uptake of C6-Pep dimer into the bilayer model in TME. Accurate protonation states of all the titratable residues of the C6-Pep dimer in TME were determined using constant pH molecular dynamics. In TME, except for two terminal Glu5 residues, all other Glu residues in the C6-Pep dimer were permanently protonated. The remaining Glu5 residues had differential pKa values, leading to the construction of four possible dimers with different fixed protonation states, and molecular dynamics was used to study their interaction with the anionic bilayer. Except for the dimer at a physiological pH, the other dimers were positively charged and could readily adsorb on the membrane surface. The free energy of insertion of these dimers in the bilayer was lower for single and double protonated Glu5-containing dimers than for the others. After the insertion of the lipopeptides into the membrane, thinning of the bilayer in the vicinity of dimers and an increase in area per lipid of the bilayer were observed for all systems, indicating destabilization of the bilayer due to this intercalation. This study shows that the anticancer lipopeptide C6-Pep utilizes the TME around a tumor cell for insertion into the membrane.

Structural and Mechanical Response of Two-Component Photoswitchable Lipid Bilayer Vesicles

Optical control of phospholipids is an attractive option for the rapid, reversible, and tunable manipulation of membrane structure and dynamics. Azo-PC, a lipid with an azobenzene group within one acyl chain, undergoes a light-induced trans-to-cis isomerization and thus arises as a powerful tool for manipulating lipid order and dynamics. Here, we report on vesicle-scale micropipette measurements and atomistic simulations to probe the elastic stretching modulus, water permeability, toughness, thickness, and membrane area upon isomerization. We investigated both dynamics and steady-state properties. In pure azo-PC membranes, we found that the molecular area in trans was 16% smaller than that in cis, the membrane's stretching modulus kA was 2.5 ± 0.3 times greater, and the water permeability PW was 3.5 ± 0.5 times smaller. We also studied mixtures of azo-PC with the miscible, unsaturated lipid DOPC. Atomistic molecular dynamics simulations show how the membrane thickness, chain order, and correlations across membrane leaflets explain the experimental data. Together, these data show how one rotating bond changes the molecular- and membrane-scale properties. These results will be useful for photopharmacology and for developing new materials whose permeability, elasticity, and toughness may be switched on demand.

Integration of an LPAR1 Antagonist into Liposomes Enhances Their Internalization and Tumor Accumulation in an Animal Model of Human Metastatic Breast Cancer

Lysophosphatidic acid receptor 1 (LPAR1) is elevated in breast cancer. The deregulation of LPAR1, including the function and level of expression, is linked to cancer initiation, progression, and metastasis. LPAR1 antagonists, AM095 or Ki16425, may be effective therapeutic molecules, yet their limited water solubility hinders in vivo delivery. In this study, we report on the synthesis of two liposomal formulations incorporating AM095 or Ki16425, embedded within the lipid bilayer, as targeted nanocarriers for metastatic breast cancer (MBC). The data show that the Ki16425 liposomal formulation exhibited a 50% increase in internalization by MBC mouse epithelial cells (4T1) and a 100% increase in tumor accumulation in a mouse model of MBC compared with that of a blank liposomal formulation (control). At the same time, normal mouse epithelial cells (EpH-4Ev) internalized the Ki16425 liposomal formulation 25% lesser than the control formulation. Molecular dynamics simulations show that the integration of AM095 or Ki16425 modified the physical and mechanical properties of the lipid bilayer, making it more flexible in these liposomal formulations compared with liposomes without drug. The incorporation of an LPAR1 antagonist within a liposomal drug delivery system represents a viable therapeutic approach for targeting the LPA-LPAR1 axis, which may hinder the progression of MBC.

Structural Phenomena in a Vesicle Membrane Obtained through an Evolution Experiment: A Study Based on MD Simulations

The chemical evolution of biomolecules was clearly affected by the overall extreme environmental conditions found on Early Earth. Periodic temperature changes inside the Earth's crust may have played a role in the emergence and survival of functional peptides embedded in vesicular compartments. In this study, all-atom molecular dynamic (MD) simulations were used to elucidate the effect of temperature on the properties of functionalized vesicle membranes. A plausible prebiotic system was selected, constituted by a model membrane bilayer from an equimolar mixture of long-chain fatty acids and fatty amines, and an octapeptide, KSPFPFAA, previously identified as an optimized functional peptide in an evolution experiment. This peptide tends to form the largest spontaneous aggregates at higher temperatures, thereby enhancing the pore-formation process and the eventual transfer of essential molecules in a prebiotic scenario. The analyses also suggest that peptide-amphiphile interactions affect the structural properties of the membrane, with a significant increase in the degree of interdigitation at the lowest temperatures under study.

The effect of metalation on antimicrobial piscidins imbedded in normal and oxidized lipid bilayers

Metalation of the N-terminal Amino Terminal Cu(ii)- and Ni(ii)-binding (ATCUN) motif may enhance the antimicrobial properties of piscidins. Molecular dynamics simulations of free and nickelated piscidins 1 and 3 (P1 and P3) were performed in 3 : 1 POPC/POPG and 2.6 : 1 : 0.4 POPC/POPG/aldo-PC bilayers (POPC, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine: POPG, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol; aldo-PC, 1-palmitoyl-2-(9'-oxo-nonanoyl)-sn-glycero-3-phosphocholine) bilayer models. Nickel(ii) binding decreases the conformation dynamics of the ATCUN motif and lowers the charge of the N-terminus to allow it to embed deeper in the bilayer without significantly changing the overall depth due to interactions of the charged half-helix of the peptide with the headgroups. Phe1⋯Ni2+ cation-π and Phe2-Phe1 CH-π interactions contribute to a small fraction of structures within the nickelated P1 simulations and may partially protect a bound metal from metal-centered chemical activity. The substitution of Phe2 for Ile2 in P3 sterically blocks conformations with cation-π interactions offering less protection to the metal. This difference between metalated P1 and P3 may indicate a mechanism by which peptide sequence can influence antimicrobial properties. Any loss of bilayer integrity due to chain reversal of the oxidized phospholipid chains of aldo-PC may be enhanced in the presence of metalated piscidins.

Characterizing the Interactions of Cell-Membrane-Disrupting Peptides with Lipid-Functionalized Single-Walled Carbon Nanotubes

Lipid-functionalized single-walled carbon nanotubes (SWNTs) have garnered significant interest for their potential use in a wide range of biomedical applications. In this work, we used molecular dynamics simulations to study the equilibrium properties of SWNTs surrounded by the phosphatidylcholine (POPC) corona phase and their interactions with three cell membrane disruptor peptides: colistin, TAT peptide, and crotamine-derived peptide. Our results show that SWNTs favor asymmetrical positioning within the POPC corona, so that one side of the SWNT, covered by the thinnest part of the corona, comes in contact with charged and polar functional groups of POPC and water. We also observed that colistin and TAT insert deeply into the POPC corona, while crotamine-derived peptide only adsorbs to the corona surface. In separate simulations, we show that three examined peptides exhibit similar insertion and adsorption behaviors when interacting with POPC bilayers, confirming that peptide-induced perturbations to POPC in conjugates and bilayers are similar in nature and magnitude. Furthermore, we observed correlations between the peptide-induced structural perturbations and the near-infrared emission of the lipid-functionalized SWNTs, which suggest that the optical signal of the conjugates transduces the morphological changes in the lipid corona. Overall, our findings indicate that lipid-functionalized SWNTs could serve as simplified cell membrane model systems for prescreening of new antimicrobial compounds that disrupt cell membranes.

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent κC modulus for bilayer bending. Our results show a strong correlation with κC values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

Molecular dynamics insight of interaction between the functionalized-carbon nanotube and cancerous cell membrane in doxorubicin delivery

BACKGROUND AND OBJECTIVE

Doxorubicin (DOX) is a known anticancer drug which is widely used in cancer therapy. Carbon nanotubes (CNTs) are among the most promising platforms for smart drug delivery applications. However, due to the toxicity and their low sulubility their application is limited and their functionalization with wide range of biomolecules are suggested. Therefore, the functionalized carbon nanotubes (f-CNT) with carboxyl (CNT-COO) and folic acid (CNT-COO-FA) were investigated as drug-carrier.

METHODS

Molecular dynamics (MD) simulation along with the Density Functional Theory (DFT) methods are being used to study the drug loading process on functionalized carbon nanotubes.

RESULTS

The results indicate that doxorubicin molecules interact more with CNT-COO-FA than CNT-COO. The embedded dipalmitoylphosphatidylcholine (DPPC) lipid bilayer with a folate receptor was considered a cancerous cell's representative model. Then the drug release from the f-CNTs near the lipid bilayer was simulated. The results showed that CNT-COO-FA with a pH and ligand-sensitive mechanism strongly interacts with cancerous cells, which led to higher drug release, in agreement with the experimental results. The conformational changes of the lipid bilayer and folate receptor during drug release were evaluated. The analysis showed that drug release from CNT-COO-FA has significantly changed lipid bilayer and receptor conformations. The obtained results were interpreted and justified by considering the molecular mechanisms which control the drug delivery in the studied systems.

CONCLUSIONS

Based on the obtained results, CNT-COO-FA has a better performance during the drug release compared to CNT-COO in delivering doxorubicin. Both pH and ligand sensitive mechanisms are found to be responsible for higher drug delivery efficiency of CNT-COO-FA. In contrast, CNT-COO can only enhance drug delivery efficiently with a pH-sensitive mechanism.

Interaction of psychedelic tryptamine derivatives with a lipid bilayer

Naturally occurring psychedelics have been used for a long time as remedies or in religious ceremonies and recreational activities. Recent studies have proven the therapeutic potential of some psychedelic compounds to safely treat a wide range of diseases such as anxiety, depression, migraine, and addiction. It is hypothesized that psychedelic compounds like tryptamines can exert their effects by two possible mechanisms: binding to the transmembrane serotonin receptor and/or modifying the properties of the neuronal membrane that can alter the conformational equilibrium and desensitize receptors. The impact of three different tryptamine class compounds with a tertiary amine (dimethyltryptamine, bufotenine, and 5-MeO-DMT) in both neutral and charged forms on a model bilayer lipid membrane are studied using all-atom MD simulations. All compounds partition into the bilayer, and change membrane properties, but to different extents. We determine the tendency of compounds to partition into the membrane by free energy calculations. Neutral tryptamines partition into the bilayer almost completely. Dimethyltryptamine and 5-MeO-DMT cross the membrane spontaneously during the simulation time, but bufotenine does not, although it has the maximum effect on the structural properties of the membrane. However, protonated compounds partition partially into the bilayer and cannot pass through the middle of the membrane during the simulation time. In this way, subtle alteration of chemical structure can play a significant role in the improvement or deterioration of partitioning of these compounds into the bilayer and their passage across the membrane.

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.

The impact of lipid polyunsaturation on the physical and mechanical properties of lipid membranes

The lipid composition of cellular membranes and the balance between the different lipid components can be impacted by aging, certain pathologies, specific diets and other factors. This is the case in a subgroup of individuals with psychiatric disorders, such as schizophrenia, where cell membranes of patients have been shown to be deprived in polyunsaturated fatty acids (PUFAs), not only in brain areas where the target receptors are expressed but also in peripheral tissues. This PUFA deprivation thus represents a biomarker of such disorders that might impact not only the interaction of antipsychotic medications with these membranes but also the activation and signaling of the targeted receptors embedded in the lipid membrane. Therefore, it is crucial to understand how PUFAs levels alterations modulate the different physical properties of membranes. In this paper, several biophysical approaches were combined (Laurdan fluorescence spectroscopy, atomic force microscopy, differential scanning calorimetry, molecular modeling) to characterize membrane properties such as fluidity, elasticity and thickness in PUFA-enriched cell membranes and lipid model systems reflecting the PUFA imbalance observed in some diseases. The impact of both the number of unsaturations and their position along the chain on the above properties was investigated. Briefly, data revealed that PUFA presence in membranes increases membrane fluidity, elasticity and flexibility and decreases its thickness and order parameter. Both the level of unsaturation and their position affect these membrane properties.

Characterizing the Interactions of Cell Membrane-Disrupting Peptides with Lipid-Functionalized Single-Walled Carbon Nanotube Systems for Antimicrobial Screening

Lipid-functionalized single-walled carbon nanotubes (SWNTs) have garnered significant interest for their potential use in a wide range of biomedical applications. In this work, we used molecular dynamics simulations to study the equilibrium properties of SWNTs surrounded by the phosphatidylcholine (POPC) corona phase, and their interactions with three cell membrane disruptor peptides: colistin, TAT peptide, and crotamine-derived peptide. Our results show that SWNTs favor asymmetrical positioning within the POPC corona, so that one side of the SWNT, covered by the thinnest part of the corona, comes in contact with charged and polar functional groups of POPC and water. We also observed that colistin and TAT insert deeply into POPC corona, while crotamine-derived peptide only adsorbs to the corona surface. Compared to crotamine-derived peptide, colistin and TAT also induce larger perturbations in the thinnest region of the corona, by allowing more water molecules to directly contact the SWNT surface. In separate simulations, we show that three examined peptides exhibit similar insertion and adsorption behaviors when interacting with POPC bilayers, confirming that peptide-induced perturbations to POPC in conjugates and bilayers are similar in nature and magnitude. Furthermore, we observed correlations between the peptide-induced structural perturbations and the near-infrared emission of the lipid-functionalized SWNTs, which suggest that the optical signal of the conjugates transduces the morphological changes in the lipid corona. Overall, our findings indicate that lipid-functionalized SWNTs could serve as simplified cell membrane model systems for pre-screening of new antimicrobial compounds that disrupt cell membranes.

Determining structure and action mechanism of LBF14 by molecular simulation

The recently discovered, membrane-active peptide LBF14 contains several non-proteinogenic amino acids and is able to transform vesicles into tubule networks. The exact membrane interaction mechanism and detailed secondary structure are yet to be determined. We performed molecular dynamics simulations of LBF14 and let it fold de novo into its ensemble of native secondary structures. Histidine protonation state effects on secondary structure were investigated. An MD simulation of the peptide with a lipid bilayer was performed. Simulation results were compared to circular dichroism and electron paramagnetic resonance data of previous studies. LBF14 contains a conserved helical section in an otherwise random structure. Helical stability is influenced by histidine protonation. The peptide localized to the polar layer of the membrane, consistent with experimental results. While the overall secondary structure is unaffected by membrane interaction, Ramachandran plot analysis yielded two distinct peptide conformations during membrane interaction. This conformational change was accompanied by residue repositioning within the membrane. LBF14 only affected the local order in the membrane, and had no measurable effect on pressure. The simulation results are consistent with the previously proposed membrane interaction mechanism of LBF14 and can additionally explain the local interaction mechanism. Communicated by Ramaswamy H. Sarma.

Computing Individual Area per Head Group Reveals Lipid Bilayer Dynamics

Lipid bilayers express a range of phases from solid-like to gel-like to liquid-like as a function of temperature and lipid surface concentration. The area occupied per lipid head group serves as one useful indicator of the bilayer phase, in conjunction with the two-dimensional radial distribution function (i.e., structure factor) within the bilayer. Typically, the area per head group is determined by dividing the bilayer area equally among all head groups. Such an approach is less satisfactory for a multicomponent set of diverse lipids. In this work, area determination is performed on a lipid-by-lipid basis by attributing to a lipid the volume that surrounds each atom. Voronoi tessellation provides this division of the interfacial region on a per-atom basis. The method is applied to a multicomponent system of water, NaCl, and 19 phospholipid types that was devised recently [Langmuir2022, 38, 9481-9499] as a computational representation of the Gram-positive Staphylococcus aureus phospholipid bilayer. Results demonstrate that lipids and water molecules occupy similar extents of area within the interfacial region; ascribing area only to head groups implicitly incorporates assumptions about head group hydration. Results further show that lipid tails provide non-negligible contributions to area on the membrane side of the bilayer-water interface. Results for minimum and maximum area of individual lipids reveal that spontaneous fluctuations displace head groups more than 10 Å from the interfacial region during an NPT simulation at 310 K, leading to a zero contribution to total area at some times. Total area fluctuations and fluctuations per individual lipid relax with a correlation time of ∼10 ns. The method complements density profile as an approach to quantify the structure and dynamics of computational lipid bilayers.

Molecular dynamics simulation of cancer cell membrane perforated by shockwave induced bubble collapse

It has been widely accepted that cancer cells are softer than their normal counterparts. This motivates us to propose, as a proof-of-concept, a method for the efficient delivery of therapeutic agents into cancer cells, while normal cells are less affected. The basic idea of this method is to use a water jet generated by the collapse of the bubble under shockwaves to perforate pores in the cell membrane. Given a combination of shockwave and bubble parameters, the cancer membrane is more susceptible to bending, stretching, and perforating than the normal membrane because the bending modulus of the cancer cell membrane is smaller than that of the normal cell membrane. Therefore, the therapeutic agent delivery into cancer cells is easier than in normal cells. Adopting two well-studied models of the normal and cancer membranes, we perform shockwave induced bubble collapse molecular dynamics simulations to investigate the difference in the response of two membranes over a range of shockwave impulse 15-30 mPa s and bubble diameter 4-10 nm. The simulation shows that the presence of bubbles is essential for generating a water jet, which is required for perforation; otherwise, pores are not formed. Given a set of shockwave impulse and bubble parameters, the pore area in the cancer membrane is always larger than that in the normal membrane. However, a too strong shockwave and/or too large bubble results in too fast disruption of membranes, and pore areas are similar between two membrane types. The pore closure time in the cancer membrane is slower than that in the normal membrane. The implications of our results for applications in real cells are discussed in some details. Our simulation may be useful for encouraging future experimental work on novel approaches for cancer treatment.

Change in membrane fluidity induced by polyphenols is highly dependent on the position and number of galloyl groups

The cell membrane fluidity was very important in adipogenesis and galloyl groups on polyphenolic structures could enhance their antiadipogenic activity. However, the effect of polyphenols on membrane fluidity and the role of galloyl groups in fluidity changes remain unclear. Therefore, the present study chose structurally different polyphenols to compare their effects on the membrane morphology and fluidity of 3T3-L1 preadipocytes, and then the reasons behind the changes of membrane fluidity induced by galloylated polyphenols were explored from structural and molecular insights using liposome model and molecular dynamic simulation technology. Our results indicated that galloylated polyphenols could significantly change 3T3-L1 cell membrane morphology and decrease membrane fluidity, while non-galloylated ones could not. The membrane interference effect of polyphenols was enhanced as the number of galloyl groups increased. Morever, the decrease in membrane fluidity induced by galloylated polyphenols was due to the disturbance of polyphenols on lipid alkyl chains in the cell membrane. Galloylated polyphenols could not only locate in the polar head, but also insert into hydrophobic center of lipid bilayer to interfere with the lipid alkyl chains arrangement, thus decreasing the membrane fluidity and showing strong affinity for the membrane. In addition, differences in position of galloyl groups in polyphenols induced distinct effect on cell membranes interactions, thus affecting the binding manner and bioactivity. The results expanded the understanding on the strong antiadipogenic activity of galloylated polyphenols through the aspect of their effects on cell membrane by both experimental and theoretically simulated ways.

Self-Uptake Mechanism of Polymyxin-Based Lipopeptide against Gram-Negative Bacterial Membrane: Role of the First Adsorbed Lipopeptide

Research in the continuously increasing threat of polymyxin-resistant multidrug-resistant Pseudomonas aeruginosa, which causes severe infection in immunocompromised patients, has resulted in the development of several polymyxin-derived cyclic lipopeptides containing l-α-γ- diamino butyric acid-like FADDI-019 (F19). In this work, F19's insertion into a minimal model of the asymmetric outer membrane of the bacterium, which contained only penta-acylated lipid A (LipA) and lacked keto-d-octulosonic acid and O-antigens, in the top leaflet and phospholipids in the bottom leaflet, was studied. F19 exhibited all of the hallmarks of the self-uptake mechanism into the asymmetric bilayer. While a single monomer of the lipopeptide did not get partitioned into the inside of the bilayer, it competitively displaced Ca²⁺ from the membrane surface, observed as a decrease in Ca²⁺ coordination number with phosphate groups (1.89 vs 1.718), resulting in membrane destabilization. This resulted in an increment of the average defect size and the probability of interplay between lipid tails and hydrophobic residues of another F19. When more than one monomer was present in the system, the first monomer remained docked on the surface, while other monomers intercalated into the bilayer interior with their hydrophobic moieties "sleeved" by lipid acyl chains. The free energy barrier for partial insertion of the lipopeptide into a bilayer in the presence of surface-docked second F19 was recorded at ∼1.3 kcal/mol using two-dimensional (2D) well-tempered metadynamics, making it a low barrier process at 300 K. This study is an attempt to demonstrate the self-uptake mechanism of F19 during intercalation process into the bilayer interior, which may help in the design of better alternates for polymyxins to work against polymyxin resistance.

Using molecular simulation to understand the skin barrier

Skin's effectiveness as a barrier to permeation of water and other chemicals rests almost entirely in the outermost layer of the epidermis, the stratum corneum (SC), which consists of layers of corneocytes surrounded by highly organized lipid lamellae. As the only continuous path through the SC, transdermal permeation necessarily involves diffusion through these lipid layers. The role of the SC as a protective barrier is supported by its exceptional lipid composition consisting of ceramides (CERs), cholesterol (CHOL), and free fatty acids (FFAs) and the complete absence of phospholipids, which are present in most biological membranes. Molecular simulation, which provides molecular level detail of lipid configurations that can be connected with barrier function, has become a popular tool for studying SC lipid systems. We review this ever-increasing body of literature with the goals of (1) enabling the experimental skin community to understand, interpret and use the information generated from the simulations, (2) providing simulation experts with a solid background in the chemistry of SC lipids including the composition, structure and organization, and barrier function, and (3) presenting a state of the art picture of the field of SC lipid simulations, highlighting the difficulties and best practices for studying these systems, to encourage the generation of robust reproducible studies in the future. This review describes molecular simulation methodology and then critically examines results derived from simulations using atomistic and then coarse-grained models.

Phospholipid Monolayer/Graphene Interfaces: Curvature Effect on Lipid Morphology and Dynamics

Phospholipids are an important class of lipids that are widely used as model platforms for the study of biological processes and interactions. These lipids can form stable interfaces with solid substrates, such as graphene, and these interfaces have potential applications in biosensing and targeted drug delivery. In this paper, we perform molecular dynamics simulations of graphene-supported lipid monolayers to characterize the lipid properties of such interfaces. We observed substantial differences between the supported monolayer and free-standing bilayer in terms of the lipid properties, such as the tail order parameters, density profiles, diffusion rates, and so on. Furthermore, we studied these interfaces on sinusoidally deformed graphene substrates to understand the effect of curvature on the supported lipids. Here, we observed that the nature of the substrate curvature, that is, concave or convex, can locally affect the lipid/substrate adhesion strength and induce structural and dynamic changes in the adsorbed lipid monolayer. Together, these results help characterize the properties of lipid/graphene interfaces and provide insights into the substrate curvature effect on these interfaces, which can enable the tuning of lipid properties for various sensor devices and drug delivery applications.

Generation and Computational Characterization of a Complex Staphylococcus aureus Lipid Bilayer

Studies indicate a crucial cell membrane role in the antibiotic resistance of Staphylococcus aureus. To simulate its membrane structure and dynamics, a complex molecular-scale computational representation of the S. aureus lipid bilayer was developed. Phospholipid types and their amounts were optimized by reverse Monte Carlo to represent characterization data from the literature, leading to 19 different phospholipid types that combine three headgroups [phosphatidylglycerol, lysyl-phosphatidylglycerol (LPG), and cardiolipin] and 10 tails, including iso- and anteiso-branched saturated chains. The averaged lipid bilayer thickness was 36.7 Å, and area per headgroup was 67.8 Ų. Phosphorus and nitrogen density profiles showed that LPG headgroups tended to be bent and oriented more parallel to the bilayer plane. The water density profile showed that small amounts reached the membrane center. Carbon density profiles indicated hydrophobic interactions for all lipids in the middle of the bilayer. Bond vector order parameters along each tail demonstrated different C-H ordering even within distinct lipids of the same type; however, all tails followed similar trends in average order parameter. These complex simulations further revealed bilayer insights beyond those attainable with monodisperse, unbranched lipids. Longer tails often extended into the opposite leaflet. Carbon at and beyond a branch showed significantly decreased ordering compared to carbon in unbranched tails; this feature arose in every branched lipid. Diverse tail lengths distributed these disordered methyl groups throughout the middle third of the bilayer. Distributions in mobility and ordering reveal diverse properties that cannot be obtained with monodisperse lipids.

Unraveling membrane properties at the organelle-level with LipidDyn

•

Lipidomics of organelles could be used to design models for molecular simulations.

•

The bottleneck is the analysis and rationalization of the data from simulations.

•

LipidDyn is an automated pipeline to streamline the analyses of lipid bilayers.

•

LipidDyn allows to collect analysis in a non-time-consuming and reproducible manner.

•

We applied LipidDyn to different case studies to illustrate its potential.

Cellular membranes are formed from different lipids in various amounts and proportions depending on the subcellular localization. The lipid composition of membranes is sensitive to changes in the cellular environment, and its alterations are linked to several diseases. Lipids not only form lipid-lipid interactions but also interact with other biomolecules, including proteins.

Molecular dynamics (MD) simulations are a powerful tool to study the properties of cellular membranes and membrane-protein interactions on different timescales and resolutions. Over the last few years, software and hardware for biomolecular simulations have been optimized to routinely run long simulations of large and complex biological systems. On the other hand, high-throughput techniques based on lipidomics provide accurate estimates of the composition of cellular membranes at the level of subcellular compartments. Lipidomic data can be analyzed to design biologically relevant models of membranes for MD simulations. Similar applications easily result in a massive amount of simulation data where the bottleneck becomes the analysis of the data. In this context, we developed LipidDyn, a Python-based pipeline to streamline the analyses of MD simulations of membranes of different compositions. Once the simulations are collected, LipidDyn provides average properties and time series for several membrane properties such as area per lipid, thickness, order parameters, diffusion motions, lipid density, and lipid enrichment/depletion. The calculations exploit parallelization, and the pipeline includes graphical outputs in a publication-ready form. We applied LipidDyn to different case studies to illustrate its potential, including membranes from cellular compartments and transmembrane protein domains. LipidDyn is available free of charge under the GNU General Public License from https://github.com/ELELAB/LipidDyn.

WIND-PVPA: Water/Ion NMR Detected PVPA to assess lipid barrier integrity in vitro through quantification of passive water- and ion transport

Water/Ion NMR Detected - Phospholipid Vesicle Permeability Assay (WIND-PVPA), is presented as a novel, straightforward and automatable method to assess lipid barrier integrity in vitro. The apparent permeability constants of water- and ions across the PVPA barriers are determined in a one-pot experiment under the influence of membrane-active guest molecules. NMR spectroscopy is used to quantify the water directly (D₂O) and the ions indirectly (complexed with EDTA) as a function of time. WIND-PVPA is demonstrated using four anti-microbial peptides, to show that membrane active molecules can be differentiated by their disruptive influence on the PVPA system. The results obtained are compared with explicit molecular dynamics simulations of lipid bilayers, AMPs, water and salt, where the motions of all individual water molecules relative to the lipid bilayer are monitored over the course of the simulations, allowing the calculation of theoretical apparent permeability constants of the corresponding single bilayer systems. Proof-of-principle is presented that WIND-PVPA can be used to evaluate the lipid barrier destabilizing effect of active guest molecules by measuring changes in passive water- and ion permeabilities upon exposure. The method is highly flexible in terms of barrier composition, choice of probes and membrane active compounds.

Amyloid β Dodecamer Disrupts the Neuronal Membrane More Strongly than the Mature Fibril: Understanding the Role of Oligomers in Neurotoxicity

The amyloid cascade hypothesis states that senile plaques, composed of amyloid β (Aβ) fibrils, play a key role in Alzheimer’s disease (AD). However, recent experiments have shown that Aβ oligomers are more toxic to neurons than highly ordered fibrils. The molecular mechanism underlying this observation remains largely unknown. One of the possible scenarios for neurotoxicity is that Aβ peptides create pores in the lipid membrane that allow Ca²⁺ ions to enter cells, resulting in a signal of cell apoptosis. Hence, one might think that oligomers are more toxic due to their higher ability to create ion channels than fibrils. In this work, we study the effect of Aβ42 dodecamer and fibrils on a neuronal membrane, which is similar to that observed in AD patients, using all-atom molecular dynamics simulations. Due to short simulation times, we cannot observe the formation of pores, but useful insight on the early events of this process has been obtained. Namely, we showed that dodecamer distorts the lipid membrane to a greater extent than fibrils, which may indicate that ion channels can be more easily formed in the presence of oligomers. Based on this result, we anticipate that oligomers are more toxic than mature fibrils, as observed experimentally. Moreover, the Aβ–membrane interaction was found to be governed by the repulsive electrostatic interaction between Aβ and the ganglioside GM1 lipid. We calculated the bending and compressibility modulus of the membrane in the absence of Aβ and obtained good agreement with the experiment. We predict that the dodecamer will increase the compressibility modulus but has little effect on the bending modulus. Due to the weak interaction with the membrane, fibrils insignificantly change the membrane elastic properties.

Active site geometry stabilization of a presenilin homolog by the lipid bilayer promotes intramembrane proteolysis

Cleavage of membrane proteins in the lipid bilayer by intramembrane proteases is crucial for health and disease. Although different lipid environments can potently modulate their activity, how this is linked to their structural dynamics is unclear. Here, we show that the carboxy-peptidase-like activity of the archaeal intramembrane protease PSH, a homolog of the Alzheimer’s disease-associated presenilin/γ-secretase is impaired in micelles and promoted in a lipid bilayer. Comparative molecular dynamics simulations revealed that important elements for substrate binding such as transmembrane domain 6a of PSH are more labile in micelles and stabilized in the lipid bilayer. Moreover, consistent with an enhanced interaction of PSH with a transition-state analog inhibitor, the bilayer promoted the formation of the enzyme’s catalytic active site geometry. Our data indicate that the lipid environment of an intramembrane protease plays a critical role in structural stabilization and active site arrangement of the enzyme-substrate complex thereby promoting intramembrane proteolysis.

Cutting proteins into pieces is a crucial process in the cell, allowing several important processes to take place, including cell differentiation (which allows cells to develop into specific types), cell death, protein quality control, or even where in the cell a protein will end up. However, the specialized proteins that carry out this task, known as proteases, can also be involved in the development of disease. For example, in the brain, a protease called γ-secretase cuts up the amyloid-β protein precursor, producing toxic forms of amyloid-β peptides that are widely believed to cause Alzheimer’s disease.

Proteases like γ-secretase carry out their role in the membrane, the layer of fats (also known as lipids) that forms the outer boundary of the cell. The environment in this area of the cell can influence the activity of proteases, but it is poorly understood how this happens.

One way to address this question would be to compare the activity of γ-secretase in the lipid environment of the membrane to its activity when it is entirely surrounded by different molecules, such as detergent molecules. Unfortunately, γ-secretase is not active when it is removed from its lipid environment by a detergent, making it difficult to perform this comparison. To overcome this issue, Feilen et al. chose to study PSH, a protease similar to γ-secretase that produces the same amyloid-β peptides but remains active in detergent.

When Feilen et al. mixed PSH with lipid molecules like those found in the membrane and amyloid-β precursor protein, PSH produced amyloid-β peptides including those that are thought to cause Alzheimer’s. However, when a detergent was substituted for the lipid molecules this led to longer amyloid-β peptides than usual, indicating that PSH was not able to cut proteins as effectively. The change in environment appeared to reduce PSH’s ability to progressively trim small segments from the peptides.

Computer modelling of the protease’s structure in lipids versus detergent supported the experimental findings: the model predicted that the areas of PSH important for recognizing and cutting other proteins would be more stable in the membrane compared to the detergent.

These results indicate that the cell membrane plays a vital role in the stability of the active regions of proteases that are cleaving in this environment. In the future, this could help to better understand how changes to the lipid molecules in the membrane may contribute to the activity of γ-secretase and its role in Alzheimer’s disease.

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.

Elastic moduli of normal and cancer cell membranes revealed by molecular dynamics simulations

Recent studies indicate that there are mechanical differences between normal cells and cancer cells. Because the cell membrane takes part in a variety of vital processes, we test the hypothesis of whether or not two fundamental alterations in the cell membrane, i.e., the overexpression of phosphatidylserine lipids in the outer leaflet and a reduction in cholesterol concentration, could cause the softening in cancer cells. Adopting ten models of normal and cancer cell membranes, we carry out 1 μs all-atom molecular dynamics simulations to compare the structural properties and elasticity properties of two membrane types. We find that the overexpression of the phosphatidylserine lipids in the outer leaflet does not significantly alter the area per lipid, the membrane thickness, the lipid order parameters and the elasticity moduli of the cancer membranes. However, a reduction in the cholesterol concentration leads to clear changes in those quantities, especially decreases in the bending, tilt and twist moduli. This implies that the reduction of cholesterol concentration in the cancer membranes could contribute to the softening of cancer cells.

Computational Study on the Microscopic Adsorption Characteristics of Linear Alkylbenzene Sulfonates with Different Chain Lengths on Anthracite Surface

In order to explore the influence of different lengths of hydrophobic carbon chains on the diffusion characteristics of surfactants on the surface of anthracite, six linear alkyl benzene sulfonates with different hydrophobic carbon chain lengths were selected (mC, m = 8, 10, 12, 14, 16, 18; m represents the numbers of carbon atoms in the hydrophobic carbon chain), and molecular dynamics (MD) simulations were adopted. Models of surfactant-anthracite, surfactant-graphite layer, and water-surfactant-anthracite were constructed. After analyzing a series of properties such as adsorption energy, diffusion coefficient, radial distribution function (RDF), and hydrophobic tail order parameters, it was found that 12C had the highest adsorption strength on the surface of anthracite; the reason was that 12C had the highest degree of aggregation near the oxygen-containing functional groups on the surface of anthracite. Further studies had found that the hydrophobic tail chain of 12C had the strongest isotropy. The study fills the gap in the systematic study of the diffusion characteristics of linear alkylbenzene sulfonates (LAS) with different chain lengths on the surface of anthracite, enriches and develops the basic theory of coal wettability, and also provides technical ideas for the design of new surfactants and new dust suppression agents.

Interactions between Rhodamine Dyes and Model Membrane Systems—Insights from Molecular Dynamics Simulations

Background: rhodamines are dyes widely used as fluorescent tags in cell imaging, probing of mitochondrial membrane potential, and as P-glycoprotein model substrates. In all these applications, detailed understanding of the interaction between rhodamines and biomembranes is fundamental. Methods: we combined atomistic molecular dynamics (MD) simulations and fluorescence spectroscopy to characterize the interaction between rhodamines 123 and B (Rh123 and RhB, respectively) and POPC bilayers. Results: while the xanthene moiety orients roughly parallel to the membrane plane in unrestrained MD simulations, variations on the relative position of the benzoic ring (below the xanthene for Rh123, above it for RhB) were observed, and related to the structure of the two dyes and their interactions with water and lipids. Subtle distinctions were found among different ionization forms of the probes. Experimentally, RhB displayed a lipid/water partition coefficient more than two orders of magnitude higher than Rh123, in agreement with free energy profiles obtained from umbrella sampling MD. Conclusions: this work provided detailed insights on the similarities and differences in the behavior of bilayer-inserted Rh123 and RhB, related to the structure of the probes. The much higher affinity of RhB for the membranes increases the local concentration and explains its higher apparent affinity for P-glycoprotein reconstituted in model membranes.

The Role of Lipids in Allosteric Modulation of Dopamine D2 Receptor-In Silico Study

The dopamine D₂ receptor, belonging to the class A G protein-coupled receptors (GPCRs), is an important drug target for several diseases, including schizophrenia and Parkinson’s disease. The D₂ receptor can be activated by the natural neurotransmitter dopamine or by synthetic ligands, which in both cases leads to the receptor coupling with a G protein. In addition to receptor modulation by orthosteric or allosteric ligands, it has been shown that lipids may affect the behaviour of membrane proteins. We constructed a model of a D₂ receptor with a long intracellular loop (ICL3) coupled with G_iα1 or G_iα2 proteins, embedded in a complex asymmetric membrane, and simulated it in complex with positive, negative or neutral allosteric ligands. In this study, we focused on the influence of ligand binding and G protein coupling on the membrane–receptor interactions. We show that there is a noticeable interplay between the cell membrane, G proteins, D₂ receptor and its modulators.

Examining the Effect of Charged Lipids on Mitochondrial Outer Membrane Dynamics Using Atomistic Simulations

The outer mitochondrial membrane (OMM) is involved in multiple cellular functions such as apoptosis, inflammation and signaling via its membrane-associated and -embedded proteins. Despite the central role of the OMM in these vital phenomena, the structure and dynamics of the membrane have regularly been investigated in silico using simple two-component models. Accordingly, the aim was to generate the realistic multi-component model of the OMM and inspect its properties using atomistic molecular dynamics (MD) simulations. All major lipid components, phosphatidylinositol (PI), phosphatidylcholine (PC), phosphatidylethanolamine (PE), and phosphatidylserine (PS), were included in the probed OMM models. Because increased levels of anionic PS lipids have potential effects on schizophrenia and, more specifically, on monoamine oxidase B enzyme activity, the effect of varying the PS concentration was explored. The MD simulations indicate that the complex membrane lipid composition (MLC) behavior is notably different from the two-component PC-PE model. The MLC changes caused relatively minor effects on the membrane structural properties such as membrane thickness or area per lipid; however, notable effects could be seen with the dynamical parameters at the water-membrane interface. Increase of PS levels appears to slow down lateral diffusion of all lipids and, in general, the presence of anionic lipids reduced hydration and slowed down the PE headgroup rotation. In addition, sodium ions could neutralize the membrane surface, when PI was the main anionic component; however, a similar effect was not seen for high PS levels. Based on these results, it is advisable for future studies on the OMM and its protein or ligand partners, especially when wanting to replicate the correct properties on the water-membrane interface, to use models that are sufficiently complex, containing anionic lipid types, PI in particular.

A method to construct the dynamic landscape of a bio-membrane with experiment and simulation

Biomolecular function is based on a complex hierarchy of molecular motions. While biophysical methods can reveal details of specific motions, a concept for the comprehensive description of molecular dynamics over a wide range of correlation times has been unattainable. Here, we report an approach to construct the dynamic landscape of biomolecules, which describes the aggregate influence of multiple motions acting on various timescales and on multiple positions in the molecule. To this end, we use 13C NMR relaxation and molecular dynamics simulation data for the characterization of fully hydrated palmitoyl-oleoyl-phosphatidylcholine bilayers. We combine dynamics detector methodology with a new frame analysis of motion that yields site-specific amplitudes of motion, separated both by type and timescale of motion. In this study, we show that this separation allows the detailed description of the dynamic landscape, which yields vast differences in motional amplitudes and correlation times depending on molecular position.

Rotational decoupling between the hydrophilic and hydrophobic regions in lipid membranes

Cells use homeostatic mechanisms to ensure an optimal composition of distinct types of lipids in cellular membranes. The hydrophilic region of biological lipid membranes is mainly composed of several types of phospholipid headgroups that interact with incoming molecules, nanoparticles, and viruses, whereas the hydrophobic region consists of a distribution of acyl chains and sterols affecting membrane fluidity/rigidity related properties and forming an environment for membrane-bound molecules such as transmembrane proteins. A fundamental open question is to what extent the motions of these regions are coupled and, consequently, how strongly the interactions of phospholipid headgroups with other molecules depend on the properties and composition of the membrane hydrophobic core. We combine advanced solid-state nuclear magnetic resonance spectroscopy with high-fidelity molecular dynamics simulations to demonstrate how the rotational dynamics of choline headgroups remain nearly unchanged (slightly faster) with incorporation of cholesterol into a phospholipid membrane, contrasting the well-known extreme slowdown of the other phospholipid segments. Notably, our results suggest a new paradigm in which phospholipid dipole headgroups interact as quasi-freely rotating flexible dipoles at the interface, independent of the properties in the hydrophobic region.

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.

Drug Meets Monolayer: Understanding the Interactions of Sterol Drugs with Models of the Lung Surfactant Monolayer Using Molecular Dynamics Simulations

The lung surfactant monolayer (LSM) is a thin layer of lipids and proteins that forms the air/water interface of the alveoli. The primary function of the LSM is to reduce the surface tension at the air/water interface during breathing. The LSM also forms the main biological barrier for any inhaled particles, including drugs, to treat lung diseases. Elucidating the mechanism by which these drugs bind to and absorb into the LSM requires a molecular-level understanding of any drug-induced changes to the morphology, structure, and phase changes of the LSM.Molecular dynamics simulations have been used extensively to study the structure and dynamics of the LSM. The monolayer is usually simulated in at least two states: the compressed state, mimicking exhalation, and the expanded state, mimicking inhalation. In this chapter, we provide detailed instructions on how to set up, run, and analyze coarse-grained MD simulations to study the concentration-dependent effect of a sterol drug on the LSM, both in the expanded and compressed state.

Phosphatidylserine Exposed Lipid Bilayer Models for Understanding Cancer Cell Selectivity of Natural Compounds: A Molecular Dynamics Simulation Study

Development of drugs that are selectively toxic to cancer cells and safe to normal cells is crucial in cancer treatment. Evaluation of membrane permeability is a key metric for successful drug development. In this study, we have used in silico molecular models of lipid bilayers to explore the effect of phosphatidylserine (PS) exposure in cancer cells on membrane permeation of natural compounds Withaferin A (Wi-A), Withanone (Wi-N), Caffeic Acid Phenethyl Ester (CAPE) and Artepillin C (ARC). Molecular dynamics simulations were performed to compute permeability coefficients. The results indicated that the exposure of PS in cancer cell membranes facilitated the permeation of Wi-A, Wi-N and CAPE through a cancer cell membrane when compared to a normal cell membrane. In the case of ARC, PS exposure did not have a notable influence on its permeability coefficient. The presented data demonstrated the potential of PS exposure-based models for studying cancer cell selectivity of drugs.

Lipids-Fluorophores Interactions Probed by Combined Nonlinear Polarized Microscopy

Studying the structural dynamics of lipid membranes requires methods that can address both microscopic and macroscopic characteristics. Fluorescence imaging is part of the most used techniques to study membrane properties in various systems from artificial membranes to cells: It benefits from a high sensitivity to local properties such as polarity and molecular orientational order, with a high spatial resolution down to the single-molecule level. The influence of embedded fluorescent lipid probes on the lipid membrane molecules is however poorly known and relies most often on molecular dynamics simulations, due to the challenges faced by experimental approaches to address the molecular-scale dimension of this question. In this work we develop an optical microscopy imaging method to probe the effect of fluorophores embedded in the membrane as lipid probes, on their lipid environment, with a lateral resolution of a few hundreds of nanometers. We combine polarized-nonlinear microscopy contrasts that can independently address the lipid probe, by polarized two-photon fluorescence, and the membrane lipids, by polarized coherent Raman scattering. Using trimethylamino derivative 1-(4-trimethylammonium-phenyl)-6-phenyl-1,3,5-hexatriene (TMA-DPH) and di-8-butyl-amino-naphthyl-ethylene-pyridinium-propyl-sulfonate (di-8-ANEPPS) as model probes, we show that both probes tend to induce an orientational disorder of their surrounding lipid CH-bonds in 1,2-dipalmitoylphosphatidylcholine (DPPC) lipids environments, while there is no noticeable effect in more disordered 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipid membranes.

Water/oil interfacial tension reduction - an interfacial entropy driven process

The interfacial tension (IFT) of a fluid-fluid interface plays an important role in a wide range of applications and processes. When low IFT is desired, surface active compounds (e.g. surfactants) can be added to the system. Numerous attempts have been made to relate changes in IFT arising from such compounds to the specific nature of the interface. However, the IFT is controlled by an interplay of factors such as temperature and molecular structure of surface-active compounds, which make it difficult to predict IFT as those conditions change. In this study, we present the results from molecular dynamics simulations revealing the specific role surfactants play in IFT. We find that, in addition to reducing direct contact between the two fluids, surfactants serve to increase the disorder at the interface (related to interfacial entropy) and consequently reduce the water/oil IFT, especially when surfactants are present at high surface density. Our results suggest that surfactants that yield more disordered interfacial films (e.g. with flexible and/or unsaturated tails) reduce the water/oil IFT more effectively than surfactants which yield highly ordered interfacial films. Our results shed light on some of the factors that control IFT and could have important practical implications in industrial applications such as the design of cosmetics, food products, and detergents.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.