Atomic-scale structure and electrostatics of anionic palmitoyloleoylphosphatidylglycerol lipid bilayers with Na+ counterions

Design and preparation of pH-sensitive cytotoxic liposomal formulations containing antitumor colchicine analogues for target release

Herein, we describe the synthesis of pH-sensitive lipophilic colchicine prodrugs for liposomal bilayer inclusion, as well as preparation and characterization of presumably stealth PEGylated liposomes with above-mentioned prodrugs. These formulations liberate strongly cytotoxic colchicinoid derivatives selectively under slightly acidic tumor-associated conditions, ensuring tumor-targeted delivery of the compounds. The design of the prodrugs is addressed to pH-triggered release of active compounds in the slight acidic media, that corresponds to tumor microenvironment, while keeping sufficient stability of the whole formulation at physiological pH. Correlations between the structure of the conjugates, their hydrolytic stability, colloidal stability, ability of the prodrug retention in the lipid bilayer are described. Several formulations were found promising for further development and in vivo investigations.

Salting-out and Competitive Adsorption of Ethanol into Lipid Bilayer Membranes: Conflicting Effects of Salts on Ethanol-Membrane Interactions Studied by Molecular Dynamics Simulations

Small amphiphilic molecules, such as ethanol, disturb the structure of lipid bilayer membranes to increase the membrane permeability, which is important for applications such as drug delivery, disinfection, and fermentation. To investigate how and the extent to which coexisting salts affect membrane disturbance, we performed molecular dynamics (MD) simulations on lipid bilayer membranes composed of zwitterionic lipids in aqueous ethanol solutions containing 0-631 mM NaCl, KCl, and KI salts. The addition of salts at low concentrations induced cationic adsorption on the lipid membrane, which competes with ethanol adsorption, thereby reducing the hydrogen bonds between ethanol and lipid molecules. This competitive adsorption mitigated the membrane disturbance and decreased the permeation of ethanol molecules into the membrane. In contrast, higher salt concentrations enhanced the membrane disturbance and permeability, which was caused by the salting-out of ethanol from the aqueous phase to the lipid bilayer. These conflicting effects appearing at different concentrations were stronger with the chloride salts than with the iodide salt. Among the two chloride salts, NaCl and KCl, the latter showed a greater enhancement in ethanol permeation at high concentrations. This seeming anti-Hofmeister salting-out behavior resulted from greater Na+ adsorption, preventing the ethanol-lipid interactions.

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.

Sensing cholesterol-induced rigidity in model membranes with time-resolved fluorescence spectroscopy and microscopy

Here, we report the characterization of cholesterol levels on membrane fluidity with a twisted intramolecular charge transfer (TICT) membrane dye, namely DI-8-ANEPPS, using fluorescence lifetime techniques such as time-correlated single photon counting (TCSPC) and fluorescence lifetime imaging microscopy (FLIM). The characterized liposomes comprised a 3 : 1 ratio of POPC and POPG, respectively, 1% DI-8-ANEPPS, and increasing cholesterol levels from 0% to 50%. Fluorescence lifetime characterization revealed that increasing the cholesterol levels from 0% to 50% increases the fluorescence lifetime of DI-8-ANEPPS from 2.36 ns to 3.65 ns, a 55% increment. Such lengthening in the fluorescence lifetime is concomitant with reduced Stokes shifts and higher quantum yield, revealing that localized excitation (LE) dominates over TICT states with increased cholesterol levels. Fluorescence anisotropy measurements revealed a less isotropic environment in the membrane upon increasing cholesterol levels, suggesting a shift from liquid-disorder (Lα) to liquid-order (LO) upon adding cholesterol. Local electrostatic and dipole characterization experiments revealed that changes in the zeta-potential (ζ-potential) and transmembrane dipole potential (Ψd) induced by changes in cholesterol levels or the POPC : POPG ratio play a minimal role in the fluorescence lifetime outcome of DI-8-ANEPPS. Instead, these results indicate that the cholesterol's effect in restricting the degree of movement of DI-8-ANEPPS dominates its photophysics over the cholesterol effect on the local dipole strength. We envision that time-resolved spectroscopy and microscopy, coupled with TICT dyes, could be a convenient tool in exploring the complex interplay between membrane lipids, sterols, and proteins and provide novel insights into membrane fluidity, organization, and function.

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.

Implicit model to capture electrostatic features of membrane environment

Membrane protein structure prediction and design are challenging due to the complexity of capturing the interactions in the lipid layer, such as those arising from electrostatics. Accurately capturing electrostatic energies in the low-dielectric membrane often requires expensive Poisson-Boltzmann calculations that are not scalable for membrane protein structure prediction and design. In this work, we have developed a fast-to-compute implicit energy function that considers the realistic characteristics of different lipid bilayers, making design calculations tractable. This method captures the impact of the lipid head group using a mean-field-based approach and uses a depth-dependent dielectric constant to characterize the membrane environment. This energy function Franklin2023 (F23) is built upon Franklin2019 (F19), which is based on experimentally derived hydrophobicity scales in the membrane bilayer. We evaluated the performance of F23 on five different tests probing (1) protein orientation in the bilayer, (2) stability, and (3) sequence recovery. Relative to F19, F23 has improved the calculation of the tilt angle of membrane proteins for 90% of WALP peptides, 15% of TM-peptides, and 25% of the adsorbed peptides. The performances for stability and design tests were equivalent for F19 and F23. The speed and calibration of the implicit model will help F23 access biophysical phenomena at long time and length scales and accelerate the membrane protein design pipeline.

Interaction mechanism of egg derived peptides RVPSL and QIGLF with dipalmitoyl phosphatidylcholine membrane: microcalorimetric and molecular dynamics simulation studies

BACKGROUND

Egg-derived peptides are becoming increasingly popular due to their biological activity and non-toxic effects. The egg-derived peptides Arg-Val-Pro-Ser-Leu (RVPSL) and Gln-Ile-Gly-Leu-Phe (QIGLF) display strong angiotensin-converting enzyme inhibitory activity and they can be taken up by intestinal epithelial cells. The interaction of the egg-derived peptides RVPSL and QIGLF with the membrane remains unclear.

RESULTS

The position and structure of the peptides in the membrane were calculated. The maximum density values of RVPSL and QIGLF were 2.27 and 1.22 nm from the center of the 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC) membrane, respectively, indicating that peptides penetrated the membrane-water interface and were embedded in the membrane. The interaction of RVPSL and QIGLF with the DPPC membrane did not affect the average area per lipid or the lipid sequence parameters. The thermodynamic parameters ΔH, ΔG, and ΔS of the interaction between the peptide RVPSL with the DPPC membrane were 17.91 kJ mol-1 , -17.63 kJ mol-1 , 187.5 J mol-1 ·k-1 , respectively. The thermodynamic parameters ΔH, ΔG, and ΔS of the interaction between peptide QIGLF with DPPC membrane were 17.10 kJ mol-1 , -17.12 kJ mol-1 , 114.8 J mol-1 ·k-1 , respectively.

CONCLUSION

The results indicated that the binding of peptides RVPSL and QIGLF to DPPC is an endothermic, spontaneous, and entropy-driven reaction. The results of the study are relevant to the problem of the low bioavailability of bioactive peptides (BP). © 2023 Society of Chemical Industry.

Effect of phospholipid liposomes on prion fragment (106-128) amyloid formation

Misfolding and aggregation of cellular prion protein (PrPc) is a major molecular process involved in the pathogenesis of prion diseases. Here, we studied the aggregation properties of a prion fragment peptide PrP(106-128). The results show that the peptide aggregates in a concentration-dependent manner in an aqueous solution and that the aggregation is sensitive to pH and the preformed amyloid seeds. Furthermore, we show that the zwitterionic POPC liposomes moderately inhibit the aggregation of PrP(106-128), whereas POPC/cholesterol (8:2) vesicles facilitate peptide aggregation likely due to the increase of the lipid packing order and membrane rigidity in the presence of cholesterol. In addition, anionic lipid vesicles of POPG and POPG/cholesterol above a certain concentration accelerate the aggregation of the peptide remarkably. The strong electrostatic interactions between the N-terminal region of the peptide and POPG may constrain the conformational plasticity of the peptide, preventing insertion of the peptide into the inner side of the membrane and thus promoting fibrillation on the membrane surface. The results suggest that the charge properties of the membrane, the composition of the liposomes, and the rigidity of lipid packing are critical in determining peptide adsorption on the membrane surface and the efficiency of the membrane in catalyzing peptide oligomeric nucleation and amyloid formation. The peptide could be used as an improved model molecule to investigate the mechanistic role of the crucial regions of PrP in aggregation in a membrane-rich environment and to screen effective inhibitors to block key interactions between these regions and membranes for preventing PrP aggregation.

Synergistic Effect of Salt and Anionic Surfactants on Interfacial Tension Reduction: Insights from Molecular Dynamics Simulations

Surfactants are commonly utilized in chemical flooding processes alongside salt to effectively decrease interfacial tension (IFT). However, the underlying microscopic mechanism for the synergistic effect of salt and surfactants on oil displacement remains ambiguous. Herein, the structure and properties of the interface between water and n-dodecane are studied by means of molecular dynamics simulations, considering three types of anionic surfactants and two types of salts. As the salt concentration (ρsalt) increases, the IFT first decreases to a minimum value, followed by a subsequent increase to higher values. The salt ions reduce the IFT only at low ρsalt due to the salt screening effect and ion bridging effect, both of which contribute to a decrease in the nearest head-to-head distance of surfactants. By incorporating salt doping, the IFTs can be reduced by at most 5%. Notably, the IFTs of different surfactants are mainly determined by the hydrogen bond interactions between oxygen atoms in the headgroup and water molecules. The presence of a greater number of oxygen atoms corresponds to lower IFT values. Specifically, for alkyl ethoxylate sulfate, the ethoxy groups play a crucial role in reducing the IFTs. This study provides valuable insights into formulating anionic surfactants that are applicable to oil recovery processes in petroleum reservoirs using saline water.

Protein Corona of Anionic Fluid-Phase Liposomes Compromises Their Integrity Rather than Uptake by Cells

Despite the undisputable role of the protein corona in the biointeractions of liposome drug carriers, the field suffers from a lack of knowledge regarding the patterns of protein deposition on lipid surfaces with different compositions. Here, we investigated the protein coronas formed on liposomes of basic compositions containing combinations of egg phosphatidylcholine (PC), palmitoyloleoyl phosphatidylglycerol (POPG), and cholesterol. Liposome-protein complexes isolated by size-exclusion chromatography were delipidated and analyzed using label-free LC-MS/MS. The addition of the anionic lipid and cholesterol both affected the relative protein abundances (and not the total bound proteins) in the coronas. Highly anionic liposomes, namely those containing 40% POPG, carried corona enriched with cationic proteins (apolipoprotein C1, beta-2-glycoprotein 1, and cathelicidins) and were the least stable in the calcein release assay. Cholesterol improved the liposome stability in the plasma. However, the differences in the corona compositions had little effect on the liposome uptake by endothelial (EA.hy926) and phagocytic cells in the culture (U937) or ex vivo (blood-derived monocytes and neutrophils). The findings emphasize that the effect of protein corona on the performance of the liposomes as drug carriers occurs through compromising particle stability rather than interfering with cellular uptake.

Implicit model to capture electrostatic features of membrane environment

Membrane protein structure prediction and design are challenging due to the complexity of capturing the interactions in the lipid layer, such as those arising from electrostatics. Accurately capturing electrostatic energies in the low-dielectric membrane often requires expensive Poisson-Boltzmann calculations that are not scalable for membrane protein structure prediction and design. In this work, we have developed a fast-to-compute implicit energy function that considers the realistic characteristics of different lipid bilayers, making design calculations tractable. This method captures the impact of the lipid head group using a mean-field-based approach and uses a depth-dependent dielectric constant to characterize the membrane environment. This energy function Franklin2023 (F23) is built upon Franklin2019 (F19), which is based on experimentally derived hydrophobicity scales in the membrane bilayer. We evaluated the performance of F23 on five different tests probing (1) protein orientation in the bilayer, (2) stability, and (3) sequence recovery. Relative to F19, F23 has improved the calculation of the tilt angle of membrane proteins for 90% of WALP peptides, 15% of TM-peptides, and 25% of the adsorbed peptides. The performances for stability and design tests were equivalent for F19 and F23. The speed and calibration of the implicit model will help F23 access biophysical phenomena at long time and length scales and accelerate the membrane protein design pipeline.

Improving HyperCEST performance by favorable xenon exchange conditions in liposomal nanocarriers

MRI reporters that combine signal enhancement from saturation transfer with hyperpolarized ¹²⁹ Xe show nanomolar detection sensitivity for in vitro studies. However, they need further improvement for accelerated CEST build-up that sufficiently dominates the intrinsic loss of hyperpolarization under in vivo conditions. This study introduces liposomes with a HyperCEST-active lipopeptide to enhance the efficiency of a well known Xe host, CrA-ma, with medium Xe exchange kinetics in aqueous environment, by two orders of magnitude. The depolarization time for constant saturation power but increasing saturation time is used as a comparative measure to rank different nanocarrier formulations. A variable cage load illustrates that the available CEST sites should be well distributed throughout the nanocarriers to avoid inefficiency from back exchange. For a liposome loading with only 2 mol% CrA-lipopeptide, the higher exchange kinetics allowed us to work even with 17-fold lower saturation power than for CrA-ma itself to achieve significant image contrast with ¹²⁹ Xe. Overall, this study illustrates the wide parameter space that is now available when incorporating CrA-labelled lipopeptides into liposomal carriers.

Effect of cholesterol on the ion-membrane interaction: Zeta potential and dynamic light scattering study

Cholesterol in a bio-membrane plays a significant role in many cellular event and is known to regulate the functional activity of protein and ion channel. In this study we report a significant effect of cholesterol on the ion-membrane interaction. We prepare large unilamellar vesicles, composed of zwitterionic lipid DOPC and anionic lipid DOPG with different cholesterol concentration. Electrostatics of anionic membranes containing cholesterol in the presence of NaCl has systematically been explored using dynamic light scattering and zeta potential. Negative zeta potential of the membrane decreases its negative value with increasing ion concentration for all cholesterol concentrations. However, zeta potential itself decreases with increasing cholesterol content even in the absence of monovalent ions. Electrostatic behaviour of the membrane is determined from well-known Gouy Chapmann model. Negative surface charge density of the membrane decreases with increasing cholesterol content. Binding constant, estimated from the electrostatic double layer theory, is found to increase significantly in the presence of cholesterol. Comparison of electrostatic parameters of the membrane in the presence and absence of cholesterol suggests that cholesterol significantly alter the electrostatic behaviour of the membrane.

Relating Molecular Dynamics Simulations to Functional Activity for Gly-Rich Membranolytic Helical Kiadin Peptides

Kiadins are in silico designed peptides with a strong similarity to diPGLa-H, a tandem sequence of PGLa-H (KIAKVALKAL) and with single, double or quadruple glycine substitutions. They were found to show high variability in their activity and selectivity against Gram-negative and Gram-positive bacteria, as well as cytotoxicity against host cells, which are influenced by the number and placing of glycine residues along the sequence. The conformational flexibility introduced by these substitutions contributes differently peptide structuring and to their interactions with the model membranes, as observed by molecular dynamics simulations. We relate these results to experimentally determined data on the structure of kiadins and their interactions with liposomes having a phospholipid membrane composition similar to simulation membrane models, as well as to their antibacterial and cytotoxic activities, and also discuss the challenges in interpreting these multiscale experiments and understanding why the presence of glycine residues in the sequence affected the antibacterial potency and toxicity towards host cells in a different manner.

Enhanced photocatalytic and antibacterial activity of acridinium-grafted g-C3N4 with broad-spectrum light absorption for antimicrobial photocatalytic therapy

As a metal-free polymeric photocatalyst, graphitic carbon nitride (g-C₃N₄) has attracted great attention owing to its high stability and low toxicity. However, g-C₃N₄ suffers from low light harvesting ability which limits its applications in antimicrobial photocatalytic therapy (APCT). Herein, acridinium (ADN)-grafted g-C₃N₄ (ADN@g-C₃N₄) nanosheets are prepared via covalent grafting of ADN to g-C₃N₄. The obtained ADN@g-C₃N₄ exhibits a narrow optical band gap (2.12 eV) and a wide optical absorption spectrum (intensity a.u. > 0.30) ranging from ultraviolet to near-infrared region. Moreover, ADN@g-C₃N₄ would produce reactive oxygen species (ROS) under light irradiation to exert effective sterilization and biofilm elimination activities against both gram-negative and gram-positive bacteria. Molecular dynamics simulation reveals that the ADN@g-C₃N₄ may move toward, tile and insert the bacterial lipid bilayer membrane through strong van der Waals and electrostatic interaction, decreasing the order parameter of the lipid while increasing the conducive of ROS migration, inducing ADN@g-C₃N₄ with improved antimicrobial and antibiofilm performance. Moreover, ADN@g-C₃N₄ could efficiently eradicate oral biofilm on artificial teeth surfaces. This work may provide a broad-spectrum light-induced photocatalytic therapy for preventing and treating dental plaque diseases and artificial teeth-related infections, showing potential applications for intractable biofilm treatment applications. An acridinium-grafted g-C₃N₄ (ADN@g-C₃N₄) with a narrow band gap and broad-spectrum light absorption was synthesized. The narrow optical band gap and improved electrostatic interaction with bacterial lipid bilayer membrane of ADN@g-C₃N₄ strengthened the ROS generation and facilitated the diffusion of ROS to bacteria surface, leading to enhanced photocatalytic and antibacterial activity against bacteria and corresponding biofilm under light irradiation. STATEMENT OF SIGNIFICANCE: An acridinium-grafted g-C₃N₄ (ADN@g-C₃N₄) with a narrow band gap and broad-spectrum light absorption was developed as an antimicrobial photocatalytic therapy agent. The ADN@g-C₃N₄ exhibited enhanced photocatalytic and antibacterial activity against bacteria and corresponding biofilm under light irradiation, showing potential applications for intractable biofilm treatment.

Molecular Modeling is an Enabling Approach to Complement and Enhance Channelopathy Research

Hundreds of human membrane proteins form channels that transport necessary ions and compounds, including drugs and metabolites, yet details of their normal function or how function is altered by genetic variants to cause diseases are often unknown. Without this knowledge, researchers are less equipped to develop approaches to diagnose and treat channelopathies. High-resolution computational approaches such as molecular modeling enable researchers to investigate channelopathy protein function, facilitate detailed hypothesis generation, and produce data that is difficult to gather experimentally. Molecular modeling can be tailored to each physiologic context that a protein may act within, some of which may currently be difficult or impossible to assay experimentally. Because many genomic variants are observed in channelopathy proteins from high-throughput sequencing studies, methods with mechanistic value are needed to interpret their effects. The eminent field of structural bioinformatics integrates techniques from multiple disciplines including molecular modeling, computational chemistry, biophysics, and biochemistry, to develop mechanistic hypotheses and enhance the information available for understanding function. Molecular modeling and simulation access 3D and time-dependent information, not currently predictable from sequence. Thus, molecular modeling is valuable for increasing the resolution with which the natural function of protein channels can be investigated, and for interpreting how genomic variants alter them to produce physiologic changes that manifest as channelopathies. © 2022 American Physiological Society. Compr Physiol 12:3141-3166, 2022.

IN SITU SOLID-STATE NMR STUDY OF ANTIMICROBIAL PEPTIDE INTERACTIONs WITH ERYTHROCYTE MEMBRANES

Antimicrobial peptides are promising therapeutic agents to mitigate the global rise of antibiotic resistance. They generally act by perturbing the bacterial cell membrane and are thus less likely to induce resistance. Because they are membrane-active molecules, it is critical to verify and understand their potential action toward eukaryotic cells to help design effective and safe drugs. In this work, we studied the interaction of two antimicrobial peptides, aurein 1.2 and caerin 1.1, with red blood cell (RBC) membranes using in situ ³¹P and ²H solid-state NMR (SS-NMR). We established a protocol to integrate up to 25% of deuterated fatty acids in the membranes of ghosts, which are obtained when hemoglobin is removed from RBCs. Fatty acid incorporation and the integrity of the lipid bilayer were confirmed by SS-NMR and fluorescence confocal microscopy. Leakage assays were performed to assess the lytic power of the antimicrobial peptides. The in situ perturbation of the ghost membranes by aurein 1.2 and caerin 1.1 revealed by ³¹P and ²H SS-NMR is consistent with membrane perturbation through a carpet mechanism for aurein 1.2, whereas caerin 1.1 acts on RBCs via pore formation. These results are compatible with fluorescence microscopy images of the ghosts. The peptides interact with eukaryotic membranes following similar mechanisms that take place in bacteria, highlighting the importance of hydrophobicity when determining such interactions. Our work bridges model membranes and in vitro studies and provides an analytical toolbox to assess drug toxicity toward eukaryotic cells.

On the Coupling between Mechanical Properties and Electrostatics in Biological Membranes

Cell membrane structure is proposed as a lipid matrix with embedded proteins, and thus, their emerging mechanical and electrostatic properties are commanded by lipid behavior and their interconnection with the included and absorbed proteins, cytoskeleton, extracellular matrix and ionic media. Structures formed by lipids are soft, dynamic and viscoelastic, and their properties depend on the lipid composition and on the general conditions, such as temperature, pH, ionic strength and electrostatic potentials. The dielectric constant of the apolar region of the lipid bilayer contrasts with that of the polar region, which also differs from the aqueous milieu, and these changes happen in the nanometer scale. Besides, an important percentage of the lipids are anionic, and the rest are dipoles or higher multipoles, and the polar regions are highly hydrated, with these water molecules forming an active part of the membrane. Therefore, electric fields (both, internal and external) affects membrane thickness, density, tension and curvature, and conversely, mechanical deformations modify membrane electrostatics. As a consequence, interfacial electrostatics appears as a highly important parameter, affecting the membrane properties in general and mechanical features in particular. In this review we focus on the electromechanical behavior of lipid and cell membranes, the physicochemical origin and the biological implications, with emphasis in signal propagation in nerve cells.

Mechanistic Insight into How PEGylation Reduces the Efficacy of pH-Sensitive Liposomes from Molecular Dynamics Simulations

Liposome-based drug delivery systems composed of DOPE stabilized with cholesteryl hemisuccinate (CHMS) have been proposed as a drug delivery mechanism with pH-triggered release as the anionic form (CHSa) is protonated (CHS) at reduced pH; PEGylation is known to decrease this pH sensitivity. In this manuscript, we set out to use molecular dynamics (MD) simulations with a model with all-atom resolution to provide insight into why incorporation of poly(ethyleneglycol) (PEG) into DOPE–CHMS liposomes reduces their pH sensitivity; we also address two additional questions: (1) How CHSa stabilizes DOPE bilayers into a lamellar conformation at a physiological pH of 7.4? and (2) how the change from CHSa to CHS at acidic pH triggers the destabilization of DOPE bilayers? We found that (A) CHSa stabilizes the DOPE lipid membrane by increasing the hydrophilicity of the bilayer surface, (B) when CHSa changes to CHS by pH reduction, DOPE bilayers are destabilized due to a reduction in bilayer hydrophilicity and a reduction in the area per lipid, and (C) PEG stabilizes DOPE bilayers into the lamellar phase, thus reducing the pH sensitivity of the liposomes by increasing the area per lipid through penetration into the bilayer, which is our main focus.

Enhanced dynamics in the anomalous melting regime of DMPG lipid membranes

Like many soft materials, lipids undergo a melting transition associated with a significant increase in their dynamics. At temperatures below the main melting transition (T_m ), all molecular and collective dynamics are suppressed, while above T_m the alkyl tail motions, lipid diffusivity, and collective membrane undulations are at least an order of magnitude faster. Here we study the collective dynamics of dimyristoylphosphatidylglycerol (DMPG, di 14:0 PG) using neutron spin echo spectroscopy throughout its anomalous phase transition that occurs over a 12 °C-20° C wide temperature window. Our results reveal that the membranes are softer and more dynamic during the phase transition than at higher temperatures corresponding to the fluid phase and provide direct experimental evidence for the predicted increase in membrane fluctuations during lipid melting. These results provide new insights into the nanoscale lipid membrane dynamics during the melting transition and demonstrate how these dynamics are coupled to changes in the membrane structure.

Characterization of lipid composition and diffusivity in OLA generated vesicles

Giant Unilamellar Vesicles (GUVs) are a versatile tool in many branches of science, including biophysics and synthetic biology. Octanol-Assisted Liposome Assembly (OLA), a recently developed microfluidic technique enables the production and testing of GUVs within a single device under highly controlled experimental conditions. It is therefore gaining significant interest as a platform for use in drug discovery, the production of artificial cells and more generally for controlled studies of the properties of lipid membranes. In this work, we expand the capabilities of the OLA technique by forming GUVs of tunable binary lipid mixtures of DOPC, DOPG and DOPE. Using fluorescence recovery after photobleaching we investigated the lateral diffusion coefficients of lipids in OLA liposomes and found the expected values in the range of 1 μm²/s for the lipid systems tested. We studied the OLA derived GUVs under a range of conditions and compared the results with electroformed vesicles. Overall, we found the lateral diffusion coefficients of lipids in vesicles obtained with OLA to be quantitatively similar to those in vesicles obtained via traditional electroformation. Our results provide a quantitative biophysical validation of the quality of OLA derived GUVs, which will facilitate the wider use of this versatile platform.

Influence of the Headgroup on the Interaction of Poly(ethylene oxide)-Poly(propylene oxide) Block Copolymers with Lipid Bilayers

The lipid headgroup plays an important role in the association of polymers with lipid bilayer membranes. Herein, we report how a glycerol headgroup versus a choline headgroup affects the interaction of poly(ethylene oxide)-b-poly(propylene oxide) (PEO-PPO) block copolymers with lipid bilayer vesicles. Unilamellar vesicles composed of phosphatidylcholine and phosphatidylglycerol at various molar ratios were used as model membranes. The interactions between the block copolymers and lipid bilayers were quantified by pulsed-field gradient nuclear magnetic resonance (PFG-NMR) based on the distinctly different mobilities of free and bound polymers. All the investigated polymer species showed significantly higher binding with 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) sodium salt (POPG) liposomes than with 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) liposomes, indicating stronger association with the glycerol headgroup compared to the choline headgroup. This effect did not become significant until the composition of mixed POPC/POPG liposomes contained more than 20 mol % POPG. A plausible explanation for the enhanced polymer binding with POPG invokes the role of hydrogen bonding between the glycerol headgroup and the ether moieties of the polymers.

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

The physical properties of lipid bilayers are sensitive to the specific type and composition of the lipids that make up the many different types of cell membranes. Studying model bilayers of representative heterogeneous compositions can provide key insights into membrane functionality. In this work, we use atomistic molecular dynamics simulations to characterize key properties in a number of bilayer membranes of varying composition. We first examine basic properties, such as lipid area, volume, and bilayer thickness, of simple, homogeneous bilayers comprising several lipid types, which are prevalent in biological membranes. Such lipids are then used in simulations of heterogeneous systems representative of bacterial, mammalian, and cancer membranes. Our analysis is especially focused on depth-dependent, transmembrane profiles; in particular, we calculate lateral pressure and dipole potential profiles, two fundamental properties which play key roles in a large number of biological functions.

Computational Modeling of Exciton-Bath Hamiltonians for Light Harvesting 2 and Light Harvesting 3 Complexes of Purple Photosynthetic Bacteria at Room Temperature

Light harvesting 2 (LH2) complex is the primary component of the photosynthetic unit of purple bacteria that is responsible for harvesting and relaying excitons. The electronic absorption line shape of LH2 contains two major bands at 800 and 850 nm wavelength regions. Under low light conditions, some species of purple bacteria replace LH2 with light harvesting 3 (LH3), a variant form with almost the same structure as the former but with distinctively different spectral features. The major difference between the absorption line shapes of LH2 and LH3 is the shift of the 850 nm band of the former to a new 820 nm region. The microscopic origin of this difference has been the subject of some theoretical/computational investigations. However, the genuine molecular level source of such a difference is not clearly understood yet. This work reports a comprehensive computational study of LH2 and LH3 complexes so as to clarify different molecular level features of LH2 and LH3 complexes and to construct simple exciton-bath models with a common form. All-atomistic molecular dynamics simulations of both LH2 and LH3 complexes provide detailed molecular level structural differences of bacteriochlorophylls (BChls) in the two complexes, in particular, in their patterns of hydrogen bonding (HB) and torsional angles of the acetyl group. Time-dependent density functional theory calculation of the excitation energies of BChls for structures sampled from the MD simulations suggests that the observed differences in the HB and torsional angles cannot fully account for the experimentally observed spectral shift of LH3. Potential sources that can explain the actual spectral shift of LH3 are discussed, and their magnitudes are assessed through fitting of experimental line shapes. These results demonstrate the feasibility of developing simple exciton-bath models for both LH2 and LH3, which can be employed for large-scale exciton quantum dynamics in their aggregates.

Driving Force for Spontaneous Perforation of Bilayers Formed by Ionic Amphiphiles in Aqueous Salt

Spontaneous perforation of amphiphilic membranes is important in both living matter and technology because of an impact on functions of biological membranes and shape transitions of self-assembling structures. Nevertheless, no definite molecular mechanism has been established so far even for simple ionic surfactant systems. We show that spontaneous perforation of a bilayer formed by an ionic amphiphile is driven by electrostatics. Creation of large pores with a concave-convex geometry of the rim is promoted by lower electrostatic free energy than that for a flat nonperforated bilayer. The opposite effect comes from the elasticity of the hydrocarbon tails of the amphiphile that prefer flat geometry of a nonperforated bilayer. The balance between electrostatics and tail deformation controls the appearance of pores; this balance is modulated by added salt that screens the electrostatic interactions. We illustrate the proposed mechanism with the aid of classical aggregation model that has been extended by including an analytical description of the electrostatic contribution for the toroidal rim of a pore. Numerical solution of the linearized Poisson-Boltzmann equation confirms the role of electrostatic forces in formation of pores. For the ionic surfactants of C_nTAB family, we predict shape transitions including bilayer perforations and formation of branched micellar networks induced by changing salinity or temperature and demonstrate the effect of surfactant's molecular parameters on these transitions.

Insight into the interactions, residue snorkeling, and membrane disordering potency of a single antimicrobial peptide into different lipid bilayers

Pardaxin, with a bend-helix-bend-helix structure, is a membrane-active antimicrobial peptide that its membrane activity depends on the lipid bilayer composition. Herein, all-atom molecular dynamics (MD) simulations were performed to provide further molecular insight into the interactions, structural dynamics, orientation behavior, and cationic residues snorkeling of pardaxin in the DMPC, DPPC, POPC, POPG, POPG/POPE (3:1), and POPG/POPE (1:3) lipid bilayers. The results showed that the C-terminal helix of the peptide was maintained in all six types of the model-bilayers and pardaxin was tilted into the DMPC, DPPC, and POPG/POPE mixed bilayers more than the POPC and POPG bilayers. As well as, the structure of zwitterionic membranes was more affected by the peptide than the anionic bilayers. Taken together, the study demonstrated that the cationic residues of pardaxin snorkeled toward the interface of lipid bilayers and all phenylalanine residues of the peptide played important roles in the peptide-membrane interactions. We hope that this work will provide a better understanding of the interactions of antimicrobial peptides with the membranes.

Effect of Curcumin Addition on the Adsorption and Transport of a Cationic Dye across DPPG-POPG Liposomes Probed by Second Harmonic Spectroscopy

The effect of addition of curcumin on the adsorption and transport characteristics of a cationic dye, LDS⁺, across negatively charged bilayers composed of POPG and DPPG lipids were investigated by the interface selective second harmonic (SH) spectroscopic technique. Curcumin induced changes in the SH electric field signal of the LDS⁺ ions (E_2ω (LDS⁺)) were observed to depend critically on the bilayer acyl chain saturation/unsaturation ratio (S/U). Following earlier works, the increase in the E_2ω (LDS⁺) signal is attributed to the release of the Na⁺ counterions present in the head group region of the bilayer by curcumin and the decay of the E_2ω (LDS⁺) signal is attributed to the bilayer intercalated state of curcumin. While the changes observed in the E_2ω (LDS⁺) signal in the presence of POPG liposomes were consistent with our earlier study ( Varshney, G. K. et al. Langmuir , 2016 , 32 , 10415 - 10421 ), they were significantly different for DPPG liposomes, following curcumin addition. While the increase in the E_2ω (LDS⁺) signal in the presence of POPG liposomes, is marginal (∼10-20%) and instantaneous (<1 s) followed by a rapid decay (completed within ∼100 s), in the presence of DPPG liposomes it was observed to increase slowly and at saturation shows a substantial increase (100-200%), following curcumin addition. When liposomes consisting of a mixture of POPG and DPPG lipids are used, curcumin induced kinetic characteristics of the E_2ω (LDS⁺) signal showed a mixture of the individual kinetic characteristics observed for the unsaturated (POPG) and saturated (DPPG) liposomes. The observed kinetic trends of the E_2ω (LDS⁺) signal following curcumin addition are explained on the basis of the relative strength of the Na⁺-POPG and Na⁺-DPPG interaction. Higher ordering of the lipid acyl chain region in DPPG liposome makes the Na⁺-DPPG interaction much stronger than the Na⁺-POPG interaction. Further, it is proposed that, in POPG-DPPG liposomes, individual domains of POPG and DPPG lipids exist at low temperature as suggested by the observed temperature dependent kinetic characteristics of the E_2ω (LDS⁺) signal following curcumin addition. These domains are dependent on the S/U ratio and phase state of the bilayer. The gel phase was observed to be more conducive for individual domain formation. Results presented in this work not only support the notion that biological activity of curcumin is associated with its bilayer altering properties, but more interestingly it provides a qualitative insight about how bilayer phase separation can be achieved by modulating the hydrophobic interactions between the lipid acyl chains.

Characterizing the structure-function relationship reveals the mode of action of a novel antimicrobial peptide, P1, from jumper ant Myrmecia pilosula

Microbial infections of antibiotic-resistant strains cause serious diseases and have a significant impact on public health worldwide, so novel antimicrobial drugs are urgently needed. Insect venoms, a rich source of bioactive components containing antimicrobial peptides (AMPs), are attractive candidates for new therapeutic agents against microbes. Recently, a novel peptide, P1, identified from the venom of the Australian jumper ant Myrmecia pilosula, showed potent antimicrobial activities against both Gram-negative and Gram-positive bacteria, but its structure-function relationship is unknown. Here, we used biochemical and biophysical techniques coupled with computational simulations to explore the mode of action of P1 interaction with dodecylphosphocholine (DPC) micelles as a model membrane system. Our circular dichroism (CD) and NMR studies revealed an amphipathic α-helical structure for P1 upon interaction with DPC micelles. A paramagnetic relaxation enhancement approach revealed that P1 orients its α-helix segment (F6-G14) into DPC micelles. In addition, the α-helix segment could be essential for membrane permeabilization and antimicrobial activity. Moreover, the arginine residues R8, R11, and R15 significantly contribute to helix formation and membrane-binding affinity. The lysine residue K19 of the C-terminus functionally guides P1 to interact with DPC micelles in the early interaction stage. Our study provides insights into the mode of action of P1, which is valuable in modifying and developing potent AMPs as antibiotic drugs.

Effect of Sodium and Chloride Binding on a Lecithin Bilayer. A Molecular Dynamics Study

The effect of ion binding on the structural, mechanical, dynamic and electrostatic properties of a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer in a 0.5 M aqueous NaCl solution is investigated using classical atomistic molecular dynamics simulation with different force-field descriptions for ion-ion and ion-lipid interactions. Most importantly, the repulsive Lennard-Jones parameters for the latter were modified, such that approximately similar binding of cations and anions to the lipid membrane is achieved. This was done to qualitatively improve the apparent ion-lipid binding constants obtained from simulations with the original force field (Berger lipids and GROMOS87 ions in combination with the SPC water model) in comparison to experimental data. Furthermore, various parameters characterizing membrane structure, elasticity, order and dynamics are analyzed. It is found that ion binding as observed in simulations involving the modified in comparison to the original force-field description leads to: (i) a smaller salt-induced change in the area per lipid, which is in closer agreement with the experiment; (ii) a decrease in the area compressibility and bilayer thickness to values comparable to a bilayer in pure water; (iii) lipid deuterium order parameters and lipid diffusion coefficients on nanosecond timescales that are very similar to the values for a membrane in pure water. In general, salt effects on the structural properties of a POPC bilayer in an aqueous sodium-chloride solution appear to be reproduced reasonably well by the new force-field description. An analysis of membrane-membrane disjoining pressure suggests that the smaller salt-induced change in area per lipid induced by the new force-field description is not due to the alteration of membrane-associated net charge, but must rather be understood as a consequence of ion-specific effects on the arrangement of lipid molecules.

Model Lung Surfactant Films: Why Composition Matters

Lung surfactant replacement therapies, Survanta and Infasurf, and two lipid-only systems both containing saturated and unsaturated phospholipids and one containing additional palmitic acid were used to study the impact of buffered saline on the surface activity, morphology, rheology, and structure of Langmuir monolayer model membranes. Isotherms and Brewster angle microscopy show that buffered saline subphases induce a film expansion, except when the cationic protein, SP-B, is present in sufficient quantities to already screen electrostatic repulsion, thus limiting the effect of changing pH and adding counterions. Grazing incidence X-ray diffraction results indicate an expansion not only of the liquid expanded phase but also an expansion of the lattice of the condensed phase. The film expansion corresponded in all cases with a significant reduction in the viscosity and elasticity of the films. The viscoelastic parameters are dominated by liquid expanded phase properties and do not appear to be dependent on the structure of the condensed phase domains in a phase separated film. The results highlight that the choice of subphase and film composition is important for meaningful interpretations of measurements using model systems.

Effect of Bilayer Partitioning of Curcumin on the Adsorption and Transport of a Cationic Dye Across POPG Liposomes Probed by Second-Harmonic Spectroscopy

The effect of Curcumin partitioning into the bilayer during the adsorption and transport of a cationic dye, LDS, across a negatively charged POPG bilayer was investigated by the interface-selective second-harmonic (SH) spectroscopic technique. The intensity of SH electric field (E_2ω) arising due to LDS adsorbed on the outer bilayer of the POPG liposome was observed to increase instantaneously (<1 s) following the addition of Curcumin. The fractional increase in the SH electric field (E^f_2ω) and the bilayer transport rates (k_T) of LDS were studied with respect to the pH of the solution and also with the Curcumin content in the lipid bilayer. Results obtained indicate that compared with the anionic form of the drug, its neutral form is more conducive of increasing the E^f_2ω of LDS. With increasing Curcumin content in the lipid bilayer, two distinct regimes could be observed in terms of E^f_2ω and k_T values of LDS. For Curcumin:Lipid (C/L) ratio ≤0.02, the E^f_2ω of LDS increases rapidly, while k_T remains unchanged; and for C/L ratio ≥0.02, the E^f_2ω values remains more or less constant, while there is a significant (∼40 times) increase followed by a modest increase in the k_T values of LDS. The observed results support an earlier two-state binding model of Curcumin with the POPG bilayer. In addition, it is further proposed that at low C/L ratio Curcumin binds to the surface of the bilayer replacing the counterions (Na⁺) bound to the lipid head groups, which changes the bilayer surface charge density, thereby causing more LDS cations to adsorb on the bilayer surface. At high C/L ratio, Curcumin intercalates within the hydrophobic domain of the bilayer, altering its hydrophobicity and inducing enhanced transport of the LDS cation. Results presented in this work provide further insights into how Curcumin alters bilayer properties when it partitions from the aqueous to the bilayer phase.

Biological Membranes in Extreme Conditions: Anionic Tetraether Lipid Membranes and Their Interactions with Sodium and Potassium

Archaea such as Sulfolobus acidocaldarius tolerate extreme temperatures and high acidity and can grow in the presence of toxic metals and low concentrations of Na⁺ or K⁺. It is believed that their unique tetraether membranes protect them from harsh environments and allow their survival under such conditions. We used molecular dynamics simulations to study membranes comprising glycerol dialkylnonitol tetraether lipids, which are the main component of S. acidocaldarius membranes, in solutions containing different concentrations of NaCl and KCl or with Na⁺ or K⁺ counterions (trace cations, 0 M). Anionic binding sites on the membranes were almost 50% occupied in the presence of counterions. The free energy of cation-phosphate complexation and the residence times of ions near the membranes were found to be both ion- and concentration-dependent. Sodium ions had more favorable interactions with the membranes and a longer residence time, whereas higher cation concentrations led to shorter ion residence times. When only counterions were present in the solutions, large residence times suggested that the membrane may function as a cation-attracting reservoir. The results suggested that the ions can be easily transferred to the cytoplasm as needed, explaining the growth curves of S. acidocaldarius under different salinities and pH.

Effect of cations on the solubilization/deposition of triclosan in sediment-water-rhamnolipid system

Cations had great influence on the self-assembly of rhamnolipid, which in turn affected the fate of triclosan. The migration of triclosan from sediment to water benefited its biodegradation but it could be transformed into more toxic compounds. To regulate the fate of triclosan and reduce environmental risks extremely, the effect of four common cations in surface water (Na(+)/K(+)/Ca(2+)/Mg(2+)) on the solubilization/deposition of triclosan in sediment-water-rhamnolipid system was investigated. The interaction among cations, triclosan and rhamnolipid was explored based on self-assembly of rhamnolipid and water solubility of triclosan in rhamnolipid solutions. Results showed that cations had little influence on the fate of triclosan in the absence of rhamnolipid. Cations, especially Ca(2+)/Mg(2+), reduced the critical micelle concentration, micellar size and zeta potential of rhamnolipid solutions. The changes in self-assembly of rhamnolipid with different cations led to the difference of residual rhamnolipid concentration in water, which was nearly invariant with 0.01 M Na(+)/K(+) while decreased significantly with 0.01 M Ca(2+)/Mg(2+). Consequently, water solubility of triclosan in rhamnolipid solutions increased with the addition of Na(+)/K(+) whereas decreased with Ca(2+)/Mg(2+). In sediment-water- rhamnolipid system, triclosan was slightly solubilized from sediment to water with Na(+)/K(+) while deposited in sediment with Ca(2+)/Mg(2+). These findings provided an alternative application of rhamnolipid for the remediation of triclosan-polluted sediment.