All-atom simulations and free-energy calculations of coiled-coil peptides with lipid bilayers: binding strength, structural transition, and effect on lipid dynamics

Coarse-Grained Molecular Dynamics of pH-Sensitive Lipids

pH-sensitive lipids represent a class of lipids that can be protonated and destabilized in acidic environments, as they become positively charged in response to low-pH conditions. They can be incorporated into lipidic nanoparticles such as liposomes, which are able to change their properties and allow specific drug delivery at the acidic conditions encountered in some pathological microenvironments. In this work, we used coarse-grained molecular-dynamic simulations to study the stability of neutral and charged lipid bilayers containing POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and various kinds of ISUCA ((F)2-(imidazol-1-yl)succinic acid)-derived lipids, which can act as pH-sensitive molecules. In order to explore such systems, we used a MARTINI-derived forcefield, previously parameterized using all-atom simulation results. We calculated the average area per lipid, the second-rank order parameter and the lipid diffusion coefficient of both lipid bilayers made of pure components and mixtures of lipids in different proportions, under neutral or acidic conditions. The results show that the use of ISUCA-derived lipids disturbs the lipid bilayer structure, with the effect being particularly marked under acidic conditions. Although more-in depth studies on these systems must be carried out, these initial results are encouraging and the lipids designed in this research could be a good basis for developing new pH-sensitive liposomes.

Molecular Dynamics Simulations of Human Beta-Defensin Type 3 Crossing Different Lipid Bilayers

Human β defensin type 3 (hBD-3) is a small cationic cysteine-rich peptide. It has a broad spectrum of antimicrobial activities. However, at high concentrations, it also shows hemolytic activity by interrupting red blood cells. To understand the selectivity of hBD-3 disrupting cell membranes, investigating the capability of hBD-3 translocating through different membranes is important. Since hBD-3 in the analogue form in which all three pairs of disulfide bonds are broken has similar antibacterial activities to the wild-type, this project investigates the structure and dynamics of an hBD-3 analogue in monomer, dimer, and tetramer forms through both zwitterionic and negatively charged lipid bilayers using molecular dynamics (MD) simulations. One tetramer structure of hBD-3 was predicted by running all-atom MD simulations on hBD-3 in water at a high concentration, which was found to be stable in water during 400 ns all-atom simulations based on root-mean-squared deviation, root-mean-squared fluctuation, buried surface area, and binding interaction energy calculations. After that, hBD-3 in different forms was placed inside different membranes, and then steered MD simulation was conducted to pull the hBD-3 out of the membrane along the z-direction to generate different configurational windows to set up umbrella-sampling (US) simulations. Because extensive sampling is important to obtain accurate free energy barriers, coarse-grained US MD simulations were performed in each window. Based on the long-term simulation result, membrane thinning was found near hBD-3 in different lipid bilayers and in different hBD-3 oligomer systems. By calculating the root-mean-squared deviation of the z-coordinate of hBD-3 molecules, rotation of the oligomer inside the bilayer and stretching of the oligomer structure along the z-direction were observed. Although reorientation of lipid heads toward the hBD-3 tetramer was observed based on the density profile calculation, the order parameter calculation shows that hBD-3 disrupts 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-l-serine (POPS) lipids more significantly and makes it less ordered than on 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipids. Calculating the free energy of hBD-3 through different lipid bilayers, it was found that generally hBD-3 encounters a lower energy barrier through negatively charged lipid membranes than the zwitterionic membrane. hBD-3 in different forms needs to overcome a lower energy barrier crossing the combined POPC+POPS bilayer through the POPS leaflet than through the POPC leaflet. Besides that, the potential of mean force result suggests that hBD-3 forms an oligomer translocating negatively charged lipid membranes at a low concentration. This study supplied new insight into the antibacterial mechanism of hBD-3 through different membranes.

Membrane Binding and Homodimerization of Atg16 Via Two Distinct Protein Regions is Essential for Autophagy in Yeast

Macroautophagy is a bulk degradation mechanism in eukaryotic cells. Efficiency of an essential step of this process in yeast, Atg8 lipidation, relies on the presence of Atg16, a subunit of the Atg12-Atg5-Atg16 complex acting as the E3-like enzyme in the ubiquitination-like reaction. A current view on the functional structure of Atg16 in the yeast S. cerevisiae comes from the two crystal structures that reveal the Atg5-interacting α-helix linked via a flexible linker to another α-helix of Atg16, which then assembles into a homodimer. This view does not explain the results of previous in vitro studies revealing Atg16-dependent deformations of membranes and liposome-binding of the Atg12-Atg5 conjugate upon addition of Atg16. Here we show that Atg16 acts as both a homodimerizing and peripheral membrane-binding polypeptide. These two characteristics are imposed by the two distinct regions that are disordered in the nascent protein. Atg16 binds to membranes in vivo via the amphipathic α-helix (amino acid residues 113-131) that has a coiled-coil-like propensity and a strong hydrophobic face for insertion into the membrane. The other protein region (residues 64-99) possesses a coiled-coil propensity, but not amphipathicity, and is dispensable for membrane anchoring of Atg16. This region acts as a Leu-zipper essential for formation of the Atg16 homodimer. Mutagenic disruption in either of these two distinct domains renders Atg16 proteins that, in contrast to wild type, completely fail to rescue the autophagy-defective phenotype of atg16Δ cells. Together, the results of this study yield a model for the molecular mechanism of Atg16 function in macroautophagy.

Design, Screening, and Testing of Non-Rational Peptide Libraries with Antimicrobial Activity: In Silico and Experimental Approaches

One of the challenges of modern biotechnology is to find new routes to mitigate the resistance to conventional antibiotics. Antimicrobial peptides (AMPs) are an alternative type of biomolecules, naturally present in a wide variety of organisms, with the capacity to overcome the current microorganism resistance threat. Here, we reviewed our recent efforts to develop a new library of non-rationally produced AMPs that relies on bacterial genome inherent diversity and compared it with rationally designed libraries. Our approach is based on a four-stage workflow process that incorporates the interplay of recent developments in four major emerging technologies: artificial intelligence, molecular dynamics, surface-display in microorganisms, and microfluidics. Implementing this framework is challenging because to obtain reliable results, the in silico algorithms to search for candidate AMPs need to overcome issues of the state-of-the-art approaches that limit the possibilities for multi-space data distribution analyses in extremely large databases. We expect to tackle this challenge by using a recently developed classification algorithm based on deep learning models that rely on convolutional layers and gated recurrent units. This will be complemented by carefully tailored molecular dynamics simulations to elucidate specific interactions with lipid bilayers. Candidate AMPs will be recombinantly-expressed on the surface of microorganisms for further screening via different droplet-based microfluidic-based strategies to identify AMPs with the desired lytic abilities. We believe that the proposed approach opens opportunities for searching and screening bioactive peptides for other applications.

Liposome Drug Delivery System across Endothelial Plasma Membrane: Role of Distance between Endothelial Cells and Blood Flow Rate

This paper discusses specific features of the interactions of small-diameter liposomes with the cytoplasmic membrane of endothelial cells using in silico methods. The movement pattern of the liposomal drug delivery system was modeled in accordance with the conditions of the near-wall layer of blood flow. Our simulation results show that the liposomes can become stuck in the intercellular gaps and even break down when the gap is reduced. Liposomes stuck in the gaps are capable of withstanding a shell deformation of ~15% with an increase in liposome energy by 26%. Critical deformation of the membrane gives an impetus to drug release from the liposome outward. We found that the liposomes moving in the near-wall layer of blood flow inevitably stick to the membrane. Liposome sticking on the membrane is accompanied by its gradual splicing with the membrane bilayer. This leads to a gradual drug release inside the cell.

Key Factors Determining Efficiency of Liquid-Liquid Extraction: Implications from Molecular Dynamics Simulations of Biphasic Behaviors of CyMe4-BTPhen and Its Am(III) Complexes

CyMe₄-BTPhen (2,9-bis(5,5,8,8-tetramethyl-5,6,7,8-tetrahydro-1,2,4-benzotriazin-3-yl)-1,10-phenanthroline, denoted as L) has been considered as a promising extractant in lanthanide(III)/actinide(III) separation. Vast endeavors in its application put forward a compelling need on the understanding of the underlying mechanism in the liquid-liquid extraction. To address the issue of its dynamics in biphasic systems, we carried out molecular dynamics (MD) simulations of L and its complexes with a heavy f-block metal ion, americium(III) (Am³⁺) in "oil"/water binary solvents. Two types of organic phases have been considered, differing in the presence of octanol in the bulk n-dodecane or not, and the distribution of the solutes and their interfacial behaviors have been investigated. Two of the key factors that determine the efficiency of a liquid-liquid extraction protocol were delineated and discussed, that is, the appropriate ligand to enhance the lipophilicity of AmL complexes and appropriate way to form ion pairs to minimize the attraction between the complexes and aqueous phase. The simulations showed that the charge states of both ligand and AmL complexes were strongly correlated with their phase behavior, and the migration of neutral species was driven by van der Waals interactions while that of charged species by electrostatic interactions, indicating stronger lipophilicity of the former than the latter. The presence of octanol facilitated the migration of the ligand from the interface to the organic phase via hydrogen bond between its polar head and the ligand or the AmL complexes and constituted a polar core in the organic phase. This work bridged the widely used liquid-liquid extraction technique in chemistry to a fundamental chemical concept, that is, minimization of hydrophilicity and maximization of lipophilicity to facilitate phase transfer from the aqueous phase to the organic phase, and is expected to improve the understanding of dynamics of ligands and their complexes with metal ions and to contribute to the development of efficient protocols for phase transfer of target species.

Interaction of SNARE Mimetic Peptides with Lipid bilayers: Effects of Secondary Structure, Bilayer Composition and Lipid Anchoring

The coiled-coil forming peptides 'K' enriched in lysine and 'E' enriched in glutamic acid have been used as a minimal SNARE mimetic system for membrane fusion. Here we describe atomistic molecular dynamics simulations to characterize the interactions of these peptides with lipid bilayers for two different compositions. For neutral phosphatidylcholine (PC)/phosphatidylethanolamine (PE) bilayers the peptides experience a strong repulsive barrier against adsorption, also observed in potential of mean force (PMF) profiles calculated with umbrella sampling. For peptide K, a minimum of -12 k_BT in the PMF provides an upper bound for the binding free energy whereas no stable membrane bound state could be observed for peptide E. In contrast, the electrostatic interactions with negatively charged phosphatidylglycerol (PG) lipids lead to fast adsorption of both peptides at the head-water interface. Experimental data using fluorescently labeled peptides confirm the stronger binding to PG containing bilayers. Lipid anchors have little effect on the peptide-bilayer interactions or peptide structure, when the peptide also binds to the bilayer in the absence of a lipid anchor. For peptide E, which does not bind to the PC bilayer without a lipid anchor, the presence of such an anchor strengthens the electrostatic interactions between the charged side chains and the zwitterionic head-groups and leads to a stabilization of the peptide's helical fold by the membrane.

Intrinsic lipid binding activity of ATG16L1 supports efficient membrane anchoring and autophagy

Membrane targeting of autophagy-related complexes is an important step that regulates their activities and prevents their aberrant engagement on non-autophagic membranes. ATG16L1 is a core autophagy protein implicated at distinct phases of autophagosome biogenesis. In this study, we dissected the recruitment of ATG16L1 to the pre-autophagosomal structure (PAS) and showed that it requires sequences within its coiled-coil domain (CCD) dispensable for homodimerisation. Structural and mutational analyses identified conserved residues within the CCD of ATG16L1 that mediate direct binding to phosphoinositides, including phosphatidylinositol 3-phosphate (PI3P). Mutating putative lipid binding residues abrogated the localisation of ATG16L1 to the PAS and inhibited LC3 lipidation. On the other hand, enhancing lipid binding of ATG16L1 by mutating negatively charged residues adjacent to the lipid binding motif also resulted in autophagy inhibition, suggesting that regulated recruitment of ATG16L1 to the PAS is required for its autophagic activity. Overall, our findings indicate that ATG16L1 harbours an intrinsic ability to bind lipids that plays an essential role during LC3 lipidation and autophagosome maturation.

Molecular Insights into the Unusual Structure of an Antifreeze Protein with a Hydrated Core

The primary driving force for protein folding is the formation of a well-packed, anhydrous core. However, recently, the crystal structure of an antifreeze protein, maxi, has been resolved where the core of the protein is filled with water, which apparently contradicts the existing notion of protein folding. Here, we have performed standard molecular dynamics (MD) simulation, replica exchange MD (REMD) simulation, and umbrella sampling using TIP4P water at various temperatures (300, 260, and 240 K) to explore the origin of this unusual structural feature. It is evident from standard MD and REMD simulations that the protein is found to be stable at 240 K in its unusual state. The core of protein has two layers of semi-clathrate water separating the methyl groups of alanine residues from different helical strands. However, with increasing temperature (260 and 300 K), the stability decreases as the core becomes dehydrated, and methyl groups of alanine are tightly packed driven by hydrophobic interactions. Calculation of the potential of mean force by an umbrella sampling technique between a pair of model hydrophobes resembling maxi protein at 240 K shows the stabilization of second solvent-separated minima (SSM), which provides a thermodynamic rationale of the unusual structural feature in terms of weakening of the hydrophobic interaction. Because the stabilization of SSMs is implicated for cold denaturation, it suggests that the maxi protein is so designed by nature where the cold denatured-like state becomes the biologically active form as it works near or below the freezing point of water.

Distinct roles of SNARE-mimicking lipopeptides during initial steps of membrane fusion

A model system for membrane fusion, inspired by SNARE proteins and based on two complementary lipopeptides CPnE4 and CPnK4, has been recently developed. It consists of cholesterol (C), a poly(ethylene glycol) linker (Pn) and either a cationic peptide K4 (KIAALKE)4 or an anionic peptide E4 (EIAALEK)4. In this paper, fluorescence spectroscopy is used to decipher distinct but complementary roles of these lipopeptides during early stages of membrane fusion. Molecular evidence is provided that different distances of E4 in CPnE4 and K4 in CPnK4 from the bilayer represent an important mechanism, which enables fusion. Whereas E4 is exposed to the bulk and solely promotes membrane binding of CPnK4, K4 loops back to the lipid-water interface where it fulfills two distinct roles: it initiates bilayer contact by binding to CPnE4 containing bilayers; and it initiates fusion by modulating the bilayer properties. The interaction between CPnE4 and CPnK4 is severely down-regulated by binding of K4 to the bilayer and possible only if the lipopeptides approach each other as constituents of different bilayers. When the complementary lipopeptides are localized in the same bilayer, hetero-coiling is disabled. These data provide crucial insights as to how fusion is initiated and highlight the importance of both peptides in this process.

Membrane-Fusogen Distance Is Critical for Efficient Coiled-Coil-Peptide-Mediated Liposome Fusion

We have developed a model system for membrane fusion that utilizes lipidated derivatives of a heterodimeric coiled-coil pair dubbed E3 (EIAALEK)3 and K3 (KIAALKE)3. In this system, peptides are conjugated to a lipid anchor via a poly(ethylene glycol) (PEG) spacer, and this contribution studies the influence of the PEG spacer length, coupled with the type of lipid anchor, on liposome-liposome fusion. The effects of these modifications on peptide secondary structure, their interactions with liposomes, and their ability to mediate fusion were studied using a variety of different content mixing experiments and CD spectroscopy. Our results demonstrate the asymmetric role of the peptides in the fusion process because alterations to the PEG spacer length affect E3 and K3 differently. We conclude that negatively charged E3 acts as a "handle" for positively charged K3 and facilitates liposome docking, the first stage of the fusion process, through coiled-coil formation. The efficacy of this E3 handle is enhanced by longer spacer lengths. K3 directs the fusion process via peptide-membrane interactions, but the length of the PEG spacer plays two competing roles: a PEG4/PEG8 spacer length is optimal for membrane destabilization; however, a PEG12 spacer increases the fusion efficiency over time by improving the peptide accessibility for successive fusion events. Both the anchor type and spacer length affect the peptide structure; a cholesterol anchor appears to enhance K3-membrane interactions and thus mediates fusion more efficiently.

Modulation of lipid membrane structural and mechanical properties by a peptidomimetic derived from reduced amide scaffold

Understanding how antimicrobial peptidomimetics interact with lipid membranes is important in battling multidrug resistant bacterial pathogens. We study the effects of a recently reported peptidomimetic on lipid bilayer structural and mechanical properties. The compound referred to as E107-3 is synthesized based on the acylated reduced amide scaffold and has been shown to exhibit good antimicrobial potency. Our vesicle leakage assay indicates that the compound increases lipid bilayer permeability. We use micropipette aspiration to explore the kinetic response of giant unilamellar vesicles (GUVs). Exposure to the compound causes the GUV protrusion length L_P to spontaneously increase and then decrease, followed by GUV rupture. Solution atomic force microscopy (AFM) is used to visualize lipid bilayer structural modulation within a nanoscopic regime. Unlike melittin, which produces pore-like structures, the peptidomimetic compound is found to induce nanoscopic heterogeneous structures. Finally, we use AFM-based force spectroscopy to study the impact of the compound on lipid bilayer mechanical properties. We find that incremental addition of the compound to planar lipid bilayers results in a moderate decrease of the bilayer puncture force F_P and a 39% decrease of the bilayer area compressibility modulus K_A. To explain our experimental data, we propose a membrane interaction model encompassing disruption of lipid chain packing and extraction of lipid molecules. The later action mode is supported by our observation of a double-bilayer structure in the presence of fusogenic calcium ions.

A Coiled-Coil Peptide Shaping Lipid Bilayers upon Fusion

A system based on two designed peptides, namely the cationic peptide K, (KIAALKE)3, and its complementary anionic counterpart called peptide E, (EIAALEK)3, has been used as a minimal model for membrane fusion, inspired by SNARE proteins. Although the fact that docking of separate vesicle populations via the formation of a dimeric E/K coiled-coil complex can be rationalized, the reasons for the peptides promoting fusion of vesicles cannot be fully explained. Therefore it is of significant interest to determine how the peptides aid in overcoming energetic barriers during lipid rearrangements leading to fusion. In this study, investigations of the peptides' interactions with neutral PC/PE/cholesterol membranes by fluorescence spectroscopy show that tryptophan-labeled K∗ binds to the membrane (KK∗ ∼6.2 103 M-1), whereas E∗ remains fully water-solvated. 15N-NMR spectroscopy, depth-dependent fluorescence quenching, CD-spectroscopy experiments, and MD simulations indicate a helix orientation of K∗ parallel to the membrane surface. Solid-state 31P-NMR of oriented lipid membranes was used to study the impact of peptide incorporation on lipid headgroup alignment. The membrane-immersed K∗ is found to locally alter the bilayer curvature, accompanied by a change of headgroup orientation relative to the membrane normal and of the lipid composition in the vicinity of the bound peptide. The NMR results were supported by molecular dynamics simulations, which showed that K reorganizes the membrane composition in its vicinity, induces positive membrane curvature, and enhances the lipid tail protrusion probability. These effects are known to be fusion relevant. The combined results support the hypothesis for a twofold role of K in the mechanism of membrane fusion: 1) to bring opposing membranes into close proximity via coiled-coil formation and 2) to destabilize both membranes thereby promoting fusion.