Browse
Submit Data
Databases
API
Help

Dataset Information

15 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

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

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

SUBMITTER: Bubnis G

PROVIDER: S-EPMC7822740 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

MD Simulations Reveal Complex Water Paths in Squalene-Hopene Cyclase: Tunnel-Obstructing Mutations Increase the Flow of Water in the Active Site.

Project description:Squalene-hopene cyclase catalyzes the cyclization of squalene to hopanoids. A previous study has identified a network of tunnels in the protein, where water molecules have been indicated to move. Blocking these tunnels by site-directed mutagenesis was found to change the activation entropy of the catalytic reaction from positive to negative with a concomitant lowering of the activation enthalpy. As a consequence, some variants are faster and others are slower than the wild type (wt) in vitro under optimal reaction conditions for the wt. In this study, molecular dynamics (MD) simulations have been performed for the wt and the variants to investigate how the mutations affect the protein structure and the water flow in the enzyme, hypothetically influencing the activation parameters. Interestingly, the tunnel-obstructing variants are associated with an increased flow of water in the active site, particularly close to the catalytic residue Asp376. MD simulations with the substrate present in the active site indicate that the distance for the rate-determining proton transfer between Asp376 and the substrate is longer in the tunnel-obstructing protein variants than in the wt. On the basis of the previous experimental results and the current MD results, we propose that the tunnel-obstructing variants, at least partly, could operate by a different catalytic mechanism, where the proton transfer may have contributions from a Grotthuss-like mechanism.

| S-EPMC6645472 | biostudies-literature

Exploring peptide-membrane interactions with coarse-grained MD simulations.

Project description:The interaction of ?-helical peptides with lipid bilayers is central to our understanding of the physicochemical principles of biological membrane organization and stability. Mutations that alter the position or orientation of an ?-helix within a membrane, or that change the probability that the ?-helix will insert into the membrane, can alter a range of membrane protein functions. We describe a comparative coarse-grained molecular dynamics simulation methodology, based on self-assembly of a lipid bilayer in the presence of an ?-helical peptide, which allows us to model membrane transmembrane helix insertion. We validate this methodology against available experimental data for synthetic model peptides (WALP23 and LS3). Simulation-based estimates of apparent free energies of insertion into a bilayer of cystic fibrosis transmembrane regulator-derived helices correlate well with published data for translocon-mediated insertion. Comparison of values of the apparent free energy of insertion from self-assembly simulations with those from coarse-grained molecular dynamics potentials of mean force for model peptides, and with translocon-mediated insertion of cystic fibrosis transmembrane regulator-derived peptides suggests a nonequilibrium model of helix insertion into bilayers.

| S-EPMC3077683 | biostudies-literature

Energetics of water permeation through fullerene membrane.

Project description:Lipid bilayer membranes are important as fundamental structures in biology and possess characteristic water-permeability, stability, and mechanical properties. Water permeation through a lipid bilayer membrane occurs readily, and more readily at higher temperature, which is largely due to an enthalpy cost of the liquid-to-gas phase transition of water. A fullerene bilayer membrane formed by dissolution of a water-soluble fullerene, Ph(5)C(60)K, has now been shown to possess properties entirely different from those of the lipid membranes. The fullerene membrane is several orders of magnitude less permeable to water than a lipid membrane, and the permeability decreases at higher temperature. Water permeation is burdened by a very large entropy loss and may be favored slightly by an enthalpy gain, which is contrary to the energetics observed for the lipid membrane. We ascribe this energetics to favorable interactions of water molecules to the surface of the fullerene molecules as they pass through the clefts of the rigid fullerene bilayer. The findings provide possibilities of membrane design in science and technology.

| S-EPMC1975687 | biostudies-literature

Water transport in human aquaporin-4: molecular dynamics (MD) simulations.

Project description:Aquaporin-4 (AQP4) is the predominant water channel in the central nervous system, where it has been reported to be involved in many pathophysiological roles including water transport. In this paper, the AQP4 tetramer was modeled from its PDB structure file, embedded in a palmitoyl-oleoyl-phosphatidyl-choline (POPC) lipid bilayer, solvated in water, then minimized and equilibrated by means of molecular dynamics simulations. Analysis of the equilibrated structure showed that the central pore along the fourfold axis of the tetramers is formed with hydrophobic amino acid residues. In particular, Phe-195, Leu-191 and Leu-75, form the narrowest part of the pore. Therefore water molecules are not expected to transport through the central pore, which was confirmed by MD simulations. Each monomer of the AQP4 tetramers forms a channel whose walls consist mostly of hydrophilic residues. There are eight water molecules in single file observed in each of the four channels, transporting through the selectivity filter containing Arg-216, His-201, Phe-77, Ala-210, and the two conserved Asn-Pro-Ala (NPA) motifs containing Asn-213 and Asn-97. By using Brownian dynamics fluctuation-dissipation-theorem (BD-FDT), the overall free-energy profile was obtained for water transporting through AQP4 for the first time, which gives a complete map of the entire channel of water permeation.

| S-EPMC3809129 | biostudies-literature

Enantiomerization of Axially Chiral Biphenyls: Polarizable MD Simulations in Water and Butylmethylether.

Project description:In this study, we investigate the influence of chiral and achiral cations on the enantiomerization of biphenylic anions in n-butylmethylether and water. In addition to the impact of the cations and solvent molecules on the free energy profile of rotation, we also explore if chirality transfer between a chiral cation and the biphenylic anion is possible, i.e., if pairing with a chiral cation can energetically favour one conformer of the anion via diastereomeric complex formation. The quantum-mechanical calculations are accompanied by polarizable MD simulations using umbrella sampling to study the impact of solvents of different polarity in more detail. We also discuss how accurate polarizable force fields for biphenylic anions can be constructed from quantum-mechanical reference data.

| S-EPMC7503397 | biostudies-literature

Simple Model of Protein Energetics To Identify Ab Initio Folding Transitions from All-Atom MD Simulations of Proteins.

Project description:A fundamental requirement to predict the native conformation, address questions of sequence design and optimization, and gain insights into the folding mechanisms of proteins lies in the definition of an unbiased reaction coordinate that reports on the folding state without the need to compare it to reference values, which might be unavailable for new (designed) sequences. Here, we introduce such a reaction coordinate, which does not depend on previous structural knowledge of the native state but relies solely on the energy partition within the protein: the spectral gap of the pair nonbonded energy matrix (ENergy Gap, ENG). This quantity can be simply calculated along unbiased MD trajectories. We show that upon folding the gap increases significantly, while its fluctuations are reduced to a minimum. This is consistently observed for a diverse set of systems and trajectories. Our approach allows one to promptly identify residues that belong to the folding core as well as residues involved in non-native contacts that need to be disrupted to guide polypeptides to the folded state. The energy gap and fluctuations criteria are then used to develop an automatic detection system which allows us to extract and analyze folding transitions from a generic MD trajectory. We speculate that our method can be used to detect conformational ensembles in dynamic and intrinsically disordered proteins, revealing potential preorganization for binding.

| S-EPMC8009504 | biostudies-literature

Molecular simulations of dodecyl-?-maltoside micelles in water: influence of the headgroup conformation and force field parameters.

Project description:This paper deals with the development and validation of new potential parameter sets, based on the CHARMM36 and GLYCAM06 force fields, to simulate micelles of the two anomeric forms (? and ?) of N-dodecyl-?-maltoside (C(12)G(2)), a surfactant widely used in the extraction and purification of membrane proteins. In this context, properties such as size, shape, internal structure, and hydration of the C(12)G(2) anomer micelles were thoroughly investigated by molecular dynamics simulations and the results compared with experiments. Additional simulations were also performed with the older CHARMM22 force field for carbohydrates (Kuttel, M.; et al. J. Comput. Chem. 2002, 23, 1236-1243). We find that our CHARMM and GLYCAM parameter sets yield similar results in the case of properties related to the micelle structure but differ for other properties such as the headgroup conformation or the micelle hydration. In agreement with experiments, our results show that for all model potentials the ?-C(12)G(2) micelles have a more pronounced ellipsoidal shape than those containing ? anomers. The computed radius of gyration is 20.2 and 25.4 Å for the ?- and ?-anomer micelles, respectively. Finally, we show that depending on the potential the water translational diffusion of the interfacial water is 7-11.5 times slower than that of bulk water due to the entrapment of the water in the micelle crevices. This retardation is independent of the headgroup in ?- or ?-anomers.

| S-EPMC3076730 | biostudies-literature

Methodologies for the analysis of instantaneous lipid diffusion in MD simulations of large membrane systems.

Project description:Interactions between lipids and membrane proteins play a key role in determining the nanoscale dynamic and structural properties of biological membranes. Molecular dynamics (MD) simulations provide a valuable tool for studying membrane models, complementing experimental approaches. It is now possible to simulate large membrane systems, such as simplified models of bacterial and viral envelope membranes. Consequently, there is a pressing need to develop tools to visualize and quantify the dynamics of these immense systems, which typically comprise millions of particles. To tackle this issue, we have developed visual and quantitative analyses of molecular positions and their velocity field using path line, vector field and streamline techniques. This allows us to highlight large, transient flow-like movements of lipids and to better understand crowding within the lipid bilayer. The current study focuses on visualization and analysis of lipid dynamics. However, the methods are flexible and can be readily applied to e.g. proteins and nanoparticles within large complex membranes. The protocols developed here are readily accessible both as a plugin for the molecular visualization program VMD and as a module for the MDAnalysis library.

| S-EPMC4208077 | biostudies-literature

Membrane association of the PTEN tumor suppressor: neutron scattering and MD simulations reveal the structure of protein-membrane complexes.

Project description:Neutron reflection (NR) from planar interfaces is an emerging technology that provides unique and otherwise inaccessible structural information on disordered molecular systems such as membrane proteins associated with fluid bilayers, thus addressing one of the remaining challenges of structural biology. Although intrinsically a low-resolution technique, using structural information from crystallography or NMR allows the construction of NR models that describe the architecture of protein-membrane complexes at high resolution. In addition, a combination of these methods with molecular dynamics (MD) simulations has the potential to reveal the dynamics of protein interactions with the bilayer in atomistic detail. We review recent advances in this area by discussing the application of these techniques to the complex formed by the PTEN phosphatase with the plasma membrane. These studies provide insights in the cellular regulation of PTEN, its interaction with PI(4,5)P2 in the inner plasma membrane and the pathway by which its substrate, PI(3,4,5)P3, accesses the PTEN catalytic site.

| S-EPMC4388770 | biostudies-literature

Common mechanism unites membrane poration by amyloid and antimicrobial peptides.

Project description:Poration of bacterial membranes by antimicrobial peptides such as magainin 2 is a significant activity performed by innate immune systems. Pore formation by soluble forms of amyloid proteins such as islet amyloid polypeptide (IAPP) is implicated in cell death in amyloidoses. Similarities in structure and poration activity of these two systems suggest a commonality of mechanism. Here, we investigate and compare the mechanisms by which these peptides induce membrane leakage and bacterial cell death through the measurement of liposome leakage kinetics and bacterial growth inhibition. For both systems, leakage occurs through the nucleation-dependent formation of stable membrane pores. Remarkably, we observe IAPP and magainin 2 to be fully cross-cooperative in the induction of leakage and inhibition of bacterial growth. The effects are dramatic, with mixtures of these peptides showing activities >100-fold greater than simple sums of the activities of individual peptides. Direct protein-protein interactions cannot be the origin of cooperativity, as IAPP and its enantiomer D-IAPP are equally cross-cooperative. We conclude that IAPP and magainin 2 induce membrane leakage and cytotoxicity through a shared, cross-cooperative, tension-induced poration mechanism.

| S-EPMC3631635 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data