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:The function of proteins is influenced not only by the atomic structure but also by the detailed structure of the solvent surrounding it. Computational studies of protein structure also critically depend on the water structure around the protein. Herein we compare the water structure obtained from molecular dynamics (MD) simulations of galectin-3 in complex with two ligands to crystallographic water molecules observed in the corresponding crystal structures. We computed MD trajectories both in a water box, which mimics a protein in solution, and in a crystallographic unit cell, which mimics a protein in a crystal. The calculations were compared to crystal structures obtained at both cryogenic and room temperature. Two types of analyses of the MD simulations were performed. First, the positions of the crystallographic water molecules were compared to peaks in the MD density after alignment of the protein in each snapshot. The results of this analysis indicate that all simulations reproduce the crystallographic water structure rather poorly. However, if we define the crystallographic water sites based on their distances to nearby protein atoms and follow these sites throughout the simulations, the MD simulations reproduce the crystallographic water sites much better. This shows that the failure of MD simulations to reproduce the water structure around proteins in crystal structures observed both in this and previous studies is caused by the problem of identifying water sites for a flexible and dynamic protein (traditionally done by overlaying the structures). Our local clustering approach solves the problem and shows that the MD simulations reasonably reproduce the water structure observed in crystals. Furthermore, analysis of the crystal MD simulations indicates a few water molecules that are close to unmodeled electron density peaks in the crystal structures, suggesting that crystal MD could be used as a complementary tool for identifying and modelling water in protein crystallography.

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