Physical properties of model biological lipid bilayers: insights from all-atom molecular dynamics simulations

Insights into Transfer of Supramolecular Doxorubicin/Congo Red Aggregates through Phospholipid Membranes

Doxorubicin (DOX) is a commonly used chemotherapeutic drug, from the anthracycline class, which is genotoxic to neoplastic cells via a DNA intercalation mechanism. It is effective and universal; however, it also causes numerous side effects. The most serious of them are cardiotoxicity and a decrease in the number of myeloid cells. For this reason, targeted DOX delivery systems are desirable, since they would allow lowering the drug dose and therefore limiting systemic side effects. Recently, synthetic dyes, in particular Congo red (CR), have been proposed as possible DOX carriers. CR is a planar molecule, built of a central biphenyl moiety and two substituted naphthalene rings, connected with diazo bonds. In water, it forms elongated ribbon-shaped supramolecular structures, which are able to selectively interact with immune complexes. In our previous studies, we have shown that CR aggregates can intercalate DOX molecules. In this way, they preclude DOX precipitation in water solutions and increase its uptake by MCF7 breast cancer cells. In the present work, we further explore the interactions between DOX, CR, and their aggregates (CR/DOX) with phospholipid membranes. In addition to neutral molecules, the protonated doxorubicin form, DXP, is also studied. Molecular dynamics simulations are employed to study the transfer of CR, DOX, DXP, and their aggregates through POPC bilayers. Interactions of CR, DOX, and CR/DOX with model monolayers are studied with Langmuir trough measurements. This study shows that CR may support the transfer of doxorubicin molecules into the bilayer. Both electrostatic and van der Waals interactions with lipids are important in this respect. The former promote the initial stages of the insertion process, the latter keep guest molecules inside the bilayer.

Multifaceted Activity of Fabimycin: Insights from Molecular Dynamics Studies on Bacterial Membrane Models

Membranes─cells' essential scaffolds─are valid molecular targets for substances with an antimicrobial effect. While certain substances, such as octenidine, have been developed to target membranes for antimicrobial purposes, the recently reported molecule, fabimycin (F2B)─a novel agent targeting drug-resistant Gram-negative bacteria─has not received adequate attention regarding its activity on membranes in the literature. The following study aims to investigate the effects of F2B on different bacterial membrane models, including simple planar bilayers and more complex bilayer systems that mimic the Escherichia coli shell equipped with double inner and outer bilayers. Our results show that F2B exhibited more pronounced interactions with bacterial membrane systems compared to the control PC system. Furthermore, we observed significant changes in local membrane property homeostasis in both the inner and outer membrane models, specifically in the case of lateral diffusion, membrane thickness, and membrane resilience (compressibility, tilt). Finally, our results showed that the effect of F2B differed in a complex system and a single membrane system. Our study provides new insights into the multifaceted activity of F2B, demonstrating its potential to disrupt bacterial membrane homeostasis, indicating that its activity extends the currently known mechanism of FabI enzyme inhibition. This disruption, coupled with the ability of F2B to penetrate the outer membrane layers, sheds new light on the behavior of this antimicrobial molecule. This highlights the importance of the interaction with the membrane, crucial in combating bacterial infections, particularly those caused by drug-resistant strains.

Elucidating Daptomycin's Antibacterial Efficacy: Insights into the Tripartite Complex with Lipid II and Phospholipids in Bacterial Septum Membrane

This study elucidated the mechanism of formation of a tripartite complex containing daptomycin (Dap), lipid II, and phospholipid phosphatidylglycerol in the bacterial septum membrane, which was previously reported as the cause of the antibacterial action of Dap against gram-positive bacteria via molecular dynamics and enhanced sampling methods. Others have suggested that this transient complex ushers in the inhibition of cell wall synthesis by obstructing the downstream polymerization and cross-linking processes involving lipid II, which is absent in the presence of cardiolipin lipid in the membrane. In this work, we observed that the complex was stabilized by Ca2+-mediated electrostatic interactions between Dap and lipid head groups, hydrophobic interaction, hydrogen bonds, and salt bridges between the lipopeptide and lipids and was associated with Dap concentration-dependent membrane depolarization, thinning of the bilayer, and increased lipid tail disorder. Residues Orn6 and Kyn13, along with the DXDG motif, made simultaneous contact with constituent lipids, hence playing a crucial role in the formation of the complex. Incorporating cardiolipin into the membrane model led to its competitively displacing lipid II away from the Dap, reducing the lifetime of the complex and the nonexistence of lipid tail disorder and membrane depolarization. No evidence of water permeation inside the membrane hydrophobic interior was noted in all of the systems studied. Additionally, it was shown that using hydrophobic contacts between Dap and lipids as collective variables for enhanced sampling gave rise to a free energy barrier for the translocation of the lipopeptide. A better understanding of Dap's antibacterial mechanism, as studied through this work, will help develop lipopeptide-based antibiotics for rising Dap-resistant bacteria.

Functional in vitro diversity of an intrinsically disordered plant protein during freeze-thawing is encoded by its structural plasticity

Intrinsically disordered late embryogenesis abundant (LEA) proteins play a central role in the tolerance of plants and other organisms to dehydration brought upon, for example, by freezing temperatures, high salt concentration, drought or desiccation, and many LEA proteins have been found to stabilize dehydration-sensitive cellular structures. Their conformational ensembles are highly sensitive to the environment, allowing them to undergo conformational changes and adopt ordered secondary and quaternary structures and to participate in formation of membraneless organelles. In an interdisciplinary approach, we discovered how the functional diversity of the Arabidopsis thaliana LEA protein COR15A found in vitro is encoded in its structural repertoire, with the stabilization of membranes being achieved at the level of secondary structure and the stabilization of enzymes accomplished by the formation of oligomeric complexes. We provide molecular details on intra- and inter-monomeric helix-helix interactions, demonstrate how oligomerization is driven by an α-helical molecular recognition feature (α-MoRF) and provide a rationale that the formation of noncanonical, loosely packed, right-handed coiled-coils might be a recurring theme for homo- and hetero-oligomerization of LEA proteins.

Unlocking bacterial defense: Exploring the potent inhibition of NorA efflux pump by coumarin derivatives in Staphylococcus aureus

The occurrence of bacterial resistance has been increasing, compromising the treatment of various infections. The high virulence of Staphylococcus aureus allows for the maintenance of the infectious process, causing many deaths and hospitalizations. The MepA and NorA efflux pumps are transporter proteins responsible for expelling antimicrobial agents such as fluoroquinolones from the bacterial cell. Coumarins are phenolic compounds that have been studied for their diverse biological actions, including against bacteria. A pharmacokinetic in silico characterization of compounds C10, C11, C13, and C14 was carried out according to the principles of Lipinski's Rule of Five, in addition to searching for similarity in ChemBL and subsequent search for publications in CAS SciFinder. All compounds were evaluated for their in vitro antibacterial and modulatory activity against standard and multidrug-resistant Gram-positive and Gram-negative strains. The effect of coumarins C9, C10, C11, C13, and C14 as efflux pump inhibitors in Staphylococcus aureus strains was evaluated using the microdilution method (MepA or NorA) and fluorimetry (NorA). The behavior of coumarins regarding the efflux pump was determined from their interaction properties with the membrane and coumarin-protein using molecular docking and molecular dynamics simulations. Only the isolated coumarin compound C13 showed antibacterial activity against standard strains of Staphylococcus aureus and Escherichia coli. However, the other tested coumarins showed modulatory capacity for fluoroquinolone and aminoglycoside antibacterials. Compounds C10, C13, and C14 were effective in reducing the MIC of both antibiotics for both multidrug-resistant strains, while C11 potentiated the effect of norfloxacin and gentamicin for Gram-positive and Gram-negative bacteria and only norfloxacin for Gram-negative. Only coumarin C14 produced synergistic effects when associated with ciprofloxacin in MepA-carrying strains. All tested coumarins have the ability to inhibit the NorA efflux pump present in Staphylococcus aureus, both in reducing the MIC and inducing increased ethidium bromide fluorescence emission in fluorimetry. The findings of this study offer an atomistic perspective on the potential of coumarins as active inhibitors of the NorA pump, highlighting their specific mode of action mainly targeting protein inhibition. In molecular docking, it was observed that coumarins are capable of interacting with various amino acid residues of the NorA pump. The simulation showed that coumarin C10 can cross the bilayer; however, the other coumarins interacted with the membrane but were unable to cross it. Coumarins demonstrated their potentiating role in the effect of norfloxacin through a dual mechanism: efflux pump inhibition through direct interaction with the protein (C9, C10, C11, and C13) and increased interaction with the membrane (C10 and C13). In the context of pharmacokinetic prediction studies, the studied structures have a suitable chemical profile for possible oral use. We suggest that coumarin derivatives may be an interesting alternative in the future for the treatment of resistant bacterial infections, with the possibility of a synergistic effect with other antibacterials, although further studies are needed to characterize their therapeutic effects and toxicity.

An update on antibacterial AlkylGuanidino Ureas: Design of new derivatives, synergism with colistin and data analysis of the whole library

Antimicrobial resistance (AMR) represents one of the most challenging global Public Health issues, with an alarmingly increasing rate of attributable mortality. This scenario highlights the urgent need for innovative medicinal strategies showing activity on resistant isolates (especially, carbapenem-resistant Gram-negative bacteria, methicillin-resistant S. aureus, and vancomycin-resistant enterococci) yielding new approaches for the treatment of bacterial infections. We previously reported AlkylGuanidino Ureas (AGUs) with broad-spectrum antibacterial activity and a putative membrane-based mechanism of action. Herein, new tetra- and mono-guanidino derivatives were designed and synthesized to expand the structure-activity relationships (SARs) and, thereby, tested on the same panel of Gram-positive and Gram-negative bacteria. The membrane-active mechanism of selected compounds was then investigated through molecular dynamics (MD) on simulated bacterial membranes. In the end, the newly synthesized series, along with the whole library of compounds (more than 70) developed in the last decade, was tested in combination with subinhibitory concentrations of the last resort antibiotic colistin to assess putative synergistic or additive effects. Moreover, all the AGUs were subjected to cheminformatic and machine learning analyses to gain a deeper knowledge of the key features required for bioactivity.

Molecular modeling of the carbohydrate corona formation on a polyvinyl chloride nanoparticle and its impact on the adhesion to lipid bilayers

The pervasive presence of nanoplastics (NPs) in the environment has gained increasing attention due to their accumulation in living organisms. These emerging contaminants inevitably interact with extracellular polymeric substances along respiratory or gastrointestinal tracts, and diverse organic coating on the surface of NPs, known as bio- or eco-corona, is formed. Although its impact on altering the NP properties and potential cell internalization has been extensively examined, studies on its role in NP partitioning in the cell membrane are elusive yet. In this work, molecular dynamics is used to investigate the formation of chitosan (CT) corona centered on a polyvinyl chloride (PVC) nanoparticle and the uptake of the resulting complex onto lipid membranes. Coarse-grained models compatible with the newly developed Martini 3.0 force field are implemented for the two polymers employing the atomistic properties as targets in the parameterization. The reliability of the coarse-grained polymer models is demonstrated by reproducing the structural properties of the PVC melt and of solvated CT strands, as well as by determining the conformation adopted by the latter at the NP surface. Results show that the spontaneous binding of CT chains of high and intermediate protonation degrees led to the formation of soft and hard corona that modulates the interaction of PVC core with model membranes. The structural changes of the corona adsorbed at the lipid-water interface enable a subsequent transfer of the NP to the center of the saturated lipid membranes and a complete or partial transition to a snorkel conformation depending on the hydrophilic/hydrophobic balance in the CT-PVC complex. Overall, the computational investigation of the coarse-grained model system provides implications for understanding how the eco-corona development influences the uptake and implicit toxicology of NPs.

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.

Computational Assessment of the Biocompatibility of Two-Dimensional g-C3N3 Toward Lipid Membranes

One of the most recent additions to the family of two-dimensional (2D) materials, graphitic C3N3 (g-C3N3), has been considered a viable contender for biomedical applications, although its potential toxicity remains elusive. We perform all-atom molecular dynamics simulations to decipher the interactions between model lipid membranes and g-C3N3 as a first step toward exploring the cytotoxicity induced at the nanoscale. We show that g-C3N3 can easily insert into the cellular membranes following a multistage mechanism consisting of simultaneous desolvation of the 2D material along with enrichment of nanomaterial-lipid interactions. Free energy calculations indicate that g-C3N3 is more stable in a membrane-bound state compared to an aqueous solution; however, the insertion of the material does not disturb the structural integrity of lipid membranes. After being inserted into a membrane, g-C3N3 is unlikely to be released into the cellular environment and is incapable of extracting lipid molecules from the membrane. The nature of interaction between the 2D material and membranes is found to be independent of the nanomaterial size. Also, the performance of g-C3N3 toward biomolecular delivery is shown to be significantly improved compared to the state-of-the-art 2D materials graphene and hexagonal boron nitride (h-BN). It is revealed that, the affinity of g-C3N3 toward lipid membranes is weaker compared to the nanotoxic graphene and h-BN, while being marginally higher than h2D-C2N, which in turn, increases the biocompatibility of the material, thereby brightening its future as a noncytotoxic material for forthcoming biomedical applications.

Stability of Cytoplasmic Membrane of Escherichia coli Bacteria in Aqueous and Ethanolic Environment

The mechanism of action of any antibacterial agent or disinfectant depends largely on their interaction with the bacterial membrane. Herein, we use the SPICA (surface property fitting coarse graining) force-field and develop a coarse-grained (CG) model for the structure of the cytoplasmic membrane of Escherichia coli (E. coli) and its interaction with water and ethanol. We elucidate the impact of different concentrations of ethanol on the cytoplasmic membrane bilayers and vesicles of E. coli using the CG molecular dynamics (CG MD) simulations. Our modeling approach first focuses on the parametrization of the required force-field for POPG lipid and its interaction with water, ethanol, and POPE lipid. Subsequently, the structural stability of the E. coli bacterial membrane in the presence of high and low concentrations of ethanol is delineated. Both flat bilayers as well as vesicles of E. coli membrane were considered for the CG MD. Our results reveal that, at low ethanol concentrations (<30 mol %), the size of the E. coli vesicles increases with discernible deformations in their shapes. Because of ethanol-induced interdigitation, thinning of the E. coli vesicular membrane is also observed. However, at higher ethanol concentrations (>30 mol %), the integrity of the vesicles is lost because of deteriorating invasion of ethanol molecules into the vesicle bilayer and significant weakening of lipid-lipid interactions. At higher ethanol concentrations (40 and 70 mol %), both the multivesicle and single-vesicle bacterial membranes exhibit a similar rupturing pattern wherein the extraction of lipids from the membrane and formation of aggregates of the component lipids are observed. These aggregates consist of polar head groups of 3-5 POPE/POPG lipids with intertwined nonpolar tails.

Overlay databank unlocks data-driven analyses of biomolecules for all

Tools based on artificial intelligence (AI) are currently revolutionising many fields, yet their applications are often limited by the lack of suitable training data in programmatically accessible format. Here we propose an effective solution to make data scattered in various locations and formats accessible for data-driven and machine learning applications using the overlay databank format. To demonstrate the practical relevance of such approach, we present the NMRlipids Databank-a community-driven, open-for-all database featuring programmatic access to quality-evaluated atom-resolution molecular dynamics simulations of cellular membranes. Cellular membrane lipid composition is implicated in diseases and controls major biological functions, but membranes are difficult to study experimentally due to their intrinsic disorder and complex phase behaviour. While MD simulations have been useful in understanding membrane systems, they require significant computational resources and often suffer from inaccuracies in model parameters. Here, we demonstrate how programmable interface for flexible implementation of data-driven and machine learning applications, and rapid access to simulation data through a graphical user interface, unlock possibilities beyond current MD simulation and experimental studies to understand cellular membranes. The proposed overlay databank concept can be further applied to other biomolecules, as well as in other fields where similar barriers hinder the AI revolution.

In-vitro antimicrobial activity of AF-DP protein and in-silico approach of cell membrane disruption

Microbial resistance against common antibiotics has become one of the most serious threats to human health. The increasing statistics on this problem show the necessity of finding a way to deal with it. In recent years, antimicrobial peptides with unique properties and the capability of targeting a wide range of pathogens, have been considered as a potential for replacing common antibiotics. A small chitin-binding protein with anticandidal activity was isolated from Moringa oleifera seeds by Neto and colleagues in 2017, which very much resembled antimicrobial peptides. In this study, the antimicrobial protein 'AF-DP' was identified and characterized. AF-DP was heterologously expressed, purified, and characterized, and its 3D structure was predicted. Six molecular dynamic simulations were performed to investigate how the protein interacts with Gram-negative inner and outer, Gram-positive, fungal, cancerous, and normal mammalian membranes. Also, its antimicrobial and anticancer activity was assessed in vitro via minimum inhibition concentration (MIC) and MTT assays, respectively. This protein with 111 amino acids and a total net charge (of 10.5) has been predicted to be mainly composed of alpha helix and random coils. Its MIC affecting the growth of Escherichia coli, Staphylococcus aureus, and Candida albicans was 30 µg/ml, 100 µg/ml, and 100 µg/ml, respectively; AF-DP showed anticancer activity against MCF-7 breast cancer cell line. Scanning electron microscopic analysis confirmed the creation of pores and scratches on the surface of the bacterial membrane. The results of this research show that AF-DP can be a candidate for the production of new drugs as an AMP with antimicrobial activity.Communicated by Ramaswamy H. Sarma.

Multimillion Atom Simulations of Di-8-ANEPPS Chromophores Embedded in a Model Plasma Membrane: Toward the Investigation of Realistic Dyed Cell Membranes

A multistep computational approach has been employed to study a multimillion all-atom dyed plasma membrane, with no less than 42 different lipid species spanning the major head groups and a variety of fatty acids, as well as cholesterol, with the objective of investigating its structure and dynamics, as well as its impact on the embedded di-8-ANEPPS dyes. The latter are commonly used as bioimaging probes and serve as local microscopes. So, they provide information on membrane morphology via their second harmonic nonlinear optical (NLO) responses, which have the advantage of being specific to interface regions and sensitive to the chromophore environment. In previous studies, this chromophore has only been studied in simpler membrane models, far from the complexity of real lipid bilayers, while, owing to the ever-increasing computational resources, multimillion lipid bilayers have been studied, giving access to the effects of its heterogeneity. First, using molecular dynamics (MD) simulations, it is found that the combination of lipids produces a more ordered and denser membrane compared to its homogeneous model counterparts, while the local environment of the embedded dyes becomes enriched in phosphatidylcholine. Subsequently, the second harmonic first hyperpolarizability of the probes was calculated at the TDDFT level on selected frames of MD, highlighting the influence of the lipid environment. Due to the complexity of the system, machine learning (ML) tools have been employed to establish relationships between the membrane structural parameters, the orientation of the probes, and their NLO responses. These ML approaches have revealed influential features, including the presence of diacylglycerol lipids close to the dye. On the whole, this work provides a first step toward understanding the cooperation, synergy, and interactions that occur in such complex guest-host environments, which have emerged as new targets for drug design and membrane lipid therapy.

Insight into the effect of ibuprofen on the permeability of the membrane: a molecular dynamic simulation study

Studying interactions between drugs and cell membranes is of great interest to designing novel drugs, optimizing drug delivery, and discerning drug mechanism action. In this study, we investigated the physical properties of the bilayer membrane model of POPC upon interaction with ibuprofen (IBU) using molecular dynamics simulations. The area per lipid (APL) was calculated to describe the effect of ibuprofen on the packing properties of the lipid bilayer. The APL was 0.58 nm2 and 0.63 nm2 for the membrane in low and high IBU respectively, and 0.57 nm2 for the membrane without IBU. Our finding showed that the mean square deviation (MSD) increased with increased ibuprofen content. In addition, the order parameter for the hydrocarbon chain of lipids increased with increased ibuprofen content. There was an increment in the transfer free energy after the head group region while it was maximum in the hydrophobic core for hydrogen peroxide (H2O2) (∼6.2 kcal.mol-1) and H2O (∼3.4 kcal.mol-1) which then decreased to respective values of (∼4.6 kcal.mol-1), and (∼2.3 kcal.mol-1) at the center of the bilayer in the presence of IBU. It seems that in the presence of ibuprofen, the free energy profile of the permeability of water and H2O2 significantly decreased. These findings show that ibuprofen significantly influences the physical properties of the bilayer by decreasing the packing and intermolecular interaction in the hydrocarbon chain region and increasing the water permeability of the bilayer. These results may provide insights into the local cytotoxic side effects of ibuprofen and its underlying molecular mechanisms.Communicated by Ramaswamy H. Sarma.

Nucleoside Analog Reverse-Transcriptase Inhibitors in Membrane Environment: Molecular Dynamics Simulations

The behavior of four drugs from the family of nucleoside analog reverse-transcriptase inhibitors (zalcitabine, stavudine, didanosine, and apricitabine) in a membrane environment was traced using molecular dynamics simulations. The simulation models included bilayers and monolayers composed of POPC and POPG phospholipids. It was demonstrated that the drugs have a higher affinity towards POPG membranes than POPC membranes due to attractive long-range electrostatic interactions. The results obtained for monolayers were consistent with those obtained for bilayers. The drugs accumulated in the phospholipid polar headgroup region. Two adsorption modes were distinguished. They differed in the degree of penetration of the hydrophilic headgroup region. Hydrogen bonds between drug molecules and phospholipid heads were responsible for adsorption. It was shown that apricitabine penetrated the hydrophilic part of the POPC and POPG membranes more effectively than the other drugs. Van der Waals interactions between S atoms and lipids were responsible for this.

Extension of the iSoLF implicit-solvent coarse-grained model for multicomponent lipid bilayers

iSoLF is a coarse-grained (CG) model for lipid molecules with the implicit-solvent approximation used in molecular dynamics (MD) simulations of biological membranes. Using the original iSoLF (iSoLFv1), MD simulations of lipid bilayers consisting of either POPC or DPPC and these bilayers, including membrane proteins, can be performed. Here, we improve the original model, explicitly treating the electrostatic interactions between different lipid molecules and adding CG particle types. As a result, the available lipid types increase to 30. To parameterize the potential functions of the new model, we performed all-atom MD simulations of each lipid at three different temperatures using the CHARMM36 force field and the modified TIP3P model. Then, we parameterized both the bonded and non-bonded interactions to fit the area per lipid and the membrane thickness of each lipid bilayer by using the multistate Boltzmann Inversion method. The final model reproduces the area per lipid and the membrane thickness of each lipid bilayer at the three temperatures. We also examined the applicability of the new model, iSoLFv2, to simulate the phase behaviors of mixtures of DOPC and DPPC at different concentrations. The simulation results with iSoLFv2 are consistent with those using Dry Martini and Martini 3, although iSoLFv2 requires much fewer computations. iSoLFv2 has been implemented in the GENESIS MD software and is publicly available.

Graph-Based Analyses of Dynamic Water-Mediated Hydrogen-Bond Networks in Phosphatidylserine: Cholesterol Membranes

Phosphatidylserine lipids are anionic molecules present in eukaryotic plasma membranes, where they have essential physiological roles. The altered distribution of phosphatidylserine in cells such as apoptotic cancer cells, which, unlike healthy cells, expose phosphatidylserine, is of direct interest for the development of biomarkers. We present here applications of a recently implemented Depth-First-Search graph algorithm to dissect the dynamics of transient water-mediated lipid clusters at the interface of a model bilayer composed of 1-palmytoyl-2-oleoyl-sn-glycero-2-phosphatidylserine (POPS) and cholesterol. Relative to a reference POPS bilayer without cholesterol, in the POPS:cholesterol bilayer there is a somewhat less frequent sampling of relatively complex and extended water-mediated hydrogen-bond networks of POPS headgroups. The analysis protocol used here is more generally applicable to other lipid:cholesterol bilayers.

Uncovering the mechanisms of cyclic peptide self-assembly in membranes with the chirality-aware MA(R/S)TINI forcefield

Cyclic peptides (CPs) formed by alternation of D- and L-amino acids (D,L-CPs) can self-assemble into nanotubes (SCPNs) by parallel or/and antiparallel stacking. Different applications have been attributed to these nanotubes, including the disruption of lipid bilayers of specific compositions and the selective transport of ions throughout membranes. Molecular dynamics (MD) simulations have significantly contributed to understand the interaction between CPs, including the structural, dynamic and transport properties of their supramolecular aggregates. The high computational cost of atomic resolution forcefields makes them impractical for simulating the self-assembly of macromolecules, so coarse-grained (CG) models might represent a more feasible solution for this purpose. However, general CG models used for the simulation of biomolecules such as the MARTINI forcefield do not explicitly consider the non-covalent interactions leading to the formation of secondary structure patterns in proteins. This becomes particularly important in the case of CPs due to the D- and L-chirality alternation in their sequence, leading to opposite orientations of the backbone polar groups on both sides of the cyclic ring plane. In order to overcome this limitation, we have extended the MARTINI forcefield to introduce chirality in each residue of the CPs. The new parametrization, which we have called MA(R/S)TINI, reproduces the expected self-assembly patterns for several CP sequences in the presence of different membrane models, explicitly considering the chirality of the CPs and with no significant extra computational cost. Our simulations provide new mechanistic information of how these systems self-assemble in presence of different lipid scenarios, showing that the CP-CP and CP-membrane interactions are sensitive to the peptide sequence chirality. This opens the door to design new bioactive CPs based on CG-MD simulations. A web-based tool for the automatic parameterization of new CP sequences using MA(R/S)TINI, among other functionalities, is under construction (see http://cyclopep.com).

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.

Studying the structural organization of non-membranous protein hemoglobin in a lipid environment after reconstitution

In our current work we have developed a supported 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer with embedded hemoglobin, reconstituted via detergent-mediated method. Microscopic studies revealed that the hemoglobin molecules could be visualized without any labelling agents. The reconstituted proteins assemble themselves as supramolecular structures to adapt to lipid bilayer environment. The nonionic detergent, n-octyl-β-d-glucoside (NOG) used for insertion of hemoglobin played an important role in formation of these structures. When concentrations of lipid, protein and detergent were raised by four folds, we observed phase separation by protein molecules within bilayer via protein-protein assembly. This phase separation process exhibited extremely slow kinetics to form large stable domains with correlation times in the order of minutes. Confocal Z-scanning images showed that these supramolecular structures generated membrane deformities. UV-Vis, Fluorescence and Circular Dichroism (CD) measurement indicated minor structural change to expose the hydrophobic regions of the protein to adjust the hydrophobic stress of the lipid environment whilst Small Angle Neutron Scattering (SANS) results indicated that the hemoglobin molecules retained their overall tetrameric form in the system. In conclusion, we state that this investigation allowed us to closely inspect some rare but noteworthy phenomena like the formation of supramolecular structures, large domain formation and membrane deformation etc.

Investigating membrane-binding properties of lipoxygenases using surface plasmon resonance

Lipoxygenases (LOXs) catalyze the oxidation of polyunsaturated fatty acids and synthesize oxylipin products that drive important cellular signaling processes in plants and animals. While there has been indirect evidence presented for the interaction of mammalian LOXs with membranes, a quantitative study of the molecular details of LOX-membrane interactions is lacking. Here, we mimicked biological membranes using surface plasmon resonance (SPR) sensor chips derivatized with 2-D planar lipophilic anchors (2D LP) to capture liposomes of varying phospholipid compositions that self-assemble into lipid bilayers on the SPR chip. The sensor chip surfaces were then used to investigate the membrane-binding properties of model LOX enzymes. SPR binding assays displayed reproducible and stable liposome capture to the sensor chip surface that allowed for the detailed characterization of LOX-membrane interactions. Our studies demonstrate a calcium-dependence for the membrane binding activities of coral 8R-LOX and human 15-LOX-2. Furthermore, our data confirm the importance of key membrane insertion loop residues in each of these LOX enzymes for membrane binding activity. Experiments utilizing model plant and human LOXs reveal differences in membrane-binding specificities. Our study establishes and validates a robust SPR-based platform using 2D LP sensor chips that allows for the detailed study of LOX-membrane interactions under different experimental conditions, including altered membrane compositions. Collectively, this investigation improves our overall understanding of LOX-membrane interaction properties, and our SPR-based approach holds potential for future use in the development of LOX-based therapeutics.

Coarse-grained modeling of polystyrene-modified CNTs and their interactions with lipid bilayers

In the present work, we describe Martini3 coarse-grained models of polystyrene and carboxyl-terminated polystyrene functionalized carbon nanotubes (CNTs) and investigate their interactions with lipid bilayers with and without cholesterol (CHOL) using molecular dynamics simulations. By changing the polystyrene chain length and grafting density at the end ring of the CNTs at two different nanotube concentrations, we observe the translocation of nanoparticles as well as changes in the lipid bilayer properties. Our results show that all developed models passively diffuse into the membranes without causing any damage to the membrane integrity, although high concentrations of CNTs induce structural and elastic changes in lipid bilayers. In the presence of CHOL, increasing CNT concentration results in decreased rates of CHOL transmembrane motions. On the other hand, CNTs are prone to lipid and polystyrene blockage, which affects their equilibrated configurations, and tilting behavior within the membranes. Hence, we demonstrate that polystyrene-functionalized CNTs are promising drug-carrier agents. However, polystyrene chain length and grafting density are important factors to consider to enhance the efficiency of drug delivery.

Cellular uptake of rose bengal is mediated by OATP1B1/1B3 transporters

Due to its fluorescent properties and high yield of singlet oxygen, rose bengal (RB) is one of the most promising photosensitizers for cancer treatment. However, the negative charge of RB molecule may significantly hamper its intracellular delivery by passive diffusion through the cell membrane. Thus, specific membrane protein transporters may be needed. The organic anion transporting polypeptides (OATPs) are a well-characterized group of membrane protein transporters, responsible for cellular uptake of a number of drugs. To our knowledge, this is the first study that evaluates cellular transport of RB mediated by the OATP transporter family. First, electrified liquid-liquid interface, together with biophysical analysis and molecular dynamics simulations were used to characterize the interaction of RB with several models of a cellular membranes. These experiments proved that RB interacts only with the membrane's surface, without spontaneously crossing the lipid bilayer. Evaluation of intracellular uptake of RB by flow cytometry and confocal microscopy showed significant differences in uptake between liver and intestinal cell line models differing in expression of OATP transporters. The use of specific pharmacological inhibitors of OATPs, together with Western blotting and in silico analysis, indicated that OATPs are crucial for cellular uptake of RB.

Lipid oxidation controls peptide self-assembly near membranes through a surface attraction mechanism

The self-assembly of peptides into supramolecular structures has been linked to neurodegenerative diseases but has also been observed in functional roles. Peptides are physiologically exposed to crowded environments of biomacromolecules, and particularly cellular membrane lipids. Previous research has shown that membranes can both accelerate and inhibit peptide self-assembly. Here, we studied the impact of membrane models that mimic cellular oxidative stress and compared this to mammalian and bacterial membranes. Using molecular dynamics simulations and experiments, we propose a model that explains how changes in peptide-membrane binding, electrostatics, and peptide secondary structure stabilization determine the nature of peptide self-assembly. We explored the influence of zwitterionic (POPC), anionic (POPG) and oxidized (PazePC) phospholipids, as well as cholesterol, and mixtures thereof, on the self-assembly kinetics of the amyloid β (1-40) peptide (Aβ40), linked to Alzheimer's disease, and the amyloid-forming antimicrobial peptide uperin 3.5 (U3.5). We show that the presence of an oxidized lipid had similar effects on peptide self-assembly as the bacterial mimetic membrane. While Aβ40 fibril formation was accelerated, U3.5 aggregation was inhibited by the same lipids at the same peptide-to-lipid ratio. We attribute these findings and peptide-specific effects to differences in peptide-membrane adsorption with U3.5 being more strongly bound to the membrane surface and stabilized in an α-helical conformation compared to Aβ40. Different peptide-to-lipid ratios resulted in different effects. We found that electrostatic interactions are a primary driving force for peptide-membrane interaction, enabling us to propose a model for predicting how cellular changes might impact peptide self-assembly in vivo.

Thermal site energy fluctuations in photosystem I: new insights from MD/QM/MM calculations

Cyanobacterial photosystem I (PSI) is one of the most efficient photosynthetic machineries found in nature. Due to the large scale and complexity of the system, the energy transfer mechanism from the antenna complex to the reaction center is still not fully understood. A central element is the accurate evaluation of the individual chlorophyll excitation energies (site energies). Such an evaluation must include a detailed treatment of site specific environmental influences on structural and electrostatic properties, but also their evolution in the temporal domain, because of the dynamic nature of the energy transfer process. In this work, we calculate the site energies of all 96 chlorophylls in a membrane-embedded model of PSI. The employed hybrid QM/MM approach using the multireference DFT/MRCI method in the QM region allows to obtain accurate site energies under explicit consideration of the natural environment. We identify energy traps and barriers in the antenna complex and discuss their implications for energy transfer to the reaction center. Going beyond previous studies, our model also accounts for the molecular dynamics of the full trimeric PSI complex. Via statistical analysis we show that the thermal fluctuations of single chlorophylls prevent the formation of a single prominent energy funnel within the antenna complex. These findings are also supported by a dipole exciton model. We conclude that energy transfer pathways may form only transiently at physiological temperatures, as thermal fluctuations overcome energy barriers. The set of site energies provided in this work sets the stage for theoretical and experimental studies on the highly efficient energy transfer mechanisms in PSI.

Fluorescent Probes cis- and trans-Parinaric Acids in Fluid and Gel Lipid Bilayers: A Molecular Dynamics Study

Fluorescence probes are indispensable tools in biochemical and biophysical membrane studies. Most of them possess extrinsic fluorophores, which often constitute a source of uncertainty and potential perturbation to the host system. In this regard, the few available intrinsically fluorescent membrane probes acquire increased importance. Among them, cis- and trans-parinaric acids (c-PnA and t-PnA, respectively) stand out as probes of membrane order and dynamics. These two compounds are long-chained fatty acids, differing solely in the configurations of two double bonds of their conjugated tetraene fluorophore. In this work, we employed all-atom and coarse-grained molecular dynamics simulations to study the behavior of c-PnA and t-PnA in lipid bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), representative of the liquid disordered and solid ordered lipid phases, respectively. All-atom simulations indicate that the two probes show similar location and orientation in the simulated systems, with the carboxylate facing the water/lipid interface and the tail spanning the membrane leaflet. The two probes establish interactions with the solvent and lipids to a similar degree in POPC. However, the almost linear t-PnA molecules have tighter lipid packing around them, especially in DPPC, where they also interact more with positively charged lipid choline groups. Probably for these reasons, while both probes show similar partition (assessed from computed free energy profiles across bilayers) to POPC, t-PnA clearly partitions more extensively than c-PnA to the gel phase. t-PnA also displays more hindered fluorophore rotation, especially in DPPC. Our results agree very well with experimental fluorescence data from the literature and allow deeper understanding of the behavior of these two reporters of membrane organization.

NS1 from Two Zika Virus Strains Differently Interact with a Membrane: Insights to Understand Their Differential Virulence

Zika virus (ZIKV) from Uganda (UG) expresses a phenotype related to fetal loss, whereas the variant from Brazil (BR) induces microcephaly in neonates. The differential virulence has a direct relation to biomolecular mechanisms that make one strain more aggressive than the other. The nonstructural protein 1 (NS1) is a key viral toxin to comprehend these viral discrepancies because of its versatility in many processes of the virus life cycle. Here, we aim to examine through coarse-grained models and molecular dynamics simulations the protein-membrane interactions for both NS1_ZIKV-UG and NS1_ZIKV-BR dimers. A first evaluation allowed us to establish that the NS1 proteins, in the membrane presence, explore new conformational spaces when compared to systems simulated without a lipid bilayer. These events derive from both differential coupling patterns and discrepant binding affinities to the membrane. The N-terminal domain, intertwined loop, and greasy finger proposed previously as binding membrane regions were also computationally confirmed by us. The anchoring sites have aromatic and ionizable residues that manage the assembly of NS1 toward the membrane, especially for the Ugandan variant. Furthermore, in the presence of the membrane, the difference in the dynamic cross-correlation of residues between the two strains is particularly pronounced in the putative epitope regions. This suggests that the protein-membrane interaction induces changes in the distal part related to putative epitopes. Taken together, these results open up new strategies for the treatment of flaviviruses that would specifically target these dynamic differences.

Molecular dynamics simulation of phosphatidylcholine membrane in low ionic strengths of sodium chloride

The one-microsecond molecular dynamics simulations of a membrane-protein complex investigate the influence of the aqueous sodium chloride solutions on the structure and dynamics of a palmitoyl-oleoyl-phosphatidylcholine bilayer membrane. The simulations were performed on five different concentrations (40, 150, 200, 300, and 400 mM) in addition to a salt-free system by using the charmm36 force field for all atoms. Four biophysical parameters, (membrane thicknesses of annular and bulk lipids, and the area per lipid of both leaflets), were computed separately. Nevertheless, the area per lipid was expressed by using the Voronoi algorithm. All time-independent analyses were carried out for the last 400 ns trajectories. Different concentrations revealed dissimilar membrane dynamics before equilibration. The biophysical properties of the membrane (thickness, area-per-lipid, and order parameter) have non-significant changes with increasing ionic strength, however, the 150 mM system had exceptional behavior. Sodium cations were dynamically penetrating the membrane forming weak coordinate bonds with single or multiple lipids. Nevertheless, the binding constant was unaffected by the cation concentration. The electrostatic and Van der Waals energies of lipid-lipid interactions were influenced by the ionic strength. On the other hand, the Fast Fourier Transform was performed to figure out the dynamics at the membrane-protein interface. The nonbonding energies of membrane-protein interactions and order parameters explained the differences in the synchronization pattern. All results were consensus with experimental and theoretical works.Communicated by Ramaswamy H. Sarma.

Study of the Rv1417 and Rv2617c Membrane Proteins and Their Interactions with Nicotine Derivatives as Potential Inhibitors of Erp Virulence-Associated Factor in Mycobacterium tuberculosis: An In Silico Approach

The increasing emergence of Mycobacterium tuberculosis (Mtb) strains resistant to traditional anti-tuberculosis drugs has alarmed health services worldwide. The search for new therapeutic targets and effective drugs that counteract the virulence and multiplication of Mtb represents a challenge for the scientific community. Several studies have considered the erp gene a possible therapeutic target in the last two decades, since its disruption negatively impacts Mtb multiplication. This gene encodes the exported repetitive protein (Erp), which is located in the cell wall of Mtb. In vitro studies have shown that the Erp protein interacts with two putative membrane proteins, Rv1417 and Rv2617c, and the impairment of their interactions can decrease Mtb replication. In this study, we present five nicotine analogs that can inhibit the formation of heterodimers and trimers between these proteins. Through DFT calculations, molecular dynamics, docking, and other advanced in silico techniques, we have analyzed the molecular complexes, and show the effect these compounds have on protein interactions. The results show that four of these analogs can be possible candidates to counteract the pathogenicity of Mtb. This study aims to combine research on the Erp protein as a therapeutic target in the search for new drugs that serve to create new therapies against tuberculosis disease.

Influence of the Nature of the Lipid Building Blocks on the Second-Order Nonlinear Optical Responses of an Embedded Di-8-ANEPPS Probe

The extensive collection of lipids found in cell membranes is justified by the fact that each lipid contributes to their overall structure, dynamics, and properties and so to the biological processes taking place within them. It also showcases that, in order to deepen our understanding of membranes, we need to have a tool to differentiate lipid bilayers of varying composition. In this work, we investigate a suite of single-component saturated glycerophospholipids varying only in their headgroup structure by analyzing the second harmonic generation (SHG) nonlinear optical (NLO) response of a probe, di-8-ANEPPS, embedded into the membranes. The seven hydrophilic heads chosen (phosphatidylcholine (PC), phosphatidylethanolamine (PE), diaglycerol (GL), phosphatidylserine (PS), phosphatidylglycerol (PG), phosphatidylinositol (PI), and phosphatidyc acid (PA)) represent all the major headgroups that are part of mammalian plasma membranes and provide an assortment of neutral, zwiterrionic, and charged species. First, molecular dynamics simulations revealed that the lipidic arrangement is strongly sensitive to the nature of the hydrophilic head and less to the variety in the hydrophobic region. Membranes exhibiting drastically opposite structural properties can be pointed out: 1,2-dihexadecanoyl-rac-glycerol (DPGL) is the thickest and most ordered and aligned system, whereas 1,2-diacyl-sn-glycero-3-phospho-(1'-sn-glycerol) (DPPG) is thinnest and least ordered and aligned system. The structural analyses are then confronted with the molecular NLO responses, β, computed at the time-dependent density functional theory (TDDFT) level. As the orientation of the chromophore is impacted by the various degrees of order within the lipid bilayers, the diagonal component of the β tensor parallel to the bilayer normal, β_ZZZ, is as well. In the end, this computational approach provides insights into the link between lipid building blocks and the NLO responses of the embedded dye.

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

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

Bombyx mori Cecropin D could trigger cancer cell apoptosis by interacting with mitochondrial cardiolipin

Cecropin D is an antimicrobial peptide from Bombyx mori displaying anticancer and pro-apoptotic activities and, together with Cecropin XJ and Cecropin A, one of the very few peptides targeting esophageal cancer. Cecropin D displays poor similarity to other cecropins but a remarkable similarity in the structure and activity spectrum with Cecropin A and Cecropin XJ, offering the possibility to highlight key motifs at the base of the biological activity. In this work we show by NMR and MD simulations that Cecropin D is partially structured in solution and stabilizes its two-helix folding upon interaction with biomimetic membranes. Simulations show that Cecropin D strongly interacts with the surface of cancer cell biomimetic bilayers where it recognises the phosphatidylserine headgroup often exposed in the outer leaflet of cancerous cells by means of specific salt bridges. Cecropin D is also able to penetrate deeply in bilayers containing cardiolipin, a phospholipid found in mitochondria, causing significant destabilization in the lipid packing which might account for its pro-apoptotic activity. In bacterial membranes, phosphatidylglycerol and phosphatidylethanolamine act synergically by electrostatically attracting cecropin D and providing access to the membrane core, respectively.

Cholesterol's inhibition effect on entering of chlorzoxazone into phosphatidylethanolamine bilayer: Relevance to cytochrome P450 drug metabolism at endoplasmic reticulum membranes

Many drugs are metabolized by cytochrome P450 (CYP) in the endoplasmic reticulum (ER) membrane. Recent studies have shown that CYP-substrate drugs reach the CYP active site after entering the lipid hydrophobic part of the ER membrane. To clarify the role of cholesterol (Chol) in the CYP-related drug metabolic process, we investigated the lipid bilayer entry of CYP-substrate drugs using a model membrane system as follows. The model membrane system comprised palmitoyl-oleoyl-phosphatidylethanolamine (POPE) and Chol. Phosphatidylethanolamine is the second major phospholipid component of ER membranes. Chlorzoxazone (CZX) was used as the CYP-substrate drug. Calorimetric measurements showed that the addition of CZX to POPE bilayers decreased the gel-liquid crystal phase transition temperature; X-ray diffraction indicated that CZX distributes into the liquid crystal phase bilayers but not practically the gel phase POPE bilayers. In the presence of Chol, dialysis and X-ray structural analyses showed that Chol inhibited CZX entry into the bilayer with an increase in Chol concentration. The Chol concentration in the ER membrane (5-10 mol%) is much lower than that in the plasma membrane (approximately 30 mol%). This fact may allow CYP-substrate drugs to enter the hydrophobic portion of the ER membrane more easily than other organelle membranes, yielding efficient drug metabolism.

Binding of His-tagged fluorophores to lipid bilayers of giant vesicles

His-tagged molecules can be attached to lipid bilayers via certain anchor lipids, a method that has been widely used for the biofunctionalization of membranes and vesicles. To observe the membrane-bound molecules, it is useful to consider His-tagged molecules that are fluorescent as well. Here, we study two such molecules, green fluorescence protein (GFP) and green-fluorescent fluorescein isothiocyanate (FITC), both of which are tagged with a chain of six histidines (6H) that bind to the anchor lipids within the bilayers. The His-tag 6H is much smaller than the GFP molecule but somewhat larger than the FITC dye. The lipid bilayers form giant unilamellar vesicles (GUVs), the behavior of which can be directly observed in the optical microscope. We apply and compare three well-established preparation methods for GUVs: electroformation on platinum wire, polyvinyl alcohol (PVA) hydrogel swelling, and electroformation on indium tin oxide (ITO) glass. Microfluidics is used to expose the GUVs to a constant fluorophore concentration in the exterior solution. The brightness of membrane-bound 6H-GFP exceeds the brightness of membrane-bound 6H-FITC, in contrast to the quantum yields of the two fluorophores in solution. In fact, 6H-FITC is observed to be strongly quenched by the anchor lipids which bind the fluorophores via Ni²⁺ ions. For both 6H-GFP and 6H-FITC, the membrane fluorescence is measured as a function of the fluorophores' molar concentration. The theoretical analysis of these data leads to the equilibrium dissociation constants K_d = 37.5 nM for 6H-GFP and K_d = 18.5 nM for 6H-FITC. We also observe a strong pH-dependence of the membrane fluorescence.

Biomembrane lipids: When physics and chemistry join to shape biological activity

Biomembranes constitute the first lines of defense of cells. While small molecules can often permeate cell walls in bacteria and plants, they are generally unable to penetrate the barrier constituted by the double layer of phospholipids, unless specific receptors or channels are present. Antimicrobial or cell-penetrating peptides are in fact highly specialized molecules able to bypass this barrier and even discriminate among different cell types. This capacity is made possible by the intrinsic properties of its phospholipids, their distribution between the internal and external leaflet, and their ability to mutually interact, modulating the membrane fluidity and the exposition of key headgroups. Although common phospholipids can be found in the membranes of most organisms, some are characteristic of specific cell types. Here, we review the properties of the most common lipids and describe how they interact with each other in biomembrane. We then discuss how their assembly in bilayers determines some key physical-chemical properties such as permeability, potential and phase status. Finally, we describe how the exposition of specific phospholipids determines the recognition of cell types by membrane-targeting molecules.

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

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

Human myelin proteolipid protein structure and lipid bilayer stacking

The myelin sheath is an essential, multilayered membrane structure that insulates axons, enabling the rapid transmission of nerve impulses. The tetraspan myelin proteolipid protein (PLP) is the most abundant protein of compact myelin in the central nervous system (CNS). The integral membrane protein PLP adheres myelin membranes together and enhances the compaction of myelin, having a fundamental role in myelin stability and axonal support. PLP is linked to severe CNS neuropathies, including inherited Pelizaeus-Merzbacher disease and spastic paraplegia type 2, as well as multiple sclerosis. Nevertheless, the structure, lipid interaction properties, and membrane organization mechanisms of PLP have remained unidentified. We expressed, purified, and structurally characterized human PLP and its shorter isoform DM20. Synchrotron radiation circular dichroism spectroscopy and small-angle X-ray and neutron scattering revealed a dimeric, α-helical conformation for both PLP and DM20 in detergent complexes, and pinpoint structural variations between the isoforms and their influence on protein function. In phosphatidylcholine membranes, reconstituted PLP and DM20 spontaneously induced formation of multilamellar myelin-like membrane assemblies. Cholesterol and sphingomyelin enhanced the membrane organization but were not crucial for membrane stacking. Electron cryomicroscopy, atomic force microscopy, and X-ray diffraction experiments for membrane-embedded PLP/DM20 illustrated effective membrane stacking and ordered organization of membrane assemblies with a repeat distance in line with CNS myelin. Our results shed light on the 3D structure of myelin PLP and DM20, their structure-function differences, as well as fundamental protein-lipid interplay in CNS compact myelin.

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.

Membrane Dipole Potential: An Emerging Approach to Explore Membrane Organization and Function

Biological membranes are complex organized molecular assemblies of lipids and proteins that provide cells and membrane-bound intracellular organelles their individual identities by morphological compartmentalization. Membrane dipole potential originates from the electrostatic potential difference within the membrane due to the nonrandom arrangement (orientation) of amphiphile and solvent (water) dipoles at the membrane interface. In this Feature Article, we will focus on the measurement of dipole potential using electrochromic fluorescent probes and highlight interesting applications. In addition, we will focus on ratiometric fluorescence microscopic imaging technique to measure dipole potential in cellular membranes, a technique that can be used to address novel problems in cell biology which are otherwise difficult to address using available approaches. We envision that membrane dipole potential could turn out to be a convenient tool in exploring the complex interplay between membrane lipids and proteins and could provide novel insights in membrane organization and function.

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.

Tailoring a Solvent-Assisted Method for Solid-Supported Hybrid Lipid-Polymer Membranes

Combining amphiphilic block copolymers and phospholipids opens new opportunities for the preparation of artificial membranes. The chemical versatility and mechanical robustness of polymers together with the fluidity and biocompatibility of lipids afford hybrid membranes with unique properties that are of great interest in the field of bioengineering. Owing to its straightforwardness, the solvent-assisted method (SA) is particularly attractive for obtaining solid-supported membranes. While the SA method was first developed for lipids and very recently extended to amphiphilic block copolymers, its potential to develop hybrid membranes has not yet been explored. Here, we tailor the SA method to prepare solid-supported polymer-lipid hybrid membranes by combining a small library of amphiphilic diblock copolymers poly(dimethyl siloxane)-poly(2-methyl-2-oxazoline) and poly(butylene oxide)-block-poly(glycidol) with phospholipids commonly found in cell membranes including 1,2-dihexadecanoyl-sn-glycero-3-phosphocholine, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine, sphingomyelin, and 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine-N-(glutaryl). The optimization of the conditions under which the SA method was applied allowed for the formation of hybrid polymer-lipid solid-supported membranes. The real-time formation and morphology of these hybrid membranes were evaluated using a combination of quartz crystal microbalance and atomic force microscopy. Depending on the type of polymer-lipid combination, significant differences in membrane coverage, formation of domains, and quality of membranes were obtained. The use of the SA method for a rapid and controlled formation of solid-supported hybrid membranes provides the basis for developing customized artificial hybrid membranes.

Algorithm to catalogue topologies of dynamic lipid hydrogen-bond networks

Lipid membrane interfaces host reactions essential for the functioning of cells. The hydrogen-bonding environment at the membrane interface is particularly important for binding of proteins, drug molecules, and ions. We present here the implementation and applications of a depth-first search algorithm that analyzes dynamic lipid interaction networks. Lipid hydrogen-bond networks sampled transiently during simulations of lipid bilayers are clustered according to main types of topologies that characterize three-dimensional arrangements of lipids connected to each other via short water bridges. We characterize the dynamics of hydrogen-bonded lipid clusters in simulations of model POPE and POPE:POPG membranes that are often used for bacterial membrane proteins, in a model of the Escherichia coli membrane with six different lipid types, and in POPS membranes. We find that all lipids sample dynamic hydrogen-bonded networks with linear, star, or circular arrangements of the lipid headgroups, and larger networks with combinations of these three types of topologies. Overall, linear lipid-water bridges tend to be short. Water-mediated lipid clusters in all membranes with PE lipids tend to be somewhat small, with about four lipids in all membranes studied here. POPS membranes allow circular arrangements of three POPS lipids to be sampled frequently, and complex arrangements of linear, star, and circular paths may also be sampled. These findings suggest a molecular picture of the membrane interface whereby lipid molecules transiently connect in clusters with somewhat small spatial extension.

Biophysical Characterization of Membrane Proteins Embedded in Nanodiscs Using Fluorescence Correlation Spectroscopy

Proteins embedded in biological membranes perform essential functions in all organisms, serving as receptors, transporters, channels, cell adhesion molecules, and other supporting cellular roles. These membrane proteins comprise ~30% of all human proteins and are the targets of ~60% of FDA-approved drugs, yet their extensive characterization using established biochemical and biophysical methods has continued to be elusive due to challenges associated with the purification of these insoluble proteins. In response, the development of nanodisc techniques, such as nanolipoprotein particles (NLPs) and styrene maleic acid polymers (SMALPs), allowed membrane proteins to be expressed and isolated in solution as part of lipid bilayer rafts with defined, consistent nanometer sizes and compositions, thus enabling solution-based measurements. Fluorescence correlation spectroscopy (FCS) is a relatively simple yet powerful optical microscopy-based technique that yields quantitative biophysical information, such as diffusion kinetics and concentrations, about individual or interacting species in solution. Here, we first summarize current nanodisc techniques and FCS fundamentals. We then provide a focused review of studies that employed FCS in combination with nanodisc technology to investigate a handful of membrane proteins, including bacteriorhodopsin, bacterial division protein ZipA, bacterial membrane insertases SecYEG and YidC, Yersinia pestis type III secretion protein YopB, yeast cell wall stress sensor Wsc1, epidermal growth factor receptor (EGFR), ABC transporters, and several G protein-coupled receptors (GPCRs).

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.

Alteration of lipid bilayer mechanics by volatile anesthetics: Insights from μs-long molecular dynamics simulations

Very few drugs in clinical practice feature the chemical diversity, narrow therapeutic window, unique route of administration, and reversible cognitive effects of volatile anesthetics. The correlation between their hydrophobicity and their potency and the increasing amount of evidence suggesting that anesthetics exert their action on transmembrane proteins, justifies the investigation of their effects on phospholipid bilayers at the molecular level, given the strong functional and structural link between transmembrane proteins and the surrounding lipid matrix. Molecular dynamics simulations of a model lipid bilayer in the presence of ethylene, desflurane, methoxyflurane, and the nonimmobilizer 1,2-dichlorohexafluorocyclobutane (also called F6 or 2N) at different concentrations highlight the structural consequences of VA partitioning in the lipid phase, with a decrease of lipid order and bilayer thickness, an increase in overall lipid lateral mobility and area-per-lipid, and a marked reduction in the mechanical stiffness of the membrane, that strongly correlates with the compounds' hydrophobicity.

•

Molecular simulations of lipid bilayer interaction with volatile anesthetics

•

Comparison of volatile anesthetics' and nonimmobilizers' effects on lipid bilayers

•

Ligand-dependent partitioning of the compounds in the lipid phase

•

Effects on bilayer thickness, stiffness, and lipid order upon ligand partitioning

Analysis of the Conformational Dynamics of the Soluble and Membrane-Bound Forms of Interleukin-1 Receptor Type-1: Insights into Linker Flexibility and Domain Rearrangements

Interleukin-1 receptor type 1 (IL-1R1) is a key player in inflammation and immune responses. This receptor regulates IL-1 activity in two forms: as a membrane-bound form and as a soluble ectodomain. The details and differences between the conformational dynamics of the membrane-bound and the soluble IL-1R1 ectodomains (ECDs) remain largely elusive. Here, we study and compare the structural dynamics of the soluble and membrane-bound IL-1R1-ECDs using molecular dynamics (MD) simulations, focusing on the flexible interdomain linker of the ECD, as well as the spatial rearrangements between the Ig-like domains of the ECD. To explore the membrane-bound conformations, a full-length IL-1R1 structural model was developed and subjected to classical equilibrium MD. Comparative analysis of multiple MD trajectories of the soluble and the membrane-bound IL-1R1-ECDs reveals that (i) as somewhat expected, the extent of the visited “open-to-closed” transitional states differs significantly between the soluble and membrane-bound forms; (ii) the soluble form presents open-closed transitions, sampling a wider rotational motion between the Ig-like domains of the ECD, visiting closed and “twisted” conformations in higher extent, whereas the membrane-bound form is characterized by more conformationally restricted states; (iii) interestingly, the backbone dihedral angles of residues Glu202, Glu203 and Asn204, located in the flexible linker, display the highest variations during the transition between discrete conformational states detected in IL-1R1, thus appearing to work as the “central wheel of a clock’s movement”. The simulations and analyses presented in this contribution offer a deeper insight into the structure and dynamics of IL-1R1, which may be explored in a drug discovery setting.

Small molecules enhance the potency of natural antimicrobial peptides

The skin-associated microbiome plays an important role in general well-being and in a variety of treatable skin conditions. In this regard, endogenous antimicrobial peptides have both a direct and indirect role in determining the composition of the microbiota. We demonstrate here that certain small molecular species can amplify the antimicrobial potency of naturally occurring antimicrobial peptides. In this study, we have used niacinamide, a form of vitamin B3 naturally found in foods and widely used in cosmetic skincare products, and two of its structural analogs, to investigate their cooperativity with the human antimicrobial peptide LL37 on the bacterium Staphylococcus aureus. We observed a clear synergistic effect of niacinamide and, to some extent, N-methylnicotinamide, whereas isonicotinamide showed no significant cooperativity with LL37. Adaptively biased molecular dynamics simulations using simplified model membrane substrates and single peptides revealed that these molecules partition into the headgroup region of an anionic bilayer used to mimic the bacterial membrane. The simulated effects on the physical properties of the simulated model membrane are well correlated with experimental activity observed in real biological assays despite the simplicity of the model. In contrast, these molecules have little effect on zwitterionic bilayers that mimic a mammalian membrane. We conclude that niacinamide and N-methylnicotinamide can therefore potentiate the activity of host peptides by modulating the physical properties of the bacterial membrane, and to a lesser extent through direct interactions with the peptide. The level of cooperativity is strongly dependent on the detailed chemistry of the additive, suggesting an opportunity to fine-tune the behavior of host peptides.

Ultrasensitive two-dimensional material-based MCF-7 cancer cell sensor driven by perturbation processes

Changes in lipid composition and structure during cell development can be markers for cell apoptosis or various diseases such as cancer. Although traditional fluorescence techniques utilising molecular probes have been studied, these methods are limited in studying these micro-changes as they require complex probe preparation and cannot be reused, making cell monitoring and detection challenging. Here, we developed a direct current (DC) resistance sensor based on two-dimensional (2D) molybdenum disulfide (MoS₂) nanosheets to enable cancer cell-specific detection dependent on micro-changes in the cancer cell membrane. Atomistic molecular dynamics (MD) simulations were used to study the interaction between 2D MoS₂ and cancer lipid bilayer systems, and revealed that previously unconsidered perturbations in the lipid bilayer can cause an increase in resistance. Under an applied DC sweep, we observed an increase in resistance when cancer cells were incubated with the nanosheets. Furthermore, a correlation was observed between the resistance and breast cancer epithelial cell (MCF-7) population, illustrating a cell population-dependent sensitivity of our method. Our method has a detection limit of ∼3 × 10³ cells, below a baseline of ∼1 × 10⁴ cells for the current state-of-the-art electrical-based biosensors using an adherent monolayer with homogenous cells. This combination of a unique 2D material and electrical resistance framework represents a promising approach for the early detection of cancerous cells and to reduce the risk of post-surgery cancer recurrence.