A Comparison Between Elastic Network Interpolation and MD Simulation of 16S Ribosomal RNA

Abstract In this paper a coarse-grained method called elastic network interpolation (ENI) is used to generate feasible transition pathways between two given conformations of the core central domain of 16S Ribosomal RNA (16S rRNA). The two given conformations are the extremes generated by a molecular dynamics (MD) simulation, which differ from each other by 10Å in root-mean-square deviation (RMSD). It takes only several hours to build an ENI pathway on a 1.5GHz Pentium with 512 MB memory, while the MD takes several weeks on high-performance multi-processor servers such as the SGI ORIGIN 2000/2100. It is shown that multiple ENI pathways capture the essential anharmonic motions of millions of timesteps in a particular MD simulation. A coarse-grained normal mode analysis (NMA) is performed on each intermediate ENI conformation, and the lowest 1% of the normal modes (representing about 40 degrees of freedom (DOF)) are used to parameterize fluctuations. This combined ENI/NMA method captures all intermediate conformations in the MD run with 1.5Å RMSD on average. In addition, if we restrict attention to the time interval of the MD run between the two extreme conformations, the RMSD between the closest ENI/NMA pathway and the MD results is about 1Å. These results may serve as a paradigm for reducedDOF dynamic simulations of large biological macromolecules as well as a method for the reduced-parameter interpretation of massive amounts of MD data.

[1]  B. Ma,et al.  Binding interactions between the core central domain of 16S rRNA and the ribosomal protein S15 determined by molecular dynamics simulations. , 2003, Nucleic acids research.

[2]  J. Williamson,et al.  Structure of the S15,S6,S18-rRNA complex: assembly of the 30S ribosome central domain. , 2000, Science.

[3]  G J Kleywegt,et al.  Phi/psi-chology: Ramachandran revisited. , 1996, Structure.

[4]  J R Williamson,et al.  A hierarchy of RNA subdomains in assembly of the central domain of the 30 S ribosomal subunit. , 2000, RNA.

[5]  M. Karplus,et al.  Harmonic dynamics of proteins: normal modes and fluctuations in bovine pancreatic trypsin inhibitor. , 1983, Proceedings of the National Academy of Sciences of the United States of America.

[6]  H. G. Petersen.  Accuracy and efficiency of the particle mesh Ewald method , 1995 .

[7]  X. Zhuang,et al.  Ligand-induced conformational changes observed in single RNA molecules. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[8]  R L Jernigan,et al.  Collective motions in HIV-1 reverse transcriptase: examination of flexibility and enzyme function. , 1999, Journal of molecular biology.

[9]  D. Ming,et al.  How to describe protein motion without amino acid sequence and atomic coordinates , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[10]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[11]  G J Kleywegt,et al.  Where freedom is given, liberties are taken. , 1995, Structure.

[12]  Bruce A. Shapiro,et al.  Molecular Dynamics Simulations of the Denaturation and Refolding of an RNA Tetraloop , 2001, Journal of biomolecular structure & dynamics.

[13]  W. Wriggers,et al.  Exploring global distortions of biological macromolecules and assemblies from low-resolution structural information and elastic network theory. , 2002, Journal of molecular biology.

[14]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[15]  R. Jernigan,et al.  Anisotropy of fluctuation dynamics of proteins with an elastic network model. , 2001, Biophysical journal.

[16]  G. Chirikjian,et al.  Elastic models of conformational transitions in macromolecules. , 2002, Journal of molecular graphics & modelling.

[17]  R. Batey,et al.  Interaction of the Bacillus stearothermophilus ribosomal protein S15 with 16 S rRNA: II. Specificity determinants of RNA-protein recognition. , 1996, Journal of molecular biology.

[18]  Peter A. Kollman,et al.  AMBER: Assisted model building with energy refinement. A general program for modeling molecules and their interactions , 1981 .

[19]  Daniel ben-Avraham,et al.  Normal modes analyses of macromolecules , 1998 .

[20]  G. Chirikjian,et al.  Efficient generation of feasible pathways for protein conformational transitions. , 2002, Biophysical journal.

[21]  W. Delano.  The PyMOL Molecular Graphics System , 2002 .

[22]  R L Jernigan,et al.  Vibrational dynamics of transfer RNAs: comparison of the free and synthetase-bound forms. , 1998, Journal of molecular biology.

[23]  G. Schulz,et al.  Movie of the structural changes during a catalytic cycle of nucleoside monophosphate kinases. , 1995, Structure.