ID code: 1AAY, PDT039 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 10 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 404 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 G02C1' 30 C03C1' 30 C03C1' 49 G04C1' 49 G04C1' 71 T05C1' 71 T05C1' 91 G06C1' 91 G06C1' 113 G07C1' 113 G07C1' 135 G08C1' 135 G08C1' 157 C09C1' 157 C09C1' 176 G10C1' 176 G10C1' 198 T11C1' 236 C53C1' 215 A52C1' 255 G54C1' 236 C53C1' 277 C55C1' 255 G54C1' 296 C56C1' 277 C55C1' 315 C57C1' 296 C56C1' 334 A58C1' 315 C57C1' 355 C59C1' 334 A58C1' 374 G60C1' 355 C59C1' 396 C61C1' 374 G60C1' 9 G02N9 31 C03N1 31 C03N1 50 G04N9 50 G04N9 72 T05N1 72 T05N1 92 G06N9 92 G06N9 114 G07N9 114 G07N9 136 G08N9 136 G08N9 158 C09N1 158 C09N1 177 G10N9 177 G10N9 199 T11N1 237 C53N1 216 A52N9 256 G54N9 237 C53N1 278 C55N1 256 G54N9 297 C56N1 278 C55N1 316 C57N1 297 C56N1 335 A58N9 316 C57N1 356 C59N1 335 A58N9 375 G60N9 356 C59N1 397 C61N1 375 G60N9 ************************************** IN ORTHONORMAL COORDINATES, HELIX AXIS IN PARAMETRIC FORM: X = 0.60786*S + -6.21428 Y = 0.79299*S + 2.31983 Z = 0.04094*S + 47.33630 IN ORIGINAL CRYSTAL COORDINATES, HELIX AXIS IN PARAMETRIC FORM: X = 0.60786*S + -6.21428 Y = 0.79299*S + 2.31983 Z = 0.04094*S + 47.33630 >>>>>> HELIX ROTATION: 32.105 DISPLACEMENT: 3.2741 STATISTICS: OVERALL STANDARD DEV.: 1.0652 SIGMA(X): 1.5884, SIGMA(Y): 1.6752, SIGMA(Z)=SIGMA(DISPLACEMENT): 0.3715, SIGMA(ROTATION): 7.143 THERE ARE 11.21 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 = 10 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.10130 -0.06329 0.99284 95.81 93.63 6.86 -0.23808 -0.20066 0.95029 103.77 101.58 18.14 -0.00932 -0.13215 0.99119 90.53 97.59 7.61 0.08558 -0.19143 0.97777 85.09 101.04 12.10 0.11024 -0.14535 0.98322 83.67 98.36 10.51 0.11083 -0.07360 0.99111 83.64 94.22 7.65 0.10877 0.05334 0.99263 83.76 86.94 6.96 0.19631 0.17790 0.96427 78.68 79.75 15.36 0.09597 0.16619 0.98141 84.49 80.43 11.06 -0.01454 0.24316 0.96988 90.83 75.93 14.10 STRAND 2 BASE NORMAL COSINES AND ANGLES COS(AX) COS(AY) COS(AZ) ANG X ANG Y ANG Z -0.23409 -0.13193 0.96322 103.54 97.58 15.59 -0.03335 -0.18065 0.98298 91.91 100.41 10.59 0.08699 -0.18221 0.97940 85.01 100.50 11.65 0.19101 -0.15772 0.96883 78.99 99.07 14.34 0.17222 0.02401 0.98477 80.08 88.62 10.01 0.19478 0.07402 0.97805 78.77 85.76 12.03 0.06845 0.28147 0.95713 86.08 73.65 16.84 0.08675 0.17920 0.97998 85.02 79.68 11.48 -0.04351 0.10051 0.99398 92.49 84.23 6.29 -0.04119 0.08280 0.99571 92.36 85.25 5.31 BASE PAIR NORMAL COSINES AND ANGLES COS(AX) COS(AY) COS(AZ) ANG X ANG Y ANG Z -0.13337 -0.08993 0.98698 97.66 95.16 9.26 -0.10023 -0.16797 0.98068 95.75 99.67 11.28 0.03578 -0.16759 0.98521 87.95 99.65 9.87 0.12515 -0.17384 0.97679 82.81 100.01 12.37 0.13498 -0.07558 0.98796 82.24 94.33 8.90 0.13133 -0.00834 0.99130 82.45 90.48 7.56 0.07644 0.14446 0.98655 85.62 81.69 9.41 0.11058 0.15304 0.98201 83.65 81.20 10.88 0.04439 0.12411 0.99127 87.46 82.87 7.57 0.01774 0.12871 0.99152 88.98 82.60 7.47 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 I= 7 PS= 90.0000 LS= 90.0000 PL= 90.0000 I= 8 PS= 90.0000 LS= 90.0000 PL= 90.0000 I= 9 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 8 9 10 I= 1 0 4 10 15 15 15 18 19 15 15 I= 2 4 0 7 12 14 16 20 22 18 18 I= 3 10 7 0 5 7 10 18 18 16 17 I= 4 15 12 5 0 5 9 18 18 17 18 I= 5 15 14 7 5 0 3 13 13 12 13 I= 6 15 16 10 9 3 0 9 9 9 10 I= 7 18 20 18 18 13 9 0 2 2 3 I= 8 19 22 18 18 13 9 2 0 4 5 I= 9 15 18 16 17 12 9 2 4 0 1 I= 10 15 18 17 18 13 10 3 5 1 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 2.09 6.53 3.48 10.61 0.00 -0.00 4.04 15.03 7.72 -9.03 0.00 -0.00 0.59 18.13 9.09-10.31 0.00 0.00 6.97 7.93 5.91 -3.57 0.00 0.00 0.99 7.55 5.18 3.80 0.00****** 6.99 9.27 6.06 -0.95 -0.00 0.00 1.65 11.99 2.59 -1.50 0.00 0.00 8.01 11.82 4.58 -9.46 0.00 0.00 -2.66 10.16 3.82 -1.54 0.00 0.00 7.19 6.94 3.14 0.10 -0.00 0.00 -3.39 6.85 7.18 1.25 0.00 0.00 5.62 10.60 13.15 -4.89 0.00 0.00 -0.21 6.96 1.89 8.55****** 0.00 12.53 11.07 -5.04 -3.18 0.00 0.00 -3.85 14.85 5.37 -2.18 0.00 0.00 1.45 11.39 6.00 -6.35 0.00 0.00 -5.91 9.32 4.59 6.22 -0.00 0.00 2.25 5.87 0.20 -1.00 0.00 0.50 -7.10 12.12 0.00 0.00 0.00 0.00 -0.69 5.26 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 4.87 0.33 4.86 -0.23 23.69 3.68 0.38 8.71 3.12 0.33 4.85 -0.14 -18.63 -1.96 8.52 -2.59 0.14 7.80 -4.08 6.65 -0.32 39.36 2.96 -0.84 10.53 4.01 -4.14 6.62 -0.24 11.91 -6.38 -10.12 -1.88 0.25 5.16 0.64 5.12 -0.46 27.42 3.54 0.47 9.03 3.94 0.74 5.09 -0.41 -1.85 -6.00 1.79 -2.17 0.87 5.70 -4.39 3.63 -0.25 34.97 3.45 -0.97 11.65 4.11 -4.49 3.45 -0.19 -3.10 -6.37 -0.06 -1.05 0.58 3.86 -1.69 3.48 -0.93 30.51 3.23 -0.62 8.79 1.36 -1.67 3.48 -0.88 6.95 -9.85 -3.16 -1.58 0.52 9.32 -1.47 9.20 -0.40 34.66 3.40 -0.07 7.55 0.47 -1.45 9.23 -0.26 0.54 -9.01 3.78 -2.01 0.09 2.03 1.80 -0.96 -0.27 28.66 3.46 0.73 7.76 5.28 1.78 -0.96 -0.22 -7.81 -12.75 4.32 -1.75 0.45 4.17 -2.65 3.22 0.18 35.18 3.38 -1.03 10.87 -0.51 -2.61 3.21 0.25 -0.00 -5.31 -3.48 -0.66 0.20 1.55 0.72 1.37 -0.48 31.24 3.45 0.90 7.01 -2.84 0.74 1.36 -0.43 -3.42 -8.16 -3.49 -1.63 0.58 - - - - - - - 5.56 -4.96 - - - - -6.45 -6.91 -0.69 0.33 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.