ID code: 1ZQI, PDE035 C ********************************************************************* C FREEHEL98.FOR 1 June 1998 C This is the Fortran program code for helix analysis program FREEHELIX. C For operating instructions, see the separate file FREEHEL98.TEX. This C program is radically changed from the earlier NEWHELIX, deleting the C old RADJ and TADJ, and adding the helix-independent parameters VALL C (total bending), VROL (roll), VTIL (tilt), VTWI (twist), VSLI (slide), C VRIS (rise) and VSHF (shift). This program can be used with DNA C helices of any degree of bending. C Richard E. Dickerson, Oxford C ********************************************************************* TITL CELL 1., 1., 1., 90., 90., 90. BRKH FPUN 0 PMIN 0 PMAX 0 BASE 7 HELX RC1' YC1' HELX RN9 YN1 BROL CYLN TRNG END 1 INPUT FRAC.COORDS. X,Y,Z AND ATOM NAME IN FORMAT: (T31,3F8.3,T20,A1,T25,2A1,T14,3A1) WILL READ 281 ATOMS: PRINT FLAG= 0 CELL CONSTANTS A,B,C,ALPHA,BETA,GAMMA ARE: 1.000 1.000 1.000 90.000 90.000 90.000 COORDINATES IN ORTHONORMAL SYSTEM: PRINT FLAG= 0 ************************************** ATOM PAIRS USED TO DETERMINE THE HELIX AXIS ARE: 8 C01C1' 27 A02C1' 27 A02C1' 48 T03C1' 48 T03C1' 68 T04C1' 68 T04C1' 88 A05C1' 88 A05C1' 109 G06C1' 109 G06C1' 131 A07C1' 169 C02C1' 149 T01C1' 188 T03C1' 169 C02C1' 208 A04C1' 188 T03C1' 229 A05C1' 208 A04C1' 250 T06C1' 229 A05C1' 270 G07C1' 250 T06C1' 9 C01N1 28 A02N9 28 A02N9 49 T03N1 49 T03N1 69 T04N1 69 T04N1 89 A05N9 89 A05N9 110 G06N9 110 G06N9 132 A07N9 170 C02N1 150 T01N1 189 T03N1 170 C02N1 209 A04N9 189 T03N1 230 A05N9 209 A04N9 251 T06N1 230 A05N9 271 G07N9 251 T06N1 ************************************** IN ORTHONORMAL COORDINATES, HELIX AXIS IN PARAMETRIC FORM: X = 0.56912*S + 20.68143 Y = 0.81147*S + -20.69543 Z = 0.13277*S + 37.83588 IN ORIGINAL CRYSTAL COORDINATES, HELIX AXIS IN PARAMETRIC FORM: X = 0.56912*S + 20.68143 Y = 0.81147*S + -20.69543 Z = 0.13277*S + 37.83588 >>>>>> HELIX ROTATION: 32.209 DISPLACEMENT: 3.4408 STATISTICS: OVERALL STANDARD DEV.: 0.6930 SIGMA(X): 0.8589, SIGMA(Y): 0.8462, SIGMA(Z)=SIGMA(DISPLACEMENT): 0.2281, SIGMA(ROTATION): 6.035 THERE ARE 11.18 RESIDUES PER TURN ***************************************************** THE FOLLOWING DIAMOND LIST HAS BEEN OUTPUT ON UNIT: 12 NORMAL END OF JOB: YOU HAVE GIVEN BIRTH TO A HELIX NUMBER BASE PAIRS = 7 1 ROLL+TILT OUTPUT, FREEHELIX98 STRAND 1 BASE NORMAL COSINES AND ANGLES COS(AX) COS(AY) COS(AZ) ANG X ANG Y ANG Z -0.21037 0.22301 0.95185 102.14 77.11 17.85 -0.13346 0.08501 0.98740 97.67 85.12 9.10 -0.15032 0.10442 0.98311 98.65 84.01 10.55 -0.13869 -0.12173 0.98283 97.97 96.99 10.63 -0.01195 -0.16247 0.98664 90.68 99.35 9.38 0.03576 -0.22098 0.97462 87.95 102.77 12.94 0.18308 -0.17992 0.96649 79.45 100.37 14.87 STRAND 2 BASE NORMAL COSINES AND ANGLES COS(AX) COS(AY) COS(AZ) ANG X ANG Y ANG Z -0.18155 0.24333 0.95280 100.46 75.92 17.67 -0.22306 -0.04416 0.97380 102.89 92.53 13.14 -0.13504 -0.16938 0.97626 97.76 99.75 12.51 0.04200 -0.19236 0.98042 87.59 101.09 11.36 0.14725 -0.22037 0.96424 81.53 102.73 15.37 0.18729 -0.05547 0.98074 79.21 93.18 11.26 0.04471 0.00316 0.99899 87.44 89.82 2.57 BASE PAIR NORMAL COSINES AND ANGLES COS(AX) COS(AY) COS(AZ) ANG X ANG Y ANG Z -0.10149 0.10365 0.98942 95.83 84.05 8.34 -0.18605 0.00223 0.98254 100.72 89.87 10.72 -0.18370 -0.02597 0.98264 100.59 91.49 10.69 -0.06531 -0.14258 0.98763 93.74 98.20 9.02 0.08635 -0.17451 0.98086 85.05 100.05 11.23 0.09351 -0.15763 0.98306 84.63 99.07 10.56 0.11287 -0.10256 0.98830 83.52 95.89 8.77 1 I= 1 PS= 90.0000 LS= 90.0000 PL= 90.0000 I= 2 PS= 90.0000 LS= 90.0000 PL= 90.0000 I= 3 PS= 90.0000 LS= 90.0000 PL= 90.0000 I= 4 PS= 90.0000 LS= 90.0000 PL= 90.0000 I= 5 PS= 90.0000 LS= 90.0000 PL= 90.0000 I= 6 PS= 90.0000 LS= 90.0000 PL= 90.0000 1 ANGLES BETWEEN NORMAL VECTORS TO BASE PAIRS J= 1 2 3 4 5 6 7 I= 1 0 7 8 14 19 18 17 I= 2 7 0 1 10 18 18 18 I= 3 8 1 0 9 17 17 17 I= 4 14 10 9 0 8 9 10 I= 5 19 18 17 8 0 1 4 I= 6 18 18 17 9 1 0 3 I= 7 17 18 17 10 4 3 0 1 ROLL+TILT OUTPUT, FREEHELIX98 STRAND 1 ROLL AND TILT ANGLES STRAND 2 ROLL AND TILT ANGLES TIP INCL ROLL TILT RADJ TADJ TIP INCL ROLL TILT RADJ TADJ -1.73 17.76 5.51 -7.21 0.00 0.00 -3.70 17.26 16.28 -4.38****** 0.00 -2.38 8.78 -1.09 0.99 0.00 0.00 6.22 11.53 7.06 -5.23 0.00 0.00 -8.50 6.19 9.99 8.37****** 0.00 6.65 10.55 7.90 -6.54 0.00 0.00 -2.86 10.23 7.65 -0.22 0.00 0.00 8.05 7.95 6.23 -0.50 0.00 0.00 -0.54 9.36 1.66 4.00****** 0.00 8.67 12.59 4.90 -8.43 0.00 0.50 -4.25 12.19 7.93 3.77 0.00 0.00 7.93 7.94 -3.96 -7.92 0.00 0.00 -3.43 14.46 0.00 0.00****** 0.00 1.48 2.10 0.00 0.00 0.00 0.00 BEST PLANE THROUGH BOTH BASES VALL VTIL VROL VSLI VTWI VRIS VSHF INCL TIP TILT ROLL SLIDE CUP PROP BUCK X DSP Y DSP 7.58 -1.46 7.44 0.37 23.33 3.55 -1.00 8.31 -0.67 2.13 7.28 0.49 3.29 1.97 -0.49 -0.44 0.20 1.62 -0.11 1.62 -0.21 35.78 3.38 -0.15 10.27 3.04 -0.18 1.61 -0.21 1.51 -8.61 2.80 -1.26 0.37 9.54 -1.26 9.45 -0.36 37.31 3.43 -0.03 10.49 -2.04 -0.42 9.56 -0.27 -6.52 -15.17 4.32 -1.16 0.29 8.90 -0.71 8.87 0.02 31.68 3.32 -0.28 8.92 1.32 -1.17 8.84 0.11 5.55 -10.91 -2.20 -1.08 0.14 1.06 -0.85 0.63 -0.22 28.94 3.51 0.37 9.97 5.11 -0.81 0.67 -0.19 -7.52 -9.21 3.35 -1.07 0.26 3.36 -1.54 2.99 0.68 29.98 3.65 -0.38 10.55 0.39 -1.65 2.91 0.73 -8.19 -12.20 -4.17 -0.50 -0.22 - - - - - - - 8.60 -1.70 - - - - -4.93 -12.37 -0.94 0.24 NOTE: Angles are calculated from 5" end to 3" end of strand 1, and signs of angles also are calculated with respect to strand 1. To examine individual strand 2 bases w.r.t. strand 2, reverse signs of Tip and Tilt. For Z-DNA, reverse signs of Incl and X Dsp. Y Dsp is correct as printed. ROLL and TILT are the simple components of base pair normals along minor and major axes of base pairs. They are the values that were calculated in NEWHEL90 and earlier. VALL, VTIL, VROL, VSLI and VTWI are the total angle betw een base pair normal vectors, and the Tilt, Roll, Slide and Twist calculated relative to a set of local axes halfway between each of the long axes, the short axes, and normal vectors for the two base pairs. They are completely independent of the choice of overall helix axis.