화학공학소재연구정보센터
Journal of the American Chemical Society, Vol.119, No.21, 4805-4825, 1997
Molecular-Dynamics Simulations Highlight the Structural Differences Among DNA-DNA, RNA-RNA, and DNA-RNA Hybrid Duplexes
Nanosecond length simulations applying the particle mesh Ewald method within AMBER 4.1 on canonical A-form and B-form geometries of d[CCAACGTTGG](2), r[CCAACGUUGG](2), and d[CCAACGTTGG]-r[CCAACGUUGG] duplexes in aqueous solution are reported. DNA duplexes only adopt a stable B-DNA geometry, in contrast to RNA duplexes which adopt both astable A-RNA and "B-RNA" geometry. The observation of a stable "B-RNA" structure is somewhat surprising and suggests significant kinetic barriers to structural conversion in RNA structures on a nanosecond time scale. The "B-RNA" can be converted to A-RNA by forcing a concerted flip in the sugar puckers from C2’-endo to C3’-endo. The A-RNA structure displays features similar to A-form crystal structures, specifically interstrand purine stacking at the central pyrimidine-purine step is observed. When started in a canonical A-form geometry, DNA:RNA hybrid duplexes converge to a structure that is characteristic of experimental solution structures; specifically, a minor groove width intermediate between A-form and B-form geometries, the RNA strand in an A-form geometry, a mixture of C2’-endo and C3’-endo sugar puckers in the DNA strand, expected distribution of backbone angles and reasonable agreement with the helicoidal parameters are observed. In all of the simulations reported, A-form geometries appear to be less flexible than B-form geometries. There are also significant differences in the patterns of hydration and counterion association between A-form and B-form duplexes. In A-RNA, sodium counterions tend to associate into "pockets" in the major groove whereas these counterions tend to associate into the minor groove in B-form structures.