Skittles - Revision history

Christoph: /* External links */

2015-03-02T20:57:22Z

‎External links

Christoph: /* See also */

2012-07-13T00:37:34Z

‎See also

Christoph at 02:35, 9 May 2012

2012-05-09T02:35:51Z

Christoph: /* SKTTLS Report */

2012-05-09T02:33:35Z

‎SKTTLS Report

Christoph at 22:26, 8 May 2012

2012-05-08T22:26:22Z

Christoph: New page: '''Skittles''' (aka '''SKTTLS''') is a program I helped to write for the analysis of results of several kinds of tests on the TLS model (SKTTLS) comparing ADPs for the atoms in bonds betwe...

2012-05-08T21:48:53Z

New page: '''Skittles''' (aka '''SKTTLS''') is a program I helped to write for the analysis of results of several kinds of tests on the TLS model (SKTTLS) comparing ADPs for the atoms in bonds betwe...

New page

'''Skittles''' (aka '''SKTTLS''') is a program I helped to write for the analysis of results of several kinds of tests on the TLS model (SKTTLS) comparing ADPs for the atoms in bonds between residues.<ref>Frank Zucker, '''P. Christoph Champ''', and Ethan A. Merritt (2010). Validation of crystallographic models containing TLS or other descriptions of anisotropy. ''Acta Cryst., D66:889-900''. {{doi|10.1107/S0907444910020421}}</ref> The code was written in C, Fortran, and Python. The same output is written to the log file after the analysis of tensors. The SKTOUT file can be used for plotting with programs such as [[gnuplot]]. It generates a "SKTTLS Report" file.

==SKTTLS Report==
Several kinds of tests on the TLS model (SKTTLS) are reported following the tensor analysis; if a SKTOUT file is given, a copy of this SKTTLS report is also written to that file. Three residuals are calculated for each bond between protein or nucleic acid residues:

;CCuij : the correlation of anisotropic ADPs from Merritt (1999)<ref>E.A.Merritt Acta Cryst., D55, 1997 (1999)</ref> (ranges from 1 down);
;rSIMU : residual of the SIMU (BFAC) restraint, rmsd of ANISOU values (ranges from 0 up); and
;rDELU : residual of the DELU (RBON) Rosenfeld rigid bond restraint (ranges from 0 up).

The 95th and 99th percentile values of these three residuals were calculated from a survey of the REFMAC refinements with segmented TLS in the PDB as of Sept. 2009. Extreme values of these residuals, e.g. CCuij below the 99th percentile, may indicate problems with the ADPs for the atoms in the bond: the TLS segment boundary was misassigned, or something went wrong during refinement, or one or both atoms have non-positive definite ADPs.

The SKTTLS report has summaries and the full table of these residuals plus a table of the distribution of anisotropies. Summaries include:

*the number of bonds for which residuals were counted, and the number of bonds with any residual beyond the 95th or 99th percentile;
*the mean, standard deviation and worst value of each residual for this structure;
*the mean and standard deviation anisotropy for all protein, nucleic acid and other non-solvent atoms in this structure;
*a list of up to 100 outliers, bonds in this structure which are beyond the 95th percentils for any residual. If there are more than 100 such outliers, this analysis should be run again after whatever caused those problems is fixed.

The table of all residuals, SKTTLS TABLE 1, and the table of anisotropy distributions, SKTTLS TABLE 2, are formatted for loggraph or xloggraph plotting. In addition, the files are formatted so that simple extracts of these tables can be plotted by programs such as gnuplot: all non-blank lines other than table data lines are prefixed with '#' and each table is bracketted by lines of the form: "### START OF SKTTLS TABLE N ###" and "### END OF SKTTLS TABLE N ###".

SKTTLS TABLE 1 contains bond identifiers, TLS group identifiers, and residual values. Bond identifiers are residue number, chain ID as text and as a number (to facilitate plotting of separate chains in gnuplot) and bond name e.g. C1-N2 or, if there are alternate atoms, C(A)1-N(A)2. TLS group IDs are TLS group number (numbered sequentially from 1 by their order in TLSIN, rather than the given group number) and a point at each segment break. The opposite of rDELU is included after the residual values so that all 3 residuals can be plotted by loggraph with minimal overlapping. The first two lines of the table contain values for the 95th and 99th percentile of each residual.

SKTTLS TABLE 2 contains anisotropy bin minimum, maximum and center, and fractions and counts for protein, nucleic acid or other non-solvent atoms. Non-positive definite ADPs are counted in the first bin, -0.05 to 0.00; isotropic ADPS are counted in the last bin, 1.00 to 1.00.

SKTTLS TABLE 3 is the list of outliers, containing bond identifiers and all three residuals values for each outlying bond. Marks (? or !) show which residuals were beyond the 95 or 99th percentile. This table is not formatted for plotting.

NOTE: This is a change to the previous behaviour (CCP4 6.1.1 and earlier) where no ANISOU records from XYZIN were written to XYZOUT.

==See also==
*[[TLSMD]]
*[[Raster3D]]

==References==
<references/>

==External links==
*[http://www.ccp4.ac.uk/html/tlsanl.html tlsanl] — analysis of TLS tensors and derived anisotropic U factors (part of [[:Category:Ccp4|CCP4]])

[[Category:Crystallography]]
[[Category:Ccp4|CCP4]]

← Older revision		Revision as of 20:57, 2 March 2015
Line 50:		Line 50:
	[[Category:Crystallography]]		[[Category:Crystallography]]
	[[Category:Ccp4\|CCP4]]		[[Category:Ccp4\|CCP4]]
		+	[[Category:Portfolio]]

← Older revision		Revision as of 02:35, 9 May 2012
Line 1:		Line 1:
−	'''Skittles''' (aka '''SKTTLS''') is a program I helped to write for the analysis of results of several kinds of tests on the TLS model (SKTTLS) comparing ADPs for the atoms in bonds between residues.<ref name="Zucker2010">Frank Zucker, '''P. Christoph Champ''', and Ethan A. Merritt (2010). Validation of crystallographic models containing TLS or other descriptions of anisotropy. ''Acta Cryst., D66:889-900''. {{doi\|10.1107/S0907444910020421}}</ref> The code was written in C, Fortran, and Python. The same output is written to the log file after the analysis of tensors. The SKTOUT file can be used for plotting with programs such as [[gnuplot]]. It generates a "SKTTLS Report" file.	+	'''Skittles''' (aka '''SKTTLS''') is a program I helped to write for the analysis of results of several kinds of tests on the TLS model (SKTTLS) comparing ADPs for the atoms in bonds between residues.<ref name="Zucker2010">Frank Zucker, '''P. Christoph Champ''', and Ethan A. Merritt (2010). "Validation of crystallographic models containing TLS or other descriptions of anisotropy". ''Acta Cryst., D66:889-900''. {{doi\|10.1107/S0907444910020421}}</ref> The code was written in C, Fortran, and [[Python]]. The same output is written to the log file after the analysis of tensors. The SKTOUT file can be used for plotting with programs such as [[gnuplot]]. It generates a "SKTTLS Report" file.

	==Description==		==Description==

@@ Line 38: / Line 38: @@
 ==See also==
 *[[TLSMD]]
 *[[Raster3D]]
 *[[Ramachandran plot]]

@@ Line 9: / Line 9: @@
 Several kinds of tests on the TLS model (SKTTLS) are reported following the tensor analysis; if a SKTOUT file is given, a copy of this SKTTLS report is also written to that file. Three residuals are calculated for each bond between protein or nucleic acid residues:
-;CCuij : the correlation of anisotropic ADPs from Merritt (1999)<ref>E.A.Merritt Acta Cryst., D55, 1997 (1999)</ref> (ranges from 1 down);
+;CCuij : the correlation of anisotropic ADPs from Merritt (1999)<ref>Merritt EA (1999). "[http://www.ncbi.nlm.nih.gov/pubmed/10666575 Comparing anisotropic displacement parameters in protein structures]". ''Acta Crystallogr D Biol Crystallogr.'', '''55'''(Pt 12):1997-2004. pmid:10666575.</ref> (ranges from 1 down);
 ;rSIMU : residual of the SIMU (BFAC) restraint, rmsd of ANISOU values (ranges from 0 up); and
 ;rDELU : residual of the DELU (RBON) Rosenfeld rigid bond restraint (ranges from 0 up).
-The 95th and 99th percentile values of these three residuals were calculated from a survey of the REFMAC refinements with segmented TLS in the PDB as of Sept. 2009. Extreme values of these residuals, e.g. CCuij below the 99th percentile, may indicate problems with the ADPs for the atoms in the bond: the TLS segment boundary was misassigned, or something went wrong during refinement, or one or both atoms have non-positive definite ADPs.
+The 95th and 99th percentile values of these three residuals were calculated from a survey of the REFMAC refinements with segmented TLS in the PDB as of September 2009. Extreme values of these residuals, e.g. CCuij below the 99th percentile, may indicate problems with the ADPs for the atoms in the bond: the TLS segment boundary was misassigned, or something went wrong during refinement, or one or both atoms have non-positive definite ADPs.
 The SKTTLS report has summaries and the full table of these residuals plus a table of the distribution of anisotropies. Summaries include:
@@ Line 22: / Line 22: @@
 *a list of up to 100 outliers, bonds in this structure which are beyond the 95th percentils for any residual. If there are more than 100 such outliers, this analysis should be run again after whatever caused those problems is fixed.
-The table of all residuals, SKTTLS TABLE 1, and the table of anisotropy distributions, SKTTLS TABLE 2, are formatted for loggraph or xloggraph plotting. In addition, the files are formatted so that simple extracts of these tables can be plotted by programs such as gnuplot: all non-blank lines other than table data lines are prefixed with '#' and each table is bracketted by lines of the form: "### START OF SKTTLS TABLE N ###" and "### END OF SKTTLS TABLE N ###".
+The table of all residuals, SKTTLS TABLE 1, and the table of anisotropy distributions, SKTTLS TABLE 2, are formatted for <code>loggraph</code> or <code>xloggraph</code> plotting. In addition, the files are formatted so that simple extracts of these tables can be plotted by programs such as [[gnuplot]]: all non-blank lines other than table data lines are prefixed with '#' and each table is bracketted by lines of the form: "<code>### START OF SKTTLS TABLE N ###</code>" and "<code>### END OF SKTTLS TABLE N ###</code>".
-SKTTLS TABLE 1 contains bond identifiers, TLS group identifiers, and residual values. Bond identifiers are residue number, chain ID as text and as a number (to facilitate plotting of separate chains in gnuplot) and bond name e.g. C1-N2 or, if there are alternate atoms, C(A)1-N(A)2. TLS group IDs are TLS group number (numbered sequentially from 1 by their order in TLSIN, rather than the given group number) and a point at each segment break. The opposite of rDELU is included after the residual values so that all 3 residuals can be plotted by loggraph with minimal overlapping. The first two lines of the table contain values for the 95th and 99th percentile of each residual.
+;SKTTLS TABLE 1: contains bond identifiers, TLS group identifiers, and residual values. Bond identifiers are residue number, chain ID as text and as a number (to facilitate plotting of separate chains in gnuplot) and bond name e.g. C1-N2 or, if there are alternate atoms, C(A)1-N(A)2. TLS group IDs are TLS group number (numbered sequentially from 1 by their order in TLSIN, rather than the given group number) and a point at each segment break. The opposite of rDELU is included after the residual values so that all 3 residuals can be plotted by loggraph with minimal overlapping. The first two lines of the table contain values for the 95th and 99th percentile of each residual.
-SKTTLS TABLE 2 contains anisotropy bin minimum, maximum and center, and fractions and counts for protein, nucleic acid or other non-solvent atoms. Non-positive definite ADPs are counted in the first bin, -0.05 to 0.00; isotropic ADPS are counted in the last bin, 1.00 to 1.00.
+;SKTTLS TABLE 2: contains anisotropy bin minimum, maximum and center, and fractions and counts for protein, nucleic acid or other non-solvent atoms. Non-positive definite ADPs are counted in the first bin, -0.05 to 0.00; isotropic ADPS are counted in the last bin, 1.00 to 1.00.
-SKTTLS TABLE 3 is the list of outliers, containing bond identifiers and all three residuals values for each outlying bond. Marks (? or !) show which residuals were beyond the 95 or 99th percentile. This table is not formatted for plotting.
+;SKTTLS TABLE 3: is the list of outliers, containing bond identifiers and all three residuals values for each outlying bond. Marks (? or !) show which residuals were beyond the 95 or 99th percentile. This table is not formatted for plotting.
-NOTE: This is a change to the previous behaviour (CCP4 6.1.1 and earlier) where no ANISOU records from XYZIN were written to XYZOUT.
+''NOTE: This is a change to the previous behaviour (CCP4 6.1.1 and earlier) where no ANISOU records from XYZIN were written to XYZOUT.''
 ==Equations==

@@ Line 1: / Line 1: @@
-'''Skittles''' (aka '''SKTTLS''') is a program I helped to write for the analysis of results of several kinds of tests on the TLS model (SKTTLS) comparing ADPs for the atoms in bonds between residues.<ref>Frank Zucker, '''P. Christoph Champ''', and Ethan A. Merritt (2010). Validation of crystallographic models containing TLS or other descriptions of anisotropy. ''Acta Cryst., D66:889-900''. {{doi|10.1107/S0907444910020421}}</ref> The code was written in C, Fortran, and Python. The same output is written to the log file after the analysis of tensors. The SKTOUT file can be used for plotting with programs such as [[gnuplot]]. It generates a "SKTTLS Report" file.
+'''Skittles''' (aka '''SKTTLS''') is a program I helped to write for the analysis of results of several kinds of tests on the TLS model (SKTTLS) comparing ADPs for the atoms in bonds between residues.<ref name="Zucker2010">Frank Zucker, '''P. Christoph Champ''', and Ethan A. Merritt (2010). Validation of crystallographic models containing TLS or other descriptions of anisotropy. ''Acta Cryst., D66:889-900''. {{doi|10.1107/S0907444910020421}}</ref> The code was written in C, Fortran, and Python. The same output is written to the log file after the analysis of tensors. The SKTOUT file can be used for plotting with programs such as [[gnuplot]]. It generates a "SKTTLS Report" file.
 ==SKTTLS Report==
@@ Line 26: / Line 31: @@
 NOTE: This is a change to the previous behaviour (CCP4 6.1.1 and earlier) where no ANISOU records from XYZIN were written to XYZOUT.
 ==See also==
 *[[TLSMD]]
 *[[Raster3D]]
 ==References==