Coarse-grained simulation studies of peptide-induced pore formation

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.

Planar aggregation of the influenza viral fusion peptide alters membrane structure and hydration, promoting poration

To infect, enveloped viruses employ spike protein, spearheaded by its amphipathic fusion peptide (FP), that upon activation extends out from the viral surface to embed into the target cellular membrane. Here we report that synthesized influenza virus FPs are membrane active, generating pores in giant unilamellar vesicles (GUV), and thus potentially explain both influenza virus' hemolytic activity and the liposome poration seen in cryo-electron tomography. Experimentally, FPs are heterogeneously distributed on the GUV at the time of poration. Consistent with this heterogeneous distribution, molecular dynamics (MD) simulations of asymmetric bilayers with different numbers of FPs in one leaflet show FP aggregation. At the center of FP aggregates, a profound change in the membrane structure results in thinning, higher water permeability, and curvature. Ultimately, a hybrid bilayer nanodomain forms with one lipidic leaflet and one peptidic leaflet. Membrane elastic theory predicts a reduced barrier to water pore formation when even a dimer of FPs thins the membrane as above, and the FPs of that dimer tilt, to continue the leaflet bending initiated by the hydrophobic mismatch between the FP dimer and the surrounding lipid.

Finite-Size Effects in Simulations of Peptide/Lipid Assembly

Abstract

Molecular dynamics simulations are an attractive tool for understanding lipid/peptide self-assembly but can be plagued by inaccuracies when the system sizes are too small. The general guidance from self-assembly simulations of homogeneous micelles is that the total number of surfactants should be three to five times greater than the equilibrium aggregate number of surfactants per micelle. Herein, the heuristic is tested on the more complicated self-assembly of lipids and amphipathic peptides using the Cooke and Martini 3 coarse-grained models. Cooke model simulations with 50 to 1000 lipids and no peptide are dominated by finite-size effects, with usually one aggregate (micelle or nanodisc) containing most of the lipids forming at each system size. Approximately 200 systems of different peptide/lipid (P/L) ratios and sizes of up to 1000 lipids yield a “finite-size phase diagram” for peptide driven self-assembly, including a coexistence region of micelles and discs. Insights from the Cooke model are applied to the assembly of dimyristoylphosphatidylcholine and the ELK-neutral peptide using the Martini 3 model. Systems of 150, 450, and 900 lipids with P/L = 1/6.25 form mixtures of lipid-rich discs that agree in size with experiment and peptide-rich micelles. Only the 150-lipid system shows finite-size effects, which arise from the long-tailed distribution of aggregate sizes. The general rule of three to five times the equilibrium aggregate size remains a practical heuristic for the Cooke and Martini 3 systems investigated here.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s00232-022-00255-9.

Enhanced translocation of amphiphilic peptides across membranes by transmembrane proteins

Cell membranes are phospholipid bilayers with a large number of embedded transmembrane proteins. Some of these proteins, such as scramblases, have properties that facilitate lipid flip-flop from one membrane leaflet to another. Scramblases and similar transmembrane proteins could also affect the translocation of other amphiphilic molecules, including cell-penetrating or antimicrobial peptides. We studied the effect of transmembrane proteins on the translocation of amphiphilic peptides through the membrane. Using two very different models, we consistently demonstrate that transmembrane proteins with a hydrophilic patch enhance the translocation of amphiphilic peptides by stabilizing the peptide in the membrane. Moreover, there is an optimum amphiphilicity because the peptide could become overstabilized in the transmembrane state, in which the peptide-protein dissociation is hampered, limiting the peptide translocation. The presence of scramblases and other proteins with similar properties could be exploited for more efficient transport into cells. The described principles could also be utilized in the design of a drug-delivery system by the addition of a translocation-enhancing peptide that would integrate into the membrane.

Stabilizing Leaflet Asymmetry under Differential Stress in a Highly Coarse-Grained Lipid Membrane Model

We present a version of the coarse-grained Cooke lipid model, modified to simulate asymmetric lipid membranes. It is inspired by a method employed by Wang et al. [ Commun. Comput. Phys. 2013, 13, 1093-1106] for artificially penalizing lipid flip-flop but copes more robustly with differential stress, at the cost of one additional bead per lipid and the concomitant increase in computational overhead. Bilayer asymmetry ultimately breaks down beyond a system size dependent critical differential stress, which can be predicted from a simple analytical model. We remeasure many important material parameters for the new model and find it to be consistent with typical fluid lipid membranes. Maintaining a stable stress asymmetry has many applications, and we give two examples: (i) connecting monolayer stress to lipid number asymmetry in order to directly measure the monolayer area modulus and (ii) finding its strain-dependent higher-order correction by monitoring the equilibrium bilayer area.

Membrane pore energetics and the pathways to membrane rupture.. [DOI: 10.1101/2020.06.29.178988] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Biological membranes owe their strength and low permeability to the phospholipid bilayers at their core. Membrane strength is determined by the energetics and dynamics of membrane pores, whose tension-dependent nucleation and growth leads to rupture. Creation of nanoscale membrane pores is central to exocytosis, trafficking and other processes fundamental to life that require breaching of secure plasma or organelle membranes, and is the basis for biotechnologies using drug delivery, delivery of genetic material for gene editing and antimicrobial peptides. A prevailing view from seminal electroporation and membrane rupture studies is that pore growth and bilayer rupture are controlled by macroscopically long-lived metastable defect states that precede fully developed pores. It was argued that defect nucleation becomes rate-limiting at high tensions, explaining the exponential tension-dependence of rupture times [E. Evans et al., Biophys. J. 85, 2342-2350 (2003)]. Here we measured membrane pore free energies and bilayer rupture using highly coarse-grained simulations that probe very long time scales. We find no evidence of metastable pore states. At lower tensions, small hydrophobic pores mature into large hydrophilic pores on the pathway to rupture, with classical tension dependence of rupture times. Above a critical tension membranes rupture directly from a small hydrophobic pore, and rupture times depend exponentially on tension. Thus, we recover the experimentally reported regimes, but the origin of the high tension exponential regime is unrelated to macroscopically long-lived pre-pore defects. It arises because hydrophilic pores cannot exist above a critical tension, leading to radically altered pore dynamics and rupture kinetics.

Saddle-curvature instability of lipid bilayer induced by amphipathic peptides: a molecular model

Amphipathic peptides that partition into lipid bilayers affect the curvature elastic properties of their host. Some of these peptides are able to shift the Gaussian modulus to positive values, thus triggering an instability with respect to the formation of saddle curvatures. To characterize the generic aspects of the underlying mechanism, we employ a molecular lipid model that accounts for the interfacial tension between the polar and apolar regions of the membrane, for interactions between the lipid headgroups, and for the energy to stretch or compress the hydrocarbon chains. Peptides are modeled as cylinders that partition into the host membrane in a parallel orientation where they diminish the space available to the lipid headgroups and chains. The penetration depth into the membrane is determined by the angular size of the peptide's hydrophilic region. We demonstrate that only peptides with a small angular size of their hydrophilic region have an intrinsic tendency to render the Gaussian modulus more positive, and we identify conditions at which the Gaussian modulus adopts a positive sign upon increasing the peptide concentration. Our model allows us to also incorporate electrostatic interactions between cationic peptides and anionic lipids on the level of the linear Debye-Hückel model. We show that electrostatic interactions tend to shift the Gaussian modulus toward more positive values. Steric and electrostatic lipid-peptide interactions jointly decrease the effective interaction strength in the headgroup region of the host membrane thus suggesting a generic mechanisms of how certain amphipathic peptides are able to induce the formation of saddle curvatures.

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.

How Melittin Inserts into Cell Membrane: Conformational Changes, Inter-Peptide Cooperation, and Disturbance on the Membrane

Antimicrobial peptides (AMPs), as a key component of the immune defense systems of organisms, are a promising solution to the serious threat of drug-resistant bacteria to public health. As one of the most representative and extensively studied AMPs, melittin has exceptional broad-spectrum activities against microorganisms, including both Gram-positive and Gram-negative bacteria. Unfortunately, the action mechanism of melittin with bacterial membranes, especially the underlying physics of peptide-induced membrane poration behaviors, is still poorly understood, which hampers efforts to develop melittin-based drugs or agents for clinical applications. In this mini-review, we focus on recent advances with respect to the membrane insertion behavior of melittin mostly from a computational aspect. Membrane insertion is a prerequisite and key step for forming transmembrane pores and bacterial killing by melittin, whose occurrence is based on overcoming a high free-energy barrier during the transition of melittin molecules from a membrane surface-binding state to a transmembrane-inserting state. Here, intriguing simulation results on such transition are highlighted from both kinetic and thermodynamic aspects. The conformational changes and inter-peptide cooperation of melittin molecules, as well as melittin-induced disturbances to membrane structure, such as deformation and lipid extraction, are regarded as key factors influencing the insertion of peptides into membranes. The associated intermediate states in peptide conformations, lipid arrangements, membrane structure, and mechanical properties during this process are specifically discussed. Finally, potential strategies for enhancing the poration ability and improving the antimicrobial performance of AMPs are included as well.

Extracting lipid vesicles from plasma membranes via self-assembly of clathrin-inspired scaffolding nanoparticles

Single-cell analysis is a new and rapidly expanding field, the goal of which is obtaining fresh information from individual cells to understand the regulatory mechanisms of cell development and diseases. Conventional approaches generally rely on the cell lysis which, however, is destructive to cells and against multiple sampling from the living cell. Here, we propose and design a scaffolding nanoparticle (NP) system that enables us to sample cytoplasmic contents without rupturing the cellular membrane, by mimicking the unusual features of clathrin. Our simulation results reveal the design principles, following which scaffolding NPs can extract lipid vesicles from plasma membranes, with both the pathway and the mechanism resembling the clathrin-mediated endocytosis, i.e. multiple NPs deposit at the membrane, assembling into cage-like structures to deform the membrane into a vesicle shape. As important design parameters, the interaction between different NPs should be properly stronger than that between each NP and the membrane to ensure the cage formation, and optimal NP concentration and the membrane surface tension are also requisite for extracting lipid vesicles. Our results provide useful guidelines for design of bio-inspired scaffolding NPs as an intelligent machine for practical use in but not limited to the single-cell analysis.

Penetration of antimicrobial peptides in a lung surfactant model

Molecular dynamics simulations were successfully performed to understand the absorption mechanism of antimicrobial peptides LL-37, CATH-2, and SMAP-29 in a lung surfactant model. The antimicrobial peptides quickly penetrate in the lung surfactant model in dozens or hundreds nanoseconds, but they electrostatically interact with the lipid polar heads during the simulation time of 2 μs. This electrostatic interaction should be the explanation for the inactivation of the antimicrobial peptides when co-administrated with lung surfactant. As they strongly interact with the lipid polar heads of the lung surfactant, there is no positive charge available on the antimicrobial peptide to attack the negatively charged bacteria membrane. In order to avoid the interaction of peptides with the lipid polar heads, sodium cholate was used to form nanoparticles which act as an absorption enhancer of all antimicrobial peptides used in this investigation. The nanoparticles of 150 molecules of sodium cholate with one peptide were inserted on the top of the lung surfactant model. The nanoparticles penetrated into the lung surfactant model, spreading the sodium cholate molecules around the lipid polar heads. The sodium cholate molecules seem to protect the peptides from the interaction with the lipid polar heads, leaving them free to be delivered to the water phase. The penetration of peptides alone or even the peptide nanoparticles with sodium cholate do not collapse the lung surfactant model, indicating to be a promisor drug delivery system to the lung. The implications of this finding are that antimicrobial peptides may only be co-administered with an absorption enhancer such as sodium cholate into lung surfactant in order to avoid inactivation of their antimicrobial activity.

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'.

Probing molecular forces in multi-component physiological membranes

Biological membranes are remarkably heterogeneous, composed of diverse lipid mixtures with distinct chemical structure and composition. By combining molecular dynamics simulations and the newly developed Lipid-Force Distribution Analysis (L-FDA), we explore force transmission in complex multi-component membrane models mimicking eukaryotic organelles. We found that the chemical-moiety based segmentation at membrane interfaces revealed a distinctive distribution of bonded and non-bonded forces in diverse membrane environment. Our molecular stress analysis could have far-reaching implications in describing the relationship between membrane mechanical properties and functional states of chemically distinct lipids.

Engineering monolayer poration for rapid exfoliation of microbial membranes

A novel mechanism of monolayer poration leading to the rapid exfoliation and lysis of microbial membranes is reported.

The spread of bacterial resistance to traditional antibiotics continues to stimulate the search for alternative antimicrobial strategies. All forms of life, from bacteria to humans, are postulated to rely on a fundamental host defense mechanism, which exploits the formation of open pores in microbial phospholipid bilayers. Here we predict that transmembrane poration is not necessary for antimicrobial activity and reveal a distinct poration mechanism that targets the outer leaflet of phospholipid bilayers. Using a combination of molecular-scale and real-time imaging, spectroscopy and spectrometry approaches, we introduce a structural motif with a universal insertion mode in reconstituted membranes and live bacteria. We demonstrate that this motif rapidly assembles into monolayer pits that coalesce during progressive membrane exfoliation, leading to bacterial cell death within minutes. The findings offer a new physical basis for designing effective antibiotics.

Lipid somersaults: Uncovering the mechanisms of protein-mediated lipid flipping

Membrane lipids diffuse rapidly in the plane of the membrane but their ability to flip spontaneously across a membrane bilayer is hampered by a significant energy barrier. Thus spontaneous flip-flop of polar lipids across membranes is very slow, even though it must occur rapidly to support diverse aspects of cellular life. Here we discuss the mechanisms by which rapid flip-flop occurs, and what role lipid flipping plays in membrane homeostasis and cell growth. We focus on conceptual aspects, highlighting mechanistic insights from biochemical and in silico experiments, and the recent, ground-breaking identification of a number of lipid scramblases.

Dimeric peptides with three different linkers self-assemble with phospholipids to form peptide nanodiscs that stabilize membrane proteins

Three dimers of the amphipathic α-helical peptide 18A have been synthesized with different interhelical linkers inserted between the two copies of 18A. The dimeric peptides were denoted 'beltides' where Beltide-1 refers to the 18A-dimer without a linker, Beltide-2 is the 18A-dimer with proline (Pro) as a linker and Beltide-3 is the 18A-dimer linked by two glycines (Gly-Gly). The self-assembly of the beltides with the phospholipid DMPC was studied with and without the incorporated membrane protein bacteriorhodopsin (bR) through a combination of coarse-grained MD simulations, size-exclusion chromatography (SEC), circular dichroism (CD) spectroscopy, small-angle scattering (SAS), static light scattering (SLS) and UV-Vis spectroscopy. For all three beltides, MD and combined small-angle X-ray and -neutron scattering were consistent with a disc structure composed by a phospholipid bilayer surrounded by a belt of peptides and with a total disc diameter of approximately 10 nm. CD confirmed that all three beltides were α-helical in the free form and with DMPC. However, as shown by SEC the different interhelical linkers clearly led to different properties of the beltides. Beltide-3, with the Gly-Gly linker, was very adaptable such that peptide nanodiscs could be formed for a broad range of different peptide to lipid stoichiometries and therefore also possible disc-sizes. On the other hand, both Beltide-2 with the Pro linker and Beltide-1 without a linker were less adaptable and would only form discs of certain peptide to lipid stoichiometries. SLS revealed that the structural stability of the formed peptide nanodiscs was also highly affected by the linkers and it was found that Beltide-1 gave more stable discs than the other two beltides. With respect to membrane protein stabilization, each of the three beltides in combination with DMPC stabilizes the seven-helix transmembrane protein bacteriorhodopsin significantly better than the detergent octyl glucoside, but no significant difference was observed between the three beltides. We conclude that adaptability, size, and structural stability can be tuned by changing the interhelical linker while maintaining the properties of the discs with respect to membrane protein stabilization.

Neutron Reflectivity as a Tool for Physics-Based Studies of Model Bacterial Membranes

The principles of neutron reflectivity and its application as a tool to provide structural information at the (sub-) molecular unit length scale from models for bacterial membranes are described. The model membranes can take the form of a monolayer for a single leaflet spread at the air/water interface, or bilayers of increasing complexity at the solid/liquid interface. Solid-supported bilayers constrain the bilayer to 2D but can be used to characterize interactions with antimicrobial peptides and benchmark high throughput lab-based techniques. Floating bilayers allow for membrane fluctuations, making the phase behaviour more representative of native membranes. Bilayers of varying levels of compositional accuracy can now be constructed, facilitating studies with aims that range from characterizing the fundamental physical interactions, through to the characterization of accurate mimetics for the inner and outer membranes of Gram-negative bacteria. Studies of the interactions of antimicrobial peptides with monolayer and bilayer models for the inner and outer membranes have revealed information about the molecular control of the outer membrane permeability, and the mode of interaction of antimicrobials with both inner and outer membranes.

Design Principles for Nanoparticles Enveloped by a Polymer-Tethered Lipid Membrane

We propose the design for a nanoparticle carrier that combines three existing motifs into a single construct: a liposome is stabilized by anchoring it to an enclosed solid core via extended polymeric tethers that are chemically grafted to the core and physisorb into the surrounding lipid membrane. Such a design would exhibit several enticing properties, among them: (i) the anchoring stabilizes the liposome against a variety of external stresses, while preserving an aqueous compartment between core and membrane; (ii) the interplay of design parameters such as polymer length or grafting density enforces strong constraints on nanoparticle size and hence ensures a high degree of uniformity; and (iii) the physical and chemical characteristics of the individual constituents equip the construct with numerous functionalities that can be exploited in many ways. However, navigating the large parameter space requires a sound prior understanding for how various design features work together, and how this impacts potential pathways for synthesizing and assembling these nanoparticles. In this paper, we examine these connections in detail, using both soft matter theory and computer simulations at all levels of resolution. We thereby derive strong constraints on the experimentally relevant parameter space, and also propose potential equilibrium and nonequilibrium pathways for nanoparticle assembly.

Biomembranes in atomistic and coarse-grained simulations

The architecture of biological membranes is tightly coupled to the localization, organization, and function of membrane proteins. The organelle-specific distribution of lipids allows for the formation of functional microdomains (also called rafts) that facilitate the segregation and aggregation of membrane proteins and thus shape their function. Molecular dynamics simulations enable to directly access the formation, structure, and dynamics of membrane microdomains at the molecular scale and the specific interactions among lipids and proteins on timescales from picoseconds to microseconds. This review focuses on the latest developments of biomembrane force fields for both atomistic and coarse-grained molecular dynamics (MD) simulations, and the different levels of coarsening of biomolecular structures. It also briefly introduces scale-bridging methods applicable to biomembrane studies, and highlights selected recent applications.

Molecular understanding of a potential functional link between antimicrobial and amyloid peptides

Antimicrobial and amyloid peptides do not share common sequences, typical secondary structures, or normal biological activity but both the classes of peptides exhibit membrane-disruption ability to induce cell toxicity. Different membrane-disruption mechanisms have been proposed for antimicrobial and amyloid peptides, individually, some of which are not exclusive to either peptide type, implying that certain common principles may govern the folding and functions of different cytolytic peptides and associated membrane disruption mechanisms. Particularly, some antimicrobial and amyloid peptides have been identified to have dual complementary amyloid and antimicrobial properties, suggesting a potential functional link between amyloid and antimicrobial peptides. Given that some similar structural and membrane-disruption characteristics exist between the two classes of peptides, this review summarizes major findings, recent advances, and future challenges related to antimicrobial and amyloid peptides and strives to illustrate the similarities, differences, and relationships in the sequences, structures, and membrane interaction modes between amyloid and antimicrobial peptides, with a special focus on direct interactions of the peptides with the membranes. We hope that this review will stimulate further research at the interface of antimicrobial and amyloid peptides - which has been studied less intensively than either type of peptides - to decipher a possible link between both amyloid pathology and antimicrobial activity, which can guide drug design and peptide engineering to influence peptide-membrane interactions important in human health and diseases.

Theoretical insight into the relationship between the structures of antimicrobial peptides and their actions on bacterial membranes

Antimicrobial peptides with diverse cationic charges, amphiphathicities, and secondary structures possess a variety of antimicrobial activities against bacteria, fungi, and other generalized targets. To illustrate the relationship between the structures of these peptide and their actions at microscopic level, we present systematic coarse-grained dissipative particle dynamics simulations of eight types of antimicrobial peptides with different secondary structures interacting with a lipid bilayer membrane. We find that the peptides use multiple mechanisms to exert their membrane-disruptive activities: A cationic charge is essential for the peptides to selectively target negatively charged bacterial membranes. This cationic charge is also responsible for promoting electroporation. A significant hydrophobic portion is necessary to disrupt the membrane through formation of a permeable pore or translocation. Alternatively, the secondary structure and the corresponding rigidity of the peptides determine the pore structure and the translocation pathway.

A coarse-grained model for peptide aggregation on a membrane surface

The aggregation of peptides on a lipid bilayer is studied using coarse-grained molecular dynamics in implicit solvent. Peptides bind to and self-assemble on the membrane surface into β-rich fibrillar aggregates, even under conditions where only disordered oligomers form in bulk solution. Relative to a solid surface, the membrane surface facilitates peptide mobility and a more complex network of morphology transitions as aggregation proceeds. Additionally, final aggregate structures realized on the membrane surface are distinct from those observed on a comparable solid surface. The aggregated fibrils alter the local structure and material properties of the lipid bilayer in their immediate vicinity but have only a modest effect on the overall bending rigidity of the bilayer.

More than the sum of its parts: coarse-grained peptide-lipid interactions from a simple cross-parametrization

Interfacial systems are at the core of fascinating phenomena in many disciplines, such as biochemistry, soft-matter physics, and food science. However, the parametrization of accurate, reliable, and consistent coarse-grained (CG) models for systems at interfaces remains a challenging endeavor. In the present work, we explore to what extent two independently developed solvent-free CG models of peptides and lipids--of different mapping schemes, parametrization methods, target functions, and validation criteria--can be combined by only tuning the cross-interactions. Our results show that the cross-parametrization can reproduce a number of structural properties of membrane peptides (for example, tilt and hydrophobic mismatch), in agreement with existing peptide-lipid CG force fields. We find encouraging results for two challenging biophysical problems: (i) membrane pore formation mediated by the cooperative action of several antimicrobial peptides, and (ii) the insertion and folding of the helix-forming peptide WALP23 in the membrane.

Simulations suggest possible novel membrane pore structure

Amphiphilic proteins and peptides can induce the formation of stable and metastable pores in membranes. Using coarse-grained simulations, we have studied the factors that affect structure of peptide-stabilized pores. Our simulations are able to reproduce the formation of the well-known barrel-stave or toroidal pores, but in addition, we find evidence for a novel "double-belt" pore structure: in this structure the peptides that coat the membrane pore are oriented parallel to the membrane plane. To check the predictions of our coarse-grained model, we have performed more detailed simulations, using the MARTINI force field. These simulations show that the double-belt structure is stable up to at least the microsecond time scale.

Self-assembling peptides form nanodiscs that stabilize membrane proteins

New methods to handle membrane bound proteins, e.g. G-protein coupled receptors (GPCRs), are highly desirable. Recently, apoliprotein A1 (ApoA1) based lipoprotein particles have emerged as a new platform for studying membrane proteins, and it has been shown that they can self-assemble in combination with phospholipids to form discoidal shaped particles that can stabilize membrane proteins. In the present study, we have investigated an ApoA1 mimetic peptide with respect to its solution structure when in complex with phospholipids. This was achieved using a powerful combination of small-angle X-ray scattering (SAXS) and small-angle neutron scattering (SANS) supported by coarse-grained molecular dynamics simulations. The detailed structure of the discs was determined in unprecedented detail and it was found that they adopt a discoidal structure very similar to the ApoA1 based nanodiscs. We furthermore show that, like the ApoA1 and derived nanodiscs, these peptide discs can accommodate and stabilize a membrane protein. Finally, we exploit their dynamic properties and show that the 18A discs may be used for transferring membrane proteins and associated phospholipids directly and gently into phospholipid nanodiscs.

Engineering antimicrobial peptides with improved antimicrobial and hemolytic activities

The rapid rise of antibiotic resistance in pathogens becomes a serious and growing threat to medicine and public health. Naturally occurring antimicrobial peptides (AMPs) are an important line of defense in the immune system against invading bacteria and microbial infection. In this work, we present a combined computational and experimental study of the biological activity and membrane interaction of the computationally designed Bac2A-based peptide library. We used the MARTINI coarse-grained molecular dynamics with adaptive biasing force method and the umbrella sampling technique to investigate the translocation of a total of 91 peptides with different amino acid substitutions through a mixed anionic POPE/POPG (3:1) bilayer and a neutral POPC bilayer, which mimic the bacterial inner membrane and the human red blood cell (hRBC) membrane, respectively. Potential of mean force (PMF, free energy profile) was obtained to measure the free energy barrier required to transfer the peptides from the bulk water phase to the water-membrane interface and to the bilayer interior. Different PMF profiles can indeed identify different membrane insertion scenarios by mapping out peptide-lipid energy landscapes, which are correlated with antimicrobial activity and hemolytic activity. Computationally designed peptides were further tested experimentally for their antimicrobial and hemolytic activities using bacteria growth inhibition assay and hemolysis assay. Comparison of PMF data with cell assay results reveals a good correlation of the peptides between predictive transmembrane activity and antimicrobial/hemolytic activity. Moreover, the most active mutants with the balanced substitutions of positively charged Arg and hydrophobic Trp residues at specific positions were discovered to achieve the improved antimicrobial activity while minimizing red blood cell lysis. Such substitutions provide more effective and cooperative interactions to distinguish the peptide interaction with different lipid bilayers. This work provides a useful computational tool to better understand the mechanism and energetics of membrane insertion of AMPs and to rationally design more effective AMPs.

Advances in the understanding of nanomaterial–biomembrane interactions and their mathematical and numerical modeling

The widespread application of nanomaterials (NMs), which has accompanied advances in nanotechnology, has increased their chances of entering an organism, for example, via the respiratory system, skin absorption or intravenous injection. Although accumulating experimental evidence has indicated the important role of NM–biomembrane interaction in these processes, the underlying mechanisms remain unclear. Computational techniques, as an alternative to experimental efforts, are effective tools to simulate complicated biological behaviors. Computer simulations can investigate NM–biomembrane interactions at the nanoscale, providing fundamental insights into dynamic processes that are challenging to experimental observation. This paper reviews the current understanding of NM–biomembrane interactions, and existing mathematical and numerical modeling methods. We highlight the advantages and limitations of each method, and also discuss the future perspectives in this field. Better understanding of NM–biomembrane interactions can benefit various fields, including nanomedicine and diagnosis.

Modeling peptide binding to anionic membrane pores

Peptide-induced pore formation in membranes can be dissected into two steps: pore formation and peptide binding to the pore. A computational method is proposed to study the second step in anionic membranes. The electrostatic potential is obtained from numerical solutions to the Poisson-Boltzmann equation and is then used in conjunction with IMM1 (implicit membrane model 1). A double charge layer model is used to incorporate the effects of the membrane dipole potential. Inhomogeneity of the charge density in the pore, characterized by explicit membrane simulations of toroidal pores, is included in the model. This approach was applied to two extensively studied peptides, magainin and melittin. In agreement with previous work, binding to toroidal pores is more favorable than binding to the flat membrane. The dependence of binding energy on anionic content exhibits different patterns for the two peptides, in correlation with the different lipid selectivity that has been observed experimentally.

Membrane sculpting by F-BAR domains studied by molecular dynamics simulations

Interplay between cellular membranes and their peripheral proteins drives many processes in eukaryotic cells. Proteins of the Bin/Amphiphysin/Rvs (BAR) domain family, in particular, play a role in cellular morphogenesis, for example curving planar membranes into tubular membranes. However, it is still unclear how F-BAR domain proteins act on membranes. Electron microscopy revealed that, in vitro, F-BAR proteins form regular lattices on cylindrically deformed membrane surfaces. Using all-atom and coarse-grained (CG) molecular dynamics simulations, we show that such lattices, indeed, induce tubes of observed radii. A 250 ns all-atom simulation reveals that F-BAR domain curves membranes via the so-called scaffolding mechanism. Plasticity of the F-BAR domain permits conformational change in response to membrane interaction, via partial unwinding of the domains 3-helix bundle structure. A CG simulation covering more than 350 µs provides a dynamic picture of membrane tubulation by lattices of F-BAR domains. A series of CG simulations identified the optimal lattice type for membrane sculpting, which matches closely the lattices seen through cryo-electron microscopy.

To generate organelles, eukaryotic cells sculpt their membranes into compartments, often employing proteins as chaperones, for example, F-BAR domains. The latter induce formation of tubular and vesicular membranes. Functional and structural studies suggest that F-BAR domains sculpt membranes through electrostatic interactions, driving the membrane to match the concave surface of the protein's banana-like shape. Cryo-electron microscopy (cryo-EM) studies provide an average static picture of how F-BAR domains form lattices on the surface of membranes to induce tube formation. Complementing the cryo-EM images, molecular dynamics simulations reported here offer a detailed, dynamic picture of membrane tubulation by a lattice of F-BAR domains and identified lattice types optimally attuned to producing high membrane curvature. The simulations reproduced also a process lasting 350 µs in which lattices of F-BAR domains form a complete tube out of an initially flat membrane. The molecular dynamics study offers, thereby, both a large-scale picture of membrane sculpting by F-BAR domain lattices as well as atomic-level dynamic information about the involvement of the individual F-BAR domain and its interactions with partner F-BAR domains and membrane in the sculpting process.

Multiscale molecular dynamics simulations of membrane proteins

The time and length scales accessible by biomolecular simulations continue to increase. This is in part due to improvements in algorithms and computing performance, but is also the result of the emergence of coarse-grained (CG) potentials, which complement and extend the information obtainable from fully detailed models. CG methods have already proven successful for a range of applications that benefit from the ability to rapidly simulate spontaneous self-assembly within a lipid membrane environment, including the insertion and/or oligomerization of a range of "toy models," transmembrane peptides, and single- and multi-domain proteins. While these simplified approaches sacrifice atomistic level detail, it is now straightforward to "reverse map" from CG to atomistic descriptions, providing a strategy to assemble membrane proteins within a lipid environment, prior to all-atom simulation. Moreover, recent developments have been made in "dual resolution" techniques, allowing different molecules in the system to be modeled with atomistic or CG resolution simultaneously.

Line-tension controlled mechanism for influenza fusion

Our molecular simulations reveal that wild-type influenza fusion peptides are able to stabilize a highly fusogenic pre-fusion structure, i.e. a peptide bundle formed by four or more trans-membrane arranged fusion peptides. We rationalize that the lipid rim around such bundle has a non-vanishing rim energy (line-tension), which is essential to (i) stabilize the initial contact point between the fusing bilayers, i.e. the stalk, and (ii) drive its subsequent evolution. Such line-tension controlled fusion event does not proceed along the hypothesized standard stalk-hemifusion pathway. In modeled influenza fusion, single point mutations in the influenza fusion peptide either completely inhibit fusion (mutants G1V and W14A) or, intriguingly, specifically arrest fusion at a hemifusion state (mutant G1S). Our simulations demonstrate that, within a line-tension controlled fusion mechanism, these known point mutations either completely inhibit fusion by impairing the peptide's ability to stabilize the required peptide bundle (G1V and W14A) or stabilize a persistent bundle that leads to a kinetically trapped hemifusion state (G1S). In addition, our results further suggest that the recently discovered leaky fusion mutant G13A, which is known to facilitate a pronounced leakage of the target membrane prior to lipid mixing, reduces the membrane integrity by forming a 'super' bundle. Our simulations offer a new interpretation for a number of experimentally observed features of the fusion reaction mediated by the prototypical fusion protein, influenza hemagglutinin, and might bring new insights into mechanisms of other viral fusion reactions.

Determining the Gaussian curvature modulus of lipid membranes in simulations

The Gaussian curvature modulus κ¯ of lipid bilayers likely contributes more than 100 kcal/mol to every cellular fission or fusion event. This huge impact on membrane remodeling energetics might be a factor that codetermines the complex lipid composition of biomembranes through tuning of κ¯. Yet, its value has been measured only for a handful of simple lipids, and no simulation has so far determined it better than a factor of two, rendering a systematic investigation of such enticing speculations impossible. Here we propose a highly accurate method to determine κ¯ in computer simulations. It relies on the interplay between curvature stress and edge tension of partially curved axisymmetric membrane disks and requires determining their closing probability. For a simplified lipid model we obtain κ¯ and its relation to the normal bending modulus κ for membranes differing both in stiffness and spontaneous lipid curvature. The elastic ratio κ¯/κ can be determined with a few percent statistical accuracy. Its value agrees with the scarce experimental data, and its change with spontaneous lipid curvature is compatible with theoretical expectations, thereby granting additional information on monolayer properties. We also show that an alternative determination of these elastic parameters based on moments of the lateral stress profile gives markedly different and unphysical values.

Biomolecule surface patterning may enhance membrane association

Under dehydration conditions, amphipathic late embryogenesis abundant proteins fold spontaneously from a random conformation into α-helical structures, and this transition is promoted by the presence of membranes. To gain insight into the thermodynamics of membrane association, we model the resulting α-helical structures as infinite rigid cylinders patterned with hydrophobic and hydrophilic stripes oriented parallel to their axis. Statistical thermodynamic calculations using single chain mean field theory show that the relative thickness of the stripes controls the free energy of interaction of the α-helices with a phospholipid bilayer, as does the bilayer structure and the depth of the equilibrium penetration of the cylinders into the bilayer. The results may suggest the optimal thickness of the stripes to mimic the association of such protein with membranes.

Modeling the self-assembly of lipids and nanotubes in solution: forming vesicles and bicelles with transmembrane nanotube channels

Via dissipative particle dynamics (DPD), we simulate the self-assembly of end-functionalized, amphiphilic nanotubes and lipids in a hydrophilic solvent. Each nanotube encompasses a hydrophobic stalk and two hydrophilic ends, which are functionalized with end-tethered chains. With a relatively low number of the nanotubes in solution, the components self-assemble into stable lipid-nanotube vesicles. As the number of nanotubes is increased, the system exhibits a vesicle-to-bicelle transition, resulting in stable hybrid bicelle. Moreover, our results reveal that the nanotubes cluster into distinct tripod-like structures within the vesicles and aggregate into a ring-like assembly within the bicelles. For both the vesicles and bicelles, the nanotubes assume trans-membrane orientations, with the tethered hairs extending into the surrounding solution or the encapsulated fluid. Thus, the hairs provide a means of regulating the transport of species through the self-assembled structures. Our findings provide guidelines for creating nanotube clusters with distinctive morphologies that might be difficult to achieve through more conventional means. The results also yield design rules for creating synthetic cell-like objects or microreactors that can exhibit biomimetic functionality.

Spontaneous buckling of lipid bilayer and vesicle budding induced by antimicrobial peptide magainin 2: a coarse-grained simulation study

Molecular mechanisms of the action of antimicrobial peptides on bacterial membranes were studied by large scale coarse-grained simulations of magainin 2-dipalmitoylphosphatidylcholine/palmitoyloleoylphosphatidylglycerol (DPPC/POPG) mixed bilayer systems with spatial extents up to 0.1 μm containing up to 1600 peptides. Equilibrium simulations exhibit disordered toroidal pores stabilized by peptides. However, when a layer of peptides is placed near the lipid head groups on one side of the bilayer only, their incorporation leads to a spontaneous buckling of the bilayer. This buckling is followed by the formation of a quasi-spherical vesicular bud connected to the bilayer by a narrow neck. The mean curvature of the budding region is consistent with what is expected based on the dependence of the area per lipid on the peptide-to-lipid ratio in equilibrium simulations. Our simulations suggest that the incorporation of antimicrobial peptides on the exterior surface of a vesicle or a bacterial cell leads to buckling and vesicle budding, presumably accompanied by nucleations of giant transient pores of sizes that are much larger than indicated by equilibrium measurements and simulations.

Cell surface engineering with polyelectrolyte multilayer thin films

Layer-by-layer assembly of polyelectrolyte multilayer (PEM) films represents a bottom-up approach for re-engineering the molecular landscape of cell surfaces with spatially continuous and molecularly uniform ultrathin films. However, fabricating PEMs on viable cells has proven challenging owing to the high cytotoxicity of polycations. Here, we report the rational engineering of a new class of PEMs with modular biological functionality and tunable physicochemical properties which have been engineered to abrogate cytotoxicity. Specifically, we have discovered a subset of cationic copolymers that undergoes a conformational change, which mitigates membrane disruption and facilitates the deposition of PEMs on cell surfaces that are tailorable in composition, reactivity, thickness, and mechanical properties. Furthermore, we demonstrate the first successful in vivo application of PEM-engineered cells, which maintained viability and function upon transplantation and were used as carriers for in vivo delivery of PEMs containing biomolecular payloads. This new class of polymeric film and the design strategies developed herein establish an enabling technology for cell transplantation and other therapies based on engineered cells.

Forming transmembrane channels using end-functionalized nanotubes

Using dissipative particle dynamics (DPD) simulations, we examine the interaction between amphiphilic nanotubes and lipid bilayer membranes. The nanotubes are represented by a hydrophobic shaft that is end-functionalized with hydrophilic groups. Nanotubes that are capped by a monolayer of hydrophilic beads or also encompass hydrophilic "hairs" on just one end of the shaft are found to spontaneously penetrate and assume a transmembrane position; the process, however, depends critically on the membrane tension. On the other hand, nanotubes that include hydrophilic hairs at both ends of the hydrophobic shaft are not observed to spontaneously self-organize into the bilayer. When the membrane is stretched to form a pore, the nanotubes with two hairy ends adsorb on the edge of the pore and become localized in the membrane, thus forming a transmembrane channel. The findings from these studies provide guidelines for creating biomimetic nanotube channels that are capable of selectively transporting molecules through the membrane in response to changes in the local environment.

A systematically coarse-grained solvent-free model for quantitative phospholipid bilayer simulations

We present an implicit solvent coarse-grained (CG) model for quantitative simulations of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers. The absence of explicit solvent enables membrane simulations on large length and time scales at moderate computational expense. Despite improved computational efficiency, the model preserves chemical specificity and quantitative accuracy. The bonded and nonbonded interactions together with the effective cohesion mimicking the hydrophobic effect were systematically tuned by matching structural and mechanical properties from experiments and all-atom bilayer simulations, such as saturated area per lipid, radial distribution functions, density and pressure profiles across the bilayer, P(2) order, etc. The CG lipid model is shown to self-assemble into a bilayer starting from a random dispersion. Its line tension and elastic properties, such as bending and stretching modulus, are semiquantitatively consistent with experiments. The effects of (i) reduced molecular friction and (ii) more efficient integration combine to an overall speed-up of 3-4 orders of magnitude compared to all-atom bilayer simulations. Our CG lipid model is especially useful for studies of large-scale phenomena in membranes that nevertheless require a fair description of chemical specificity, e.g., membrane patches interacting with movable and transformable membrane proteins and peptides.

Membrane poration by antimicrobial peptides combining atomistic and coarse-grained descriptions

Antimicrobial peptides (AMPs) comprise a large family of peptides that include small cationic peptides, such as magainins, which permeabilize lipid membranes. Previous atomistic level simulations of magainin-H2 peptides show that they act by forming toroidal transmembrane pores. However, due to the atomistic level of description, these simulations were necessarily limited to small system sizes and sub-microsecond time scales. Here, we study the long-time relaxation properties of these pores by evolving the systems using a coarse-grain (CG) description. The disordered nature and the topology of the atomistic pores are maintained at the CG level. The peptides sample different orientations but at any given time, only a few peptides insert into the pore. Key states observed at the CG level are subsequently back-transformed to the atomistic level using a resolution-transformation protocol. The configurations sampled at the CG level are stable in the atomistic simulation. The effect of helicity on pore stability is investigated at the CG level and we find that partial helicity is required to form stable pores. We also show that the current CG scheme can be used to study spontaneous poration by magainin-H2 peptides. Overall, our simulations provide a multi-scale view of a fundamental biophysical membrane process involving a complex interplay between peptides and lipids.

Cause and effect of melittin-induced pore formation: a computational approach

Melittin embedded in a palmitoyl oleyl phosphatidylcholine bilayer at a high peptide/lipid ratio (1:30) was simulated in the presence of explicit water and ions. The simulation results indicate the incipience of an ion-permeable water pore through collective membrane perturbation by bound peptides. The positively charged residues of melittin not only act as "anchors" but also disrupt the membrane, leading to cell lysis. A detailed analysis of the lipid tail order parameter profile depicts localized membrane perturbation. The lipids in the vicinity of the aqueous cavity adopt a tilted conformation, which allows local bilayer thinning. The prepore thus formed can be considered as the melittin-induced structural defects in the bilayer membrane. Because of the strong cationic nature, the melittin-induced prepore exhibits selectivity toward anions over cations. As Cl(-) ions entered into the prepore, they are electrostatically entrapped by positively charged residues located at its wall. The confined motion of the Cl(-) ions in the membrane interior is obvious from calculated diffusion coefficients. Moreover, reorientation of the local lipids occurs in such a way that few lipid heads along with peptide helices can line the surface of the penetrating aqueous phase. The flipping of lipids argued in favor of melittin-induced toroidal pore over a barrel-stave mechanism. Thus, our result provides atomistic level details of the mechanism of membrane disruption by antimicrobial peptide melittin.

Molecular interactions between magainin 2 and model membranes in situ

In this paper, we investigated the molecular interactions of magainin 2 with model cell membranes using sum frequency generation (SFG) vibrational spectroscopy and attenuated total reflectance-Fourier transform infrared spectroscopy (ATR-FTIR). Symmetric 1-palmitoyl-2-oleoyl-sn-glycero-3-[Phospho-rac-(1-glycerol)] (POPG) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers, which model the bacterial and mammalian cell membranes, respectively, were used in the studies. It was observed by SFG that magainin 2 orients relatively parallel to the POPG lipid bilayer surface at low solution concentrations, around 200 nM. When increasing the magainin 2 concentration to 800 nM, both SFG and ATR-FTIR results indicate that magainin 2 molecules insert into the POPG bilayer and adopt a transmembrane orientation with an angle of about 20 degrees from the POPG bilayer normal. For the POPC bilayer, even at a much higher peptide concentration of 2.0 microM, no ATR-FTIR signal was detected. For this concentration on POPC, SFG studies indicated that magainin 2 molecules adopt an orientation nearly parallel to the bilayer surface, with an orientation angle of about 75 degrees from the surface normal. This shows that SFG has a much better detection limit than ATR-FTIR and can therefore be applied to study interfacial molecules with a much lower surface coverage. This magainin 2 orientation study and further investigation of the lipid bilayer SFG signals support the proposed toroidal pore model for the antimicrobial activity of magainin 2.

Lipids on the move: simulations of membrane pores, domains, stalks and curves

In this review we describe the state-of-the-art of computer simulation studies of lipid membranes. We focus on collective lipid-lipid and lipid-protein interactions that trigger deformations of the natural lamellar membrane state, showing that many important biological processes including self-aggregation of membrane components into domains, the formation of non-lamellar phases, and membrane poration and curving, are now amenable to detailed simulation studies.