Biophys J, August 2001, p. 630-642, Vol. 81, No. 2

Molecular Dynamics Simulation of the Human U2B" Protein Complex with U2 snRNA Hairpin IV in Aqueous Solution

Eppley Institute, University of Nebraska Medical Center, Omaha, Nebraska 68198-6805 USA

A 2200-ps molecular dynamics (MD) simulation of the U2 snRNA hairpin IV/U2B" complex was performed in aqueous solution using the particle mesh Ewald method to consider long-range electrostatic interactions. To investigate the interaction and recognition process between the RNA and protein, the free energy contributions resulting from individual amino acids of the protein component of the RNA/protein complex were calculated using the recently developed glycine-scanning method. The results revealed that the loop region of the U2 snRNA hairpin IV interacted mainly with three regions of the U2B" protein: 1) 1-helix A, 2) 2-3, and 3) 4-helix C. U2 snRNA hairpin IV bound U2B" in a similar orientation as that previously described for U1 snRNA with the U1A' protein; however, the details of the interaction differed in several aspects. In particular, 1-helix A and 4-helix C in U2B" were not observed to interact with RNA in the U1A' protein complex. Most of the polar and charged residues in the interacting regions had larger mutant free energies than the nonpolar residues, indicating that electrostatic interactions were important for stabilizing the RNA/protein complex. The interaction was further stabilized by a network of hydrogen bonds and salt bridges formed between RNA and protein that was maintained throughout the MD trajectory. In addition to the direct interactions between RNA and the protein, solvent-mediated interactions also contributed significantly to complex stability. A detailed analysis of the ordered water molecules in the hydration of the RNA/protein complex revealed that bridged water molecules reside at the interface of RNA and protein as long as 2100 ps in the 2200-ps trajectory. At least 20 bridged water molecules, on average, contributed to the instantaneous stability of the RNA/protein complex. The stabilizing interaction energy due to bridging water molecules was obtained from ab initio Hartree-Fock and density functional theory calculations.

Biophys J, August 2001, p. 630-642, Vol. 81, No. 2
© 2001 by the Biophysical Society   0006-3495/01/08/630/13  $2.00

This article has been cited by other articles:

				M. V. Krasovska, J. Sefcikova, K. Reblova, B. Schneider, N. G. Walter, and J. Sponer Cations and Hydration in Catalytic RNA: Molecular Dynamics of the Hepatitis Delta Virus Ribozyme Biophys. J., July 15, 2006; 91(2): 626 - 638. [Abstract] [Full Text] [PDF]

				V. Cojocaru, R. Klement, and T. M. Jovin Loss of G-A base pairs is insufficient for achieving a large opening of U4 snRNA K-turn motif Nucleic Acids Res., June 13, 2005; 33(10): 3435 - 3446. [Abstract] [Full Text] [PDF]

				F. Razga, J. Koca, J. Sponer, and N. B. Leontis Hinge-Like Motions in RNA Kink-Turns: The Role of the Second A-Minor Motif and Nominally Unpaired Bases Biophys. J., May 1, 2005; 88(5): 3466 - 3485. [Abstract] [Full Text] [PDF]

				K. Reblova, N. Spackova, J. Koca, N. B. Leontis, and J. Sponer Long-Residency Hydration, Cation Binding, and Dynamics of Loop E/Helix IV rRNA-L25 Protein Complex Biophys. J., November 1, 2004; 87(5): 3397 - 3412. [Abstract] [Full Text] [PDF]