Dataset Information

Molecular dynamics simulations of membrane proteins under asymmetric ionic concentrations.

ABSTRACT: 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.

SUBMITTER: Khalili-Araghi F

PROVIDER: S-EPMC3787774 | biostudies-literature | 2013 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Molecular dynamics simulations of membrane proteins under asymmetric ionic concentrations.

Khalili-Araghi Fatemeh F Ziervogel Brigitte B Gumbart James C JC Roux Benoît B

The Journal of general physiology 20131001 4

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 ...[more]

PMID: 24081985

Similar Datasets

Project description:The membrane-bound component F(0), which is a major component of the F(0)F(1)-ATP synthase, works as a rotary motor and plays a central role in driving the F(1) component to transform chemiosmotic energy into ATP synthesis. We conducted molecular dynamics simulations of b(2)-free F(0) in a 1-palmitoyl-2-oleoyl-phosphatidylcholine lipid bilayer for tens of nanoseconds with two different protonation states of the cAsp-61 residue at the interface of the a-c complex in the absence of electric fields and under electric fields of +/-0.03 V/nm across the membrane. To our surprise, we observed that the upper half of the N-terminal helix of the c(1) subunit rotated about its axis clockwise by 30 degrees . An energetic analysis revealed that the electrostatic repulsion between this N-terminal helix and subunit c(12) was a major contributor to the observed rotation. A correlation map analysis indicated that the correlated motions of residues in the interface of the a-c complex were significantly reduced by external electric fields. The deuterium order parameter (S(CD)) profile calculated by averaging all the lipids in the F(0)-bound bilayer was not very different from that of the pure bilayer system, in agreement with recent (2)H solid-state NMR experiments. However, by delineating the lipid properties according to their vicinity to F(0), we found that the S(CD) profiles of different lipid shells were prominently different. Lipids close to F(0) formed a more ordered structure. Similarly, the lateral diffusion of lipids on the membrane surface also followed a shell-dependent behavior. The lipids in the proximity of F(0) exhibited very significantly reduced diffusional motion. The numerical value of S(CD) was anticorrelated with that of the diffusion coefficient, i.e., the more ordered lipid structures led to slower lipid diffusion. Our findings will help elucidate the dynamics of F(0) depending on the protonation state and electric field, and may also shed some light on the interactions between the motor F(0) and its surrounding lipids under physiological conditions, which could help to rationalize its extraordinary energy conversion efficiency.

Dataset Information

Molecular dynamics simulations of membrane proteins under asymmetric ionic concentrations.

Publications

Molecular dynamics simulations of membrane proteins under asymmetric ionic concentrations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets