Molecular dynamics simulations of a double unit cell in a protein crystal: Volume relaxation at constant pressure and correlation of motions between the two unit cells

Eight molecular dynamics simulations of a double crystal unit cell of ubiquitin were performed to investigate the effects of simulating at constant pressure and of simulating two unit cells compared to a single unit cell. To examine the influence of different simulation conditions, the constant‐pressure and constant‐volume simulations were each performed with and without counterions and using two different treatments of the long‐range electrostatic interactions (lattice‐sum and reaction‐field methods). The constant‐pressure simulations were analyzed in terms of unit cell deformation and accompanying protein deformations. Energetic and structural properties of the proteins in the simulations of the double unit cell were compared to the results of previously reported one‐unit‐cell simulations. Correlation between the two unit cells was also investigated based on relative translational and rotational movements of the proteins and on dipole fluctuations. The box in the constant‐pressure simulations is found to deform slowly to reach convergence only after 5–10 ns. This deformation does not result from a distortion in the structure of the proteins but rather from changes in protein packing within the unit cell. The results of the double‐unit‐cell simulations are closely similar to the results of the single‐unit‐cell simulations, and little motional correlation is found between the two unit cells. Proteins 2002;48:327–340. © 2002 Wiley‐Liss, Inc.

[1]  Wilfred F. van Gunsteren,et al.  A generalized reaction field method for molecular dynamics simulations , 1995 .

[2]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[3]  Andrew E. Torda,et al.  The GROMOS biomolecular simulation program package , 1999 .

[4]  Gerhard Hummer,et al.  Simulation and Theory of Electrostatic Interactions in Solution: Computational Chemistry, Biophysics and Aqueous Solutions, Santa Fe, New Mexico, U. S. A., 23-25 June 1999 , 1999 .

[5]  E. Pietsch,et al.  Gmelin's Handbuch Der Anorganischen Chemie , 1956 .

[6]  W. V. van Gunsteren,et al.  On the similarity of properties in solution or in the crystalline state: A molecular dynamics study of hen lysozyme , 2000, Journal of biomolecular NMR.

[7]  Jeremy C. Smith,et al.  X-ray diffuse scattering and rigid-body motion in crystalline lysozyme probed by molecular dynamics simulation. , 1998, Journal of molecular biology.

[8]  H. Berendsen,et al.  Simulations of Proteins in Water a , 1986, Annals of the New York Academy of Sciences.

[9]  M. Marchi,et al.  SIMULATION OF A PROTEIN CRYSTAL AT CONSTANT PRESSURE , 1997 .

[10]  D. Nguyen,et al.  On achieving better than 1-A accuracy in a simulation of a large protein: Streptomyces griseus protease A. , 1993, Proceedings of the National Academy of Sciences of the United States of America.

[11]  MD simulation of subtilisin BPN′ in a crystal environment , 1992, Proteins.

[12]  H. Berendsen,et al.  Computer simulation of the dynamics of hydrated protein crystals and its comparison with x-ray data. , 1983, Proceedings of the National Academy of Sciences of the United States of America.

[13]  J. Doucet,et al.  Molecular dynamics studied by analysis of the X-ray diffuse scattering from lysozyme crystals , 1987, Nature.

[14]  J. Clarage,et al.  Liquid-like movements in crystalline insulin , 1988, Nature.

[15]  T. Darden,et al.  Atomic-level accuracy in simulations of large protein crystals. , 1994, Proceedings of the National Academy of Sciences of the United States of America.

[16]  T. Morozova,et al.  Viscoelasticity of protein crystal as a probe of the mechanical properties of a protein molecule. Hen egg-white lysozyme. , 1982, Journal of molecular biology.

[17]  W. F. Gunsteren,et al.  A comparison of the structure and dynamics of avian pancreatic polypeptide hormone in solution and in the crystal , 2004, European Biophysics Journal.

[18]  W. V. van Gunsteren,et al.  Comparison of different schemes to treat long‐range electrostatic interactions in molecular dynamics simulations of a protein crystal , 2001, Proteins.

[19]  Thomas Thüne,et al.  Thermal diffuse X-ray scattering and its contribution to understanding protein dynamics. , 1995, Progress in biophysics and molecular biology.

[20]  J. Doucet,et al.  Diffuse scattering in protein crystallography , 1995, Quarterly Reviews of Biophysics.

[21]  S. Yoneda.  A further implementation of the rotational symmetry boundary conditions for calculations of P4(3)2(1)2 symmetry crystals. , 1997, Journal of molecular graphics & modelling.

[22]  Wilfred F. van Gunsteren,et al.  A molecular dynamics simulation study of chloroform , 1994 .

[23]  Bernard R. Brooks,et al.  Removal of pressure and free energy artifacts in charged periodic systems via net charge corrections to the Ewald potential , 1998 .

[24]  Wilfred F. van Gunsteren,et al.  Determination of force field parameters for molecular simulation by molecular simulation: An application of the weak-coupling method , 1995 .

[25]  Philippe H. Hünenberger,et al.  Optimal charge-shaping functions for the particle–particle—particle–mesh (P3M) method for computing electrostatic interactions in molecular simulations , 2000 .

[26]  Philippe H. Hünenberger,et al.  Calculation of the group-based pressure in molecular simulations. I. A general formulation including Ewald and particle-particle-particle-mesh electrostatics , 2002 .

[27]  P. Argos,et al.  Knowledge‐based protein secondary structure assignment , 1995, Proteins.

[28]  H. Berendsen,et al.  COMPUTER-SIMULATION OF MOLECULAR-DYNAMICS - METHODOLOGY, APPLICATIONS, AND PERSPECTIVES IN CHEMISTRY , 1990 .

[29]  Wilfred F. van Gunsteren,et al.  GROMOS Force Field , 2002 .

[30]  Baldomero Oliva,et al.  Calculation of the group-based pressure in molecular simulations. II. Numerical tests and application to liquid water , 2002 .

[31]  J. Benoit,et al.  Molecular rigid-body displacements in a tetragonal lysozyme crystal confirmed by X-ray diffuse scattering. , 1996, Acta crystallographica. Section D, Biological crystallography.

[32]  Thermal motion in protein crystals estimated using laser‐generated ultrasound and Young's modulus measurements , 1990 .

[33]  G J Williams,et al.  The Protein Data Bank: a computer-based archival file for macromolecular structures. , 1978, Archives of biochemistry and biophysics.

[34]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[35]  D. Moss,et al.  Diffuse X-ray scattering from macromolecular crystals using synchrotron radiation , 1995 .

[36]  Wilfred F. van Gunsteren,et al.  Consistent dielectric properties of the simple point charge and extended simple point charge water models at 277 and 300 K , 1994 .

[37]  Philippe H. Hünenberger,et al.  Lattice-sum methods for computing electrostatic interactions in molecular simulations , 1999 .

[38]  C. Bugg,et al.  Structure of ubiquitin refined at 1.8 A resolution. , 1987, Journal of molecular biology.

[39]  M. Karplus,et al.  Protein dynamics in solution and in a crystalline environment: a molecular dynamics study. , 1982, Biochemistry.

[40]  Molecular dynamics simulation of despentapeptide insulin in a crystalline environment. , 1988 .

[41]  M. Deserno,et al.  HOW TO MESH UP EWALD SUMS. II. AN ACCURATE ERROR ESTIMATE FOR THE PARTICLE-PARTICLE-PARTICLE-MESH ALGORITHM , 1998, cond-mat/9807100.

[42]  J Moult,et al.  Molecular dynamics study of the structure and dynamics of a protein molecule in a crystalline ionic environment, Streptomyces griseus protease A. , 1990, Biochemistry.

[43]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .