Calculate global helical-axis, curvatures and tangents from local helical axis.
It can be used to calculate global helical-axis and further curvatures and tangents along it.
The local helical axis coordinates should be present in input HDF5 file. It can be stored in HDF5 file using ‘dnaMD saveasH5’ command with ‘-ha/–helical-axis’ option.
At first global helical axis is calculated by spline interpolation using local helical axis. This method makes the local helical axis smooth and can be used as global helical axis.
After determining the global helical axis, “-ctan/–curv-tangent” option can be used to calculate curvature and tangents at each helical axis position (each base-step). Curvature and tangents data will be stored in same HDF5 file.
This HDF5 file can be used later to extract and calculate data such as bending angle, overall curvature etc.
usage: dnaMD smoothenAxis [-h] [-tbp total-bp-number] [-fbp 1] [-io inout.h5] [-bs bp/s-start-number] [-be b/ps-end-number] [-ctan] [-s 500.0] [-sp 3] [-fp 6] [-cta 20] [-ap helical_axis.pdb] [-apo both] [-scp 1.0]
-h, --help show this help message and exit -tbp total-bp-number, --total-bp total-bp-number Total number of basepair in DNA/RNA. It is an essential input. -fbp 1, --first-bp 1 Basepair number of first base-pair. Usually it is one. Therefore, if this option is not provided, base-pair numbering will start from one. In rare cases, base-pair numbering might start with other number. In those cases, use this option to start numbering of basepair from other number than one. -io inout.h5, --in-out-put inout.h5 Name of input-output file. It is a HDF5 format file generated from saveAsH5 tool. It should contain local helical axis coordinates generated from do_x3dna. It is also the output file as global helical axis, curvature and tangents will be stored/saved in same file. -bs bp/s-start-number, --bp-start bp/s-start-number First BP/BPS of considered DNA segment. If it is not given, first base-step will be considered. -be b/ps-end-number, --bp-end b/ps-end-number Last BP/BPS of considered DNA segment. If it is not given, last base-step will be considered. -ctan, --curv-tangent Calculate curvature and tangents. If this options is used, curvature and tangents will be calculated from global helical axis and stored in same HDF5 file. If this option is not used, curvature and tangents will not be calculated. The tangents can be used later to calculate bending angle. -s 500.0, --smooth 500.0 A smoothing condition for spline interpolation. For more details, see below in the link about "s = None", which is passed into "scipy.interpolate.splprep()" function: http://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.interpolate.splprep.html#scipy.interpolate.splprep NOTE: => Lower value may lead to an artifact of local sharp kink in the smoothed axis. => Higher value may lead to the calculation of wrong helical axis. -sp 3, --spline 3 Degree of spline. For more details, see below in the link about "k = 3", which is passed into "scipy.interpolate.splprep()" function: http://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.interpolate.splprep.html#scipy.interpolate.splprep Its value should be from one to five. -fp 6, --fill-point 6 Number of interpolated points between two adjacent helical-axis coordinates. Its value should be from one to ten. -cta 20, --cut-off-angle 20 Cut-off bending angle to define sharp kink in fitted 3D curve or axis. If angle in fitted curve is larger than this cut-off, refitting will be performed after deleting few of the original helical axis positions. If after this deletions, bending angle will not reduce below cut-off angle, value of ""smooth"" will be increased by 100 and entire cycle of fitting-refitting will be performed. When, value of ""smooth"" increases to more than 10000 during this fitting-refitting cycles, fitting process will be stopped with a warning message. If fitting process is not successful (angle > cutoff-angle), a mask will be added for the respective frame. This mask can be used later to discard these frames during subsequent analysis. -ap helical_axis.pdb, --axis-pdb helical_axis.pdb Output helical-axis as PDB file. This option enable the output of helical axis to a PDB file. Either local or global or both axis can be included in output file. Also, scaled curvature can be included in b-factor field of PDB file. -apo both, --axis-pdb-option both Option for helical-axis PDB file. Following three options can be used to modify the output in PDB file: * 'global': Only global smoothed helical axis will be present. * 'local' : Only local helical axis calculated from do_x3dna will be present. * 'both' : Both global and local helical axis will be present. -scp 1.0, --scale-curv-pdb 1.0 Store curvature in PDB file This option writes the curvature values of smoothed global helical axis in B-factor column of PDB file. This value further scale the calculated curvature value in PDB file. For example: '-scp 10' option will multiply all curvature values with 10 and write it to PDB file's B-factor column.