.\" Automatically generated by Pod::Man 4.14 (Pod::Simple 3.40) .\" .\" Standard preamble: .\" ======================================================================== .de Sp \" Vertical space (when we can't use .PP) .if t .sp .5v .if n .sp .. .de Vb \" Begin verbatim text .ft CW .nf .ne \\$1 .. .de Ve \" End verbatim text .ft R .fi .. .\" Set up some character translations and predefined strings. \*(-- will .\" give an unbreakable dash, \*(PI will give pi, \*(L" will give a left .\" double quote, and \*(R" will give a right double quote. \*(C+ will .\" give a nicer C++. Capital omega is used to do unbreakable dashes and .\" therefore won't be available. \*(C` and \*(C' expand to `' in nroff, .\" nothing in troff, for use with C<>. .tr \(*W- .ds C+ C\v'-.1v'\h'-1p'\s-2+\h'-1p'+\s0\v'.1v'\h'-1p' .ie n \{\ . ds -- \(*W- . ds PI pi . if (\n(.H=4u)&(1m=24u) .ds -- \(*W\h'-12u'\(*W\h'-12u'-\" diablo 10 pitch . if (\n(.H=4u)&(1m=20u) .ds -- \(*W\h'-12u'\(*W\h'-8u'-\" diablo 12 pitch . ds L" "" . ds R" "" . ds C` "" . ds C' "" 'br\} .el\{\ . ds -- \|\(em\| . ds PI \(*p . ds L" `` . ds R" '' . ds C` . ds C' 'br\} .\" .\" Escape single quotes in literal strings from groff's Unicode transform. .ie \n(.g .ds Aq \(aq .el .ds Aq ' .\" .\" If the F register is >0, we'll generate index entries on stderr for .\" titles (.TH), headers (.SH), subsections (.SS), items (.Ip), and index .\" entries marked with X<> in POD. Of course, you'll have to process the .\" output yourself in some meaningful fashion. .\" .\" Avoid warning from groff about undefined register 'F'. .de IX .. .nr rF 0 .if \n(.g .if rF .nr rF 1 .if (\n(rF:(\n(.g==0)) \{\ . if \nF \{\ . de IX . tm Index:\\$1\t\\n%\t"\\$2" .. . if !\nF==2 \{\ . nr % 0 . nr F 2 . \} . \} .\} .rr rF .\" .\" Accent mark definitions (@(#)ms.acc 1.5 88/02/08 SMI; from UCB 4.2). .\" Fear. Run. Save yourself. No user-serviceable parts. . \" fudge factors for nroff and troff .if n \{\ . ds #H 0 . ds #V .8m . ds #F .3m . ds #[ \f1 . ds #] \fP .\} .if t \{\ . ds #H ((1u-(\\\\n(.fu%2u))*.13m) . ds #V .6m . ds #F 0 . ds #[ \& . ds #] \& .\} . \" simple accents for nroff and troff .if n \{\ . ds ' \& . ds ` \& . ds ^ \& . ds , \& . ds ~ ~ . ds / .\} .if t \{\ . ds ' \\k:\h'-(\\n(.wu*8/10-\*(#H)'\'\h"|\\n:u" . ds ` \\k:\h'-(\\n(.wu*8/10-\*(#H)'\`\h'|\\n:u' . ds ^ \\k:\h'-(\\n(.wu*10/11-\*(#H)'^\h'|\\n:u' . ds , \\k:\h'-(\\n(.wu*8/10)',\h'|\\n:u' . ds ~ \\k:\h'-(\\n(.wu-\*(#H-.1m)'~\h'|\\n:u' . ds / \\k:\h'-(\\n(.wu*8/10-\*(#H)'\z\(sl\h'|\\n:u' .\} . \" troff and (daisy-wheel) nroff accents .ds : \\k:\h'-(\\n(.wu*8/10-\*(#H+.1m+\*(#F)'\v'-\*(#V'\z.\h'.2m+\*(#F'.\h'|\\n:u'\v'\*(#V' .ds 8 \h'\*(#H'\(*b\h'-\*(#H' .ds o \\k:\h'-(\\n(.wu+\w'\(de'u-\*(#H)/2u'\v'-.3n'\*(#[\z\(de\v'.3n'\h'|\\n:u'\*(#] .ds d- \h'\*(#H'\(pd\h'-\w'~'u'\v'-.25m'\f2\(hy\fP\v'.25m'\h'-\*(#H' .ds D- D\\k:\h'-\w'D'u'\v'-.11m'\z\(hy\v'.11m'\h'|\\n:u' .ds th \*(#[\v'.3m'\s+1I\s-1\v'-.3m'\h'-(\w'I'u*2/3)'\s-1o\s+1\*(#] .ds Th \*(#[\s+2I\s-2\h'-\w'I'u*3/5'\v'-.3m'o\v'.3m'\*(#] .ds ae a\h'-(\w'a'u*4/10)'e .ds Ae A\h'-(\w'A'u*4/10)'E . \" corrections for vroff .if v .ds ~ \\k:\h'-(\\n(.wu*9/10-\*(#H)'\s-2\u~\d\s+2\h'|\\n:u' .if v .ds ^ \\k:\h'-(\\n(.wu*10/11-\*(#H)'\v'-.4m'^\v'.4m'\h'|\\n:u' . \" for low resolution devices (crt and lpr) .if \n(.H>23 .if \n(.V>19 \ \{\ . ds : e . ds 8 ss . ds o a . ds d- d\h'-1'\(ga . ds D- D\h'-1'\(hy . ds th \o'bp' . ds Th \o'LP' . ds ae ae . ds Ae AE .\} .rm #[ #] #H #V #F C .\" ======================================================================== .\" .IX Title "PROFBVAL 1" .TH PROFBVAL 1 "2022-01-18" "1.0.22" "User Commands" .\" For nroff, turn off justification. Always turn off hyphenation; it makes .\" way too many mistakes in technical documents. .if n .ad l .nh .SH "NAME" profbval \- predicts flexibile/rigid residues from sequence .SH "SYNOPSIS" .IX Header "SYNOPSIS" profbval <\s-1FASTA_FILE\s0> <\s-1RDBPROF_FILE\s0> <\s-1HSSP_FILE\s0> <\s-1WINDOW\s0> <\s-1DEBUG\s0> .SH "DESCRIPTION" .IX Header "DESCRIPTION" PROFbval is a neural-network based trained on backbone B\-value data from X\-ray structure. PROFbval was trained on a sequence unique set of high-resolution protein structures from the \s-1PDB.\s0 .PP The mobility of a given residue on the protein surface is related to its functional role. A common measure of atom mobility in proteins is B\-value data from x\-ray crystallography structures. PROFbval is a method predicting backbone B\-values from amino-acid sequence. PROFbval can be useful for both protein structure and function predictions. For instance, a biologist can locate potentially antigenic determinants by identifying the most flexible residues on the protein surface. Additionally, a crystallographer can locate residues that potentially have high experimental B\-values. .SS "Conversion of PSI-BLAST alignment to \s-1HSSP\s0 format" .IX Subsection "Conversion of PSI-BLAST alignment to HSSP format" The most up-to-date procedure can be found at . .IP "1. Convert \s-1BLAST\s0 output to a Single Alignment Format (\s-1SAF\s0):" 4 .IX Item "1. Convert BLAST output to a Single Alignment Format (SAF):" .Vb 2 \& /usr/share/librg\-utils\-perl/blast2saf.pl fasta= maxAli=3000 eSaf=1 \e \& saf= .Ve .IP "2. Convert \s-1SAF\s0 format to \s-1HSSP:\s0" 4 .IX Item "2. Convert SAF format to HSSP:" .Vb 2 \& /usr/share/librg\-utils\-perl/copf.pl formatIn=saf formatOut=hssp \e \& fileOut= exeConvertSeq=convert_seq .Ve .IP "3. Filter results to 80% redundancy:" 4 .IX Item "3. Filter results to 80% redundancy:" .Vb 1 \& /usr/share/librg\-utils\-perl/hssp_filter.pl red=80 fileOut= .Ve .SS "Output format" .IX Subsection "Output format" This description applies to the default output format. .IP "number" 4 .IX Item "number" residue number .IP "residue" 4 .IX Item "residue" residue type .IP "raw" 4 .IX Item "raw" raw value of the different between the two output nodes .IP "Bnorm" 4 .IX Item "Bnorm" predicted normalized B\-value the values are normalized so in a given sequence the average value is 0 and the standard deviation is 1. The higher the value the more flexible a residue is predicted to be. .IP "Non-strict mode (\s-1NS\s0)" 4 .IX Item "Non-strict mode (NS)" indicates on a flexible residue. According to the \s-1NS\s0 mode most of the residues are flexible; hence, a residue on the surface that is predicted to be rigid is likely to have a functional role. .IP "Strict mode (S)" 4 .IX Item "Strict mode (S)" indicates on a flexible residue. According to the S mode about a third of the residues are flexible, therefore, a stretch of residues that is predicted to be flexible might be important for the protein function. .SH "REFERENCES" .IX Header "REFERENCES" .IP "Schlessinger A, Yachdav G, & Rost B. PROFbval: predict flexible and rigid residues in proteins. Bioinformatics. 2006 Apr 1;22(7):891\-3." 4 .IX Item "Schlessinger A, Yachdav G, & Rost B. PROFbval: predict flexible and rigid residues in proteins. Bioinformatics. 2006 Apr 1;22(7):891-3." .PD 0 .IP "Schlessinger A, Rost B. Protein flexibility and rigidity predicted from sequence. Proteins. 2005 Oct 1;61(1):115\-26." 4 .IX Item "Schlessinger A, Rost B. Protein flexibility and rigidity predicted from sequence. Proteins. 2005 Oct 1;61(1):115-26." .PD .SH "OPTIONS" .IX Header "OPTIONS" .IP "\s-1FASTA_FILE\s0" 4 .IX Item "FASTA_FILE" File containing protein amino-acid sequence in fasta format. .IP "\s-1RDBPROF_FILE\s0" 4 .IX Item "RDBPROF_FILE" Secondary structure and solvent accessibility prediction by \s-1PROF\s0 in rdb format. .IP "\s-1HSSP_FILE\s0" 4 .IX Item "HSSP_FILE" PSI-BLAST alignment profile file converted to \s-1HSSP\s0 format. .IP "\s-1OUTPUT_FILE\s0" 4 .IX Item "OUTPUT_FILE" The name of the final PROFbavl output file. You may give multiple output files. You then probably want to give a list of output modes as well, one for each output file. Separate entries with ',' (comma is not allowed in the file name). .IP "\s-1WINDOW\s0" 4 .IX Item "WINDOW" This is the desired smoothing window for the output. The method was trained on window of 9. .IP "\s-1OUTPUT_MODE\s0" 4 .IX Item "OUTPUT_MODE" PROFbval can create output files in different formats for different purposes. The default output mode for the package is 6. You may give multiple output modes. You then probably want to give a list of output files as well, one for each output mode. Separate entries with ','. .Sp Modes: '\-1', '0', '1', '2', '3', '4', '5', '6', 'snap', 'snapfun' .RS 4 .IP "5" 4 .IX Item "5" For \fBmetadisorder\fR\|(1) .IP "snap" 4 .IX Item "snap" .PD 0 .IP "snapfun" 4 .IX Item "snapfun" .PD For \fBsnapfun\fR\|(1) .RE .RS 4 .RE .IP "\s-1DEBUG\s0" 4 .IX Item "DEBUG" Default: 0. Enable with 1. .SH "EXAMPLE" .IX Header "EXAMPLE" .Vb 1 \& profbval /usr/share/doc/profbval/examples/cad23.f /usr/share/doc/profbval/examples/cad23\-fil.rdbProf /usr/share/doc/profbval/examples/cad23\-fil.hssp cad23.profbval 9 6 .Ve .SH "ENVIRONMENT" .IX Header "ENVIRONMENT" .IP "\s-1PROFBVAL_ROOT\s0" 4 .IX Item "PROFBVAL_ROOT" The directory used to look up \fI./scr/createDataFile.pl\fR, \fI./scr/PROFbval.pl\fR and \fI./nn_files/jct.in\fR. If unset \fI/usr/share/profbval\fR is used. .SH "FILES" .IX Header "FILES" .IP "\fI*.profbval\fR" 4 .IX Item "*.profbval" default output file extension .IP "\fI/usr/share/doc/profbval/examples\fR" 4 .IX Item "/usr/share/doc/profbval/examples" default precomputed input files directory .SH "NOTES" .IX Header "NOTES" .IP "1. It is recommended to create the profiles using 3 iteration of PSI-BLAST against big database" 4 .IX Item "1. It is recommended to create the profiles using 3 iteration of PSI-BLAST against big database" .PD 0 .IP "2. It is also recommended to filter the hssp files using hssp_filter.pl from the Prof package using the following command: perl hssp_filter.pl hssp_file red=80" 4 .IX Item "2. It is also recommended to filter the hssp files using hssp_filter.pl from the Prof package using the following command: perl hssp_filter.pl hssp_file red=80" .PD .SH "AUTHOR" .IX Header "AUTHOR" .IP "A. Schlessinger " 4 .IX Item "A. Schlessinger " .SH "SEE ALSO" .IX Header "SEE ALSO" .PD 0 .IP "Main website" 4 .IX Item "Main website" .PD