Characterization of the internal motions of Escherichia coli ribonuclease HI by a combination of 15N-NMR relaxation analysis and molecular dynamics simulation: examination of dynamic models.

The backbone dynamics of Escherichia coli ribonuclease HI (RNase HI) in the picosecond to nanosecond time scale were characterized by a combination of measurements of 15N-NMR relaxation (T1, T2, and NOE), analyzed by a model-free approach, and molecular dynamics (MD) simulation in water. The MD simulations in water were carried out with long-range Coulomb interactions to avoid the artificial fluctuation caused by the cutoff approximation. The model-free analysis of the 15N-NMR relaxation indicated that RNase HI has a rotational correlation time of 10.9 ns at 27 degrees C. The generalized order parameter (S2) for the internal motions varied from 0.15 to 1.0, with an average value of 0.85, which is much larger than that of the RNase H domain of HIV-1 reverse transcriptase (0.78). Large internal motions (small order parameters) were observed in the N-terminal region (Leu2-Lys3), the loop between beta-strands A and B (Cys13-Gly15), the turn between alpha-helix I and beta-strand D (Glu61, His62), the loop between beta-strand D and alpha-helix II (Asp70-Tyr71), the loop between alpha-helices III and IV (Ala93-Lys96), the loop between beta-strand E and alpha-helix V (Gly123-His127), and the C-terminal region (Gln152-Val155). The effective correlation time observed in these regions varied from 0.45 ns (Glu61, Lys96) to 2.2 ns (Leu14). The order parameters calculated from the MD agreed well with those from the NMR experiment, with a few exceptions. The distributions of most of the backbone N-H vectors obtained by MD are approximately consistent with the diffusion-in-a-cone model. These distributions, however, were elliptic, with a long axis perpendicular to the plane defined by the N-H and N-C alpha vectors. Distributions supporting the axial fluctuation model or the jump-between-two-cones model were also observed in the MD simulation.

[1]  Kenneth Levenberg.  A METHOD FOR THE SOLUTION OF CERTAIN NON – LINEAR PROBLEMS IN LEAST SQUARES , 1944 .

[2]  L. Kay,et al.  Pulse sequences for removal of the effects of cross correlation between dipolar and chemical-shift anisotropy relaxation mechanisms on the measurement of heteronuclear T1 and T2 values in proteins , 1992 .

[3]  J. Keeler,et al.  The solution structure and dynamics of the DNA-binding domain of HMG-D from Drosophila melanogaster. , 1994, Structure.

[4]  G. Wagner,et al.  Measurement of 13C spin-spin relaxation times by two-dimensional heteronuclear 1H13C correlation spectroscopy , 1989 .

[5]  A. Gronenborn,et al.  Analysis of the backbone dynamics of interleukin-1.beta. using two-dimensional inverse detected heteronuclear nitrogen-15-proton NMR spectroscopy , 1990 .

[6]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[7]  M L Connolly,et al.  The molecular surface package. , 1993, Journal of molecular graphics.

[8]  Minoru Saito,et al.  Molecular dynamics simulations of proteins in water without the truncation of long-range Coulomb interactions , 1992 .

[9]  D. Fushman,et al.  Determination of the backbone mobility of ribonuclease T1 and its 2'GMP complex using molecular dynamics simulations and NMR relaxation data. , 1994, Journal of biomolecular structure & dynamics.

[10]  Y. Satow,et al.  Structure of ribonuclease H phased at 2 A resolution by MAD analysis of the selenomethionyl protein. , 1990, Science.

[11]  M. Oobatake,et al.  pH-dependent thermostabilization of Escherichia coli ribonuclease HI by histidine to alanine substitutions. , 1993, Journal of biotechnology.

[12]  U. Singh,et al.  A NEW FORCE FIELD FOR MOLECULAR MECHANICAL SIMULATION OF NUCLEIC ACIDS AND PROTEINS , 1984 .

[13]  D. Woessner.  Spin Relaxation Processes in a Two‐Proton System Undergoing Anisotropic Reorientation , 1962 .

[14]  Haruki Nakamura,et al.  Presto(protein Engineering Simulator): A Vectorized Molecular Mechanics Program for Biopolymers , 1992, Comput. Chem..

[15]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[16]  L. Nilsson,et al.  A comparison of 15N NMR relaxation measurements with a molecular dynamics simulation: Backbone dynamics of the glucocorticoid receptor DNA‐binding domain , 1993, Proteins.

[17]  S. Kanaya,et al.  DNA sequence of the gene coding for Escherichia coli ribonuclease H. , 1983, The Journal of biological chemistry.

[18]  K. Morikawa,et al.  Three-dimensional structure of ribonuclease H from E. coli , 1990, Nature.

[19]  D. Fushman,et al.  Backbone dynamics of ribonuclease T1 and its complex with 2′GMP studied by two-dimensional heteronuclear NMR spectroscopy , 1994, Journal of biomolecular NMR.

[20]  M. Yoshida,et al.  Complete assignments of magnetic resonances of ribonuclease H from Escherichia coli by double- and triple-resonance 2D and 3D NMR spectroscopies. , 1993, Biochemistry.

[21]  K. Sharp,et al.  How much is a stabilizing bond worth? , 1994, Trends in biochemical sciences.

[22]  L. Kay,et al.  Backbone dynamics of proteins as studied by 15N inverse detected heteronuclear NMR spectroscopy: application to staphylococcal nuclease. , 1989, Biochemistry.

[23]  M Ikehara,et al.  Importance of the positive charge cluster in Escherichia coli ribonuclease HI for the effective binding of the substrate. , 1991, The Journal of biological chemistry.

[24]  M Itaya,et al.  A novel strategy for stabilization of Escherichia coli ribonuclease HI involving a screen for an intragenic suppressor of carboxyl-terminal deletions. , 1994, The Journal of biological chemistry.

[25]  A. Szabó,et al.  Effect of librational motion on fluorescence depolarization and nuclear magnetic resonance relaxation in macromolecules and membranes. , 1980, Biophysical journal.

[26]  H. Nakamura,et al.  Assignments of backbone 1H, 13C, and 15N resonances and secondary structure of ribonuclease H from Escherichia coli by heteronuclear three-dimensional NMR spectroscopy. , 1991, Biochemistry.

[27]  K. Kinosita,et al.  A theory of fluorescence polarization decay in membranes. , 1977, Biophysical journal.

[28]  A. Szabó,et al.  Theory of NMR relaxation in macromolecules: Restricted diffusion and jump models for multiple internal rotations in amino acid side chains , 1978 .

[29]  R. Levy,et al.  Protein dynamics and NMR relaxation: comparison of simulations with experiment , 1982, Nature.

[30]  A. Szabó,et al.  Model-free approach to the interpretation of nuclear magnetic resonance relaxation in macromolecules. 1. Theory and range of validity , 1982 .

[31]  H. Carr,et al.  The Principles of Nuclear Magnetism , 1961 .

[32]  B. Brooks,et al.  A 500 ps molecular dynamics simulation study of interleukin-1 beta in water. Correlation with nuclear magnetic resonance spectroscopy and crystallography. , 1992, Journal of molecular biology.

[33]  J. H. R. Clarke,et al.  A comparison of constant energy, constant temperature and constant pressure ensembles in molecular dynamics simulations of atomic liquids , 1984 .

[34]  G. Wagner,et al.  Measurement of 13C relaxation times in proteins by two-dimensional heteronuclear 1H-13C correlation spectroscopy , 1988 .

[35]  K Morikawa,et al.  Structural details of ribonuclease H from Escherichia coli as refined to an atomic resolution. , 1992, Journal of molecular biology.

[36]  Minoru Saito,et al.  Molecular dynamics simulations of proteins in solution: Artifacts caused by the cutoff approximation , 1994 .

[37]  D. Torchia,et al.  Determination of 15N chemical shift tensor via 15N-2H dipolar coupling in Boc-glycylglycyl[15N]glycine benzyl ester , 1988 .

[38]  K Wüthrich,et al.  Observation of internal motility of proteins by nuclear magnetic resonance in solution. , 1986, Methods in enzymology.

[39]  Ad Bax,et al.  Rapid recording of 2D NMR spectra without phase cycling. Application to the study of hydrogen exchange in proteins , 1989 .

[40]  P. Wright,et al.  NMR Order Parameters of Biomolecules: A New Analytical Representation and Application to the Gaussian Axial Fluctuation Model , 1994 .

[41]  P Argos,et al.  Intramolecular cavities in globular proteins. , 1994, Protein engineering.

[42]  S. Kanaya,et al.  Role of histidine 124 in the catalytic function of ribonuclease HI from Escherichia coli. , 1993, The Journal of biological chemistry.

[43]  R. Powers,et al.  Analysis of the backbone dynamics of the ribonuclease H domain of the human immunodeficiency virus reverse transcriptase using 15N relaxation measurements. , 1992, Biochemistry.

[44]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[45]  A. J. Shaka,et al.  Computer-optimized decoupling scheme for wideband applications and low-level operation , 1985 .

[46]  M. L. Connolly.  Solvent-accessible surfaces of proteins and nucleic acids. , 1983, Science.

[47]  P. Wright,et al.  Intramolecular motions of a zinc finger DNA-binding domain from Xfin characterized by proton-detected natural abundance carbon-13 heteronuclear NMR spectroscopy , 1991 .

[48]  Comparison of the 13C relaxation times and proton scalar couplings of BPTI with values predicted by molecular dynamics. , 1994, Journal of magnetic resonance. Series B.