Membrane pore formation in atomistic and coarse-grained simulations

Free Energy of Membrane Pore Formation and Stability from Molecular Dynamics Simulations

Understanding the molecular mechanisms of pore formation is crucial for elucidating fundamental biological processes and developing therapeutic strategies, such as the design of drug delivery systems and antimicrobial agents. Although experimental methods can provide valuable information, they often lack the temporal and spatial resolution necessary to fully capture the dynamic stages of pore formation. In this study, we present two novel collective variables (CVs) designed to characterize membrane pore behavior, particularly its energetics, through molecular dynamics (MD) simulations. The first CV─termed Full-Path─effectively tracks both the nucleation and expansion phases of pore formation. The second CV─called Rapid─is tailored to accurately assess pore expansion in the limit of large pores, providing quick and reliable method for evaluating membrane line tension under various conditions. Our results clearly demonstrate that the line tension predictions from both our CVs are in excellent agreement. Moreover, these predictions align qualitatively with available experimental data. Specifically, they reflect higher line tension of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membranes containing 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS) lipids compared to pure POPC, the decrease in line tension of POPC vesicles as the 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) content increases, and higher line tension when ionic concentration is increased. Notably, these experimental trends are accurately captured only by the all-atom CHARMM36 and prosECCo75 force fields. In contrast, the all-atom Slipids force field, along with the coarse-grained Martini 2.2, Martini 2.2 polarizable, and Martini 3 models, show varying degrees of agreement with experiments. Our developed CVs can be adapted to various MD simulation engines for studying pore formation, with potential implications in membrane biophysics. They are also applicable to simulations involving external agents, offering an efficient alternative to existing methodologies.

Structural Dynamics of Neutral Amino Acid Transporter SLC6A19 in Simple and Complex Lipid Bilayers

B0AT1 (SLC6A19) is a major sodium-coupled neutral amino acid transporter that relies on angiotensin converting enzyme 2 (ACE2) or collectrin for membrane trafficking. Despite its significant role in disorders associated with amino acid metabolism, there is a deficit of comprehensive structure-function understanding of B0AT1 in lipid environment. Herein, we have employed molecular dynamics (MD) simulations to explore the architectural characteristics of B0AT1 in two distinct environments: a simplified POPC bilayer and a complex lipid system replicating the native membrane composition. Notably, our B0AT1 analysis in terms of structural stability and regions of maximum flexibility shows consistency in both the systems with enhanced structural features in the case of complex lipid system. Our findings suggest that diacylglycerol phospholipids significantly alter the pore radius, hydrophobic index, and surface charge distribution of B0AT1, thereby affecting the flexibility of transmembrane helices TM7, TM12, and loop connecting TM7-TM8, crucial for ACE2-B0AT1 interaction. Pro41, Ser190, Arg214, Arg240, Ser413, Pro414, Cys463, and Val582 are among the most prominent lipid binding residues that might influence B0AT1 functionality. We also perceive notable lipid mediated deviation in the degree of tilt and loss of helicity in TM1 and TM6 which might affect the substrate binding sites S1 and S2 in B0AT1. Considerably, destabilization in the structure of B0AT1 in lipid environment was evident upon mutation in TM domain, associated with Hartnup disorder through various structure-based protein stability tools. Our two-tiered approach allowed us to validate the use of POPC as a baseline for initial analyses of SLC transporters. Altogether, our all-atoms MD study provides a platform for future investigations into the structure-function mechanism of B0AT1 in realistic lipid mimetic bilayers and offers a framework for developing new therapeutic agents targeting this transporter.

Remote loading in liposome: a review of current strategies and recent developments

Liposomes have gained prominence as nanocarriers in drug delivery, and the number of products in the market is increasing steadily, particularly in cancer therapeutics. Remote loading of drugs in liposomes is a significant step in the translation and commercialization of the first liposomal product. Low drug loading and drug leakage from liposomes is a translational hurdle that was effectively circumvented by the remote loading process. Remote loading or active loading could load nearly 100% of the drug, which was not possible with the passive loading procedure. A major drawback of conventional remote loading is that only a very small percentage of the drugs are amenable to this method. Therefore, methods for drug loading are still a problem for several drugs. The loading of multiple drugs in liposomes to improve the efficacy and safety of nanomedicine has gained prominence recently with the introduction of a marketed formulation (Vyxeos) that improves overall survival in acute myeloid leukemia. Different strategies for modifying the remote loading process to overcome the drawbacks of the conventional method are discussed here. The review aims to discuss the latest developments in remote loading technology and its implications in liposomal drug delivery.

Uncoated gold nanoparticles create fewer and less localized defects in model prokaryotic than in model eukaryotic lipid membranes

The rise of the populations of antibiotic resistant bacteria represents an increasing threat to human health. In addition to the synthesis of new antibiotics, which is an extremely expensive and time-consuming process, one of the ways to combat bacterial infections is the use of gold nanoparticles (Au NPs) as the vehicles for targeted delivery of therapeutic drugs. Since such a strategy requires the investigation of the effect of Au NPs (with and without drugs) on both bacterial and human cells, we investigated how the presence of coating-free Au NPs affects the physicochemical properties of lipid membranes that model prokaryotic (PRO) and eukaryotic (EU) cells. PRO/EU systems prepared as multilamellar liposomes (MLVs) and hybrid structures (HSs) from 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphatidylglycerol (DPPG)/1,2-dipalmitoyl-sn-glycero-3-phosphoserine (DPPS) in the absence (MLVs)/presence (HSs) of differently distributed Au NPs (sizes ∼20 nm) reported stabilization of the gel phase of PRO systems in comparison with EU one (DSC data of PRO/EU were Tm(MLVs) ≈ 41.8 °C/42.0 °C, Tm¯ (HSs) ≈ 43.1 °C/42.4 °C, whereas UV-Vis response Tm(MLVs) ≈ 41.5 °C/42.0 °C, Tm¯ (HSs) ≈ 42.9 °C/41.1 °C). Vibrational spectroscopic data unraveled a substantial impact of Au NPs on the non-polar part of lipid bilayers, emphasizing the increase of kink and gauche conformers of the hydrocarbon chain. By interpreting the latter as Au NPs-induced defects, which exert the greatest effect when Au NPs are found exclusively outside the lipid membrane, these findings suggested that Au NPs reduced the compactness of EU-based lipid bilayers much more than in analogous PRO systems. Since the uncoated Au NPs manifested adverse effects when applied as antimicrobials, the results obtained in this work contribute towards recognizing AuNP functionalization as a strategy in tuning and reversing this effect.

Influence of phospholipid head and tail molecular structures on cell membrane mechanical response under tension

Biological cell membranes are primarily comprised of a diverse lipid bilayer with multiple phospholipid (lipid) types, each of which is comprised of a hydrophilic headgroup and two hydrophobic hydrocarbon tails. The lipid type determines the molecular structure of head and tail groups, which can affect membrane mechanics at nanoscale and subsequently cell viability under mechanical loading. Hence, using molecular dynamics simulations, the current study investigated seven membrane phospholipids and the effect of their structural differences on physical deformation, mechanoporation damage, and mechanical failure of the membranes under tension. The inspected phospholipids showed similar yield stresses and strains, as well as pore evolution and damage, but significantly different failure strains. In general, failure occurred at a lower strain for lipids with a larger equilibrium area per lipid. The obtained results suggest that larger headgroup structure, greater degree of unsaturation, and tail-length asymmetry influenced the phospholipids' ability to pack against each other, increased the fluidity and equilibrium area per lipid of the membrane, and resulted in lower failure strain. Overall, this study provides insights on how different phospholipid structures affect membrane physical responses at the molecular level and serves as a reference for future studies of more complex membrane systems with intricate biophysical properties.

Chain-Length Dependence of Peptide-Lipid Bilayer Interaction Strength and Binding Kinetics: A Combined Theoretical and Experimental Approach

Physical interactions between polypeptide chains and lipid membranes underlie critical cellular processes. Yet, despite fundamental importance, key mechanistic aspects of these interactions remain elusive. Bulk experiments have revealed a linear relationship between free energy and peptide chain length in a model system, but does this linearity extend to the interaction strength and to the kinetics of lipid binding? To address these questions, we utilized a combination of coarse-grained molecular dynamics (CG MD) simulations, analytical modeling, and atomic force microscopy (AFM)-based single molecule force spectroscopy. Following previous bulk experiments, we focused on interactions between short hydrophobic peptides (WLn, n = 1, ..., 5) with 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) bilayers, a simple system that probes peptide primary structure effects. Potentials of mean force extracted from CG MD recapitulated the linearity of free energy with the chain length. Simulation results were quantitatively connected to bulk biochemical experiments via a single scaling factor of order unity, corroborating the methodology. Additionally, CG MD revealed an increase in the distance to the transition state, a result that weakens the dependence of the dissociation force on the peptide chain length. AFM experiments elucidated rupture force distributions and, through modeling, intrinsic dissociation rates. Taken together, the analysis indicates a rupture force plateau in the WLn-POPC system, suggesting that the final rupture event involves the last 2 or 3 residues. In contrast, the linear dependence on chain length was preserved in the intrinsic dissociation rate. This study advances the understanding of peptide-lipid interactions and provides potentially useful insights for the design of peptides with tailored membrane-interacting properties.

Effect of Force Field Resolution on Membrane Mechanical Response and Mechanoporation Damage under Deformation Simulations

Damage induced by transient disruption and mechanoporation in an intact cell membrane is a vital nanoscale biomechanical mechanism that critically affects cell viability. To complement experimental studies of mechanical membrane damage and disruption, molecular dynamics (MD) simulations have been performed at different force field resolutions, each of which follows different parameterization strategies and thus may influence the properties and dynamics of membrane systems. Therefore, the current study performed tensile deformation MD simulations of bilayer membranes using all-atom (AA), united-atom (UA), and coarse-grained Martini (CG-M) models to investigate how the damage biomechanics differs across atomistic and coarse-grained (CG) simulations. The mechanical response and mechanoporation damage were qualitatively similar but quantitatively different in the three models, including some progressive changes based on the coarse-graining level. The membranes yielded and reached ultimate strength at similar strains; however, the coarser systems exhibited lower average yield stresses and failure strains. The average failure strain in the UA model was approximately 7% lower than the AA, and the CG-M was 20% lower than UA and 27% lower than AA. The CG systems also nucleated a higher number of pores and larger pores, which resulted in higher damage during the deformation process. Overall, the study provides insight on the impact of force field-a critical factor in modeling biomolecular systems and their interactions-in inspecting membrane mechanosensitive responses and serves as a reference for justifying the appropriate force field for future studies of more complex membranes and more diverse biomolecular assemblies.

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.

Coarse-grained molecular simulation of extracellular vesicle squeezing for drug loading

In recent years, extracellular vesicles have become promising carriers as next-generation drug delivery platforms. Effective loading of exogenous cargos without compromising the extracellular vesicle membrane is a major challenge. Rapid squeezing through nanofluidic channels is a widely used approach to load exogenous cargoes into the EV through the nanopores generated temporarily on the membrane. However, the exact mechanism and dynamics of nanopore opening, as well as cargo loading through nanopores during the squeezing process remains unknown and it is impossible to visualize or quantify it experimentally due to the small size of the EV and the fast transient process. This paper developed a systemic algorithm to simulate nanopore formation and predict drug loading during extracellular vesicle (EV) squeezing by leveraging the power of coarse-grain (CG) molecular dynamics simulations with fluid dynamics. The EV CG beads are coupled with implicit the fluctuating lattice Boltzmann solvent. The effects of EV properties and various squeezing test parameters, such as EV size, flow velocity, channel width, and length, on pore formation and drug loading efficiency are analyzed. Based on the simulation results, a phase diagram is provided as a design guide for nanochannel geometry and squeezing velocity to generate pores on the membrane without damaging the EV. This method can be utilized to optimize the nanofluidic device configuration and flow setup to obtain desired drug loading into EVs.

Interaction of Uperin Peptides with Model Membranes: Molecular Dynamics Study

The interaction of antimicrobial and amyloid peptides with cell membranes is a critical step in their activities. Peptides of the uperin family obtained from the skin secretion of Australian amphibians demonstrate antimicrobial and amyloidogenic properties. All-atomic molecular dynamics and an umbrella sampling approach were used to study the interaction of uperins with model bacterial membrane. Two stable configurations of peptides were found. In the bound state, the peptides in helical form were located right under the head group region in parallel orientation with respect to the bilayer surface. Stable transmembrane configuration was observed for wild-type uperin and its alanine mutant in both alpha-helical and extended unstructured forms. The potential of mean force characterized the process of peptide binding from water to the lipid bilayer and its insertion into the membrane, and revealed that the transition of uperins from the bound state to the transmembrane position was accompanied by the rotation of peptides and passes through the energy barrier of 4-5 kcal/mol. Uperins have a weak effect on membrane properties.

Different effects of two Poloxamers (L61 and F68) on the conductance of bilayer lipid membranes

The integral conductance of planar lipid bilayer membranes in the presence of two Poloxamers (Pluronics) L61 and F68 with the same lengths of hydrophobic poly(propylene oxide) blocks and the different lengths of hydrophilic poly(ethylene oxide) blocks increases with an increase in the concentration of both Pluronics; however, the shape of the conductance-concentration curves is super linear for L61 and sublinear for F68. In the presence of both Pluronics, rare discrete current jumps are observed against the background of continuous current. At high concentrations, the I-V curves of membranes with both L61 and F68 became nonlinear at sufficiently low voltages but differed significantly. At voltages greater than 50 mV, the conductance of membranes with L61 increased sharply and quantized jumps were observed toward higher conductance, which could be interpreted as the appearance of additional pores. On the contrary, the conductance of membranes with F68 decreased and quantized jumps to lower conductance were observed, which could be interpreted as blocking of already existing pores. We attributed the differences in the conductance-concentration and I-V curves of these two Pluronics to their different effects on the dynamics of membrane hydration and, accordingly, on the probability of formation of conducting pores.

The Hidden Intricacies of Aquaporins: Remarkable Details in a Common Structural Scaffold

Evolution turned aquaporins (AQPs) into the most efficient facilitators of passive water flow through cell membranes at no expense of solute discrimination. In spite of a plethora of solved AQP structures, many structural details remain hidden. Here, by combining extensive sequence- and structural-based analysis of a unique set of 20 non-redundant high-resolution structures and molecular dynamics simulations of four representatives, key aspects of AQP stability, gating, selectivity, pore geometry, and oligomerization, with a potential impact on channel functionality, are identified. The general view of AQPs possessing a continuous open water pore is challenged and it is depicted that AQPs' selectivity is not exclusively shaped by pore-lining residues but also by the relative arrangement of transmembrane helices. Moreover, this analysis reveals that hydrophobic interactions constitute the main determinant of protein thermal stability. Finally, a numbering scheme of the conserved AQP scaffold is established, facilitating direct comparison of, for example, disease-causing mutations and prediction of potential structural consequences. Additionally, the results pave the way for the design of optimized AQP water channels to be utilized in biotechnological applications.

Propidium uptake and ATP release in A549 cells share similar transport mechanisms

The lipid bilayer of eukaryotic cells' plasma membrane is almost impermeable to small ions and large polar molecules, but its miniscule basal permeability in intact cells is poorly characterized. This report describes the intrinsic membrane permeability of A549 cells toward the charged molecules propidium (Pr²⁺) and ATP^4-. Under isotonic conditions, we detected with quantitative fluorescence microscopy, a continuous low-rate uptake of Pr (∼150 × 10^-21 moles (zmol)/h/cell, [Pr]_o = 150 μM, 32°C). It was stimulated transiently but strongly by 66% hypotonic cell swelling reaching an influx amplitude of ∼1500 (zmol/h)/cell. The progressive Pr uptake with increasing [Pr]_o (30, 150, and 750 μM) suggested a permeation mechanism by simple diffusion. We quantified separately ATP release with custom wide-field-of-view chemiluminescence imaging. The strong proportionality between ATP efflux and Pr²⁺ influx during hypotonic challenge, and the absence of stimulation of transmembrane transport following 300% hypertonic shock, indicated that ATP and Pr travel the same conductive pathway. The fluorescence images revealed a homogeneously distributed intracellular uptake of Pr not consistent with high-conductance channels expressed at low density on the plasma membrane. We hypothesized that the pathway consists of transiently formed water pores evenly spread across the plasma membrane. The abolition of cell swelling-induced Pr uptake with 500 μM gadolinium, a known modulator of membrane fluidity, supported the involvement of water pores whose formation depends on the membrane fluidity. Our study suggests an alternative model of a direct permeation of ATP (and other molecules) through the phospholipid bilayer, which may have important physiological implications.

The Role of Key Amino Acids in the Antimicrobial Mechanism of a Bacteriocin Model Revealed by Molecular Simulations

The AS-48 bacteriocin is a potent antimicrobial polypeptide with enhanced stability due to its circular sequence of peptidic bonds. The mechanism of biological action is still not well understood in spite of both the elucidation of the molecular structure some years ago and several experiments performed that yielded valuable information about the AS-48 bacterial membrane poration activity. In this work, we present a computational study at an atomistic scale to analyze the membrane disruption mechanism. The process is based on the two-stage model: (1) peptide binding to the bilayer surface and (2) membrane poration due to the surface tension exerted by the peptide. Indeed, the induced membrane tension mechanism is able to explain stable formation of pores leading to membrane disruption. The atomistic detail obtained from the simulations allows one to envisage the contribution of the different amino acids during the poration process. Clustering of cationic residues and hydrophobic interactions between peptide and lipids seem to be essential ingredients in the process. GLU amino acids have shown to enhance the membrane disrupting ability of the bacteriocin. TRP24–TRP24 interactions make also an important contribution in the initial stages of the poration mechanism. The detailed atomistic information obtained from the simulations can serve to better understand bacteriocin structural characteristics to design more potent antimicrobial therapies.

Genetic, cellular, and structural characterization of the membrane potential-dependent cell-penetrating peptide translocation pore

Cell-penetrating peptides (CPPs) allow intracellular delivery of bioactive cargo molecules. The mechanisms allowing CPPs to enter cells are ill-defined. Using a CRISPR/Cas9-based screening, we discovered that KCNQ5, KCNN4, and KCNK5 potassium channels positively modulate cationic CPP direct translocation into cells by decreasing the transmembrane potential (V_m). These findings provide the first unbiased genetic validation of the role of V_m in CPP translocation in cells. In silico modeling and live cell experiments indicate that CPPs, by bringing positive charges on the outer surface of the plasma membrane, decrease the V_m to very low values (–150 mV or less), a situation we have coined megapolarization that then triggers formation of water pores used by CPPs to enter cells. Megapolarization lowers the free energy barrier associated with CPP membrane translocation. Using dyes of varying dimensions in CPP co-entry experiments, the diameter of the water pores in living cells was estimated to be 2 (–5) nm, in accordance with the structural characteristics of the pores predicted by in silico modeling. Pharmacological manipulation to lower transmembrane potential boosted CPP cellular internalization in zebrafish and mouse models. Besides identifying the first proteins that regulate CPP translocation, this work characterized key mechanistic steps used by CPPs to cross cellular membranes. This opens the ground for strategies aimed at improving the ability of cells to capture CPP-linked cargos in vitro and in vivo.

Before a drug can have its desired effect, it must reach its target tissue or organ, and enter its cells. This is not easy because cells are surrounded by the plasma membrane, a fat-based barrier that separates the cell from its external environment. The plasma membrane contains proteins that act as channels, shuttling specific molecules in and out of the cell, and it also holds charge, with its inside surface being more negatively charged than its outside surface.

Cell-penetrating peptides are short sequences of amino acids (the building blocks that form proteins) that carry positive charges. These positive charges allow them to cross the membrane easily, but it is not well understood how.

To find out how cell-penetrating peptides cross the membrane, Trofimenko et al. attached them to dyes of different sizes. This revealed that the cell-penetrating peptides enter the cell through temporary holes called water pores, which measure about two nanometres across. The water pores form when the membrane becomes ‘megapolarized’, this is, when the difference in charge between the inside and the outside of the membrane becomes greater than normal. This can happen when the negative charge on the inside surface or the positive charge on the outer surface of the membrane increase. Megapolarization depends on potassium channels, which transport positive potassium ions outside the cell, making the outside of the membrane positive. When cell-penetrating peptides arrive at the outer surface of the cell near potassium channels, they make it even more positive. This increases the charge difference between the inside and the outside of the cell, allowing water pores to form. Once the peptides pass through the pores, the charge difference between the inside and the outside of the cell membrane dissipates, and the pores collapse.

Drug developers are experimenting with attaching cell-penetrating peptides to drugs to help them get inside their target cells. Currently there are several experimental medications of this kind in clinical trials. Understanding how these peptides gain entry, and what size of molecule they could carry with them, provides solid ground for further drug development.

Virus-inspired designs of antimicrobial nanocapsules

Antimicrobial resistance is becoming a serious burden for drug design. The challenges are in finding novel approaches for effectively targeting a number of different bacterial strains, and in delivering these to the site of action. We propose here a novel approach that exploits the assembly of antimicrobial peptidic units in nanocapsules that can penetrate and rupture the bacterial membrane. Additionally, the chemical versatility of the designed units can be tailored to specific targets and to the delivery of genetic material in the cell. The proposed design exploits a β-annulus (sequence ITHVGGVGGSIMAPVAVSRQLVGS) triskelion unit from the Tomato Bushy Stunt Virus, able to self assemble in solution, and functionalised with antimicrobial sequences to form dodecahedral antimicrobial nanocapsules. The stability and the activity of the antimicrobial β-annulus capsule is measured by molecular dynamics simulations in water and in the presence of model membranes.

In silico study of PEI-PEG-squalene-dsDNA polyplex formation: the delicate role of the PEG length in the binding of PEI to DNA

Biocompatible hydrophilic polyethylene glycol (PEG) is widely used in biomedical applications, such as drug or gene delivery, tissue engineering or as an antifouling component in biomedical devices. Experimental studies have shown that the size of PEG can weaken polycation-polyanion interactions, like those between branched polyethyleneimine (b-PEI) and DNA in gene carriers, but details of its cause and underlying interactions on the atomic scale are still not clear. To better understand the interaction mechanisms in the formation of polyplexes between b-PEI-PEG based carriers and DNA, we have used a combination of in silico tools and experiments on three multicomponent systems differing in PEG MW. Using the PEI-PEG-squalene-dsDNA systems of the same size, both in the all-atom MD simulations and in experimental in-gel electrophoresis measurements, we found that the binding between DNA and the vectors is highly influenced by the size of PEG, with the binding efficiency increasing with a shorter PEG length. The mechanism of how PEG interferes with the binding between PEI and DNA is explained using a two-step MD simulation protocol that showed that the DNA-vector interactions are influenced by the PEG length due to the hydrogen bond formation between PEI and PEG. Although computationally demanding we find it important to study molecular systems of the same size both in silico and in a laboratory and to simulate the behaviour of the carrier prior to the addition of bioactive molecules to understand the molecular mechanisms involved in the formation of the polyplex.

Molecular Mechanism of Ultrasound-Induced Structural Defects in Liposomes: A Nonequilibrium Molecular Dynamics Simulation Study

The use of ultrasound in combination with liposomes is a promising approach to improve drug delivery. To achieve an optimal drug release rate, it is important to understand how ultrasound induces pathways on the liposome surface where drugs can be released from the liposome. To this end, we carry out large-scale ultrasound-induced molecular dynamics simulations for three single lipid component liposomes formed from the commonly used phospholipids: 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), 1,2-dipalmitoylphosphatidylcholine (DPPC), or phosphatidylcholine (POPC). The results show that ultrasound induces the detachment of two leaflets of the DOPC surface, suggesting that the drug release pathway may be through the low lipid packing areas on the stretched surface. In contrast, ultrasound induces pore formation on the surface of DPPC and DOPC, where drugs could escape from the liposomes. While the leaflet detachment and transient pore formation are the mechanisms of DOPC and DPPC, respectively, in both liquid-ordered and liquid-disordered phases, the leaflet detachment mechanism is switched to the transient pore formation mechanism on going from the liquid-ordered phase to the liquid-disordered phase in the POPC liposome. By adding 30% mol cholesterol, the leaflet detachment mechanism is observed in all liposomes. We found that the molecular origin that determines a mechanism is the competition between the intraleaflet and interleaflet interacting energy of lipids. The connection to experimental and theoretical modeling is discussed in some detail.

Mutually Beneficial Combination of Molecular Dynamics Computer Simulations and Scattering Experiments

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

Effect of short-chain chlorinated paraffins (SCCPs) on lipid membranes: Combination of molecular dynamics and membrane damage experiments

In recent years, more attention has been paid to the biological effects of short-chain chlorinated paraffin (SCCP). Studies have shown that SCCPs exposure could cause metabolic damage and lipid metabolic damage. In the present work, based on E. coli membrane damage experiments and molecular dynamics (MD) simulation, the effects of SCCPs on the membrane structure and membrane properties were studied to explore the possible toxic damage effects of SCCPs on cell membrane. Experiments results showed that SCCPs had a significant inhibitory effect on E. coli. The E. coli cell membrane of the bacteria was broken and the macromolecules of the cell flowed out when exposed to SCCPs. SCCPs would lead to the decrease and depolarization of cell membrane potential, and then affect the integrity and permeability of cell membrane. The further molecular dynamic simulation revealed that SCCP molecules can easily enter the lipid DPPC membranes from the aqueous phase and tended to aggregate inside bilayer stably. The bound of SCCPs could lead to significant variations in DPPC bilayer with a less dense, more disorder and rougher layer, which thus made the damage of cell membrane. In a word, although the overall toxicity of SCCPs to cell was relatively weak, the damage to the cell membrane may be one of the mechanisms of its toxicity. MAIN FINDING OF THE WORK: The exposure of SCCPs could cause structural change of cell membrane in E. coli, which verified the damage to the cell membrane may be one of the mechanisms of its toxicity.

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.

Modification of lipid membrane compressibility induced by an electric field

Changes in membrane deformation and compressibility, induced by an external electric field, are investigated using coarse-grained martini force field simulations in a salt-free environment. We observe changes in the area of the membrane above a critical electric field. Below this value, the membrane compressibility modulus is found to decrease monotonically. For higher electric fields, the membrane projected area remains constant while the net interfacial area increases, with the corresponding compressibility moduli, show the opposite behavior. We find that the mechanical parameters, surface tension and bending modulus, of a freely floating membrane in the absence of explicit ions, are unaffected by the presence of the electric field. We believe these results have a bearing on our understanding of the electroformation of uncharged lipids in a salt-free environment.

Spontaneous Membrane Nanodomain Formation in the Absence or Presence of the Neurotransmitter Serotonin

Detailed knowledge on the formation of biomembrane domains, their structure, composition, and physical characteristics is scarce. Despite its frequently discussed importance in signaling, e.g., in obtaining localized non-homogeneous receptor compositions in the plasma membrane, the nanometer size as well as the dynamic and transient nature of domains impede their experimental characterization. In turn, atomistic molecular dynamics (MD) simulations combine both, high spatial and high temporal resolution. Here, using microsecond atomistic MD simulations, we characterize the spontaneous and unbiased formation of nano-domains in a plasma membrane model containing phosphatidylcholine (POPC), palmitoyl-sphingomyelin (PSM), and cholesterol (Chol) in the presence or absence of the neurotransmitter serotonin at different temperatures. In the ternary mixture, highly ordered and highly disordered domains of similar composition coexist at 303 K. The distinction of domains by lipid acyl chain order gets lost at lower temperatures of 298 and 294 K, suggesting a phase transition at ambient temperature. By comparison of domain ordering and composition, we demonstrate how the domain-specific binding of the neurotransmitter serotonin results in a modified domain lipid composition and a substantial downward shift of the phase transition temperature. Our simulations thus suggest a novel mode of action of neurotransmitters possibly of importance in neuronal signal transmission.

Sequential Water and Headgroup Merger: Membrane Poration Paths and Energetics from MD Simulations

Membrane topology changes such as poration, stalk formation, and hemifusion rupture are essential to cellular function, but their molecular details, energetics, and kinetics are still not fully understood. Here, we present a unified energetic and mechanistic picture of metastable pore defects in tensionless lipid membranes. We used an exhaustive committor analysis to test and select optimal reaction coordinates and also to determine the nucleation mechanism. These reaction coordinates were used to calculate free-energy landscapes that capture the full process and end states. The identified barriers agree with the committor analysis. To enable sufficient sampling of the complete transition path for our molecular dynamics simulations, we developed a “gizmo” potential biasing scheme. The simulations suggest that the essential step in the nucleation is the initial merger of lipid headgroups at the nascent pore center. To facilitate this event, an indentation pathway is energetically preferred to a hydrophobic defect. Continuous water columns that span the indentation were determined to be on-path transients that precede the nucleation barrier. This study gives a quantitative description of the nucleation mechanism and energetics of small metastable pores and illustrates a systematic approach to uncover the mechanisms of diverse cellular membrane remodeling processes.

Aggregation State of Synergistic Antimicrobial Peptides

By integrating various simulation and experimental techniques, we discovered that antimicrobial peptides (AMPs) may achieve synergy at an optimal concentration and ratio, which can be caused by aggregation of the synergistic peptides. On multiple time and length scales, our studies obtain novel evidence of how peptide coaggregation in solution can affect the disruption of membranes by synergistic AMPs. Our findings provide crucial details about the complex molecular origins of AMP synergy, which will help guide the future development of synergistic AMPs as well as applications of anti-infective peptide cocktail therapies.

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.

The Effect of KcsA Channel on Lipid Bilayer Electroporation Induced by Picosecond Pulse Trains

Membrane proteins are the major component of plasma membranes, and they play crucial roles in all organisms. To understand the influence of the presence of KcsA channel on cell membrane electroporation induced by picosecond pulse trains (psPT), in this paper, the electroporation of KcsA membrane protein system and bare lipid bilayer system (POPC) with the applied psPT are simulated using molecular dynamics (MD) method. First, we find that the average pore formation time of the KcsA system is longer than the bare system with the applied psPT. In the KcsA system, water protrusions appear more slowly. Then, the system size effects of psPT in the MD simulations are investigated. When the system size decreases, the average pore formation time of small KcsA membrane protein system is shorter than the bare system with the applied psPT. It is found that the psPT makes the protein fluctuation of small system increase greatly; meanwhile the instability of protein disturbs the water and then affects the water protrusion appearance time. Furthermore, it shows that the protein fluctuation of constant electric field is smaller than that of psPT and no field, and protein fluctuation increases with the psPT repetition frequency increasing.

Effect of helical kink in antimicrobial peptides on membrane pore formation

Every cell is protected by a semipermeable membrane. Peptides with the right properties, for example Antimicrobial peptides (AMPs), can disrupt this protective barrier by formation of leaky pores. Unfortunately, matching peptide properties with their ability to selectively form pores in bacterial membranes remains elusive. In particular, the proline/glycine kink in helical peptides was reported to both increase and decrease antimicrobial activity. We used computer simulations and fluorescence experiments to show that a kink in helices affects the formation of membrane pores by stabilizing toroidal pores but disrupting barrel-stave pores. The position of the proline/glycine kink in the sequence further controls the specific structure of toroidal pore. Moreover, we demonstrate that two helical peptides can form a kink-like connection with similar behavior as one long helical peptide with a kink. The provided molecular-level insight can be utilized for design and modification of pore-forming antibacterial peptides or toxins.

Chemical manipulations to facilitate membrane blebbing and vesicle shedding on the cellular cortex

OBJECTIVES

Most attention has been focused on physiologically generated membrane blebs on the cellular cortex, whereas artificial membrane blebs induced by chemicals are studied to a lesser extent.

RESULTS

We found that exposure of HeLa human cervical cancer cells to paraformaldehyde (PFA), followed by incubation in phosphate-buffered saline (PBS) efficiently induced large membrane blebs on the cellular cortex. Intriguingly, sequential exposure of the PFA-treated cells to PBS containing dimethyl sulfoxide (DMSO) facilitated shedding of the blebs from the cellular cortex, yielding a high quantity of large extracellular vesicles in the supernatant, which was applicable to assess the potentials of compounds and proteins as membrane influencers. Similar effects of PFA and DMSO were detected on the cellular cortex of other human, mouse, and fish cells.

CONCLUSIONS

Our procedure to facilitate membrane blebbing and vesicle shedding by chemicals may be practical for the manipulation of membrane dynamics and the development of vesicle-inspired technologies using a wide variety of cell types.

Computer simulations of protein-membrane systems

The interactions between proteins and membranes play critical roles in signal transduction, cell motility, and transport, and they are involved in many types of diseases. Molecular dynamics (MD) simulations have greatly contributed to our understanding of protein-membrane interactions, promoted by a dramatic development of MD-related software, increasingly accurate force fields, and available computer power. In this chapter, we present available methods for studying protein-membrane systems with MD simulations, including an overview about the various all-atom and coarse-grained force fields for lipids, and useful software for membrane simulation setup and analysis. A large set of case studies is discussed.

Increase in the current variance in bilayer lipid membranes near phase transition as a result of the occurrence of hydrophobic defects

Most researchers associate the increase in the permeability of lipid bilayers of artificial and biological membranes observed in various experiments with the formation of hypothetical hydrophobic and hydrophilic pores. Although the existence of hydrophobic defects, as the first stage of the formation of a hydrophilic pore, was hypothesized decades ago from electroporation experiments, the difficulty of describing this stage is determined by the lack of experimental data confirming the existence or at least associated with hydrophobic pores. We explored the increase in the current variance through the lipid membrane, observed when approaching the phase transition from the side of high temperatures, and have associated it with capacitive currents arising in response to the formation of hydrophobic pores. Assuming that the number of hydrophobic pores in a membrane follows a Poisson distribution, and thus, the mean number of hydrophobic pores is equal to the variance of that number, we used the measurements of the membrane current variance to evaluate the number of hydrophobic pores. Analysis of experimental data within this model allows us to estimate the number of hydrophobic pores at or above the phase transition and shows that the number of hydrophobic pores in a membrane close to the phase transition increased 20 times compared to the number of hydrophobic pores existing in the membrane far from the melting transition.

Modelling lipid systems in fluid with Lattice Boltzmann Molecular Dynamics simulations and hydrodynamics

In this work we present the coupling between Dry Martini, an efficient implicit solvent coarse-grained model for lipids, and the Lattice Boltzmann Molecular Dynamics (LBMD) simulation technique in order to include naturally hydrodynamic interactions in implicit solvent simulations of lipid systems. After validating the implementation of the model, we explored several systems where the action of a perturbing fluid plays an important role. Namely, we investigated the role of an external shear flow on the dynamics of a vesicle, the dynamics of substrate release under shear, and inquired the dynamics of proteins and substrates confined inside the core of a vesicle. Our methodology enables future exploration of a large variety of biological entities and processes involving lipid systems at the mesoscopic scale where hydrodynamics plays an essential role, e.g. by modulating the migration of proteins in the proximity of membranes, the dynamics of vesicle-based drug delivery systems, or, more generally, the behaviour of proteins in cellular compartments.

Chirality-Dependent Adsorption between Amphipathic Peptide and POPC Membrane

The interactions between chiral molecules and cell membranes have attracted more and more attention in recent decades, due to their importance in molecular science and medical applications. It is observed that some peptides composed of different chiral amino acids may have distinct interactions with a membrane. How does the membrane exhibit a selective behavior related to the chirality of the peptides? Microscopically, the interactions between the peptides and the membrane are poorly understood. In this work, we study the interactions between an amphipathic peptide (C6) and POPC membrane with simulations. The kinetics and thermodynamics of peptide enantiomers during the adsorption to the membrane are characterized with direct simulations and umbrella sampling. It is observed that there are slow kinetics for the peptide composed of D-type amino acids. Along the observed pathways, the free energy landscapes are determined with umbrella sampling techniques. A free-energy barrier for the peptide composed of D-amino acids is observed, which is consistent with the kinetic observations. The results indicate the concurrent adsorption and rotation of the peptide helix. The local interactions between the peptides and the membrane are examined in detail, including the contact interactions between the peptides and the membrane, and the distributions of the lipids around the peptide. There are observable differences of the local interactions for the cases related to different peptide enantiomers. These results further demonstrate the importance of the rotation of peptide helix during the adsorption. More interestingly, all these kinetic differences between peptide enantiomers can be explained based on the conformations of the residue Trp and interactions between Trp and lipid molecules. These results give us a molecular understanding of the mechanism of the chirality-dependent peptide-membrane interactions, and may provide clues to designing systems which are sensitive to the chirality of membranes.

Fat SIRAH: Coarse-Grained Phospholipids To Explore Membrane-Protein Dynamics

The capability to handle highly heterogeneous molecular assemblies in a consistent manner is among the greatest challenges faced when deriving simulation parameters. This is particularly the case for coarse-grained (CG) simulations in which chemical functional groups are lumped into effective interaction centers for which transferability between different chemical environments is not guaranteed. Here, we introduce the parametrization of a set of CG phospholipids compatible with the latest version of the SIRAH force field for proteins. The newly introduced lipid species include different acylic chain lengths and partial unsaturation, as well as polar and acidic head groups that show a very good reproduction of structural membrane determinants, such as areas per lipid, thickness, order parameter, etc., and their dependence with temperature. Simulation of membrane proteins showed unprecedented accuracy in the unbiased description of the thickness-dependent membrane-protein orientation in systems where this information is experimentally available (namely, the SarcoEndoplasmic Reticulum Calcium-SERCA-pump and its regulator Phospholamban). The interactions that lead to this faithful reproduction can be traced down to the single amino acid-lipid interaction level and show full agreement with biochemical data present in the literature. Finally, the present parametrization is implemented in the GROMACS and AMBER simulation packages facilitating its use by a wide portion of the biocomputing community.

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.

Interaction mechanism between the focused ultrasound and lipid membrane at the molecular level

Focused ultrasound (FUS) has a wide range of medical applications. Nowadays, the diagnostic and therapeutic ultrasound procedures are routinely used; effects of ultrasound on biological systems at the molecular level are, however, not fully understood. Experimental results on the interaction of the cell membrane, a simplest but important system component, with ultrasound are controversial. Molecular dynamics (MD) simulations could provide valuable insights, but there is no single study on the mechanism of the FUS induced structural changes in cell membranes. With this in mind, we develop a simple method to include FUS into a standard MD simulation. Adopting the 1,2-dioleoyl-sn-glycero-3-phosphocholine lipid membrane as a representative model described by the MARTINI coarse-grained force field, and using experimental values of the ultrasound frequency and intensity, we show that the heat and bubble cavitation are not the primary direct mechanisms that cause structural changes in the membrane. The spatial pressure gradients between the focused and free regions and between the parallel and perpendicular directions to the membrane are the origin of the mechanism. These gradients force lipids to move out of the focused region, forming a lipid flow along the membrane diagonal. Lipids in the free region move in the opposite direction due to the conservation of the total momentum. These opposite motions create wrinkles along the membrane diagonal at low FUS intensities and tear up the membrane at high FUS intensities. Once the membrane is torn up, it is not easy to reform. The implication of our findings in the FUS-induced drug delivery is discussed in some detail.

Coupling of Membrane Nanodomain Formation and Enhanced Electroporation near Phase Transition

Biological cells are enveloped by a heterogeneous lipid bilayer that prevents the uncontrolled exchange of substances between the cell interior and its environment. In particular, membranes act as a continuous barrier for salt and macromolecules to ensure proper physiological functions within the cell. However, it has been shown that membrane permeability strongly depends on temperature and, for phospholipid bilayers, displays a maximum at the transition between the gel and fluid phase. Here, extensive molecular dynamics simulations of dipalmitoylphosphatidylcholine bilayers were employed to characterize the membrane structure and dynamics close to phase transition, as well as its stability with respect to an external electric field. Atomistic simulations revealed the dynamic appearance and disappearance of spatially related nanometer-sized thick ordered and thin interdigitating domains in a fluid-like bilayer close to the phase transition temperature (Tm). These structures likely represent metastable precursors of the ripple phase that vanished at increased temperatures. Similarly, a two-phase bilayer with coexisting gel and fluid domains featured a thickness minimum at the interface because of splaying and interdigitating lipids. For all systems, application of an external electric field revealed a reduced bilayer stability with respect to pore formation for temperatures close to Tm. Pore formation occurred exclusively in thin interdigitating membrane nanodomains. These findings provide a link between the increased membrane permeability and the structural heterogeneity close to phase transition.

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

Membrane-mediated interaction of amphipathic peptides can be described by a one-dimensional approach

Amphipathic alpha-helical peptides, among other peripheral components of plasma membranes, are promising antimicrobial agents. Partial incorporation of a peptide into a lipid monolayer causes elastic deformations. Deformations induced by two peptides distant from each other are independent; when peptides get closer, interference between the deformations causes effective lateral interaction. We quantified the energy of membrane deformations for arbitrary configuration of two amphipathic peptides on the membrane surface. The global minimum of the deformation energy proved to be achieved when two parallel peptides are in registry at the distance of about 6 nm between the axes of peptides. The energy calculated in the unidimensional approach provides a good approximation for the dependence of the energy of peptides being in the registered configuration upon the distance between them, valid for a broad range of peptide lengths. The effective interactional length of peptides for the unidimensional approach is close to their actual length. If two parallel peptides are shifted along their axes with respect to each other, the interaction energy is also well approximated by the unidimensional potential, within the projection of one peptide onto the other. In the case when the axes of alpha helices cross at a substantial angle, the main contribution to peptide interactions comes from their edges: the effective length of peptides for the unidimensional approach is almost equal to the characteristic length of decay of deformations. Based on the results we obtained it can be concluded that interaction of membrane inclusions is quite adequately described by the potential calculated in the unidimensional approach.

Membrane Electroporation and Electropermeabilization: Mechanisms and Models

Exposure of biological cells to high-voltage, short-duration electric pulses causes a transient increase in their plasma membrane permeability, allowing transmembrane transport of otherwise impermeant molecules. In recent years, large steps were made in the understanding of underlying events. Formation of aqueous pores in the lipid bilayer is now a widely recognized mechanism, but evidence is growing that changes to individual membrane lipids and proteins also contribute, substantiating the need for terminological distinction between electroporation and electropermeabilization. We first revisit experimental evidence for electrically induced membrane permeability, its correlation with transmembrane voltage, and continuum models of electropermeabilization that disregard the molecular-level structure and events. We then present insights from molecular-level modeling, particularly atomistic simulations that enhance understanding of pore formation, and evidence of chemical modifications of membrane lipids and functional modulation of membrane proteins affecting membrane permeability. Finally, we discuss the remaining challenges to our full understanding of electroporation and electropermeabilization.

Terahertz Electric Field-Induced Membrane Electroporation by Molecular Dynamics Simulations

In this paper, the membrane electroporation induced by the terahertz electric field is simulated by means of the molecular dynamics method. The influences of the waveform and frequency of the applied terahertz electric field on the electroporation and the unique features of the process of the electroporation with the applied terahertz electric field are given. It shows that whether the electroporation can happen depends on the waveform of the applied terahertz electric field when the magnitude is not large enough. No pore appears if the terahertz electric field direction periodically reverses, and dipole moments of the interfacial water and the bulk water keep reversing. The nm-scale single pore forms with the applied terahertz trapezoidal electric field. It is found that the average pore formation time is strongly influenced by the terahertz electric field frequency. An abnormal variation region that shows decline exists on the correlation curve of the average pore formation time and the trapezoidal electric field frequency, whereas the overall trend of the curve is increasing. The decrease of the water oriented polarization degree results in the increase of the electroporation time, and the abnormal variation region appearance may be related to the drastic change of average water hydrogen bond number that is resulted from the resonance of water hydrogen bond network and the applied electric field. Compared to the nanosecond electric pulse and constant electric field, the numbers of the water protrusions and the water bridges are smaller and the pore formation time is relatively longer with the applied terahertz electric field.

Computational studies of peptide-induced membrane pore formation

A variety of peptides induce pores in biological membranes; the most common ones are naturally produced antimicrobial peptides (AMPs), which are small, usually cationic, and defend diverse organisms against biological threats. Because it is not possible to observe these pores directly on a molecular scale, the structure of AMP-induced pores and the exact sequence of steps leading to their formation remain uncertain. Hence, these questions have been investigated via molecular modelling. In this article, we review computational studies of AMP pore formation using all-atom, coarse-grained, and implicit solvent models; evaluate the results obtained and suggest future research directions to further elucidate the pore formation mechanism of AMPs.This article is part of the themed issue 'Membrane pores: from structure and assembly, to medicine and technology'.

Unperturbed hydrocarbon chains and liquid phase bilayer lipid chains: a computer simulation study

In this work, the properties of saturated and unsaturated fatty acid acyl chains 16:0, 18:0, 18:1(n-9)cis, 18:2(n-6)cis, 18:3(n-3)cis, 18:4(n-3)cis, 18:5(n-3)cis, 20:4(n-6)cis, 20:5(n-3)cis and 22:6(n-3)cis in a bilayer liquid crystalline state and similar hydrocarbon chains (with CH[Formula: see text] terminal groups instead of C=O groups) in the unperturbed state characterised by a lack of long-range interaction were investigated. The unperturbed hydrocarbon chains were modelled by Monte Carlo simulations at temperature [Formula: see text] K; sixteen fully hydrated homogeneous liquid crystalline phosphatidylcholine bilayers containing these chains were studied by molecular dynamics simulations at the same temperature. To eliminate effects of the simulation parameters, the molecular dynamics and Monte Carlo simulations were carried out using the same structural data and force field coefficients. From these computer simulations, the average distances between terminal carbon atoms of the chains (end-to-end distances) were calculated and compared. The trends in the end-to-end distances obtained for the unperturbed chains were found to be qualitatively similar to those obtained for the same lipid chains in the bilayers. So, for understanding of a number of processes in biological membranes (e.g., changes in fatty acid composition caused by environmental changes such as temperature and pressure), it is possible to use, at least as a first approximation, the relationships between the structure and properties for unperturbed or isolated hydrocarbon chains.

Membrane Active Antimicrobial Peptides: Translating Mechanistic Insights to Design

Antimicrobial peptides (AMPs) are promising next generation antibiotics that hold great potential for combating bacterial resistance. AMPs can be both bacteriostatic and bactericidal, induce rapid killing and display a lower propensity to develop resistance than do conventional antibiotics. Despite significant progress in the past 30 years, no peptide antibiotic has reached the clinic yet. Poor understanding of the action mechanisms and lack of rational design principles have been the two major obstacles that have slowed progress. Technological developments are now enabling multidisciplinary approaches including molecular dynamics simulations combined with biophysics and microbiology toward providing valuable insights into the interactions of AMPs with membranes at atomic level. This has led to increasingly robust models of the mechanisms of action of AMPs and has begun to contribute meaningfully toward the discovery of new AMPs. This review discusses the detailed action mechanisms that have been put forward, with detailed atomistic insights into how the AMPs interact with bacterial membranes. The review further discusses how this knowledge is exploited toward developing design principles for novel AMPs. Finally, the current status, associated challenges, and future directions for the development of AMP therapeutics are discussed.

Membrane negative curvature induced by a hybrid peptide from pediocin PA-1 and plantaricin 149 as revealed by atomistic molecular dynamics simulations

Antimicrobial peptides (AMPs) are cationic peptides that kill bacteria with a broad spectrum of action, low toxicity to mammalian cells and exceptionally low rates of bacterial resistance. These features have led to considerable efforts in developing AMPs as an alternative antibacterial therapy. In vitro studies have shown that AMPs interfere with membrane bilayer integrity via several possible mechanisms, which are not entirely understood. We have performed the synthesis, membrane lysis measurements, and biophysical characterization of a novel hybrid peptide. These measurements show that PA-Pln149 does not form nanopores, but instead promotes membrane rupture. It causes fast rupture of the bacterial model membrane (POPG-rich) at concentrations 100-fold lower than that required for the disruption of mammalian model membranes (POPC-rich). Atomistic molecular dynamics (MD) simulations were performed for single and multiple copies of PA-Pln149 in the presence of mixed and pure POPC/POPG bilayers to investigate the concentration-dependent membrane disruption by the hybrid peptide. These simulations reproduced the experimental trend and provided a potential mechanism of action for PA-Pln149. It shows that the PA-Pln149 does not form nanopores, but instead promotes membrane destabilization through peptide aggregation and induction of membrane negative curvature with the collapse of the lamellar arrangement. The sequence of events depicted for PA-Pln149 may offer insights into the mechanism of action of AMPs previously shown to induce negative deformation of membrane curvature and often associated with peptide translocation via non-bilayer intermediate structures.