F. Lankas
Institute of Organic Chemistry and Biochemistry
Flemingovo nam. 2
166 10 Praha 6
Czech Republic


F. Lankas, T.E. Cheatham III, N. Spackova, P. Hobza, J. Langowski, J. Sponer
Biophys. J. 2002 82:2592-2609
Critical effect of the N2 amino group on structure, dynamics, and elasticity of DNA polypurine tracts
Unrestrained 5-20-ns explicit-solvent molecular dynamics simulations using the Cornell et al. force field have been carried out for d[GCG(N)11GCG]2 (N, purine base) considering guanine·cytosine (G·C), adenine·thymine (A·T), inosine·5-methyl-cytosine (I·mC), and 2-amino-adenine·thymine (D·T) basepairs. The simulations unambiguously show that the structure and elasticity of N-tracts is primarily determined by the presence of the aminogroup in the minor groove. Simulated A-, I-, and AI-tracts show almost identical structures, with high propeller twist and minor groove narrowing. G- and D-tracts have small propeller twisting and are partly shifted toward the A-form. The elastic properties also differ between the two groups. The sequence-dependent electrostatic component of base stacking seems to play a minor role. Our conclusions are entirely consistent with available experimental data. Nevertheless, the propeller twist and helical twist in the simulated A-tract appear to be underestimated compared to crystallographic studies. To obtain further insight into the possible force field deficiencies, additional multiple simulations have been made for d(A)10, systematically comparing four major force fields currently used in DNA simulations and utilizing B and A-DNA forms as the starting structure. This comparison shows that the conclusions of the present work are not influenced by the force field choice.


This is a contributed parameter set for D-2-AMINO-ADENOSINE with 5'-phosphate group and 3'-O(minus) group, to be used with the AMBER parmbsc0 force field.

The following bonded parameters were missing in the parmbsc0 force field:
CA-NC-CA corresponding to the angle C6-N1-C2
NC-CA-NC corresponding to the angle N1-C2-N3

They were replaced by analogous parmbsc0 parameters:
CA-NC-CQ in analogy with the adenine angle C6-N1-C2
NA-CA-NC in analogy with the guanine angle N1-C2-N3

The charges were computed using a standard RESP procedure.
The parameters were obtained by Filip Lankas and Tomas Drsata.

[ APD Home ] [ Submit Parameters ] [ Amber Home ] [ UK Mirror ] [ Pharmacy Home ] [ Contact ]

Last modified 17 Sep 2012