H-bond stability in the tRNA(Asp) anticodon hairpin: 3 ns of multiple molecular dynamics simulations

Multiple molecular dynamics trajectories of the solvated and neutralized 17-residue tRNA(Asp) anticodon hairpin were generated for a total of 3 ns. Explicit treatment of all long-ranged electrostatic interactions by the particle mesh Ewald algorithm, as implemented in the AMBER MD software package, effected a degree of structural stabilization not previously achieved by use of a long 16-A solvent interaction truncation scheme. The increased stability of this multiple molecular dynamics set was appropriate for an in-depth analysis of the six 500-ps-long trajectories and allowed the characterization of a number of key structural interactions. The dynamical behavior of the standard Watson-Crick base pairs, the noncanonical G30-U40 "wobble" base pair, and the psi 32-C38 pseudo-base pair is presented as well as that of two C--H... O hydrogen bonds found to contribute to the array of tertiary interactions that stabilize the seven-nucleotide native loop conformation. The least mobile residue in the loop is U33, which forms the U-turn motif and which participates in several hydrogen-bonding interactions, whereas the most mobile residue is the apical residue G34 at the wobble position, a factor undoubtedly important in its biological function. The set of multiple molecular dynamics trajectories obtained does not converge on a 500-ps time scale to a unique dynamical model but instead describes an ensemble of structural microstates accessible to the system under the present simulation protocol, which is the result of local structural heterogeneity rather than of global conformational changes.