Minimum energy path to membrane pore formation and rupture.
Ontology highlight
ABSTRACT:
SUBMITTER: Ting CL
PROVIDER: S-EPMC3225126 | biostudies-literature | 2011 Apr
REPOSITORIES: biostudies-literature
Publications
Minimum energy path to membrane pore formation and rupture.
Ting Christina L CL Appelö Daniel D Wang Zhen-Gang ZG
Physical review letters 20110418 16
We combine dynamic self-consistent field theory with the string method to calculate the minimum energy path to membrane pore formation and rupture. In the regime where nucleation can occur on experimentally relevant time scales, the structure of the critical nucleus is between a solvophilic stalk and a locally thinned membrane. Classical nucleation theory fails to capture these molecular details and significantly overestimates the free energy barrier. Our results suggest that thermally nucleated ...[more]
Similar Datasets
Project description:Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme.
| S-EPMC3369945 | biostudies-literature
Project description:The productive fusion pore in membrane fusion is generally thought to be toroidally shaped. Theoretical studies and recent experiments suggest that its formation, in some scenarios, may be preceded by an initial pore formed near the rim of the extended hemifusion diaphragm (HD), a rim-pore. This rim-pore is characterized by a nontoroidal shape that changes with size. To determine this shape as well as the free energy along the pathway of rim-pore expansion, we derived a simple analytical free energy model. We argue that dilation of HD material via expansion of a rim-pore is favored over a regular, circular pore. Further, the expanding rim-pore faces a free energy barrier that linearly increases with HD size. In contrast, the tension required to expand the rim-pore decreases with HD size. Pore flickering, followed by sudden opening, occurs when the tension in the HD competes with the line energy of the rim-pore, and the rim-pore reaches its equilibrium size before reaching the critical pore size. The experimental observation of flickering and closing fusion pores (kiss-and-run) is very well explained by the observed behavior of rim-pores. Finally, the free energy landscape of rim-pore expansion/HD dilation may very well explain why some cellular fusion reactions, in their attempt to minimize energetic costs, progress via alternative formation and dilation of microscopic hemifusion intermediates.
| S-EPMC4241460 | biostudies-literature
Project description:Daptomycin is a calcium-dependent lipodepsipeptide antibiotic clinically used to treat serious infections caused by Gram-positive pathogens. Its precise mode of action is somewhat controversial; the biggest issue is daptomycin pore formation, which we directly investigated here. We first performed a screening experiment using propidium iodide (PI) entry to Bacillus subtilis cells and chose the optimum and therapeutically relevant conditions (10 µg/ml daptomycin and 1.25 mM CaCl2) for the subsequent analyses. Using conductance measurements on planar lipid bilayers, we show that daptomycin forms nonuniform oligomeric pores with conductance ranging from 120 pS to 14 nS. The smallest conductance unit is probably a dimer; however, tetramers and pentamers occur in the membrane most frequently. Moreover, daptomycin pore-forming activity is exponentially dependent on the applied membrane voltage. We further analyzed the membrane-permeabilizing activity in B. subtilis cells using fluorescence methods [PI and DiSC3(5)]. Daptomycin most rapidly permeabilizes cells with high initial membrane potential and dissipates it within a few minutes. Low initial membrane potential hinders daptomycin pore formation.
| S-EPMC6325215 | biostudies-literature
Project description:Gasdermin-D (GSDMD), a member of the gasdermin protein family, mediates pyroptosis in human and murine cells. Cleaved by inflammatory caspases, GSDMD inserts its N-terminal domain (GSDMDNterm) into cellular membranes and assembles large oligomeric complexes permeabilizing the membrane. So far, the mechanisms of GSDMDNterm insertion, oligomerization, and pore formation are poorly understood. Here, we apply high-resolution (≤ 2 nm) atomic force microscopy (AFM) to describe how GSDMDNterm inserts and assembles in membranes. We observe GSDMDNterm inserting into a variety of lipid compositions, among which phosphatidylinositide (PI(4,5)P2) increases and cholesterol reduces insertion. Once inserted, GSDMDNterm assembles arc-, slit-, and ring-shaped oligomers, each of which being able to form transmembrane pores. This assembly and pore formation process is independent on whether GSDMD has been cleaved by caspase-1, caspase-4, or caspase-5. Using time-lapse AFM, we monitor how GSDMDNterm assembles into arc-shaped oligomers that can transform into larger slit-shaped and finally into stable ring-shaped oligomers. Our observations translate into a mechanistic model of GSDMDNterm transmembrane pore assembly, which is likely shared within the gasdermin protein family.
| S-EPMC6043855 | biostudies-literature
Project description:The U6 RNA internal stem-loop (U6 ISL) is a highly conserved domain of the spliceosome that is important for pre-mRNA splicing. The U6 ISL contains an internal loop that is in equilibrium between two conformations controlled by the protonation state of an adenine (pK(a)=6.5). Lower pH favors formation of a protonated C-A(+) wobble pair and base flipping of the adjacent uracil. Higher pH favors stacking of the uracil and allows an essential metal ion to bind at this position. Here, we define the minimal-energy path for this conformational transition. To do this, we solved the U6 ISL structure at higher pH (8.0) in order to eliminate interference from the low-pH conformer. This structure reveals disruption of the protonated C-A(+) pair and formation of a new C-U pair, which explains the preference for a stacked uracil at higher pH. Next, we used nudged elastic band molecular dynamics simulations to calculate the minimum-energy path between the two conformations. Our results indicate that the C-U pair is dynamic, which allows formation of the more stable C-A(+) pair upon adenine protonation. After formation of the C-A(+) pair, the unpaired uracil follows a minor-groove base-flipping pathway. Molecular dynamics simulations suggest that the extrahelical uracil is stabilized by contacts with the adjacent helix.
| S-EPMC2799254 | biostudies-literature
Project description:The increasing interest in novel foldamer constructs demands an accurate computational treatment on an extensive timescale. However, it is still a challenge to derive a force field (FF) that can reproduce the experimentally known fold while also allowing the spontaneous exploration of other structures. Here, aiming at a realistic reproduction of backbone torsional barriers, the relevant proper dihedrals of acyclic β2-, β3- and β2,3-amino acids were added to the CHARMM FF and optimized using a novel, self-consistent iterative procedure based on quantum chemical relaxed scans. The new FF was validated by molecular dynamics simulations on three acyclic peptides. While they resided most of the time in their preferred fold (>80 % in helices and >50 % in hairpin), they also visited other conformations. Owing to the CHARMM36m-consistent parametrization, the proposed extension is suitable for exploring new foldamer structures and assemblies, and their interactions with diverse biomolecules.
| S-EPMC6686720 | biostudies-literature
Project description:A quantum mechanical/molecular mechanical minimum free energy path (QM/MM-MFEP) method was developed to calculate the redox free energies of large systems in solution with greatly enhanced efficiency for conformation sampling. The QM/MM-MFEP method describes the thermodynamics of a system on the potential of mean force surface of the solute degrees of freedom. The molecular dynamics (MD) sampling is only carried out with the QM subsystem fixed. It thus avoids "on-the-fly" QM calculations and thus overcomes the high computational cost in the direct QM/MM MD sampling. In the applications to two metal complexes in aqueous solution, the new QM/MM-MFEP method yielded redox free energies in good agreement with those calculated from the direct QM/MM MD method. Two larger biologically important redox molecules, lumichrome and riboflavin, were further investigated to demonstrate the efficiency of the method. The enhanced efficiency and uncompromised accuracy are especially significant for biochemical systems. The QM/MM-MFEP method thus provides an efficient approach to free energy simulation of complex electron transfer reactions.
| S-EPMC2736613 | biostudies-other
Project description:Many cationic antimicrobial peptides (AMPs) target the unique lipid composition of the prokaryotic cell membrane. However, the micromolar activities common for these peptides are considered weak in comparison to nisin, which follows a targeted, pore-forming mode of action. Here we show that AMPs can be modified with a high-affinity targeting module, which enables membrane permeabilization at low concentration. Magainin 2 and a truncated peptide analog were conjugated to vancomycin using click chemistry, and could be directed towards specific membrane embedded receptors both in model membrane systems and whole cells. Compared with untargeted vesicles, a gain in permeabilization efficacy of two orders of magnitude was reached with large unilamellar vesicles that included lipid II, the target of vancomycin. The truncated vancomycin-peptide conjugate showed an increased activity against vancomycin resistant Enterococci, whereas the full-length conjugate was more active against a targeted eukaryotic cell model: lipid II containing erythrocytes. This study highlights that AMPs can be made more selective and more potent against biological membranes that contain structures that can be targeted.
| S-EPMC3387250 | biostudies-literature
Project description:Protegrin is an antimicrobial peptide with a β-hairpin structure stabilized by a pair of disulfide bonds. It has been extensively studied by solid-state NMR and computational methods. Here we use implicit membrane models to examine the binding of monomers on the surface and in the interior of the membrane, the energetics of dimerization, the binding to membrane pores, and the stability of different membrane barrel structures in pores. Our results challenge a number of conclusions based on previous experimental and theoretical work. The burial of monomers into the membrane interior is found to be unfavorable for any membrane thickness. Because of its imperfect amphipathicity, protegrin binds weakly, at most, on the surface of zwitterionic membranes. However, it binds more favorably onto toroidal pores. Anionic charge on the membrane facilitates the binding due to electrostatic interactions. Solid-state NMR results have suggested a parallel NCCN association of monomers in dimers and association of dimers to form octameric or decameric β-barrels. We find that this structure is not energetically plausible for binding to bilayers, because in this configuration the hydrophobic sides of two monomers point in opposite directions. In contrast, the antiparallel NCCN and especially the parallel NCNC octamers are stable and exhibit a favorable binding energy to the pore. The results of 100-ns simulations in explicit bilayers corroborate the higher stability of the parallel NCNC barrel compared with the parallel NCCN barrel. The ability to form pores in zwitterionic membranes provides a rationalization for the peptide's cytotoxicity. The discrepancies between our results and experiment are discussed, and new experiments are proposed to resolve them and to test the validity of the models.
| S-EPMC3566446 | biostudies-literature
Project description:Bax protein plays a key role in mitochondrial membrane permeabilization and cytochrome c release upon apoptosis. Our recent data have indicated that the 20-residue C-terminal peptide of Bax (BaxC-KK; VTIFVAGVLTASLTIWKKMG), when expressed intracellularly, translocates to the mitochondria and exerts lethal effect on cancer cells. Moreover, the BaxC-KK peptide, as well as two mutants where the two lysines are replaced with glutamate (BaxC-EE) or leucine (BaxC-LL), have been shown to form relatively large pores in lipid membranes, composed of up to eight peptide molecules per pore. Here the pore structure is analyzed by polarized Fourier transform infrared, circular dichroism, and fluorescence experiments on the peptides reconstituted in phospholipid membranes. The peptides assume an ?/?-type secondary structure within membranes. Both ?-strands and ?-helices are significantly (by 30-60 deg) tilted relative to the membrane normal. The tryptophan residue embeds into zwitterionic membranes at 8-9 Å from the membrane center. The membrane anionic charge causes a deeper insertion of tryptophan for BaxC-KK and BaxC-LL but not for BaxC-EE. Combined with the pore stoichiometry determined earlier, these structural constraints allow construction of a model of the pore where eight peptide molecules form an "?/?-ring" structure within the membrane. These results identify a strong membranotropic activity of Bax C-terminus and propose a new mechanism by which peptides can efficiently perforate cell membranes. Knowledge on the pore forming mechanism of the peptide may facilitate development of peptide-based therapies to kill cancer or other detrimental cells such as bacteria or fungi.
| S-EPMC4537061 | biostudies-literature