Shedding light on the structural properties of lipid bilayers using molecular dynamics simulation: a review study

Lateral diffusion of ions near membrane surface

Biological membranes isolate living cells from their environment, while allowing selective molecular transport between the inner and outer realms. For example, Na+ and K+ permeability through ionic channels contributes to neural conduction. Whether the ionic currents arise directly from cations in the bulk, or from the interface, is currently unclear. There are only scant results concerning lateral diffusion of ions on aquated membrane surfaces (and strong belief that this occurs through binding to a diffusing lipid). We performed classical molecular dynamics (MD) simulations of monovalent ions, Na+, K+, and Cl-, near the surface of the zwitterionic palmitoyl-oleoyl-phosphatidylcholine (POPC) membrane. Realistic force-fields for lipids (Amber's Lipid17 and Lipid21) and water (TIP4P-Ew) are tested for the mass and charge densities and the electrostatic potential across the membrane. These calculations reveal that the chloride can bind to the choline moiety through an intervening water molecule by forming a CH⋯OH hydrogen bond, while cations bind to both the phosphatic and carbonyl oxygens of phosphatidylcholine moieties. Upon transitioning from the bulk to the interface, a cation sheds some of its hydration water, which are replaced by headgroup atoms. Notably, an interfacial cation can bind 1-4 headgroup atoms, which is a key to understanding its surface hopping mechanism. We find that cation binding to three headgroup atoms immobilizes it, while binding to four energizes it. Consequently, the lateral cation diffusion rate is only 15-25 times slower than in the bulk, and 4-5 times faster than lipid self-diffusion. K+ diffusion is notably more anomalous than Na+, switching from sub- to super-diffusion after about 2 ns.

Molecular simulation on the interaction between trehalose and asymmetric lipid bilayer mimicking the membrane of human red blood cells

Trehalose is widely acknowledged for its ability to stabilize plasma membranes during dehydration. However, the exact mechanism by which trehalose interacts with lipid bilayers remains presently unclear. In this study, we conducted atomistic molecular dynamic simulations on asymmetric model bilayers that mimic the membrane of human red blood cells at various trehalose and water contents. We considered three different hydration levels mimicking the full hydration to desiccation scenarios. Results indicate that the asymmetric distribution of lipids did not significantly influence the computed structural characteristics at full and low hydration. At dehydration, however, the order parameter obtained from the symmetric bilayer is significantly higher compared to those obtained from asymmetric ones. Analysis of hydrogen bonds revealed that the protective ability of trehalose is well described by the water replacement hypothesis at full and low hydration, while at dehydration other interaction mechanisms associated with trehalose exclusion from the bilayer may involve. In addition, we found that trehalose exclusion is not attributed to sugar saturation but rather to the reduction in hydration levels. It can be concluded that the protective effect of trehalose is not only related to the hydration level of the bilayer, but also closely tied to the asymmetric distribution of lipids within each leaflet.

Comparative analysis of whole cell-derived vesicular delivery systems for photodynamic therapy of extrahepatic cholangiocarcinoma

This first-in-its-class proof-of-concept study explored the use of bionanovesicles for the delivery of photosensitizer into cultured cholangiocarcinoma cells and subsequent treatment by photodynamic therapy (PDT). Two types of bionanovesicles were prepared: cellular vesicles (CVs) were fabricated by sonication-mediated nanosizing of cholangiocarcinoma (TFK-1) cells, whereas cell membrane vesicles (CMVs) were produced by TFK-1 cell and organelle membrane isolation and subsequent nanovesicularization by sonication. The bionanovesicles were loaded with zinc phthalocyanine (ZnPC). The CVs and CMVs were characterized (size, polydispersity index, zeta potential, stability, ZnPC encapsulation efficiency, spectral properties) and assayed for tumor (TFK-1) cell association and uptake (flow cytometry, confocal microscopy), intracellular ZnPC distribution (confocal microscopy), dark toxicity (MTS assay), and PDT efficacy (MTS assay). The mean ± SD diameter, polydispersity index, and zeta potential were 134 ± 1 nm, -16.1 ± 0.9, and 0.220 ± 0.013, respectively, for CVs and 172 ± 3 nm, -16.4 ± 1.1, and 0.167 ± 0.022, respectively, for CMVs. Cold storage for 1 wk and incorporation of ZnPC increased bionanovesicular diameter slightly but size remained within the recommended range for in vivo application (136-220 nm). ZnPC was incorporated into CVs and CMVs at an optimal photosensitizer:lipid molar ratio of 0.006 and 0.01, respectively. Both bionanovesicles were avidly taken up by TFK-1 cells, resulting in homogenous intracellular ZnPC dispersion. Photosensitization of TFK-1 cells did not cause dark toxicity, while illumination at 671 nm (35.3 J/cm2) produced LC50 values of 1.11 μM (CVs) and 0.51 μM (CMVs) at 24 h post-PDT, which is superior to most LC50 values generated in tumor cells photosensitized with liposomal ZnPC. In conclusion, CVs and CMVs constitute a potent photosensitizer platform with no inherent cytotoxicity and high PDT efficacy in vitro.

Modulation of Phospholipase A2 Membrane Activity by Anti-inflammatory Drugs

The phospholipase A2 (PLA2) superfamily consists of lipolytic enzymes that hydrolyze specific cell membrane phospholipids and have long been considered a central hub of biosynthetic pathways, where their lipid metabolites exert a variety of physiological roles. A misregulated PLA2 activity is associated with mainly inflammatory-derived pathologies and thus has shown relevant therapeutic potential. Many natural and synthetic anti-inflammatory drugs (AIDs) have been proposed as direct modulators of PLA2 activity. However, despite the specific chemical properties that these drugs share in common, little is known about the indirect modulation able to finely tune membrane structural changes at the precise lipid-binding site. Here, we use a novel experimental strategy based on differential scanning calorimetry to systematically study the structural properties of lipid membrane systems during PLA2 cleavage and under the influence of several AIDs. For a better understanding of the AIDs-membrane interaction, we present a comprehensive and comparative set of molecular dynamics (MD) simulations. Our thermodynamic results clearly demonstrate that PLA2 cleavage is hindered by those AIDs that significantly reduce the lipid membrane cooperativity, while the rest of the AIDs oppositely tend to catalyze PLA2 activity to different extents. On the other hand, our MD simulations support experimental results by providing atomistic details on the binding, insertion, and dynamics of each AID on a pure lipid system; the drug efficacy to impact membrane cooperativity is related to the lipid order perturbation. This work suggests a membrane-based mechanism of action for diverse AIDs against PLA2 activity and provides relevant clues that must be considered in its modulation.

The antibacterial properties of branched peptides based on poly(l-arginine): In vitro antibacterial evaluation and molecular dynamic simulations

The emergence of bacterial strains resistant to antibiotics is a major issue in the medical field. Antimicrobial peptides are widely studied as they do not generate as much resistant bacterial strains as conventional antibiotics and present a broad range of activity. Among them, the homopolypeptide poly(l-arginine) presents promising antibacterial properties, especially in the perspective of its use in biomaterials. Linear poly(l-arginine) has been extensively studied but the impact of its 3D structure remains unknown. In this study, the antibacterial properties of newly synthesized branched poly(l-arginine) peptides, belonging to the family of multiple antigenic peptides, are evaluated. First, in vitro activities of the peptides shows that branched poly(l-arginine) is more efficient than linear poly(l-arginine) containing the same number of arginine residues. Surprisingly, peptides with more arms and more residues are not the most effective. To better understand these unexpected results, interactions between these peptides and the membranes of Gram positive and Gram negative bacteria are simulated thanks to molecular dynamic. It is observed that the bacterial membrane is more distorted by the branched structure than by the linear one and by peptides containing smaller arms. This mechanism of action is in full agreement with in vitro results and suggest that our simulations form a robust model to evaluate peptide efficiency towards pathogenic bacteria.

The intricate link between membrane lipid structure and composition and membrane structural properties in bacterial membranes

It is now evident that the cell manipulates lipid composition to regulate different processes such as membrane protein insertion, assembly and function. Moreover, changes in membrane structure and properties, lipid homeostasis during growth and differentiation with associated changes in cell size and shape, and responses to external stress have been related to drug resistance across mammalian species and a range of microorganisms. While it is well known that the biomembrane is a fluid self-assembled nanostructure, the link between the lipid components and the structural properties of the lipid bilayer are not well understood. This perspective aims to address this topic with a view to a more detailed understanding of the factors that regulate bilayer structure and flexibility. We describe a selection of recent studies that address the dynamic nature of bacterial lipid diversity and membrane properties in response to stress conditions. This emerging area has important implications for a broad range of cellular processes and may open new avenues of drug design for selective cell targeting.

Effect of Glycone Diversity on the Interaction of Triterpenoid Saponins and Lipid Bilayers

Triterpenoid saponins are organic compounds widely available in the plant kingdom. These molecules have received extensive attention due to their antibacterial activity against both Gram-negative and Gram-positive bacteria. Recent studies identified the antibacterial activity of saponins closely relates to their interaction with bacterial membrane lipids; however, molecular details of this interaction remain unclear. Increased understanding of the mechanisms to disrupt bacterial lipid bilayers can help to mitigate development of antibiotic resistance. Here, we examined the effect of chemical structure and deprotonation states of saponin on its interaction with a bacterial membrane model using molecular dynamics simulations. We run multiple simulations with a ternary lipid mixture of POPE/POPG/DPPG (80/15/5 mol %) and different saponin molecules. While all saponin structures can permanently bind the membrane, their location and orientation inside the bilayer depend on the sugar chains attached to their backbone. Similarly, cluster formation and stability also depend on the chemical structure of the saponin molecule. Deprotonation site affects interactions with the bilayer by modulating hydrophilicity of the molecules. At the low concentrations simulated in this work, there is no statistically significant change in the membrane properties upon saponin(s) binding, but the molecules do preferentially partition to POPE lipid environment.

Computational prediction of analog compounds of the membrane protein MCHR1 antagonists ALB-127158 and KRX-104130

Obesity, which is already pervasive throughout the world, endangers public health by raising the prevalence of metabolic disorders and making their treatment more difficult. The development of drugs to treat obesity is a focus of effort. Melanin concentrated hormone receptor 1 (MCHR1) is the target of some of these therapeutic possibilities since as increased levels of melanin concentrated hormone have been found in obesity models. Known MCHR1 antagonists include BMS-830216, GW-856464, NGD-4715, ALB-127158, and AMG 076, but many have failed phase-I clinical studies. As a potential treatment for cardiotoxicity, KRX-104130 has only recently been identified. As MCH system is potentially effective target for treatment of obesity, in silico research into interaction between MCHR1 and its antagonists at molecular level was the primary goal of this study. Analogues ALB-127158 and KRX-104130 were screened among the RealEnamine library. The complexes obtained by molecular docking were embedded in mimics brain-cell membrane and simulated for 540 ns, and then MM-GBSA were calculated with MMPBSA.py. With all these computational studies, similar or different aspects of selected analogous compounds to ALB-127158 and KRX-104130 were investigated. The specificity of this study was that it analyzed MCHR1 protein as embedded in membrane. It was concluded that KRX-104130's analogue Z1922310273 and ALB-127158's analogue PV-002757495233 did not cause a difference in terms of phospholipid membrane properties. In addition, all ligands remained stable in putative binding site. It has been suggested that PV-002757495233 and Z1922310273 compounds can be evaluated as MCHR1 antagonists when all these outputs are considered in melting pots.

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.

Understanding the Mechanical Properties of Ultradeformable Liposomes Using Molecular Dynamics Simulations

Improving drug delivery efficiency to solid tumor sites is a central challenge in anticancer therapeutic research. Our previous experimental study (Guo et al., Nat. Commun. 2018, 9, 130) showed that soft, elastic liposomes had increased uptake and accumulation in cancer cells and tumors in vitro and in vivo respectively, relative to rigid particles. As a first step toward understanding how liposomes' molecular structure and composition modulates their elasticity, we performed all-atom and coarse-grained classical molecular dynamics (MD) simulations of lipid bilayers formed by mixing a long-tailed unsaturated phospholipid with a short-tailed saturated lipid with the same headgroup. The former types of phospholipids considered were 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dipalmitoleoyl-sn-glycero-3-phosphocholine (termed here DPMPC). The shorter saturated lipids examined were 1,2-diheptanoyl-sn-glycero-3-phosphocholine (DHPC), 1,2-didecanoyl-sn-glycero-3-phosphocholine (DDPC), 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). Several lipid concentrations and surface tensions were considered. Our results show that DOPC or DPMPC systems having 25-35 mol % of the shortest lipids DHPC or DDPC are the least rigid, having area compressibility moduli KA that are ∼10% smaller than the values observed in pure DOPC or DPMPC bilayers. These results agree with experimental measurements of the stretching modulus and lysis tension in liposomes with the same compositions. These mixed systems also have lower areas per lipid and form more uneven x-y interfaces with water, the tails of both primary and secondary lipids are more disordered, and the terminal methyl groups in the tails of the long lipid DOPC or DPMPC wriggle more in the vertical direction, compared to pure DOPC or DPMPC bilayers or their mixtures with the longer saturated lipid DLPC or DMPC. These observations confirm our hypothesis that adding increasing concentrations of the short unsaturated lipid DHPC or DDPC to DOPC or DPMPC bilayers alters lipid packing and thus makes the resulting liposomes more elastic and less rigid. No formation of lipid nanodomains was noted in our simulations, and no clear trends were observed in the lateral diffusivities of the lipids as the concentration, type of secondary lipid, and surface tension were varied.

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.

Concentration-dependent mechanism of the binding behavior of ibuprofen to the cell membrane: A molecular dynamic simulation study

Ibuprofen is a commonly used drug for treating headaches, pain, and fever. The lipid bilayer is the primary and most important interface for drugs to interact with biological systems. However, the molecular interactions between ibuprofen and the cell membrane are not well understood. Our findings suggest that the interactions between ibuprofen and the bilayer involve multiple steps and depend on the concentration of the drug. At low concentrations of ibuprofen, it can bind to the surface of the lipid bilayer. The electrostatic and vdW energies of IBU-lipid at 0 ns of the simulation were -22.5 ± 3.2 and -5.9 ± 1.2 kj.mol-1 Fig. 2. In the following, the vdW energy of the IBU-lipid was increased by around -134.6 ± 3.7 kj.mol-1 whereas the electrostatic energy of the IBU-lipid was significantly decreased. This binding is facilitated by electrostatic and vdW interactions between ibuprofen and the head group of lipids. In the second step, ibuprofen is inserted into the lipid bilayer and positioned at the interface between the bilayer and the aqueous phase. In high concentrations of ibuprofen, it moved to the central region of the lipid bilayer. At this concentration, the physical and structural properties of the cell membrane change significantly. Results from the radial distribution function analysis indicate that at low concentrations, ibuprofen molecules are situated close to the head groups of phosphate groups. However, at high concentrations of ibuprofen, these molecules move to the inner side of the lipid bilayer. In addition, our findings indicate that at low concentrations of ibuprofen, these molecules did not significantly alter the physical properties of the cell membrane. In contrast, at high concentrations of ibuprofen, the physical parameters of the hydrocarbon tails, such as thickness, fluidity, and order, changed dramatically. APL parameter for POPC membrane increased slightly to 0.60 and 0.63 nm2 in the presence of low and high concentrations of ibuprofen molecules. The three-step interaction between ibuprofen and the lipid bilayer involves several events, such as the movement of ibuprofen molecules towards the central region of the lipid bilayer and the deformation and alteration of the structural and stability properties of the cell membrane. These effects are observed only at high concentrations of ibuprofen. It appears that the side effects of ibuprofen overdose are related to changes in the properties of the cell membrane and, subsequently, the function of membrane-anchored target proteins.

Martini 3 Coarse-Grained Force Field for Cholesterol

Cholesterol plays a crucial role in biomembranes by regulating various properties, such as fluidity, rigidity, permeability, and organization of lipid bilayers. The latest version of the Martini model, Martini 3, offers significant improvements in interaction balance, molecular packing, and inclusion of new bead types and sizes. However, the release of the new model resulted in the need to reparameterize many core molecules, including cholesterol. Here, we describe the development and validation of a Martini 3 cholesterol model, addressing issues related to its bonded setup, shape, volume, and hydrophobicity. The proposed model mitigates some limitations of its Martini 2 predecessor while maintaining or improving the overall behavior.

Molecular simulation of the rheological properties and shear thinning principles of supramolecular drilling fluids at different burial depths

In order to investigate the rheological properties and shear thinning principles of supramolecular drilling fluids, the salt-responsive supramolecular ionomer polymers with different components were designed and the change in shear viscosity of supramolecular polymer drilling fluid system with shear rate was studied using the molecular dynamics simulation method. The result indicated that the ionic supramolecular polymer drilling fluid system exhibits better self-assembly performance than the nonionic acrylamide drilling fluid system. Moreover, the drilling fluid system exhibits the best rheological properties and self-assembly performance when the feeding ratios of the three monomers in the two polymers are m : n : o = 5 : 90 : 5 and m : n : o = 30 : 40 : 30, respectively. The shear viscosity recovery rate of the #3 ionic supramolecular polymer drilling fluid system at different burial depths (1-5 km) is >87%, where the shear viscosity is mainly determined at ambient pressure. The shear thinning phenomenon of the supramolecular polymer drilling fluid system occurs because of the combined effect of the polymer molecular orientation and entanglement structure. When the shear rate is above a critical value, the polymer molecules are oriented along the flow field direction, decreasing the shear viscosity. However, when the shear rate is very high, the entanglement structure of the molecules is opened and the mesh structure of the fluids is disrupted, decreasing the shear viscosity of the drilling fluid.

Potential of liposomes and lipid membranes for the separation of β-blockers by capillary electromigration and liquid chromatographic techniques

β-Blockers belong to a frequently used class of drugs primarily used to treat heart and circulatory conditions. Here we describe the use of lipid vesicles and liposomes as cell membrane biomimicking models in capillary electromigration (CE) and liquid chromatography (LC) techniques for the investigation of interactions between lipid membranes and β-blockers. In addition to liposomes, the use of commercial intravenous lipid emulsions, and their interactions with β-blockers are also discussed. Different CE and LC instrumental techniques designed for these purposes are introduced. Other methodologies for studying interactions between β-blockers and lipid membranes are also briefly discussed, and the different methodologies are compared. The aim is to give the reader a good overview on the status of the use of liposomes and lipids in CE and LC for studying β-blocker interactions.

Effects of Nitro-Oxidative Stress on Biomolecules: Part 1-Non-Reactive Molecular Dynamics Simulations

Plasma medicine, or the biomedical application of cold atmospheric plasma (CAP), is an expanding field within plasma research. CAP has demonstrated remarkable versatility in diverse biological applications, including cancer treatment, wound healing, microorganism inactivation, and skin disease therapy. However, the precise mechanisms underlying the effects of CAP remain incompletely understood. The therapeutic effects of CAP are largely attributed to the generation of reactive oxygen and nitrogen species (RONS), which play a crucial role in the biological responses induced by CAP. Specifically, RONS produced during CAP treatment have the ability to chemically modify cell membranes and membrane proteins, causing nitro-oxidative stress, thereby leading to changes in membrane permeability and disruption of cellular processes. To gain atomic-level insights into these interactions, non-reactive molecular dynamics (MD) simulations have emerged as a valuable tool. These simulations facilitate the examination of larger-scale system dynamics, including protein-protein and protein-membrane interactions. In this comprehensive review, we focus on the applications of non-reactive MD simulations in studying the effects of CAP on cellular components and interactions at the atomic level, providing a detailed overview of the potential of CAP in medicine. We also review the results of other MD studies that are not related to plasma medicine but explore the effects of nitro-oxidative stress on cellular components and are therefore important for a broader understanding of the underlying processes.

Building Asymmetric Lipid Bilayers for Molecular Dynamics Simulations: What Methods Exist and How to Choose One?

The compositional asymmetry of biological membranes has attracted significant attention over the last decade. Harboring more differences from symmetric membranes than previously appreciated, asymmetric bilayers have proven quite challenging to study with familiar concepts and techniques, leaving many unanswered questions about the reach of the asymmetry effects. One particular area of active research is the computational investigation of composition- and number-asymmetric lipid bilayers with molecular dynamics (MD) simulations. Offering a high level of detail into the organization and properties of the simulated systems, MD has emerged as an indispensable tool in the study of membrane asymmetry. However, the realization that results depend heavily on the protocol used for constructing the asymmetric bilayer models has sparked an ongoing debate about how to choose the most appropriate approach. Here we discuss the underlying source of the discrepant results and review the existing methods for creating asymmetric bilayers for MD simulations. Considering the available data, we argue that each method is well suited for specific applications and hence there is no single best approach. Instead, the choice of a construction protocol-and consequently, its perceived accuracy-must be based primarily on the scientific question that the simulations are designed to address.

Curvature Matters: Modeling Calcium Binding to Neutral and Anionic Phospholipid Bilayers

In this work, the influence of membrane curvature on the Ca²⁺ binding to phospholipid bilayers is investigated by means of molecular dynamics simulations. In particular, we compared Ca²⁺ binding to flat, elastically buckled, or uniformly bent zwitterionic and anionic phospholipid bilayers. We demonstrate that Ca²⁺ ions bind preferably to the concave membrane surfaces in both types of bilayers. We also show that the membrane curvature leads to pronounced changes in Ca²⁺ binding including differences in free ion concentrations, lipid coordination distributions, and the patterns of ion binding to different chemical groups of lipids. Moreover, these effects differ substantially for the concave and convex membrane monolayers. Comparison between force fields with either full or scaled charges indicates that charge scaling results in reduction of the Ca²⁺ binding to curved phosphatidylserine bilayers, while for phosphatidylcholine membranes, calcium binds only weakly for both force fields.

APL@voro-interactive visualization and analysis of cell membrane simulations

SUMMARY

Molecular dynamics (MD) simulations of cell membranes allow for a better understanding of complex processes such as changing membrane dynamics, lipid rafts and the incorporation/passing of macromolecules into/through membranes. To explore and understand cell membrane compositions, dynamics and processes, visual analytics can help to interpret MD simulation data. APL@Voro is a software for the interactive visualization and analysis of cell membrane simulations. Here, we present the new APL@Voro, which has been continuously developed since its initial release in 2013. We discuss newly implemented algorithms, methodologies and features, such as the interactive comparison of related simulations and methods to assign lipids to either the upper or lower leaflet.

AVAILABILITY AND IMPLEMENTATION

The current open-source version of APL@Voro can be downloaded from http://aplvoro.com.

The effect of rhamnolipids on fungal membrane models as described by their interactions with phospholipids and sterols: An in silico study

Introduction: Rhamnolipids (RLs) are secondary metabolites naturally produced by bacteria of the genera Pseudomonas and Burkholderia with biosurfactant properties. A specific interest raised from their potential as biocontrol agents for crop culture protection in regard to direct antifungal and elicitor activities. As for other amphiphilic compounds, a direct interaction with membrane lipids has been suggested as the key feature for the perception and subsequent activity of RLs. Methods: Molecular Dynamics (MD) simulations are used in this work to provide an atomistic description of their interactions with different membranous lipids and focusing on their antifungal properties. Results and discussion: Our results suggest the insertion of RLs into the modelled bilayers just below the plane drawn by lipid phosphate groups, a placement that is effective in promoting significant membrane fluidification of the hydrophobic core. This localization is promoted by the formation of ionic bonds between the carboxylate group of RLs and the amino group of the phosphatidylethanolamine (PE) or phosphatidylserine (PS) headgroups. Moreover, RL acyl chains adhere to the ergosterol structure, forming a significantly higher number of van der Waals contact with respect to what is observed for phospholipid acyl chains. All these interactions might be essential for the membranotropic-driven biological actions of RLs.

Effects of the RNA-Polymerase Inhibitors Remdesivir and Favipiravir on the Structure of Lipid Bilayers-An MD Study

The structure and dynamics of membranes are crucial to ensure the proper functioning of cells. There are some compounds used in therapeutics that show nonspecific interactions with membranes in addition to their specific molecular target. Among them, two compounds recently used in therapeutics against COVID-19, remdesivir and favipiravir, were subjected to molecular dynamics simulation assays. In these, we demonstrated that the compounds can spontaneously bind to model lipid membranes in the presence or absence of cholesterol. These findings correlate with the corresponding experimental results recently reported by our group. In conclusion, insertion of the compounds into the membrane is observed, with a mean position close to the phospholipid head groups.

The dynamical Matryoshka model: 3. Diffusive nature of the atomic motions contained in a new dynamical model for deciphering local lipid dynamics

In accompanying papers [Bicout et al., BioRxiv https://doi.org/10.1101/2021.09.21.461198 (2021); Cissé et al., BioRxiv https://doi.org/10.1101/2022.03.30.486370 (2022)], a new model called Matryoshka model has been proposed to describe the geometry of atomic motions in phospholipid molecules in bilayers and multilamellar vesicles based on their quasielastic neutron scattering (QENS) spectra. Here, in order to characterize the relaxational aspects of this model, the energy widths of the QENS spectra of the samples were analyzed first in a model-free way. The spectra were decomposed into three Lorentzian functions, which are classified as slow, intermediate, and fast motions depending on their widths. The analysis provides the diffusion coefficients, residence times, and geometrical parameters for the three classes of motions. The results corroborate the parameter values such as the amplitudes and the mobile fractions of atomic motions obtained by the application of the Matryoshka model to the same samples. Since the current analysis was carried out independently of the development of the Matryoshka model, the present results enhance the validity of the model. The model will serve as a powerful tool to decipher the dynamics of lipid molecules not only in model systems, but also in more complex systems such as mixtures of different kinds of lipids or natural cell membranes.

Dissecting the mechanisms of environment sensitivity of smart probes for quantitative assessment of membrane properties

The plasma membrane, as a highly complex cell organelle, serves as a crucial platform for a multitude of cellular processes. Its collective biophysical properties are largely determined by the structural diversity of the different lipid species it accommodates. Therefore, a detailed investigation of biophysical properties of the plasma membrane is of utmost importance for a comprehensive understanding of biological processes occurring therein. During the past two decades, several environment-sensitive probes have been developed and become popular tools to investigate membrane properties. Although these probes are assumed to report on membrane order in similar ways, their individual mechanisms remain to be elucidated. In this study, using model membrane systems, we characterized the probes Pro12A, NR12S and NR12A in depth and examined their sensitivity to parameters with potential biological implications, such as the degree of lipid saturation, double bond position and configuration (cis versus trans), phospholipid headgroup and cholesterol content. Applying spectral imaging together with atomistic molecular dynamics simulations and time-dependent fluorescent shift analyses, we unravelled individual sensitivities of these probes to different biophysical properties, their distinct localizations and specific relaxation processes in membranes. Overall, Pro12A, NR12S and NR12A serve together as a toolbox with a wide range of applications allowing to select the most appropriate probe for each specific research question.

Copper Phosphonate Lamella Intermediates Control the Shape of Colloidal Copper Nanocrystals

Understanding the structure and behavior of intermediates in chemical reactions is the key to developing greater control over the reaction outcome. This principle is particularly important in the synthesis of metal nanocrystals (NCs), where the reduction, nucleation, and growth of the reaction intermediates will determine the final size and shape of the product. The shape of metal NCs plays a major role in determining their catalytic, photochemical, and electronic properties and, thus, the potential applications of the material. In this work, we demonstrate that layered coordination polymers, called lamellae, are reaction intermediates in Cu NC synthesis. Importantly, we discover that the lamella structure can be fine-tuned using organic ligands of different lengths and that these structural changes control the shape of the final NC. Specifically, we show that short-chain phosphonate ligands generate lamellae that are stable enough at the reaction temperature to facilitate the growth of Cu nuclei into anisotropic Cu NCs, being primarily triangular plates. In contrast, lamellae formed from long-chain ligands lose their structure and form spherical Cu NCs. The synthetic approach presented here provides a versatile tool for the future development of metal NCs, including other anisotropic structures.

Impact of Surface Polarity on Lipid Assembly under Spatial Confinement

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

Ultralong recovery time in nanosecond electroporation systems enabled by orientational-disordering processes

The growing importance of applications based on molecular medicine and genetic engineering is driving the need to develop high-performance electroporation technologies. The electroporation phenomenon involves disruption of the cell for increasing membrane permeability. Although there is a multitude of research focused on exploring new electroporation techniques, the engineering of programming schemes suitable for these electroporation methods remains a challenge. Nanosecond stimulations could be promising candidates for these techniques owing to their ability to generate a wide range of biological responses. Here we control the membrane permeabilization of cancer cells using different numbers of electric-field pulses through orientational disordering effects. We then report our exploration of a few-volt nanosecond alternating-current (AC) stimulation method with an increased number of pulses for developing electroporation systems. A recovery time of ∼720 min was achieved, which is above the average of ∼76 min for existing electroporation methods using medium cell populations, as well as a previously unreported increased conductance with an increase in the number of pulses using weak bias amplitudes. All-atom molecular dynamics (MD) simulations reveal the orientation-disordering-facilitated increase in the degree of permeabilization. These findings highlight the potential of few-volt nanosecond AC-stimulation with an increased number of pulse strategies for the development of next-generation low-power electroporation systems.

Unravelling hierarchical levels of structure in lipid membranes

•

A set of hierarchical levels of structure can be proposed for lipid bilayers.

•

The composition of a lipid bilayer is identified as a fluid version of primary structure.

•

The interaction between leaflets is taken as the secondary structure of lipid bilayers.

•

A method to identify membrane domains and their interaction is proposed as a third level of structure.

•

The highly specific lipid perturbation around embedded macromolecules is taken as the quaternary level of structure.

In analogy with the hierarchical levels typically used to describe the structure of nucleic acids or proteins and keeping in mind that lipid bilayers are not just mere envelopers for biological material but directly responsible for many important functions of life, it is discussed here how membrane models can also be interpreted in terms of different hierarchies in their structure. Namely, lipid composition, interaction between leaflets, existence and interaction of domains arising from the coordinate behavior of lipids and their properties, plus the manifest and specific perturbation of the lipid organization around macromolecules embedded in a membrane are hereby used to define the primary, secondary, tertiary and quaternary structures, respectively. Molecular Dynamics simulations are used to illustrate this proposal. Alternative levels of organization and methods to define domains can be proposed but the final aim is to highlight the paradigm arising from this description which is expected to have significant consequences on deciphering the underlying factors governing membranes and their interactions with other molecules.

Revealing local molecular distribution, orientation, phase separation, and formation of domains in artificial lipid layers: Towards comprehensive characterization of biological membranes

Lipids, together with molecules such as DNA and proteins, are one of the most relevant systems responsible for the existence of life. Selected lipids are able to assembly into various organized structures, such as lipid membranes. The unique properties of lipid membranes determine their complex functions, not only to separate biological environments, but also to participate in regulatory functions, absorption of nutrients, cell-cell communication, endocytosis, cell signaling, and many others. Despite numerous scientific efforts, still little is known about the reason underlying the variability within lipid membranes, and its biochemical significance. In this review, we discuss the structural complexity of lipid membranes, as well as the importance to simplify studied systems in order to understand phenomena occurring in natural, complex membranes. Such systems require a model interface to be analyzed. Therefore, here we focused on analytical studies of artificial systems at various interfaces. The molecular structure of lipid membranes, specifically the nanometric thickens of molecular bilayer, limits in a major extent the choice of highly sensitive methods suitable to study such structures. Therefore, we focused on methods that combine high sensitivity, and/or chemical selectivity, and/or nanometric spatial resolution, such as atomic force microscopy, nanospectroscopy (tip-enhanced Raman spectroscopy, infrared nanospectroscopy), phase modulation infrared reflection-absorption spectroscopy, sum-frequency generation spectroscopy. We summarized experimental and theoretical approaches providing information about molecular structure and composition, lipid spatial distribution (phase separation), organization (domain shape, molecular orientation) of lipid membranes, and real-time visualization of the influence of various molecules (proteins, drugs) on their integrity. An integral part of this review discusses the latest achievements in the field of lipid layer-based biosensors.

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

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

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

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

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.

Ceramide-1-phosphate transfer protein promotes sphingolipid reorientation needed for binding during membrane interaction

Lipid transfer proteins acquire and release their lipid cargoes by interacting transiently with source and destination biomembranes. In the GlycoLipid Transfer Protein (GLTP) superfamily, the two-layer all-α-helical GLTP-fold defines proteins that specifically target sphingolipids (SLs) containing either sugar or phosphate headgroups via their conserved but evolutionarily-modified SL recognitions centers. Despite comprehensive structural insights provided by X-ray crystallography, the conformational dynamics associated with membrane interaction and SL uptake/release by GLTP superfamily members have remained unknown. Herein, we report insights gained from molecular dynamics (MD) simulations into the conformational dynamics that enable ceramide-1-phosphate transfer proteins (CPTPs) to acquire and deliver ceramide-1-phosphate (C1P) during interaction with 1-palmitoyl-2-oleoyl phosphatidylcholine bilayers. The focus on CPTP reflects this protein's involvement in regulating pro-inflammatory eicosanoid production and autophagy-dependent inflammasome assembly that drives interleukin (IL-1β and IL-18) production and release by surveillance cells. We found that membrane penetration by CPTP involved α-6 helix and the α-2 helix N-terminal region, was confined to one bilayer leaflet, and was relatively shallow. Large-scale dynamic conformational changes were minimal for CPTP during membrane interaction or C1P uptake except for the α-3/α-4 helices connecting loop, which is located near the membrane interface and interacts with certain phosphoinositide headgroups. Apart from functioning as a shallow membrane-docking element, α-6 helix was found to adeptly reorient membrane lipids to help guide C1P hydrocarbon chain insertion into the interior hydrophobic pocket of the SL binding site.These findings support a proposed 'hydrocarbon chain-first' mechanism for C1P uptake, in contrast to the 'lipid polar headgroup-first' uptake used by most lipid-transfer proteins.

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.

How Do Urea and Trimethylamine N-Oxide Influence the Dehydration-Induced Phase Transition of a Lipid Membrane?

Living organisms are often exposed to extreme dehydration, which is detrimental to the structure and function of the cell membrane. The lipid membrane undergoes fluid-to-gel phase transition due to dehydration and thus loses fluidity and functionality. To protect the fluid phase of the bilayer these organisms adopt several strategies. Enhanced production of small polar organic solutes (also called osmolytes) is one such strategy. Urea and trimethylamine N-oxide (TMAO) are two osmolytes found in different organisms combating osmotic stress. Previous experiments have found that both these osmolytes have strong effects on lipid membrane under different hydration conditions. Urea prevents the dehydration-induced phase transition of the lipid membrane by directly interacting with the lipids, while TMAO does not inhibit the phase transition. To provide atomistic insights, we have carried out all-atom molecular dynamics (MD) simulation of a lipid membrane under varying hydration levels and studied the effect of these osmolytes on different structural and dynamic properties of the membrane. This study suggests that urea significantly inhibits the dehydration-induced fluid-to-gel phase transition by strongly interacting with the lipid membrane via hydrogen bonds, which balances the reduced lipid hydration due to the decreasing water content. In contrast, TMAO is excluded from the membrane surface due to unfavorable interaction with the lipids. This induces further dehydration of the lipids which reinforces the fluid-to-gel phase transition. We have also studied the counteractive role of TMAO on the effect of urea on lipid membrane when both the osmolytes are present. TMAO draws some urea molecules out of the membrane and thereby reduces the effect of urea on the lipid membrane at lower hydration levels. This is similar to the counteraction of urea's deleterious effects on protein by TMAO. All these observations are consistent with the experimental results and thus provide deep molecular insights into the role of these osmolytes in protecting the fluid phase of the membrane, the key survival strategy against osmotic-stress-induced dehydration.

Iron coordination to pyochelin siderophore influences dynamics of FptA receptor from Pseudomonas aeruginosa: a molecular dynamics simulation study

FptA is a TonB-dependent transporter that permits the high affinity binding and transport of Fe(III)-pyochelin complex across the outer membrane of Pseudomonas aeruginosa. Molecular dynamics simulations were employed to FptA receptor and its complexes with pyochelin, and co-crystallized Fe(III)-pyochelin-ethanediol and Fe(III)-pyochelin-water embedded in dilauroyl phosphatidyl choline bilayer for the evaluation of their structural and dynamical properties. The evaluation of properties of the receptor bound to pyochelin molecule and Fe(III)-pyochelin complexes helped to figure out the iron coordination effect on the receptor properties. Moreover, comparison of these four simulation systems revealed further information on the dynamical changes occurred in extracellular loops, in particular loop-7 corresponding to the missing amino acid residues including the close-by loop-8 that was largely affected by the metal coordination to pyochelin. The binding of iron to pyochelin molecule affected the overall structure of the receptor therefore, evaluation fo the gyration radii and hydrogen bonding were evaluated as well as analysis of the pore size were also carried out to understand the effect of metal coordination on the dynamics of the helices which form a kind of translocation channel to transport the siderophore across the FptA protein into the periplasmic space. The properties of each component of the molecular systems were therefore observed to be perturbed by the incorporation of iron to the pyochelin molecule thus demonstrating that the bacteria use its receptor to abstract and transport iron from extracellular environment for its survival and that was made possible to understand at the molecular level through successful implementation of molecular dynamics simulations.

Modulation of Phospholipid Bilayer Properties by Simvastatin

Simvastatin (Zocor) is one of the most prescribed drugs for reducing high cholesterol. Although simvastatin is ingested in its inactive lactone form, it is converted to its active dihydroxyheptanoate form by carboxylesterases in the liver. The dihydroxyheptanoate form can also be converted back to its original lactone form. Unfortunately, some of the side effects associated with the intake of simvastatin and other lipophilic statins at higher doses include statin-associated myopathy (SAM) and, in more severe cases, kidney failure. While the cause of SAM is unknown, it is hypothesized that these side effects are dependent on the localization of statins in lipid bilayers and their impact on bilayer properties. In this work, we carry out all-atom molecular dynamics simulations on both the lactone and dihydroxyheptanoate forms of simvastatin (termed "SN" and "SA", respectively) with a pure 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer and a POPC/cholesterol (30 mol %) binary mixture as membrane models. Additional simulations were carried out with multiple simvastatin molecules to mimic in vitro conditions that produced pleiotropic effects. Both SN and SA spontaneously diffused into the lipid bilayer, and a longer simulation time of 4 μs was needed for the complete incorporation of multiple SAs into the bilayer. By constructing potential mean force and electron density profiles, we find that SN localizes deeper within the hydrophobic interior of the bilayer and that SA has a greater tendency to form hydrogen-bonding interactions with neighboring water molecules and lipid headgroups. For the pure POPC bilayer, both SN and SA increase membrane order, while membrane fluidity increases for the POPC/cholesterol bilayer.

On the Structure and Flip-flop of Free Docosahexaenoic Acid in a Model Human Brain Membrane

Among the omega-3 fatty acids, docosahexaenoic acid (DHA, sn22:6) is particularly vital in human brain cell membranes. There is considerable interest in DHA because low-level DHA has been associated with declined cognitive function and poor visual acuity. In this work, atomistic molecular dynamics simulations were used to investigate the effects of free protonated DHA (DHAP) in molar fractions of 0, 17, 30, and 38% in a realistic model of a healthy brain cell membrane comprising 26 lipid types. Numerous flip-flop events of DHAP were observed and categorized as successful or aborted. Novel use of the machine learning technique, density-based spatial clustering of applications with noise (DBSCAN), effectively identified flip-flop events by way of clustering. Our data show that increasing amounts of DHAP in the membrane disorder the bilayer packing, fluidize the membrane, and increase the rates of successful flip-flop from k = 0.2 μs^-1 (17% DHAP) to 0.8 μs^-1 (30% DHAP) and to 1.3 μs^-1 (38% DHAP). In addition, we also provided a detailed understanding of the flip-flop mechanism of DHAP across this complex membrane. Interestingly, we noted the role of hydrogen bonds in two distinct coordinated flip-flop phenomena between two DHAP molecules: double flip-flop and assisted flip-flop. Understanding the effects of various concentrations of DHAP on the dynamics within a lipid membrane and the resulting structural properties of the membrane are important when considering the use of DHAP as a dietary supplement or as a potential therapeutic.

Mutually Beneficial Combination of Molecular Dynamics Computer Simulations and Scattering Experiments

We showcase the combination of experimental neutron scattering data and molecular dynamics (MD) simulations for exemplary phospholipid membrane systems. Neutron and X-ray reflectometry and small-angle scattering measurements are determined by the scattering length density profile in real space, but it is not usually possible to retrieve this profile unambiguously from the data alone. MD simulations predict these density profiles, but they require experimental control. Both issues can be addressed simultaneously by cross-validating scattering data and MD results. The strengths and weaknesses of each technique are discussed in detail with the aim of optimizing the opportunities provided by this combination.

Membrane-coated 3D architectures for bottom-up synthetic biology

One of the great challenges of bottom-up synthetic biology is to recreate the cellular geometry and surface functionality required for biological reactions. Of particular interest are lipid membrane interfaces where many protein functions take place. However, cellular 3D geometries are often complex, and custom-shaping stable lipid membranes on relevant spatial scales in the micrometer range has been hard to accomplish reproducibly. Here, we use two-photon direct laser writing to 3D print microenvironments with length scales relevant to cellular processes and reactions. We formed lipid bilayers on the surfaces of these printed structures, and we evaluated multiple combinatorial scenarios, where physiologically relevant membrane compositions were generated on several different polymer surfaces. Functional dynamic protein systems were reconstituted in vitro and their self-organization was observed in response to the 3D geometry. This method proves very useful to template biological membranes with an additional spatial dimension, and thus allows a better understanding of protein function in relation to the complex morphology of cells and organelles.

Estimating localization of various statins within a POPC bilayer

As a class of drugs prescribed to heart disease patients, statins are among the most popular prescription drugs in the world. Over the years, statins have been shown to have beneficial effects on patients via pathways independent of their effect on cholesterol. These pleiotropic effects vary across the different statins, and a growing hypothesis is that they are related to the localization of the statins in and their effect on the membrane. In this study, we use molecular dynamics (MD) simulations with the CHARMM36 all-atom force field to investigate the localization of statins (atorvastatin, cerivastatin, lovastatin, and pravastatin) in a POPC bilayer and how they affect the acyl chain order parameters (S_CD), surface area per lipid (APL), and thicknesses of the bilayer. The data obtained from 500 ns simulations suggests that lovastatin is localized deepest in the membrane, mostly interacting with the hydrophobic core, cerivastatin is slightly closer to the bilayer/solvent interface than lovastatin and interacts with the headgroups via its dihydroxy acid group, and pravastatin is found closest to the bilayer/solvent interface, its hydrophobic rings interacting mostly with the region around the acyl's carbonyl and its dihydroxy acid interacting with the solvent and the headgroups. Consistent binding of atorvastatin to the bilayer is not observed during our simulation due to self-aggregation. The statins differentially alter the S_CD and APL and most of the bilayer thicknesses, but these effects are modest. Overall, as expected, the localization of statins seems to follow their hydrophilicity, and given previous data showing the relationship between statins' hydrophobicity and pleiotropic effects, one would expect statins that localize and interact with different regions of the membrane to have different effects. This research provides some important insight into statin localization in a simplified model of a cellular membrane.

Exploring the putative mechanism of allosteric modulations by mixed-action kappa/mu opioid receptor bitopic modulators

The modulation and selectivity mechanisms of seven mixed-action kappa opioid receptor (KOR)/mu opioid receptor (MOR) bitopic modulators were explored. Molecular modeling results indicated that the 'message' moiety of seven bitopic modulators shared the same binding mode with the orthosteric site of the KOR and MOR, whereas the 'address' moiety bound with different subdomains of the allosteric site of the KOR and MOR. The 'address' moiety of seven bitopic modulators bound to different subdomains of the allosteric site of the KOR and MOR may exhibit distinguishable allosteric modulations to the binding affinity and/or efficacy of the 'message' moiety. Moreover, the 3-hydroxy group on the phenolic moiety of the seven bitopic modulators induced selectivity to the KOR over the MOR.

Effect of TMAO on the Structure and Phase Transition of Lipid Membranes: Potential Role of TMAO in Stabilizing Cell Membranes under Osmotic Stress

Extremophiles adopt strategies to deal with different environmental stresses, some of which are severely damaging to their cell membrane. To combat high osmotic stress, deep-sea organisms synthesize osmolytes, small polar organic molecules, like trimethylamine-N-oxide (TMAO), and incorporate them in the cell. TMAO is known to protect cells from high osmotic or hydrostatic pressure. Several experimental and simulation studies have revealed the roles of such osmolytes on stabilizing proteins. In contrast, the effect of osmolytes on the lipid membrane is poorly understood and broadly debated. A recent experiment has found strong evidence of the possible role of TMAO in stabilizing lipid membranes. Using the molecular dynamics (MD) simulation technique, we have demonstrated the effect of TMAO on two saturated fully hydrated lipid membranes in their fluid and gel phases. We have captured the impact of TMAO's concentration on the membrane's structural properties along with the fluid/gel phase transition temperatures. On increasing the concentration of TMAO, we see a substantial increase in the packing density of the membrane (estimated by area, thickness, and volume) and enhancement in the orientational order of lipid molecules. Having repulsive interaction with the lipid head group, the TMAO molecules are expelled away from the membrane surface, which induces dehydration of the lipid head groups, increasing the packing density. The addition of TMAO also increases the fluid/gel phase transition temperature of the membrane. All of these results are in close agreement with the experimental observations. This study, therefore, provides a molecular-level understanding of how TMAO can influence the cell membrane of deep-sea organisms and help in combating the osmotic stress condition.

Tentative Peptide‒Lipid Bilayer Models Elucidating Molecular Behaviors and Interactions Driving Passive Cellular Uptake of Collagen-Derived Small Peptides

Collagen contains hydroxyproline (Hyp), which is a unique amino acid. Three collagen-derived small peptides (Gly-Pro-Hyp, Pro-Hyp, and Gly-Hyp) interacting across a lipid bilayer (POPC model membrane) for cellular uptakes of these collagen-derived small peptides were studied using accelerated molecular dynamics simulation. The ligands were investigated for their binding modes, hydrogen bonds in each coordinate frame, and mean square displacement (MSD) in the Z direction. The lipid bilayers were evaluated for mass and electron density profiles of the lipid molecules, surface area of the head groups, and root mean square deviation (RMSD). The simulation results show that hydrogen bonding between the small collagen peptides and plasma membrane plays a significant role in their internalization. The translocation of the small collagen peptides across the cell membranes was shown. Pro-Hyp laterally condensed the membrane, resulting in an increase in the bilayer thickness and rigidity. Perception regarding molecular behaviors of collagen-derived peptides within the cell membrane, including their interactions, provides the novel design of specific bioactive collagen peptides for their applications.

Cholesterol content in the membrane promotes key lipid-protein interactions in a pentameric serotonin-gated ion channel

Pentameric ligand-gated ion channels (pLGICs), embedded in the lipid membranes of nerve cells, mediate fast synaptic transmission and are major pharmaceutical targets. Because of their complexity and the limited knowledge of their structure, their working mechanisms have still to be fully unraveled at the molecular level. Over the past few years, evidence that the lipid membrane may modulate the function of membrane proteins, including pLGICs, has emerged. Here, we investigate, by means of molecular dynamics simulations, the behavior of the lipid membrane at the interface with the 5-HT_3A receptor (5-HT_3AR), a representative pLGIC which is the target of nausea-suppressant drugs, in a nonconductive state. Three lipid compositions are studied, spanning different concentrations of the phospholipids, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine, and of cholesterol, hence a range of viscosities. A variety of lipid interactions and persistent binding events to different parts of the receptor are revealed in the investigated models, providing snapshots of the dynamical environment at the membrane-receptor interface. Some of these events result in lipid intercalation within the transmembrane domain, and others reach out to protein key sections for signal transmission and receptor activation, such as the Cys-loop and the M2-M3 loop. In particular, phospholipids, with their long hydrophobic tails, play an important role in these interactions, potentially providing a bridge between these two structures. A higher cholesterol content appears to promote lipid persistent binding to the receptor.

Binding selectivity of inhibitors toward the first over the second bromodomain of BRD4: theoretical insights from free energy calculations and multiple short molecular dynamics simulations

Bromodomain-containing protein 4 (BRD4) plays an important role in mediating gene transcription involved in cancers and non-cancer diseases such as acute heart failure and inflammatory diseases. In this work, multiple short molecular dynamics (MSMD) simulations are integrated with a molecular mechanics generalized Born surface area (MM-GBSA) approach to decipher binding selectivity of three inhibitors 8NS, 82Y, and 837 toward two domains BD1 and BD2 of BRD4. The results demonstrate that the enthalpy effects play critical roles in selectivity identification of inhibitors toward BD1 and BD2, determining that 8NS has better selectivity toward BD2 than BD1, while 82Y and 837 more favorably bind to BD1 than BD2. A residue-based free-energy decomposition method was used to calculate an inhibitor-residue interaction spectrum and unveil contributions of separate residues to binding selectivity. The results identify six common residues, containing (P82, P375), (V87, V380), (L92, L385), (L94, L387), (N140, N433), and (I146, V439) individually belonging to (BD1, BD2) of BRD4, and yield a considerable binding difference of inhibitors to BD1 and BD2, suggesting that these residues play key roles in binding selectivity of inhibitors toward BD1 and BD2 of BRD4. Therefore, these results provide useful dynamics information and a structure affinity relationship for the development of highly selective inhibitors targeting BD1 and BD2 of BRD4.

Design, Screening, and Testing of Non-Rational Peptide Libraries with Antimicrobial Activity: In Silico and Experimental Approaches

One of the challenges of modern biotechnology is to find new routes to mitigate the resistance to conventional antibiotics. Antimicrobial peptides (AMPs) are an alternative type of biomolecules, naturally present in a wide variety of organisms, with the capacity to overcome the current microorganism resistance threat. Here, we reviewed our recent efforts to develop a new library of non-rationally produced AMPs that relies on bacterial genome inherent diversity and compared it with rationally designed libraries. Our approach is based on a four-stage workflow process that incorporates the interplay of recent developments in four major emerging technologies: artificial intelligence, molecular dynamics, surface-display in microorganisms, and microfluidics. Implementing this framework is challenging because to obtain reliable results, the in silico algorithms to search for candidate AMPs need to overcome issues of the state-of-the-art approaches that limit the possibilities for multi-space data distribution analyses in extremely large databases. We expect to tackle this challenge by using a recently developed classification algorithm based on deep learning models that rely on convolutional layers and gated recurrent units. This will be complemented by carefully tailored molecular dynamics simulations to elucidate specific interactions with lipid bilayers. Candidate AMPs will be recombinantly-expressed on the surface of microorganisms for further screening via different droplet-based microfluidic-based strategies to identify AMPs with the desired lytic abilities. We believe that the proposed approach opens opportunities for searching and screening bioactive peptides for other applications.

Nano-mechanical characterization of asymmetric DLPC/DSPC supported lipid bilayers

Asymmetric distribution of lipid molecules in the inner and outer leaflets of the plasma membrane is a common occurrence in the membrane formation. Such asymmetric arrangement is a crucial parameter to manipulate the properties of the cell membrane. It controls signal transduction, endocytosis, exocytosis in the cells. The artificial membrane is often used to study the lateral and transverse arrangement of the lipid molecules in place of the cell membrane. Nano-mechanical characterization of the model membrane helps to understand the mechanical stability of the lipid bilayer. The stability is sensitive to the variations in the lipid composition and their local organization. In this article, we present both topographical and nano-mechanical properties of lipid bilayer characterized by atomic force microscopy (AFM). The results show that the asymmetric lipid bilayer formation is an intrinsic character. We have selected a bi-component fluid-gel phase 1,2-dilauroyl-sn-glycero-3-phosphocholine:1,2-disteroyl-sn-glycero-3-phosphocholine (DLPC: DSPC) system for our studies. We have observed domain formation and phase separation in the bilayer by increasing the composition of the gel phase DSPC. In force spectroscopy studies, we determine the mechanical strength of the bilayer for unique mixtures of DLPC: DSPC by measuring the breakthrough force. These results also show the effect of asymmetry in the lipid bilayer. Besides AFM studies, we have implemented a coarse-grained (CG) molecular dynamics (MD) simulation using the gromacs package at room temperature and 1 bar pressure. The results from the simulation study have been compared with AFM study. It was found that the simulation studies corroborated the findings from AFM such as an increase in the bilayer thickness, change in the phase state, asymmetric and symmetric domain formation in the lipid bilayer.

Subdomain cryo-EM structure of nodaviral replication protein A crown complex provides mechanistic insights into RNA genome replication

For positive-strand RNA [(+)RNA] viruses, the major target for antiviral therapies is genomic RNA replication, which occurs at poorly understood membrane-bound viral RNA replication complexes. Recent cryoelectron microscopy (cryo-EM) of nodavirus RNA replication complexes revealed that the viral double-stranded RNA replication template is coiled inside a 30- to 90-nm invagination of the outer mitochondrial membrane, whose necked aperture to the cytoplasm is gated by a 12-fold symmetric, 35-nm diameter "crown" complex that contains multifunctional viral RNA replication protein A. Here we report optimizing cryo-EM tomography and image processing to improve crown resolution from 33 to 8.5 Å. This resolves the crown into 12 distinct vertical segments, each with 3 major subdomains: A membrane-connected basal lobe and an apical lobe that together comprise the ∼19-nm-diameter central turret, and a leg emerging from the basal lobe that connects to the membrane at ∼35-nm diameter. Despite widely varying replication vesicle diameters, the resulting two rings of membrane interaction sites constrain the vesicle neck to a highly uniform shape. Labeling protein A with a His-tag that binds 5-nm Ni-nanogold allowed cryo-EM tomography mapping of the C terminus of protein A to the apical lobe, which correlates well with the predicted structure of the C-proximal polymerase domain of protein A. These and other results indicate that the crown contains 12 copies of protein A arranged basally to apically in an N-to-C orientation. Moreover, the apical polymerase localization has significant mechanistic implications for template RNA recruitment and (-) and (+)RNA synthesis.