The importance of membrane defects-lessons from simulations

A Coronin 1-Derived Peptide Inhibits Membrane Fusion by Modulating Membrane Organization and Dynamics

Membrane fusion is considered the first step in the entry of enveloped viruses into the host cell. Several targeted strategies have been implemented to block viral entry by limiting the fusion protein to form a six-helix bundle, which is a prerequisite for fusion. Nonetheless, the development of broad-spectrum fusion inhibitors is essential to combat emerging and re-emerging viral infections. TG-23, a coronin 1, a tryptophan-aspartate-rich phagosomal protein-derived peptide, demonstrated inhibition of fusion between small unilamellar vesicles (SUVs) by modulating the membrane's physical properties. However, its inhibitory efficacy reduces with an increasing concentration of membrane cholesterol. The present work aims to develop a fusion inhibitor whose efficacy would be unaltered in the presence of membrane cholesterol. A stretch of the tryptophan-aspartic acid-containing peptide with a similar secondary structure and hydrophobicity profile of TG-23 from coronin 1 was synthesized, and its ability to inhibit SUV-SUV fusion with varying concentrations of membrane cholesterol was evaluated. Our results demonstrate that the GG-21 peptide inhibits fusion irrespective of the cholesterol content of the membrane. We have further evaluated the peptide-induced change in the membrane organization and dynamics utilizing arrays of steady-state and time-resolved fluorescence measurements and correlated these results with their effect on fusion. Interestingly, GG-21 displays inhibitory efficacy in a wide variety of lipid compositions despite having a secondary structure and physical properties similar to those of TG-23. Overall, our results advocate that the secondary structure and physical properties of the peptide may not be sufficient to predict its inhibitory efficacy.

Impact of the Unstirred Water Layer on the Permeation of Small-Molecule Drugs

Over the last two decades, numerous molecular dynamics (MD) simulation-based investigations have attempted to predict the membrane permeability to small-molecule drugs as indicators of their bioavailability, a majority of which utilize the inhomogeneous solubility diffusion (ISD) model. However, MD-based membrane permeability is routinely 3-4 orders of magnitude larger than the values measured with the intestinal perfusion technique. There have been contentious discussions on the sources of the large discrepancies, and the two indisputable, potentially dominant ones are the fixed protonation state of the permeant and the neglect of the unstirred water layer (UWL). Employing six small-molecule drugs of different biopharmaceutical classification system classes, the current MD study relies on the ISD model but introduces the (de)protonation of the permeant by characterizing the permeation free energy of both neutral and charged states. In addition, the role of the UWL as a potential resistance against permeation is explored. The new MD protocol closely mimics the nature of small-molecule permeation and yields estimates that agree well with in vivo intestinal permeability.

pSPICA Force Field Extended for Proteins and Peptides

Many coarse-grained (CG) molecular dynamics (MD) studies have been performed to investigate biological processes involving proteins and lipids. CG force fields (FFs) in these MD studies often use implicit or nonpolar water models to reduce computational costs. CG-MD using water models cannot properly describe electrostatic screening effects owing to the hydration of ionic segments and thus cannot appropriately describe molecular events involving water channels and pores through lipid membranes. To overcome this issue, we developed a protein model in the pSPICA FF, in which a polar CG water model showing the proper dielectric response was adopted. The developed CG model greatly improved the transfer free energy profiles of charged side chain analogues across the lipid membrane. Application studies on melittin-induced membrane pores and mechanosensitive channels in lipid membranes demonstrated that CG-MDs using the pSPICA FF correctly reproduced the structure and stability of the pores and channels. Furthermore, the adsorption behavior of the highly charged nona-arginine peptides on lipid membranes changed with salt concentration, indicating the pSPICA FF is also useful for simulating protein adsorption on membrane surfaces.

Phospholipids are imported into mitochondria by VDAC, a dimeric beta barrel scramblase

Mitochondria are double-membrane-bounded organelles that depend critically on phospholipids supplied by the endoplasmic reticulum. These lipids must cross the outer membrane to support mitochondrial function, but how they do this is unclear. We identify the Voltage Dependent Anion Channel (VDAC), an abundant outer membrane protein, as a scramblase-type lipid transporter that catalyzes lipid entry. On reconstitution into membrane vesicles, dimers of human VDAC1 and VDAC2 catalyze rapid transbilayer translocation of phospholipids by a mechanism that is unrelated to their channel activity. Coarse-grained molecular dynamics simulations of VDAC1 reveal that lipid scrambling occurs at a specific dimer interface where polar residues induce large water defects and bilayer thinning. The rate of phospholipid import into yeast mitochondria is an order of magnitude lower in the absence of VDAC homologs, indicating that VDACs provide the main pathway for lipid entry. Thus, VDAC isoforms, members of a superfamily of beta barrel proteins, moonlight as a class of phospholipid scramblases - distinct from alpha-helical scramblase proteins - that act to import lipids into mitochondria.

Peptide Flexibility and the Hydrophobic Moment are Determinants to Evaluate the Clinical Potential of Magainins

Using a flexibility prediction algorithm and in silico structural modeling, we have calculated the intrinsic flexibility of several magainin derivatives. In the case of magainin-2 (Mag-2) and magainin H2 (MAG-H2) we have found that MAG-2 is more flexible than its hydrophobic analog, Mag-H2. This affects the degree of bending of both peptides, with a kink around two central residues (R10, R11), whereas, in Mag-H2, W10 stiffens the peptide. Moreover, this increases the hydrophobic moment of Mag-H2, which could explain its propensity to form pores in POPC model membranes, which exhibit near-to-zero spontaneous curvatures. Likewise, the protective effect described in DOPC membranes for this peptide regarding its facilitation in pore formation would be related to the propensity of this lipid to form membranes with negative spontaneous curvature. The flexibility of another magainin analog (MSI-78) is even greater than that of Mag-2. This facilitates the peptide to present a kind of hinge around the central F12 as well as a C-terminal end prone to be disordered. Such characteristics are key to understanding the broad-spectrum antimicrobial actions exhibited by this peptide. These data reinforce the hypothesis on the determinant role of spontaneous membrane curvature, intrinsic peptide flexibility, and specific hydrophobic moment in assessing the bioactivity of membrane-active antimicrobial peptides.

SARS-CoV-2 fusion peptide sculpting of a membrane with insertion of charged and polar groups

The fusion peptide of SARS-CoV-2 spike is essential for infection. How this charged and hydrophobic domain occupies and affects membranes needs clarification. Its depth in zwitterionic, bilayered micelles at pH 5 (resembling late endosomes) was measured by paramagnetic NMR relaxation enhancements used to bias molecular dynamics simulations. Asp830 inserted deeply, along with Lys825 or Lys835. Protonation of Asp830 appeared to enhance agreement of simulated and NMR-measured depths. While the fusion peptide occupied a leaflet of the DMPC bilayer, the opposite leaflet invaginated with influx of water and choline head groups in around Asp830 and bilayer-inserted polar side chains. NMR-detected hydrogen exchange found corroborating hydration of the backbone of Thr827-Phe833 inserted deeply in bicelles. Pinching of the membrane at the inserted charge and the intramembrane hydration of polar groups agree with theory. Formation of corridors of hydrated, inward-turned head groups was accompanied by flip-flop of head groups. Potential roles of the defects are discussed.

Molecular Dynamics Simulations of Drug-Conjugated Cell-Penetrating Peptides

Cell-penetrating peptides (CPPs) are small peptides capable of translocating through biological membranes carrying various attached cargo into cells and even into the nucleus. They may also participate in transcellular transport. Our in silico study intends to model several peptides and their conjugates. We have selected three CPPs with a linear backbone, including penetratin, a naturally occurring oligopeptide; two of its modified sequence analogues (6,14-Phe-penetratin and dodeca-penetratin); and three natural CPPs with a cyclic backbone: Kalata B1, the Sunflower trypsin inhibitor 1 (SFT1), and Momordica cochinchinensis trypsin inhibitor II (MCoTI-II). We have also built conjugates with the small-molecule drug compounds doxorubicin, zidovudine, and rasagiline for each peptide. Molecular dynamics (MD) simulations were carried out with explicit membrane models. The analysis of the trajectories showed that the interaction of penetratin with the membrane led to spectacular rearrangements in the secondary structure of the peptide, while cyclic peptides remained unchanged due to their high conformational stability. Membrane-peptide and membrane-conjugate interactions have been identified and compared. Taking into account well-known examples from the literature, our simulations demonstrated the utility of computational methods for CPP complexes, and they may contribute to a better understanding of the mechanism of penetration, which could serve as the basis for delivering conjugated drug molecules to their intracellular targets.

Interaction of Near-Infrared (NIR)-Light Responsive Probes with Lipid Membranes: A Combined Simulation and Experimental Study

Cancer is considered a major societal challenge for the next decade worldwide. Developing strategies for simultaneous diagnosis and treatment has been considered a promising tool for fighting cancer. For this, the development of nanomaterials incorporating prototypic near-infrared (NIR)-light responsive probes, such as heptamethine cyanines, has been showing very promising results. The heptamethine cyanine-incorporating nanomaterials can be used for a tumor's visualization and, upon interaction with NIR light, can also produce a photothermal/photodynamic effect with a high spatio-temporal resolution and minimal side effects, leading to an improved therapeutic outcome. In this work, we studied the interaction of 12 NIR-light responsive probes with lipid membrane models by molecular dynamics simulations. We performed a detailed characterization of the location, orientation, and local perturbation effects of these molecules on the lipid bilayer. Based on this information, the probes were divided into two groups, predicting a lower and higher perturbation of the lipid bilayer. From each group, one molecule was selected for testing in a membrane leakage assay. The experimental data validate the hypothesis that molecules with charged substituents, which function as two polar anchors for the aqueous phase while spanning the membrane thickness, are more likely to disturb the membrane by the formation of defects and pores, increasing the membrane leakage. The obtained results are expected to contribute to the selection of the most suitable molecules for the desired application or eventually guiding the design of probe modifications for achieving an optimal interaction with tumor cell membranes.

Simultaneous inhibition of endocytic recycling and lysosomal fusion sensitizes cells and tissues to oligonucleotide therapeutics

Inefficient endosomal escape remains the primary barrier to the broad application of oligonucleotide therapeutics. Liver uptake after systemic administration is sufficiently robust that a therapeutic effect can be achieved but targeting extrahepatic tissues remains challenging. Prior attempts to improve oligonucleotide activity using small molecules that increase the leakiness of endosomes have failed due to unacceptable toxicity. Here, we show that the well-tolerated and orally bioavailable synthetic sphingolipid analog, SH-BC-893, increases the activity of antisense oligonucleotides (ASOs) and small interfering RNAs (siRNAs) up to 200-fold in vitro without permeabilizing endosomes. SH-BC-893 treatment trapped endocytosed oligonucleotides within extra-lysosomal compartments thought to be more permeable due to frequent membrane fission and fusion events. Simultaneous disruption of ARF6-dependent endocytic recycling and PIKfyve-dependent lysosomal fusion was necessary and sufficient for SH-BC-893 to increase non-lysosomal oligonucleotide levels and enhance their activity. In mice, oral administration of SH-BC-893 increased ASO potency in the liver by 15-fold without toxicity. More importantly, SH-BC-893 enabled target RNA knockdown in the CNS and lungs of mice treated subcutaneously with cholesterol-functionalized duplexed oligonucleotides or unmodified ASOs, respectively. Together, these results establish the feasibility of using a small molecule that disrupts endolysosomal trafficking to improve the activity of oligonucleotides in extrahepatic tissues.

Surface-Functionalized Polystyrene Nanoparticles Alter the Transmembrane Potential via Ion-Selective Pores Maintaining Global Bilayer Integrity

Although nanoplastics have well-known toxic effects toward the environment and living organisms, their molecular toxicity mechanisms, including the nature of nanoparticle-cell membrane interactions, are still under investigation. Here, we employ dynamic light scattering, quartz crystal microbalance with dissipation monitoring, and electrophysiology to investigate the interaction between polystyrene nanoparticles (PS NPs) and phospholipid membranes. Our results show that PS NPs adsorb onto lipid bilayers creating soft inhomogeneous films that include disordered defects. PS NPs form an integral part of the generated channels so that the surface functionalization and charge of the NP determine the pore conductive properties. The large difference in size between the NP diameter and the lipid bilayer thickness (∼60 vs ∼5 nm) suggests a particular and complex lipid-NP assembly that is able to maintain overall membrane integrity. In view of this, we suggest that NP-induced toxicity in cells could operate in more subtle ways than membrane disintegration, such as inducing lipid reorganization and transmembrane ionic fluxes that disrupt the membrane potential.

PEO-b-PBD Diblock Copolymers Induce Packing Defects in Lipid/Hybrid Membranes and Improve Insertion Rates of Natively Folded Peptides

Hybrid membranes assembled from biological lipids and synthetic polymers are a promising scaffold for the reconstitution and utilization of membrane proteins. Recent observations indicate that inclusion of small fractions of polymer in lipid membranes can improve protein folding and function, but the exact structural and physical changes a given polymer sequence imparts on a membrane often remain unclear. Here, we use all-atom molecular dynamics simulations to study the structure of hybrid membranes assembled from DOPC phospholipids and PEO-b-PBD diblock copolymers. We verified our computational model using new and existing experimental data and obtained a detailed picture of the polymer conformations in the lipid membrane that we can relate to changes in membrane elastic properties. We find that inclusion of low polymer fractions induces transient packing defects into the membrane. These packing defects act as insertion sites for two model peptides, and in this way, small amounts of polymer content in lipid membranes can lead to large increases in peptide insertion rates. Additionally, we report the peptide conformational space in both pure lipid and hybrid membranes. Both membranes support similar alpha helical peptide structures, exemplifying the biocompatibility of hybrid membranes.

Dynamics of the Acinetobacter baumannii inner membrane under exogenous polyunsaturated fatty acid stress

Exogenous polyunsaturated fatty acids (PUFAs) are readily incorporated into the synthesis pathways of A. baumannii membrane phospholipids, where they contribute to reduced bacterial fitness and increased antimicrobial susceptibility. Here we examine the impact of PUFA membrane modification on membrane organisation and biophysical properties using coarse grained MARTINI simulations of chemically representative membrane models developed from mass-spectrometry datasets of an untreated, arachidonic acid (AA) treated and docosahexaenoic acid (DHA) treated A. baumannii membranes. Enzymatic integration of AA or DHA into phospholipids of the A. baumannii membrane resulted in modulation of membrane biophysical properties. Membrane thickness decreased slightly following PUFA treatment, concomitant with changes in the lateral area per lipid of each lipid headgroup class. PUFA treatment resulted in a decrease in membrane ordering and an increase in lipid lateral diffusion. Changes in lateral membrane organisation were observed in the PUFA treated membranes, with a concurrent increase in ordered cardiolipin domains and disordered PUFA-containing domains. Notably, separation between ordered and disordered domains was enhanced and was more pronounced for DHA relative to AA, providing a possible mechanism for greater antimicrobial action of DHA relative to AA observed experimentally. Furthermore, the membrane active antimicrobial, pentamidine, preferentially adsorbs to cardiolipin domains of the A. baumannii model membranes. This interaction, and membrane penetration of pentamidine, was enhanced following PUFA treatment. Cumulatively, this work explores the wide-ranging effects of PUFA incorporation on the A. baumannii membrane and provides a molecular basis for bacterial inner membrane disruption by PUFAs.

Is Lipid Specificity Key to the Potential Antiviral Activity of Mouthwash Reagent Chlorhexidine against SARS-CoV-2?

Chlorhexidine (CHX), a popular antibacterial drug, is widely used for oral health. Emerging pieces of evidence suggest that commercially available chlorhexidine mouthwash formulations are effective in suppressing the spread of SARS-CoV-2, possibly through destabilization of the viral lipid envelope. CHX is known for its membrane-active properties; however, the molecular mechanism revealing how it damages the viral lipid envelope is yet to be understood. Here we used extensive conventional and umbrella sampling simulations to quantify the effects of CHX on model membranes mimicking the composition of the SARS-CoV-2 outer lipid membrane as well as the host plasma membrane. Our results show that the lipid composition and physical properties of the membrane play an important role in binding and insertion, with CHX binding favorably to the viral membrane over the plasma membrane. Among the simulated lipids, CHX preferentially binds to anionic lipids, PS and PI, which are more concentrated in the viral membrane. The deeper and stable binding of CHX to the viral membrane results in more pronounced swelling of the membrane laterally with a thinning of the bilayer. The overall free energies of pore formation are strongly reduced for the viral membrane compared to the plasma membrane; however, CHX has a larger concentration-dependent effect on free energies of pore formation in the plasma membrane than the viral membrane. The results indicate that CHX is less toxic to the human plasma membrane at low concentrations. Our simulations reveal that CHX facilitates pore formation by the combination of thinning the membrane and accumulation at the water defect. This study provides insights into the mechanism underlying the anti-SARS-CoV-2 potency of CHX, supporting its potential for application as an effective and safe oral rinse agent for preventing viral transmission.

TMEM16 scramblases thin the membrane to enable lipid scrambling

TMEM16 scramblases dissipate the plasma membrane lipid asymmetry to activate multiple eukaryotic cellular pathways. Scrambling was proposed to occur with lipid headgroups moving between leaflets through a membrane-spanning hydrophilic groove. Direct information on lipid-groove interactions is lacking. We report the 2.3 Å resolution cryogenic electron microscopy structure of the nanodisc-reconstituted Ca2+-bound afTMEM16 scramblase showing how rearrangement of individual lipids at the open pathway results in pronounced membrane thinning. Only the groove's intracellular vestibule contacts lipids, and mutagenesis suggests scrambling does not require specific protein-lipid interactions with the extracellular vestibule. We find scrambling can occur outside a closed groove in thinner membranes and is inhibited in thicker membranes, despite an open pathway. Our results show afTMEM16 thins the membrane to enable scrambling and that an open hydrophilic pathway is not a structural requirement to allow rapid transbilayer movement of lipids. This mechanism could be extended to other scramblases lacking a hydrophilic groove.

Modulating intracellular pathways to improve non-viral delivery of RNA therapeutics

RNA therapeutics (e.g. siRNA, oligonucleotides, mRNA, etc.) show great potential for the treatment of a myriad of diseases. However, to reach their site of action in the cytosol or nucleus of target cells, multiple intra- and extracellular barriers have to be surmounted. Several non-viral delivery systems, such as nanoparticles and conjugates, have been successfully developed to meet this requirement. Unfortunately, despite these clear advances, state-of-the-art delivery agents still suffer from relatively low intracellular delivery efficiencies. Notably, our current understanding of the intracellular delivery process is largely oversimplified. Gaining mechanistic insight into how RNA formulations are processed by cells will fuel rational design of the next generation of delivery carriers. In addition, identifying which intracellular pathways contribute to productive RNA delivery could provide opportunities to boost the delivery performance of existing nanoformulations. In this review, we discuss both established as well as emerging techniques that can be used to assess the impact of different intracellular barriers on RNA transfection performance. Next, we highlight how several modulators, including small molecules but also genetic perturbation technologies, can boost RNA delivery by intervening at differing stages of the intracellular delivery process, such as cellular uptake, intracellular trafficking, endosomal escape, autophagy and exocytosis.

Cardiolipin prevents pore formation in phosphatidylglycerol bacterial membrane models

Several antimicrobial peptides, including magainin and the human cathelicidin LL-37, act by forming pores in bacterial membranes. Bacteria such as Staphylococcus aureus modify their membrane's cardiolipin composition to resist such types of perturbations that compromise their membrane stability. Here, we used molecular dynamic simulations to quantify the role of cardiolipin on the formation of pores in simple bacterial-like membrane models composed of phosphatidylglycerol and cardiolipin mixtures. Cardiolipin modified the structure and ordering of the lipid bilayer, making it less susceptible to mechanical changes. Accordingly, the free-energy barrier for the formation of a transmembrane pore and its kinetic instability augmented by increasing the cardiolipin concentration. This is attributed to the unfavorable positioning of cardiolipin near the formed pore, due to its small polar head and bulky hydrophobic body. Overall, our study demonstrates how cardiolipin prevents membrane-pore formation and this constitutes a plausible mechanism used by bacteria to act against stress perturbations and, thereby, gain resistance to antimicrobial agents.

The Transporter-Mediated Cellular Uptake and Efflux of Pharmaceutical Drugs and Biotechnology Products: How and Why Phospholipid Bilayer Transport Is Negligible in Real Biomembranes

Over the years, my colleagues and I have come to realise that the likelihood of pharmaceutical drugs being able to diffuse through whatever unhindered phospholipid bilayer may exist in intact biological membranes in vivo is vanishingly low. This is because (i) most real biomembranes are mostly protein, not lipid, (ii) unlike purely lipid bilayers that can form transient aqueous channels, the high concentrations of proteins serve to stop such activity, (iii) natural evolution long ago selected against transport methods that just let any undesirable products enter a cell, (iv) transporters have now been identified for all kinds of molecules (even water) that were once thought not to require them, (v) many experiments show a massive variation in the uptake of drugs between different cells, tissues, and organisms, that cannot be explained if lipid bilayer transport is significant or if efflux were the only differentiator, and (vi) many experiments that manipulate the expression level of individual transporters as an independent variable demonstrate their role in drug and nutrient uptake (including in cytotoxicity or adverse drug reactions). This makes such transporters valuable both as a means of targeting drugs (not least anti-infectives) to selected cells or tissues and also as drug targets. The same considerations apply to the exploitation of substrate uptake and product efflux transporters in biotechnology. We are also beginning to recognise that transporters are more promiscuous, and antiporter activity is much more widespread, than had been realised, and that such processes are adaptive (i.e., were selected by natural evolution). The purpose of the present review is to summarise the above, and to rehearse and update readers on recent developments. These developments lead us to retain and indeed to strengthen our contention that for transmembrane pharmaceutical drug transport "phospholipid bilayer transport is negligible".

Phospholipid Membrane Formation Templated by Coacervate Droplets

We report the formation of coacervate-supported phospholipid membranes by hydrating a dried lipid film in the presence of coacervate droplets. Coacervate-supported membranes were characterized by fluorescence imaging, polarization, fluorescence recovery after photobleaching of labeled lipids, lipid quenching experiments, and solute uptake experiments. Our findings are consistent with the presence of lipid membranes around the coacervates, with many droplets fully coated by what appear to be continuous lipid bilayers. In contrast to traditional giant lipid vesicles formed by gentle hydration in the absence of coacervates, the coacervate-templated membrane vesicles are more uniform in size, shape, and apparent lamellarity. Due to their fully coacervate model cytoplasm, these simple artificial cells are macromolecularly crowded and can be easily pre-loaded with high concentrations of proteins or nucleic acids. Within the same population, in addition to coacervate droplets having intact lipid membrane coatings, other coacervate droplets are coated with membranes having defects or pores that permit solute entry, and some are coated with multilayered membranes. Membranes surrounding protein-based coacervate droplets provided protection from a protease added to the external solution. The simplicity of producing artificial cells having a coacervate model cytoplasm surrounded by a model membrane is at the same time interesting as a potential mechanism for prebiotic protocell formation and appealing for biotechnology. We anticipate that such structures could serve as a new type of model system for understanding interactions between intracellular phases and cell or organelle membranes, which are implicated in a growing number of processes ranging from neurotransmission to signaling.

Recent Experimental Developments in Studying Passive Membrane Transport of Drug Molecules

The ability to measure the passive membrane permeation of drug-like molecules is of fundamental biological and pharmaceutical importance. Of significance, passive diffusion across the cellular membranes plays an effective role in the delivery of many pharmaceutical agents to intracellular targets. Hence, approaches for quantitative measurement of membrane permeability have been the topics of research for decades, resulting in sophisticated biomimetic systems coupled with advanced techniques. In this review, recent developments in experimental approaches along with theoretical models for quantitative and real-time analysis of membrane transport of drug-like molecules through mimetic and living cell membranes are discussed. The focus is on time-resolved fluorescence-based, surface plasmon resonance, and second-harmonic light scattering approaches. The current understanding of how properties of the membrane and permeant affect the permeation process is discussed.

Tuning the transdermal transport by application of external continuous electric field: a coarse-grained molecular dynamics study

The control of skin permeability to specific substances (e.g., medications, vitamins, and nutrients) through stratum corneum is a challenge. Iontophoresis is an option in spite of the lack of a detailed understanding of the underlying molecular mechanism. In the present work, the simulations concerning application of an external continuous electric field to stratum corneum, in a range of low intensity (0-24 mV nm^-1), were carried out using the coarse-grained molecular dynamics approach. Using a set of random seed replicas of the starting configuration, we observed that in the range of electric field intensity of 22-23 mV nm^-1, water-rich lipid vesicles were formed in 20% of cases. Pores appeared in the remaining 80%. We argue that lipids undergo fast re-orientations under electric field inducing mechanical instability, which originates the pores. We presented a simple electrostatic model to interpret the results where the mismatch between electrical permittivities of the membrane and external media and the gradient of the local electric field in the membrane surface, govern the time scales and electric fields for vesicle formation. Our results indicate that just 10% difference between electrical permittivities of the membrane and external media decreases 1/6 the minimal time required for vesicle formation. The minimal electric field required decreases 10 times. The control and tunning of formation of biologically compatible vesicles, capable of transporting substances under low-intensity electric fields, has a promising application in fields such as drug therapy and dermo-cosmetics allowing the use of hydrophilic substances in dermal applications.

Joint Reaction Coordinate for Computing the Free-Energy Landscape of Pore Nucleation and Pore Expansion in Lipid Membranes

Topological transitions of membranes, such as pore formation or membrane fusion, play key roles in biology, biotechnology, and in medical applications. Calculating the related free-energy landscapes has been complicated by the fact that such processes involve a sequence of transitions along highly distinct directions in conformational space, making it difficult to define good reaction coordinates (RCs) for the overall process. In this study, a new RC capable of driving both pore nucleation and pore expansion in lipid membranes is presented. The potential of mean force (PMF) along the RC computed with molecular dynamics simulations provides a comprehensive view on the free-energy landscape of pore formation, including a barrier for pore nucleation; the size, free energy, and metastability of the open pore; and the energetic cost for further pore expansion against the line tension of the pore rim. The RC is illustrated by quantifying the effects of (i) simulation system size and (ii) the addition of dimethyl sulfoxide on the free-energy landscape of pore formation. PMF calculations along the RC provide mechanistic and energetic understanding of pore formation, hence they will be useful to rationalize the effects of membrane-active peptides, electric fields, and membrane composition on transmembrane pores.

Molecular Dynamics Simulations Based on Polarizable Models Show that Ion Permeation Interconverts between Different Mechanisms as a Function of Membrane Thickness

Different mechanisms have been proposed to explain the permeation of charged compounds through lipid membranes. Overall, it is expected that an ion-induced defect permeation mechanism, where substantial membrane deformations accompany ion movement, should be dominant in thin membranes but that a solubility-diffusion mechanism, where ions partition into the membrane core with large associated dehydration energy costs, becomes dominant in thicker membranes. However, while this physical picture is intuitively reasonable, capturing the interconversion between these two permeation mechanisms in molecular dynamics (MD) simulations based on atomic models is challenging. In particular, simulations relying on nonpolarizable force fields are artificially unfavorable to the solubility-diffusion mechanism, as induced polarization of the nonpolar hydrocarbon is ignored, causing overestimated free energy costs for charged molecules to enter into this region of the membrane. In this study, all-atom MD simulations based on nonpolarizable and polarizable force fields are used to quantitatively characterize the permeation process for the arginine side chain analog methyl-guanidinium through bilayer membranes of mono-unsaturated phosphatidylcholine lipids with and without cholesterol, resulting in thicknesses spanning from ∼24 to ∼42 Å. With simulations based on a nonpolarizable force field, ion translocation can take place solely through an ion-induced defect mechanism, with free energy barriers increasing linearly from 14 to 40 kcal/mol, depending on the thickness. However, with simulations based on a polarizable force field, ion translocation is predominantly dominated by an ion-induced defect mechanism in thin membranes, which progressively converts to a solubility-diffusion mechanism as the membranes get thicker. The transition between the two mechanisms occurs at a thickness of ∼29 Å, with lipid tails of 22 or more carbon atoms. This situation appears to represent the upper limit for ion-induced defect permeation within the current polarizable models. Beyond this thickness, it becomes energetically preferable for the ion to dehydrate and partition into the membrane core-a phenomenon that cannot be captured using the nonpolarizable models. Induced electronic polarizability therefore leads not just to a shift in permeation energetics but to an interconversion between two strikingly different physical mechanisms. The result highlights the importance of induced polarizability in modeling lipid membranes.

An electroporation strategy to synthesize the membrane-coated nanoparticles for enhanced anti-inflammation therapy in bone infection

The cell membrane-coated nanoparticles (MNPs) showed great potential in treating infectious disease due to their superior biofunctions in improving biocompatibility of nanoparticles and neutralization of pathogen or toxins. However, bone infection is accompanied with severe inflammation and bone loss, which also requires anti-inflammatory and osteoconductive treatment. The conventional membrane coating method has to undergo ultrasonication and extrusion procedures, which reduces the functionality of cell membrane and limits the choice of nanoparticles. In this study, we proposed an electroporation-based membrane coating strategy to facilitate the synthesis of MNPs to tackle those problems. Methods: Magnetic composite nanoparticles with osteoconductive Ca3(PO4)2 and bactericidal TiO2 were assembled into macrophages through phagocytosis and then collected to expose in electric field for obtaining macrophage membrane-coating nanoparticles. By using molecular dynamics simulation and materials characterizations, the cell membrane coating efficiency was confirmed. The in vitro anti-bacterial and anti-inflammatory abilities were tested by bacteria culturing and immune cells activation. Then drug-resistant bacteria induced bone infection model was established to verify its in vivo therapeutic effects. Results: The coated membrane prepared through electroporation reserved the integrality of membrane structure and right-sidedness, with more functional proteins. Those led to the superior properties of recognition and adsorption with bacteria, toxins and inflammatory cytokines. Owing to the benefits of electroporation, the MNPs exhibited significant better antibacterial and anti-inflammatory abilities for enhancing the tissue repair process. Conclusion: This study provides a novel self-assembly cell membrane coating strategy by electroporation to construct multifunctional membrane-coating nanoparticles for bone infection treatment. This strategy not only improves the functions of coated membrane, but is also proved to be universal for varies nanoparticles or cells, indicating a great potential for future applications in the bioengineering field.

Breakdown of classical paradigms in relation to membrane structure and functions

Updates of the mosaic fluid membrane model implicitly sustain the paradigms that bilayers are closed systems conserving a state of fluidity and behaving as a dielectric slab. All of them are a consequence of disregarding water as part of the membrane structure and its essential role in the thermodynamics and kinetics of membrane response to bioeffectors. A correlation of the thermodynamic properties with the structural features of water makes possible to introduce the lipid membrane as a responsive structure due to the relaxation of water rearrangements in the kinetics of bioeffectors' interactions. This analysis concludes that the lipid membranes are open systems and, according to thermodynamic of irreversible formalism, bilayers and monolayers can be reasonable compared under controlled conditions. The inclusion of water in the complex structure makes feasible to reconsider the concept of dielectric slab and fluidity.

Cooperative Effects of an Antifungal Moiety and DMSO on Pore Formation over Lipid Membranes Revealed by Free Energy Calculations

Itraconazole is a triazole drug widely used in the treatment of fungal infections, and it is in clinical trials for treatment of several cancers. However, the drug suffers from poor solubility, while experiments have shown that itraconazole delivery in liposome nanocarriers improves both circulation half-life and tissue distribution. The drug release mechanism from the nanocarrier is still unknown, and it depends on several factors including membrane stability against defect formation. In this work, we used molecular dynamics simulations and potential of mean force (PMF) calculations to quantify the influence of itraconazole on pore formation over lipid membranes, and we compared the effect by itraconazole with a pore-stabilizing effect by the organic solvent dimethyl sulfoxide (DMSO). According to the PMFs, both itraconazole and DMSO greatly reduce the free energy of pore formation, by up to ∼20 kJ mol^-1. However, whereas large concentrations of itraconazole of 8 mol % (relative to lipid) were required, only small concentrations of a few mole % DMSO (relative to water) were sufficient to stabilize pores. In addition, itraconazole and DMSO facilitate pore formation by different mechanisms. Whereas itraconazole predominantly aids the formation of a partial defect with a locally thinned membrane, DMSO mainly stabilizes a transmembrane water needle by shielding it from the hydrophobic core. Notably, the two distinct mechanisms act cooperatively upon adding both itraconazole and DMSO to the membrane, as revealed by an additional reduction of the pore free energy. Overall, our simulations reveal molecular mechanisms and free energies of membrane pore formation by small molecules. We suggest that the stabilization of a locally thinned membrane as well as the shielding of a transmembrane water needle from the hydrophobic membrane core may be a general mechanism by which amphiphilic molecules facilitate pore formation over lipid membranes at sufficient concentrations.

Cationic Antiseptics Facilitate Pore Formation in Model Bacterial Membranes

Antiseptics are an essential line of defense against bacterial and viral infections in modern medical practice. Many of them are supposed to act on microbial membranes. However, the detailed mechanisms of their action are still elusive. Here, we utilized coarse-grained molecular dynamics simulations to investigate interactions of different types of cationic antiseptics (CAs) with a model bacterial membrane. The simulations revealed qualitatively distinct patterns of dynamic and structural alterations of membrane induced by different types of antiseptics although none of them caused disintegration or solubilization of the bilayer even at the highest explored concentration. At the same time, the adsorption of antiseptics rendered membranes more vulnerable to poration under exposure to the external electric field. We further discuss the possible relation of the enhanced pore formation induced by CAs to their cytotoxic action.

Temperature-Dependent Dynamics of Giant Vesicles Composed of Hydrolysable Lipids Having an Amide Linkage

Various amphiphiles including surfactants and lipids have been designed and synthesized to improve and create new functionalities. In particular, the emergence of cell-like behaviors of giant vesicles (GVs) composed of synthetic lipids has drawn much attention in the development of chemical models for cells. The aim of this study was to measure temperature-dependent morphological changes of GVs induced by fragmentation and subsequent growth using hydrolysable cationic lipids having an amide linkage. Results from differential scanning calorimetry, fluorescence spectroscopy using an environment-responsive probe, and confocal Raman microscopy showed that the dynamics observed were due to changes in the vesicle membrane, including variation in the lipid composition, induced by thermal stimulation.

Spontaneous and Stress-Induced Pore Formation in Membranes: Theory, Experiments and Simulations

The large plasticity, dynamics and adaptability of biological membranes allow different modes of intrinsic and inducible permeability. These phenomena are of physiological importance for a number of natural functions related to cell death and can also be manipulated artificially for practical purposes like gene transfer, drug delivery, prevention of infections or anticancer therapy. For these advances to develop in a controllable and specific way, we need a sufficient understanding of the membrane permeability phenomena. Since the formulation of early concepts of pore formation, there has been an enormous effort to describe membrane permeability by using theory, simulations and experiments. A major breakthrough has come recently through theoretical developments that allow building continuous trajectories of pore formation both in the absence and presence of stress conditions. The new model provides a coherent quantitative view of membrane permeabilization, useful to test the impact of known lipid properties, make predictions and postulate specific pore intermediates that can be studied by simulations. For example, this theory predicts unprecedented dependencies of the line tension on the pore radius and on applied lateral tension which explain previous puzzling results. In parallel, important concepts have also come from molecular dynamics simulations, of which the role of water for membrane permeabilization is of special interest. These advances open new challenges and perspectives for future progress in the study of membrane permeability, as experiments and simulations will need to test the theoretical predictions, while theory achieves new refinements that provide a physical ground for observations.

Heterogeneous Dielectric Implicit Membrane Model for the Calculation of MMPBSA Binding Free Energies

Membrane-bound protein receptors are a primary biological drug target, but the computational analysis of membrane proteins has been limited. In order to improve molecular mechanics Poisson-Boltzmann surface area (MMPBSA) binding free energy calculations for membrane protein-ligand systems, we have optimized a new heterogeneous dielectric implicit membrane model, with respect to free energy simulations in explicit membrane and explicit water, and implemented it into the Amber software suite. This new model supersedes our previous uniform, single dielectric implicit membrane model by allowing the dielectric constant to vary with depth within the membrane. We calculated MMPBSA binding free energies for the human purinergic platelet receptor (P2Y₁₂R) and two of the muscarinic acetylcholine receptors (M2R and M3R) bound to various antagonist ligands using both membrane models, and we found that the heterogeneous dielectric membrane model has a stronger correlation with experimental binding affinities compared to the older model under otherwise identical conditions. This improved membrane model increases the utility of MMPBSA calculations for the rational design and improvement of future drug candidates.

Transport of a graphene nanosheet sandwiched inside cell membranes

The transport of nanoparticles at bio-nano interfaces is essential for many cellular responses and biomedical applications. How two-dimensional nanomaterials, such as graphene and transition-metal dichalcogenides, diffuse along the cell membrane is, however, unknown, posing an urgent and important issue to promote their applications in the biomedical area. Here, we show that the transport of graphene oxides (GOs) sandwiched inside cell membranes varies from Brownian to Lévy and even directional dynamics. Specifically, experiments evidence sandwiched graphene-cell membrane superstructures in different cells. Combined simulations and analysis identify a sandwiched GO-induced pore in cell membrane leaflets, spanning unstable, metastable, and stable states. An analytical model that rationalizes the regimes of these membrane-pore states fits simulations quantitatively, resulting in a mechanistic interpretation of the emergence of Lévy and directional dynamics. We finally demonstrate the applicability of sandwiched GOs in enhanced efficiency of membrane-specific drug delivery. Our findings inform approaches to programming intramembrane transport of two-dimensional nanomaterials toward advantageous biomedical applications.

Sugar-based bactericides targeting phosphatidylethanolamine-enriched membranes

Anthrax is an infectious disease caused by Bacillus anthracis, a bioterrorism agent that develops resistance to clinically used antibiotics. Therefore, alternative mechanisms of action remain a challenge. Herein, we disclose deoxy glycosides responsible for specific carbohydrate-phospholipid interactions, causing phosphatidylethanolamine lamellar-to-inverted hexagonal phase transition and acting over B. anthracis and Bacillus cereus as potent and selective bactericides. Biological studies of the synthesized compound series differing in the anomeric atom, glycone configuration and deoxygenation pattern show that the latter is indeed a key modulator of efficacy and selectivity. Biomolecular simulations show no tendency to pore formation, whereas differential metabolomics and genomics rule out proteins as targets. Complete bacteria cell death in 10 min and cellular envelope disruption corroborate an effect over lipid polymorphism. Biophysical approaches show monolayer and bilayer reorganization with fast and high permeabilizing activity toward phosphatidylethanolamine membranes. Absence of bacterial resistance further supports this mechanism, triggering innovation on membrane-targeting antimicrobials.

Bacillus anthracis causes the infectious disease anthrax. Here, the authors synthesized deoxy glycosides that are effective against B. anthracis and related bacteria and found that these amphiphilic compounds kill bacteria via an unusual mechanism of action.

Intracellular Delivery by Membrane Disruption: Mechanisms, Strategies, and Concepts

Intracellular delivery is a key step in biological research and has enabled decades of biomedical discoveries. It is also becoming increasingly important in industrial and medical applications ranging from biomanufacture to cell-based therapies. Here, we review techniques for membrane disruption-based intracellular delivery from 1911 until the present. These methods achieve rapid, direct, and universal delivery of almost any cargo molecule or material that can be dispersed in solution. We start by covering the motivations for intracellular delivery and the challenges associated with the different cargo types-small molecules, proteins/peptides, nucleic acids, synthetic nanomaterials, and large cargo. The review then presents a broad comparison of delivery strategies followed by an analysis of membrane disruption mechanisms and the biology of the cell response. We cover mechanical, electrical, thermal, optical, and chemical strategies of membrane disruption with a particular emphasis on their applications and challenges to implementation. Throughout, we highlight specific mechanisms of membrane disruption and suggest areas in need of further experimentation. We hope the concepts discussed in our review inspire scientists and engineers with further ideas to improve intracellular delivery.

Elucidating mysteries of phase-segregated membranes: mobile-lipid recruitment facilitates pores' passage to the fluid phase

Phase segregation of multicomponent lipid bilayers leads to, under phase-coexistence conditions, domain formation, featuring delimitation by essentially one-dimensional borders. (Micro-)phase segregation of bilayers is proposed to influence the physiological behaviour of cell membranes and provides the driving force for lipid-raft formation. Experiments show a maximum in the electrical-conductivity of membranes at the phase-transition point, which has been conjectured to arise from border-nucleated transmembrane-conducting defects or pores. However, recent electroporation experiments on phase-segregated bilayers demonstrate electro-pore detection in the liquid disordered phase (Ld), wherein they diffuse over macroscopic periods without absorption into the liquid ordered phase (Lo). Here, we scrutinise transmembrane-pore formation via molecular dynamics simulations on a multicomponent phase-segregated bilayer. We find that pores created in Lo domains always migrate spontaneously to the Ld phase, via 'recruitment' of unsaturated lipids to the pore's rim to transport the pore to the fluid phase under a large stress-field driving force. Once in Ld domains, pores migrate towards their centre, never returning or pinning to Lo. These findings are explained by thermodynamics. By comparing the free-energy cost for creating pores in the bulk of Ld and Lo membranes, and in the phase-segregated system, we show that it is always more energetically tractable to create pores in Ld domains, independent of the pore size.

Intracellular Trafficking and Endosomal Release of Oligonucleotides: What We Know and What We Don't

Understanding the cellular uptake and intracellular trafficking of oligonucleotides provides an important basic underpinning for the developing field of oligonucleotide-based therapeutics. Whether delivered as "free" oligonucleotides, as ligand-oligonucleotide conjugates, or in association with various nanocarriers, all forms of oligonucleotide enter cells by endocytosis and are initially ensconced within membrane-limited vesicles. Accordingly, the locus and extent of release to the cytosol and nucleus are key determinants of the pharmacological actions of oligonucleotides. A number of recent studies have explored the intracellular trafficking of various forms of oligonucleotides and their release from endomembrane compartments. These studies reveal a surprising convergence on an early-intermediate compartment in the trafficking pathway as the key locus of release for oligonucleotides administered in "free" form as well as those delivered with lipid complexes. Thus, oligonucleotide release from multivesicular bodies or from late endosomes seems to be the crucial endogenous process for attaining pharmacological effects. This intrinsic process of oligonucleotide release may be amplified by delivery agents such as lipid complexes or small molecule enhancers.

In situ examination of a charged amino acid-induced structural change in lipid bilayers by sum frequency generation vibrational spectroscopy

The interactions between amino acids (AAs) and membranes represent various short-range and long-range interactions for biological phenomena; however, they are still poorly understood. In this study, we used cationic lysine and arginine as AA models, and systematically investigated the interactions between charged AAs and lipid bilayers using sum frequency generation vibrational spectroscopy (SFG-VS) in situ and in real time. The AA-induced dynamic structural changes of the lipid bilayer were experimentally monitored using the spectral features of CD₂, CD₃, the lipid head phosphate, and carbonyl groups in real time. Time-dependent SFG changes in the structure of the lipid bilayer provide direct evidence for the different interactions of lysine and arginine with the membrane. It was found that the discrepancy between lysine and arginine in binding with the lipid bilayer is due to the nature of the terminal functional groups. Arginine exhibits a more drastic impact on the membrane than lysine. SFG responses of the acyl chains, phosphate groups, and carbonyl groups provide evidence that the interaction between AAs and the membrane most likely follows an electrostatics and hydrogen bond-induced defect model. This work presents an exemplary method for comprehensive investigations of interactions between membranes and other functionally significant substances.

The Road from Host-Defense Peptides to a New Generation of Antimicrobial Drugs

Host-defense peptides, also called antimicrobial peptides (AMPs), whose protective action has been used by animals for millions of years, fulfill many requirements of the pharmaceutical industry, such as: (1) broad spectrum of activity; (2) unlike classic antibiotics, they induce very little resistance; (3) they act synergically with conventional antibiotics; (4) they neutralize endotoxins and are active in animal models. However, it is considered that many natural peptides are not suitable for drug development due to stability and biodisponibility problems, or high production costs. This review describes the efforts to overcome these problems and develop new antimicrobial drugs from these peptides or inspired by them. The discovery process of natural AMPs is discussed, as well as the development of synthetic analogs with improved pharmacological properties. The production of these compounds at acceptable costs, using different chemical and biotechnological methods, is also commented. Once these challenges are overcome, a new generation of versatile, potent and long-lasting antimicrobial drugs is expected.

Membrane-Active Properties of an Amphitropic Peptide from the CyaA Toxin Translocation Region

The adenylate cyclase toxin CyaA is involved in the early stages of infection by Bordetella pertussis, the causative agent of whooping cough. CyaA intoxicates target cells by a direct translocation of its catalytic domain (AC) across the plasma membrane and produces supraphysiological levels of cAMP, leading to cell death. The molecular process of AC translocation remains largely unknown, however. We have previously shown that deletion of residues 375–485 of CyaA selectively abrogates AC translocation into eukaryotic cells. We further identified within this “translocation region” (TR), P454 (residues 454–484), a peptide that exhibits membrane-active properties, i.e., is able to bind and permeabilize lipid vesicles. Here, we analyze various sequences from CyaA predicted to be amphipatic and show that although several of these peptides can bind membranes and adopt a helical conformation, only the P454 peptide is able to permeabilize membranes. We further characterize the contributions of the two arginine residues of P454 to membrane partitioning and permeabilization by analyzing the peptide variants in which these residues are substituted by different amino acids (e.g., A, K, Q, and E). Our data shows that both arginine residues significantly contribute, although diversely, to the membrane-active properties of P454, i.e., interactions with both neutral and anionic lipids, helix formation in membranes, and disruption of lipid bilayer integrity. These results are discussed in the context of the translocation process of the full-length CyaA toxin.

Membrane Permeability of Fatty Acyl Compounds Studied via Molecular Simulation

Interest in fatty acid-derived products as fuel and chemical precursors has grown substantially. Microbes can be genetically engineered to produce fatty acid-derived products that are able to cross host membranes and can be extracted into an applied organic overlay. This process is thought to be passive, with a rate dependent on the chemistry of the crossing compound. The relationship between the chemical composition and the energetics and kinetics of product accumulation within the overlay is not well understood. Through biased and unbiased molecular simulation, we compute the membrane permeability coefficients from production to extraction for different fatty acyl products, including fatty acids, fatty alcohols, fatty aldehydes, alkanes, and alkenes. These simulations identify specific interactions that accelerate the transit of aldehydes across the membrane bilayer relative to other oxidized products, specifically the lack of hydrogen bonds to the surrounding membrane environment. However, since extraction from the outer membrane leaflet into the organic phase is found to be rate limiting for the entire process, we find that fatty alcohols and fatty aldehydes would both manifest similar fluxes into a dodecane overlay under equivalent conditions, outpacing the accumulation of acids or alkanes into the organic phase. Since aldehydes are known to be highly reactive as well as toxic in high quantities, the findings suggest that indeed fatty alcohols are the optimal long-tail fatty acyl product for extraction.

Probing a Continuous Polar Defect: A Reaction Coordinate for Pore Formation in Lipid Membranes

Various biophysical processes involve the formation of aqueous pores over lipid membranes, including processes of membrane fusion, antimicrobial peptide activity, lipid flip-flop, and membrane permeation. Reliable and efficient free-energy calculations of pore formation using molecular dynamics simulations remained challenging due to the lack of good reaction coordinates (RCs) for pore formation. We present a new RC for pore formation that probes the formation and rupture of a continuous polar defect over the membrane. Potential of mean force (PMF) calculations along the new RC rapidly converge and exhibit no hysteresis between pore-opening and pore-closing pathways, in contrast to calculations based on previous RCs. We show that restraints along the new RC may restrain the system tightly to the transition state of pore formation, rationalizing the absence of hysteresis. We observe that the PMF of pore formation in a tension-free membrane of dimyristoylphosphatidylcholine (DMPC) reveals a free-energy barrier for pore nucleation, confirming a long-hypothesized metastable prepore state. We test the influence of the lipid force field, the cutoff distance used for Lennard-Jones interactions, and the lateral membrane size on the free energies of pore formation. In contrast to PMF calculations based on previous RCs, we find that such parameters have a relatively small influence on the free energies of pore nucleation. However, the metastability of the open pore in DMPC may depend on such parameters. The RC has been implemented into an extension of the GROMACS simulation software. The new RC allows for reliable and computationally efficient free-energy calculations of pore formation in lipid membranes.

Reconciling Differences between Lipid Transfer in Free-Standing and Solid Supported Membranes: A Time-Resolved Small-Angle Neutron Scattering Study

Maintaining compositional lipid gradients across membranes in animal cells is essential to biological function, but what is the energetic cost to maintain these differences? It has long been recognized that studying the passive movement of lipids in membranes can provide insight into this toll. Confusingly the reported values of inter- and, particularly, intra-lipid transport rates of lipids in membranes show significant differences. To overcome this difficulty, biases introduced by experimental approaches have to be identified. The present study addresses the difference in the reported intramembrane transport rates of dimyristoylphosphatidylcholine (DMPC) on flat solid supports (fast flipping) and in curved free-standing membranes (slow flipping). Two possible scenarios are potentially at play: one is the difference in curvature of the membranes studied and the other the presence (or not) of the support. Using DMPC vesicles and DMPC supported membranes on silica nanoparticles of different radii, we found that an increase in curvature (from a diameter of 30 nm to a diameter of 100 nm) does not change the rates significantly, differing only by factors of order ∼1. Additionally, we found that the exchange rates of DMPC in supported membranes are similar to the ones in vesicles. And as previously reported, we found that the activation energies for exchange on free-standing and supported membranes are similar (84 and 78 kJ/mol, respectively). However, DMPC's flip-flop rates increase significantly when in a supported membrane, surpassing the exchange rates and no longer limiting the exchange process. Although the presence of holes or cracks in supported membranes explains the occurrence of fast lipid flip-flop in many studies, in defect-free supported membranes we find that fast flip-flop is driven by the surface's induced disorder of the bilayer's acyl chain packing as evidenced from their broad melting temperature behavior.

Computational Models of the Gastrointestinal Environment. 1. The Effect of Digestion on the Phase Behavior of Intestinal Fluids

Improved models of the gastrointestinal environment have great potential to assist the complex process of drug formulation. Molecular dynamics (MD) is a powerful method for investigating phase behavior at a molecular level. In this study we use multiple MD simulations to calculate phase diagrams for bile before and after digestion. In these computational models, undigested bile is represented by mixtures of palmitoyl-oleoylphosphatidylcholine (POPC), sodium glycodeoxycholate (GDX), and water. Digested bile is modeled using a 1:1 mixture of oleic acid and palmitoylphosphatidylcholine (lysophosphatidylcholine, LPC), GDX, and water. The computational phase diagrams of undigested and digested bile are compared, and we describe the typical intermolecular interactions that occur between phospholipids and bile salts. The diffusion coefficients measured from MD simulation are compared to experimental diffusion data measured by DOSY-NMR, where we observe good qualitative agreement. In an additional set of simulations, the effect of different ionization states of oleic acid on micelle formation is investigated.

Generalized Langevin Methods for Calculating Transmembrane Diffusivity

The membrane permeability coefficient of a solute can be estimated using the solubility-diffusion model. This model requires the diffusivity profile (D(z)) of the solute as it moves along the transmembrane axis, z. The generalized Langevin equation provides one strategy for calculating position-dependent diffusivity from straightforward molecular dynamics simulations where the solute is restrained to a series of positions on the z-coordinate by a harmonic potential. The diffusivity of the solute is calculated from its correlation functions, which are related to the friction experienced by the solute. Roux and Hummer have derived expressions for the diffusion coefficient from the velocity autocorrelation function (VACF) and position autocorrelation function (PACF), respectively. In this work, these methods are validated by calculating the diffusivity of H₂O and O₂ in homogeneous liquids. These methods are then used to calculate transmembrane diffusivity profiles. The VACF method is less sensitive to thermostat forces and has incrementally lower errors but is more sensitive to the spring constant of the harmonic restraint. For the permeation of a solute through a lipid bilayer, the diffusion coefficients calculated using these methods provided significantly different results. Long-lived correlations of the restrained solute due to inhomogeneities in the bilayer can result in spuriously low diffusivity when using the PACF method. The method based on the VACF does not have this issue and predicts higher rates of diffusion inside the bilayer.

Effect of an Anionic Lipid on the Barotropic Behavior of a Ternary Bicellar Mixture

Dispersions of lipid mixtures comprising long- and short-chain phospholipids (bicellar mixtures) can form small isotropically reorienting particles (bilayered micelles), magnetically orientable stuctures, or unorientable lamellar structures. Application of hydrostatic pressure can also induce interdigitation of the long-chain lipid components. In this work, variable-pressure ²H NMR was used to study the effect of head group charge on the barotropic behavior of bicellar mixtures. Observations at pressures up to 152 MPa and temperatures up to 64 °C were combined with earlier observations at lower pressure and lower temperature to obtain a pressure-temperature phase diagram for DMPC-d₅₄/DMPG/DHPC (3:1:1). In this phase diagram, a region corresponding to small, isotropically reorienting particles at lower temperature and higher pressure is separated from a region corresponding to unorientable lamellar organization, at higher temperature and lower pressure, by a band in which the magnetically orientable phase is stable below ∼100 MPa and in which an interdigitated gel phase is stable above ∼120 MPa. From ∼46 to ∼52 °C, the dispersion transforms directly from the unorientable lamellar to isotropically reorienting particle phases upon isothermal pressurization. The extent to which this behavior reflects the presence of anionic lipid in the long-chain fraction of this mixture is illustrated by comparison with spectral series obtained during isothermal pressurization of DMPC-d₅₄/DHPC (4:1) and DMPC-d₅₄/DMPG/DHPC (2.7:1.3:1) at selected temperatures. These observations show how electrostatic interactions at a bilayer surface can affect the balance between hydrophobic and hydrophilic interactions that is reflected by a dispersion's barotropic phase behavior.

CHARMM-GUI HMMM Builder for Membrane Simulations with the Highly Mobile Membrane-Mimetic Model

Slow diffusion of the lipids in conventional all-atom simulations of membrane systems makes it difficult to sample large rearrangements of lipids and protein-lipid interactions. Recently, Tajkhorshid and co-workers developed the highly mobile membrane-mimetic (HMMM) model with accelerated lipid motion by replacing the lipid tails with small organic molecules. The HMMM model provides accelerated lipid diffusion by one to two orders of magnitude, and is particularly useful in studying membrane-protein associations. However, building an HMMM simulation system is not easy, as it requires sophisticated treatment of the lipid tails. In this study, we have developed CHARMM-GUI HMMM Builder (http://www.charmm-gui.org/input/hmmm) to provide users with ready-to-go input files for simulating HMMM membrane systems with/without proteins. Various lipid-only and protein-lipid systems are simulated to validate the qualities of the systems generated by HMMM Builder with focus on the basic properties and advantages of the HMMM model. HMMM Builder supports all lipid types available in CHARMM-GUI and also provides a module to convert back and forth between an HMMM membrane and a full-length membrane. We expect HMMM Builder to be a useful tool in studying membrane systems with enhanced lipid diffusion.

Thermodynamics of cell-penetrating HIV1 TAT peptide insertion into PC/PS/CHOL model bilayers through transmembrane pores: the roles of cholesterol and anionic lipids

Efficient delivery of pharmaceutically active molecules across cellular membranes using cell penetrating peptides (CPPs), such as the cationic human immunodeficiency virus-1 trans-acting activator of transcription peptide (HIV-1 TAT), continues to attract scientific attention in drug design and disease treatment. Experimental results show that the TAT peptide is not only capable of directly penetrating the biological membrane in a passive manner, but also forming physical, membrane-spanning pores that may facilitate transport. Experiments further show that anionic lipids accelerate peptide permeation within a range of mole percentage composition. In this work, we explored the structures and translocation thermodynamics of the cationic TAT peptide across a series of DPPC/DPPS model membranes with the presence of 0-30 mol% cholesterol. We computed the potentials of the mean force by using umbrella sampling molecular dynamics simulations coupled to the Martini coarse-grained force field. We systematically investigated the roles of cholesterol and anionic lipids (membrane surface charge) in TAT peptide translocation. In qualitative agreement with experimental findings, the barrier heights were significantly reduced in the presence of anionic lipids. A toroidal hydrophilic pore was strongly suggested by membrane structure analysis. Cholesterol stabilizes the liquid-ordered (Lo) phase of membranes and increases the elastic stiffness of bilayers. Consequently, it hinders transmembrane pore formation and thus modulates solute permeability, since the liquid-ordered phase suppresses reorientation of the lipid molecules on simulation time scales. Though cholesterol contributes marginally to the total free energy associated with peptide permeation, the coordination of cholesterol to the peptide weakens more favorable peptide-lipid interactions. The addition of the anionic lipid DPPS to the neutral DPPC bilayer leads to the emergence and further enhancement of an interfacially stable state of the peptide due to the favorable peptide-anionic lipid interactions. Translocation free energy barriers decrease in lockstep with increasing DPPS composition in the model bilayers simulated. Finally, we investigated the size of hydrophilic pores emerging in our simulations, as well as the qualitative mobility of the peptide on the membrane surface.