Methods: Atomic Coordinate Files.
The 19 base-pair DNA segment in 1D66.PDB was chosen because the helix is nearly straight. The protein and water were deleted from the PDB file which speeded up the rotation frames/second in Chime by about 6-fold. The end base pairs were remarked-out so as to produce 3' hydroxyl's and 5' phosphates on the ends, leaving 17 base pairs in 1D66-PWZ.PDB. All hydrogens are missing since the atomic coordinates in 1D66.PDB were determined by X-ray diffraction at a resolution (2.7 Å) unable to resolve hydrogen atoms.
The AT and GC base pairs are residues A27 T12 and C33 G6 taken from 1D66.PDB. The free program WebLab (which like Chime is also derived in part from RasMol's source code) from Molecular Simulations Inc. was used to extract each base pair into a separate atomic coordinate file. WebLab allows atom-by-atom mouse picking as well as selecting by named residue. Each base pair was saved in a MOL file because this, unlike WebLab's PDB-saving, preserved double bond information. Some manual correction of bond orders was required. WebLab was then used to add hydrogen atoms (a few inappropriately added were removed). In the case of the GC, the hydrogens involved in hydrogen bonding were not well oriented towards their bonded oxygens/nitrogens, so WebLab's bond torsion capability was used to position these hydrogens in a more favorable orientation. Both real (from H to N/O) and simplified hydrogen bonds (from N to N/O) were added manually to the bond list.
(Another approach would be to use theoretically-generated atomic coordinate files. These can be generated for any specified sequence, in A- or B-form DNA or A-form RNA, at Pittsburgh Supercomputer Center Web Tools. Although I did not know about this resource in 1996 when I implemented my nonlinear tutorial on DNA, I like the "real world" irregularities in the empirical 1D66 structure.)
Acknowledgements.
Thanks to
MDLI
for donating
Chime
to the world. Thanks to
Tim Maffett
at MDLI for making script-enhanced Chime available during
its development. Thanks to
Roger Sayle
who pioneered the molecular
structure rendering methods in RasMol, and whose source code was used in
the development of Chime.
Professor Anne Simon suggested the page on ends and antiparallelism.
Dr. Frieda Reichsman (
moleculesinmotion.com) made several good suggestions I adopted in
version 2.8, and caught several bugs I missed!
The unique
"MoveTo" utility server for Chime by
Todd A. Carlson and Clark Wells
was invaluable
in scripting the zooms from the double helix to single base pairs, and back.
References.
The structure 1D66.PDB, from which the DNA coordinates
were excerpted, was determined by:
This page is maintained by
Eric Martz.