home software sectionSeveral years ago, the European Molecular Biology Organization published guidelines defining how nucleic acid duplex structures can be examined. The guidelines defined three rotational and three translational parameters between the bases comprising a base pair as well as from one base pair to the next. There are five programs readily distributed in order to perform the calculations. We are releasing a sixth which was created to be able to handle not only the calculations of normal duplexes, but also comparably examine mispaired bases, hoogstein base pairs, looped out bases, a single base adjacent to a base pair, intercalated drugs for which a coordinate frame has been defined, etc.
The nucleic acid structure analysis program is written in Fortran and runs in a UNIX or VMS environment (please specify preference). The program calculates all of the rotational and translational parameters for complementary base pairs, neighboring base pairs in Cartesian and helical coordinate frames, and base to base Cartesian and helical parameters along each strand. A simple interactive user interface allows for one to specify what file is examined and which parameters to calculate. The program was designed so that the user needs to spend a minimal time reading the documentation in order for the program to run. Full disclosure of the mathematics has been made and published so that the user can readily understand what their parameters mean.
The mathematics are unique for this type of calculation and aviod many of the problems previously encountered. In particular, the calculations are performed by a single rotation not sequential rotations avoiding the use of a midway coordinate frame for rotational parameter calculations. This ensures that the magnitude of the rotational parameters is strand free and direction free. Because a single rotation is used, the equations for each rotational parameter are equivalent, therefore, a 5 degree rotation about the X axis is equivalent to a 5 degree rotation about either the Y or the Z axis. This is not always true of other mathematical formulations presently being used. The mathematics used for calculating the complementary base parameters is identical to the math used to calculate the neighboring base/base pair parameters except that different axes are involved. This is not generally the case of the other available programs. In addition, since the calculations are bases upon local considerations, the value of the parameter depends on only the bases involved and are independent of adjacent bases. This allows for more accurate comparison of parameters from one structure to another.
Copies of the code can be acquired by e-mailing me at:
Please send your comments and your suggestions to Yves Epelboin, epelboin@lmcp.jussieu.fr .