Dataset Information

Electrophoresis in protein crystal: nonequilibrium molecular dynamics simulations.

ABSTRACT: Electrophoresis of a mixture of NaCl and CaCl2 in a lysozyme crystal is investigated using nonequilibrium molecular dynamics (MD) simulations. Upon exposure to an electric field, the stability of lysozyme is found to decrease slightly. This finding is demonstrated by increases in the root mean-square deviations of the heavy atoms of lysozyme, in the solvent-accessible surface area of hydrophobic residues, and in the number of hydrogen bonds between lysozyme and water. The solvent-accessible surface area of hydrophilic residues changes marginally, and the number of hydrogen bonds between lysozyme molecules decreases. Water molecules tend to align preferentially parallel to the electric field, and the dipole moment along the pore axis increases linearly with increasing field strength. Two pronounced layered structures are observed for Na+ and Ca2+ in the vicinity of protein surface, but only one enriched layer is observed for Cl-. The number distributions of all ions are nearly independent of the electric field. The water coordination numbers of all ions are smaller in the crystal than in aqueous bulk solution; however, the reverse is found for the Cl- coordination numbers of cations. Both the water and the Cl- coordination numbers are insensitive to the electric field. Ion diffusivities in the crystal are approximately 2 orders of magnitude smaller than those in aqueous bulk solution. The drift velocities of ions increase proportionally to the electric field, particularly at high strengths, and depend on ionic charge and coordination with oppositely charged ions. Electrical current exhibits a linear relationship with the field strength. The zero-field electrical conductivity is estimated to be 0.56 S/m, which is very close to 0.61 S/m as predicted by the Nernst-Einstein equation.

SUBMITTER: Hu Z

PROVIDER: S-EPMC2567938 | biostudies-literature | 2008 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Electrophoresis in protein crystal: nonequilibrium molecular dynamics simulations.

Hu Zhongqiao Z Jiang Jianwen J

Biophysical journal 20080718 9

Electrophoresis of a mixture of NaCl and CaCl2 in a lysozyme crystal is investigated using nonequilibrium molecular dynamics (MD) simulations. Upon exposure to an electric field, the stability of lysozyme is found to decrease slightly. This finding is demonstrated by increases in the root mean-square deviations of the heavy atoms of lysozyme, in the solvent-accessible surface area of hydrophobic residues, and in the number of hydrogen bonds between lysozyme and water. The solvent-accessible surf ...[more]

PMID: 18641079

Similar Datasets

Project description:We present results for the steady state shear rheology of non-concatenated, unentangled and marginally entangled ring poly(ethylene oxide) (PEO) melts from detailed, atomistic nonequilibrium molecular dynamics (NEMD) simulations, and compare them to the behavior of the corresponding linear melts. The applied flow field spans a wide range of shear rates, from the linear (Newtonian) to the highly non-linear (described by a power law) regime. For all melts studied, rings are found to exhibit shear thinning but to a lesser degree compared to linear counterparts, mostly due to their reduced deformability and stronger resistance to alignment in the direction of flow. These features are attributed to the more compact structure of ring molecules compared to linear chains; the latter are capable of adopting wider and more open conformations even under shear due to the freedom provided by the free ends. Similar to linear melts, rings also exhibit a first and a second normal stress coefficient; the latter is negative. The ratio of the magnitude of the two coefficients remains practically constant with shear rate and is systematically higher than the corresponding one for linear melts. Emphasis was also given to the statistics of terminal (re-orientational) relaxation times which we computed by analyzing all chains in the simulated systems one by one; it was demonstrated that long time dynamics are strongly heterogeneous both for rings and (especially) linears. Repeating the analysis under flow conditions, and as expected, we found that the applied flow field significantly suppresses dynamic heterogeneity, especially for high shear rates well beyond the Newtonian plateau. Finally, a detailed geometrical analysis revealed that the average population of ring-ring threading events in the longest melt studied here (the PEO-5k ring) remains practically unaffected by the imposed flow rate even at strong shear rates, except for multi-threadings which disappear. To further analyze this peculiar and rather unexpected effect, we computed the corresponding survival times and penetration lengths, and found that the overwhelming majority of threadings under shear are extremely weak constraints, as they are characterized by very small penetration lengths, thus also by short survival times. They are expected therefore to play only a minor (if any) role on chain dynamics.

Dataset Information

Electrophoresis in protein crystal: nonequilibrium molecular dynamics simulations.

Publications

Electrophoresis in protein crystal: nonequilibrium molecular dynamics simulations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets