Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

The implementation of slab geometry for membrane-channel molecular dynamics simulations.

ABSTRACT: Slab geometric boundary conditions are applied in the molecular dynamics simulation of a simple membrane-channel system. The results of the simulation were compared to those of an analogous system using normal three-dimensional periodic boundary conditions. Analysis of the dynamics and electrostatics of the system show that slab geometric periodicity eliminates the artificial bulk water orientational polarization that is present while using normal three-dimensional periodicity. Furthermore, even though the water occupancy and volume of our simple channel is the same when using either method, the electrostatic properties are considerably different when using slab geometry. In particular, the orientational polarization of water is seen to be different in the interior of the channel. This gives rise to a markedly different electric field within the channel. We discuss the implications of slab geometry for the future simulation of this type of system and for the study of channel transport properties.

SUBMITTER: Bostick D

PROVIDER: S-EPMC1303069 | biostudies-literature | 2003 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

The implementation of slab geometry for membrane-channel molecular dynamics simulations.

Bostick David D Berkowitz Max L ML

Biophysical journal 20030701 1

Slab geometric boundary conditions are applied in the molecular dynamics simulation of a simple membrane-channel system. The results of the simulation were compared to those of an analogous system using normal three-dimensional periodic boundary conditions. Analysis of the dynamics and electrostatics of the system show that slab geometric periodicity eliminates the artificial bulk water orientational polarization that is present while using normal three-dimensional periodicity. Furthermore, even ...[more]

PMID: 12829468

Similar Datasets

Molecular dynamics simulations of Piezo1 channel opening by increases in membrane tension.

Project description:Piezo1 is a mechanosensitive channel involved in many cellular functions and responsible for sensing shear stress and pressure forces in cells. Piezo1 has a unique trilobed topology with a curved membrane region in the closed state. It has been suggested that upon activation Piezo1 adopts a flattened conformation, but the molecular and structural changes underpinning the Piezo1 gating and opening mechanisms and how the channel senses forces in the membrane remain elusive. Here, we used molecular dynamics simulations to reveal the structural rearrangements that occur when Piezo1 moves from a closed to an open state in response to increased mechanical tension applied to a model membrane. We find that membrane stretching causes Piezo1 to flatten and expand its blade region, resulting in tilting and lateral movement of the pore-lining transmembrane helices 37 and 38. This is associated with the opening of the channel and movement of lipids out of the pore region. Our results reveal that because of the rather loose packing of Piezo1 pore region, movement of the lipids outside the pore region is critical for the opening of the pore. Our simulations also suggest synchronous flattening of the Piezo1 blades during Piezo1 activation. The flattened structure lifts the C-terminal extracellular domain up, exposing it more to the extracellular space. Our studies support the idea that it is the blade region of Piezo1 that senses tension in the membrane because pore opening failed in the absence of the blades. Additionally, our simulations reveal that upon opening, water molecules occupy lateral fenestrations in the cytosolic region of Piezo1, which might be likely paths for ion permeation. Our results provide a model for how mechanical force opens the Piezo1 channel and thus how it might couple mechanical force to biological response.

| S-EPMC8105709 | biostudies-literature

Implementation of force distribution analysis for molecular dynamics simulations.

Project description:BackgroundThe way mechanical stress is distributed inside and propagated by proteins and other biopolymers largely defines their function. Yet, determining the network of interactions propagating internal strain remains a challenge for both, experiment and theory. Based on molecular dynamics simulations, we developed force distribution analysis (FDA), a method that allows visualizing strain propagation in macromolecules.ResultsTo be immediately applicable to a wide range of systems, FDA was implemented as an extension to Gromacs, a commonly used package for molecular simulations. The FDA code comes with an easy-to-use command line interface and can directly be applied to every system built using Gromacs. We provide an additional R-package providing functions for advanced statistical analysis and presentation of the FDA data.ConclusionsUsing FDA, we were able to explain the origin of mechanical robustness in immunoglobulin domains and silk fibers. By elucidating propagation of internal strain upon ligand binding, we previously also successfully revealed the functionality of a stiff allosteric protein. FDA thus has the potential to be a valuable tool in the investigation and rational design of mechanical properties in proteins and nano-materials.

| S-EPMC3094238 | biostudies-literature

Self-assembly of a simple membrane protein: coarse-grained molecular dynamics simulations of the influenza M2 channel.

Project description:The transmembrane (TM) domain of the M2 channel protein from influenza A is a homotetrameric bundle of alpha-helices and provides a model system for computational approaches to self-assembly of membrane proteins. Coarse-grained molecular dynamics (CG-MD) simulations have been used to explore partitioning into a membrane of M2 TM helices during bilayer self-assembly from lipids. CG-MD is also used to explore tetramerization of preinserted M2 TM helices. The M2 helix monomer adopts a membrane spanning orientation in a lipid (DPPC) bilayer. Multiple extended CG-MD simulations (5 x 5 micros) were used to study the tetramerization of inserted M2 helices. The resultant tetramers were evaluated in terms of the most populated conformations and the dynamics of their interconversion. This analysis reveals that the M2 tetramer has 2x rotationally symmetrical packing of the helices. The helices form a left-handed bundle, with a helix tilt angle of approximately 16 degrees. The M2 helix bundle generated by CG-MD was converted to an atomistic model. Simulations of this model reveal that the bundle's stability depends on the assumed protonation state of the H37 side chains. These simulations alongside comparison with recent x-ray (3BKD) and NMR (2RLF) structures of the M2 bundle suggest that the model yielded by CG-MD may correspond to a closed state of the channel.

| S-EPMC2553146 | biostudies-literature

Sidekick for Membrane Simulations: Automated Ensemble Molecular Dynamics Simulations of Transmembrane Helices.

Project description:The interactions of transmembrane (TM) α-helices with the phospholipid membrane and with one another are central to understanding the structure and stability of integral membrane proteins. These interactions may be analyzed via coarse grained molecular dynamics (CGMD) simulations. To obtain statistically meaningful analysis of TM helix interactions, large (N ca. 100) ensembles of CGMD simulations are needed. To facilitate the running and analysis of such ensembles of simulations, we have developed Sidekick, an automated pipeline software for performing high throughput CGMD simulations of α-helical peptides in lipid bilayer membranes. Through an end-to-end approach, which takes as input a helix sequence and outputs analytical metrics derived from CGMD simulations, we are able to predict the orientation and likelihood of insertion into a lipid bilayer of a given helix of a family of helix sequences. We illustrate this software via analyses of insertion into a membrane of short hydrophobic TM helices containing a single cationic arginine residue positioned at different positions along the length of the helix. From analyses of these ensembles of simulations, we estimate apparent energy barriers to insertion which are comparable to experimentally determined values. In a second application, we use CGMD simulations to examine the self-assembly of dimers of TM helices from the ErbB1 receptor tyrosine kinase and analyze the numbers of simulation repeats necessary to obtain convergence of simple descriptors of the mode of packing of the two helices within a dimer. Our approach offers a proof-of-principle platform for the further employment of automation in large ensemble CGMD simulations of membrane proteins.

| S-EPMC4871227 | biostudies-literature

Structure refinement of membrane proteins via molecular dynamics simulations.

Project description:A refinement protocol based on physics-based techniques established for water soluble proteins is tested for membrane protein structures. Initial structures were generated by homology modeling and sampled via molecular dynamics simulations in explicit lipid bilayer and aqueous solvent systems. Snapshots from the simulations were selected based on scoring with either knowledge-based or implicit membrane-based scoring functions and averaged to obtain refined models. The protocol resulted in consistent and significant refinement of the membrane protein structures similar to the performance of refinement methods for soluble proteins. Refinement success was similar between sampling in the presence of lipid bilayers and aqueous solvent but the presence of lipid bilayers may benefit the improvement of lipid-facing residues. Scoring with knowledge-based functions (DFIRE and RWplus) was found to be as good as scoring using implicit membrane-based scoring functions suggesting that differences in internal packing is more important than orientations relative to the membrane during the refinement of membrane protein homology models.

| S-EPMC6013386 | biostudies-literature

OmpT: molecular dynamics simulations of an outer membrane enzyme.

Project description:Five molecular dynamics simulations (total duration >25 ns) have been performed on the Escherichia coli outer membrane protease OmpT embedded in a dimyristoylphosphatidylcholine lipid bilayer. Globally the protein is conformationally stable. Some degree of tilt of the beta-barrel is observed relative to the bilayer plane. The greatest degree of conformational flexibility is seen in the extracellular loops. A complex network of fluctuating H-bonds is formed between the active site residues, such that the Asp210-His212 interaction is maintained throughout, whereas His212 and Asp83 are often bridged by a water molecule. This supports a catalytic mechanism whereby Asp83 and His212 bind a water molecule that attacks the peptide carbonyl. A configuration yielded by docking calculations of OmpT simulation snapshots and a model substrate peptide Ala-Arg-Arg-Ala was used as the starting point for an extended Huckel calculation on the docked peptide. These placed the lowest unoccupied molecular orbital mainly on the carbon atom of the central C=O in the scissile peptide bond, thus favoring attack on the central peptide by the water held by residues Asp83 and His212. The trajectories of water molecules reveal exchange of waters between the intracellular face of the membrane and the interior of the barrel but no exchange at the extracellular mouth. This suggests that the pore-like region in the center of OmpT may enable access of water to the active site from below. The simulations appear to reveal the presence of specific lipid interaction sites on the surface of the OmpT barrel. This reveals the ability of extended MD simulations to provide meaningful information on protein-lipid interactions.

| S-EPMC1304768 | biostudies-literature

Formation Mechanism of Ion Channel in Channelrhodopsin-2: Molecular Dynamics Simulation and Steering Molecular Dynamics Simulations.

Project description:Channelrhodopsin-2 (ChR2) is a light-activated and non-selective cationic channel protein that can be easily expressed in specific neurons to control neuronal activity by light. Although ChR2 has been extensively used as an optogenetic tool in neuroscience research, the molecular mechanism of cation channel formation following retinal photoisomerization in ChR2 is not well understood. In this paper, studies of the closed and opened state ChR2 structures are presented. The formation of the cationic channel is elucidated in atomic detail using molecular dynamics simulations on the all-trans-retinal (ChR2-trans) configuration of ChR2 and its isomerization products, 13-cis-retinal (ChR2-cis) configuration, respectively. Photoisomerization of the retinal-chromophore causes the destruction of interactions among the crucial residues (e.g., E90, E82, N258, and R268) around the channel and the extended H-bond network mediated by numerous water molecules, which opens the pore. Steering molecular dynamics (SMD) simulations show that the electrostatic interactions at the binding sites in intracellular gate (ICG) and central gate (CG) can influence the transmembrane transport of Na+ in ChR2-cis obviously. Potential of mean force (PMF) constructed by SMD and umbrella sampling also found the existing energy wells at these two binding sites during the transportation of Na+. These wells partly hinder the penetration of Na+ into cytoplasm through the ion channel. This investigation provides a theoretical insight on the formation mechanism of ion channels and the mechanism of ion permeation.

| S-EPMC6695816 | biostudies-literature

Determination of membrane-insertion free energies by molecular dynamics simulations.

Project description:The accurate prediction of membrane-insertion probability for arbitrary protein sequences is a critical challenge to identifying membrane proteins and determining their folded structures. Although algorithms based on sequence statistics have had moderate success, a complete understanding of the energetic factors that drive the insertion of membrane proteins is essential to thoroughly meeting this challenge. In the last few years, numerous attempts to define a free-energy scale for amino-acid insertion have been made, yet disagreement between most experimental and theoretical scales persists. However, for a recently resolved water-to-bilayer scale, it is found that molecular dynamics simulations that carefully mimic the conditions of the experiment can reproduce experimental free energies, even when using the same force field as previous computational studies that were cited as evidence of this disagreement. Therefore, it is suggested that experimental and simulation-based scales can both be accurate and that discrepancies stem from disparities in the microscopic processes being considered rather than methodological errors. Furthermore, these disparities make the development of a single universally applicable membrane-insertion free energy scale difficult.

| S-EPMC3283816 | biostudies-literature

Molecular dynamics simulations of membrane proteins under asymmetric ionic concentrations.

Project description:A computational method is developed to allow molecular dynamics simulations of biomembrane systems under realistic ionic gradients and asymmetric salt concentrations while maintaining the conventional periodic boundary conditions required to minimize finite-size effects in an all-atom explicit solvent representation. The method, which consists of introducing a nonperiodic energy step acting on the ionic species at the edge of the simulation cell, is first tested with illustrative applications to a simple membrane slab model and a phospholipid membrane bilayer. The nonperiodic energy-step method is then used to calculate the reversal potential of the bacterial porin OmpF, a large cation-specific ?-barrel channel, by simulating the I-V curve under an asymmetric 10:1 KCl concentration gradient. The calculated reversal potential of 28.6 mV is found to be in excellent agreement with the values of 26-27 mV measured from lipid bilayer experiments, thereby demonstrating that the method allows realistic simulations of nonequilibrium membrane transport with quantitative accuracy. As a final example, the pore domain of Kv1.2, a highly selective voltage-activated K(+) channel, is simulated in a lipid bilayer under conditions that recreate, for the first time, the physiological K(+) and Na(+) concentration gradients and the electrostatic potential difference of living cells.

| S-EPMC3787774 | biostudies-literature

Modeling of annexin A2-Membrane interactions by molecular dynamics simulations.

Project description:The annexins are a family of Ca2+-regulated phospholipid binding proteins that are involved in membrane domain organization and membrane trafficking. Although they are widely studied and crystal structures are available for several soluble annexins their mode of membrane association has never been studied at the molecular level. Here we obtained molecular information on the annexin-membrane interaction that could serve as paradigm for the peripheral membrane association of cytosolic proteins by Molecular Dynamics simulations. We analyzed systems containing the monomeric annexin A2 (AnxA2), a membrane with negatively charged phosphatidylserine (POPS) lipids as well as Ca2+ ions. On the atomic level we identify the AnxA2 orientations and the respective residues which display the strongest interaction with Ca2+ ions and the membrane. The simulation results fully agree with earlier experimental findings concerning the positioning of bound Ca2+ ions. Furthermore, we identify for the first time a significant interaction between lysine residues of the protein and POPS lipids that occurs independently of Ca2+ suggesting that AnxA2-membrane interactions can also occur in a low Ca2+ environment. Finally, by varying Ca2+ concentrations and lipid composition in our simulations we observe a calcium-induced negative curvature of the membrane as well as an AnxA2-induced lipid ordering.

| S-EPMC5609761 | 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