Fluctuation and cross-correlation analysis of protein motions observed in nanosecond molecular dynamics simulations.

Nanosecond molecular dynamics simulations of bovine pancreatic trypsin inhibitor and lysozyme in water are analyzed in terms of backbone atomic positional fluctuations and dynamical cross-correlations. It is found that although the molecular systems are stable, B-factors calculated over a time period as long as 500 ps are not representative for the motions within the proteins. This is especially true for the most mobile residues. On a nanosecond time-scale, the B-factors calculated from the simulations of the proteins in solution are considerably larger than those obtained by structure refinement of the proteins in crystals, based on X-ray data. The time evolution of the atomic fluctuations shows that for large portions of the proteins under study, atomic positional fluctuations are not yet converged after a nanosecond. Cross-correlations do not converge faster than the fluctuations themselves. Most display very erratic behavior if the sampling covers less than about 200 ps. It is also shown that inclusion of mobile atoms into the procedure used to remove rigid-body motion from the simulation can lead to spurious correlations between the motions of the atoms at the surface of the protein.

[1]  B. Lee,et al.  The interpretation of protein structures: estimation of static accessibility. , 1971, Journal of molecular biology.

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

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

[4]  Hans Frauenfelder,et al.  Temperature-dependent X-ray diffraction as a probe of protein structural dynamics , 1979, Nature.

[5]  M. Karplus,et al.  Method for estimating the configurational entropy of macromolecules , 1981 .

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

[7]  M Karplus,et al.  Time dependence of atomic fluctuations in proteins: analysis of local and collective motions in bovine pancreatic trypsin inhibitor. , 1982, Biochemistry.

[8]  Functional significance of flexibility in proteins , 1983, Biopolymers.

[9]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

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

[11]  M. Levitt,et al.  Protein normal-mode dynamics: trypsin inhibitor, crambin, ribonuclease and lysozyme. , 1985, Journal of molecular biology.

[12]  M Karplus,et al.  Active site dynamics of ribonuclease. , 1985, Proceedings of the National Academy of Sciences of the United States of America.

[13]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[14]  SRLSQ refinement of triclinic lysozyme , 1987 .

[15]  A Wlodawer,et al.  Comparison of two highly refined structures of bovine pancreatic trypsin inhibitor. , 1987, Journal of molecular biology.

[16]  R. Sharon,et al.  Accurate simulation of protein dynamics in solution. , 1988, Proceedings of the National Academy of Sciences of the United States of America.

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

[18]  M. Karplus,et al.  Proteins: A Theoretical Perspective of Dynamics, Structure, and Thermodynamics , 1988 .

[19]  A. Cooper.  Dynamics of Proteins and Nucleic Acids , 1988 .

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

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

[22]  Y. Komeiji,et al.  Molecular dynamics simulation of trp-aporepressor in a solvent. , 1991, Protein engineering.

[23]  D. Caspar,et al.  Plasticity of crystalline proteins , 1991 .

[24]  M. Karplus,et al.  Collective motions in proteins: A covariance analysis of atomic fluctuations in molecular dynamics and normal mode simulations , 1991, Proteins.

[25]  S. Swaminathan,et al.  Investigation of domain structure in proteins via molecular dynamics simulation: application to HIV-1 protease dimer , 1991 .

[26]  D. Nguyen,et al.  Conformational Variability of Insulin: a Molecular Dynamics Analysis , 1991 .

[27]  K. P. Murphy,et al.  Molecular basis of co-operativity in protein folding. , 1992, Journal of molecular biology.

[28]  R M Sweet,et al.  Correlations of atomic movements in lysozyme crystals , 1992, Proteins.

[29]  K. P. Murphy,et al.  Molecular basis of co-operativity in protein folding. III. Structural identification of cooperative folding units and folding intermediates. , 1992, Journal of molecular biology.

[30]  S. Swaminathan,et al.  Molecular dynamics of HIV‐1 protease , 1992, Proteins.

[31]  N Go,et al.  Normal mode refinement: crystallographic refinement of protein dynamic structure. II. Application to human lysozyme. , 1992, Journal of molecular biology.

[32]  N Go,et al.  Normal mode refinement: crystallographic refinement of protein dynamic structure. I. Theory and test by simulated diffraction data. , 1992, Journal of molecular biology.

[33]  N. Go,et al.  Effect of solvent on collective motions in globular protein. , 1993, Journal of molecular biology.

[34]  M. Levitt,et al.  Protein unfolding pathways explored through molecular dynamics simulations. , 1993, Journal of molecular biology.

[35]  A. Elofsson,et al.  How consistent are molecular dynamics simulations? Comparing structure and dynamics in reduced and oxidized Escherichia coli thioredoxin. , 1993, Journal of molecular biology.

[36]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[37]  J. Schlitter.  Estimation of absolute and relative entropies of macromolecules using the covariance matrix , 1993 .

[38]  Wilfred F. van Gunsteren,et al.  CONVERGENCE PROPERTIES OF FREE ENERGY CALCULATIONS : ALPHA -CYCLODEXTRIN COMPLEXES AS A CASE STUDY , 1994 .

[39]  O Jardetzky,et al.  Protein dynamics. , 1994, FEBS letters.

[40]  D. Case.  Normal mode analysis of protein dynamics , 1994 .

[41]  Y. Komeiji,et al.  Molecular dynamics simulations of trp apo‐and holorepressors: Domain structure and ligand–protein interaction , 1994, Proteins.

[42]  A. Mark,et al.  Investigation of shape variations in the antibody binding site by molecular dynamics computer simulation. , 1994, Journal of molecular biology.

[43]  L. Ptaszek,et al.  Molecular dynamics studies of the human CD4 protein , 1994, Biopolymers.

[44]  C. Dobson,et al.  Comparison of MD simulations and NMR experiments for hen lysozyme. Analysis of local fluctuations, cooperative motions, and global changes. , 1995, Biochemistry.

[45]  Alan E. Mark,et al.  Computer simulation of protein motion , 1995 .

[46]  C. Schiffer,et al.  Time-averaging crystallographic refinement: possibilities and limitations using alpha-cyclodextrin as a test system. , 1995, Acta crystallographica. Section D, Biological crystallography.