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: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: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: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: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: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: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: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
Project description:Medin, a 50-amino-acid cleavage product of the milk fat globule-EGF factor 8 protein, is one of the most common forms of localized amyloid found in the vasculature of individuals older than 50 years. Medin induces endothelial dysfunction and vascular inflammation, yet despite its prevalence in the human aorta and multiple arterial beds, little is known about the nature of its pathology. Medin oligomers have been implicated in the pathology of aortic aneurysm, aortic dissection, and more recently, vascular dementia. Recent in vitro biomechanical measurements found increased oligomer levels in aneurysm patients with altered aortic wall integrity. Our results suggest an oligomer-mediated toxicity mechanism for medin pathology. Using lipid bilayer electrophysiology, we show that medin oligomers induce ionic membrane permeability by pore formation. Pore activity was primarily observed for preaggregated medin species from the growth-phase and rarely for lag-phase species. Atomic force microscopy (AFM) imaging of medin aggregates at different stages of aggregation revealed the gradual formation of flat domains resembling the morphology of supported lipid bilayers. Transmission electron microscopy images showed the coexistence of compact oligomers, largely consistent with the AFM data, and larger protofibrillar structures. Circular dichroism spectroscopy revealed the presence of largely disordered species and suggested the presence of ?-sheets. This observation and the significantly lower thioflavin T fluorescence emitted by medin aggregates compared to amyloid-? fibrils, along with the absence of amyloid fibers in the AFM and transmission electron microscopy images, suggest that medin aggregation into pores follows a nonamyloidogenic pathway. In silico modeling by molecular dynamics simulations provides atomic-level structural detail of medin pores with the CNpNC barrel topology and diameters comparable to values estimated from experimental pore conductances.
| S-EPMC7264854 | biostudies-literature