User:Eric Martz/Sandbox 11

From Proteopedia

< User:Eric Martz(Difference between revisions)
Jump to: navigation, search
(→TM-Align)
Current revision (22:56, 28 September 2012) (edit) (undo)
 
(32 intermediate revisions not shown.)
Line 1: Line 1:
-
Proposed page title: <big>'''Structural Alignment Tools'''</big>
+
<big><b>Proposed title: Knots in Proteins: Methods</b></big>
-
==Structural Alignment Servers==
+
This page describes technical details of the methods used for the morph at [[Knots in proteins]].
-
Alphabetical, by server name:
+
-
===CE===
+
The [[PDB file]] that displays the animated morph of the protein knot being "shrunk" is [[Image:1yve i knot morph.pdb]]. This file contains 51 models. Model 1 is a slightly smoothed backbone trace of the knotted segment of chain I in [[1yve]], namely 234 residues, sequence numbers 312-545. Model 51 is a simplified knot, which is the result of shrinking the backbone while keeping the ends fixed in space. These models were kindly provided by William R. Taylor ([http://mathbio.nimr.mrc.ac.uk/ Division of Mathematical Biology, National Institute for Medical Research, the Ridgeway, Mill Hill, London UK]). These models represent the residues with alpha carbon atoms only, and all residues are labeled GLY.
-
The ''Combinatorial Extension'' structural alignment server.
+
-
*Server: [http://cl.sdsc.edu/ CE Home Page].
+
-
*Publication (1998)<ref>PMID: 9796821</ref>
+
-
*'''N.B. Database of structure neighbors has not been updated since 2004. Java applet for viewing results is not working in Sept. 2010. Finding structure neighbors from the entire PDB database ("ALL") appears to have been broken since 2001.'''
+
-
*Rigid alignment: ONLY (according to FATCAT<ref name="fatcat" />)
+
-
*Structure-based sequence alignment: YES.
+
-
===DALI===
+
For the [http://www.umass.edu/microbio/chime/knots/ website developed in 2000], the coordinates of these models were reduced (using the MS-DOS [http://www.umass.edu/microbio/rasmol/pdbtools.htm PDB Tools] program shrink.exe) so the visualization plug-in MDL Chime would automatically connect the alpha carbon atoms (some being unrealistically distant from each other) with a backbone trace. However, Jmol did not connect consecutive alpha carbon atoms due to the mixture of unrealistically small and large distances between them. Therefore, the following script was employed to generate bonds (''wireframe'') in Jmol between consecutive alpha carbon atoms:
-
*Rigid alignment: ONLY (according to FATCAT<ref name="fatcat" />)
+
<pre>select all # all 51 models
 +
wireframe off
 +
backbone off
 +
spacefill off
 +
connect delete # deletes all the automatically generated inter-atomic bonds
 +
# now connect consecutive alpha carbons. model 1 contains 234 alpha carbons. model 51 contains 14.
 +
for (i=1; i<=234; i++)
 +
{
 +
ii = i + 1;
 +
connect (atomno=i) (atomno=ii);
 +
}
 +
wireframe 0.05 # display the new connections
 +
color group</pre>
-
===FATCAT===
+
The morph for the comparison example that has no knot, [[1tph]] chain "1", is in [[Image:1tph 1 noknot morph.pdb]]. Similar methods were employed, adjusting the maximum residue count in the above script from 234 to 242.
-
 
+
-
*Server: [http://fatcat.burnham.org/ fatcat.burnham.org] Flexible structure AlignmenT by Chaining AFPs (Aligned Fragment Pairs) with Twists (FATCAT)
+
-
*Publication (2003)<ref name="fatcat">PMID: 14534198</ref> "... the FATCAT algorithm achieves more accurate structure alignments than current methods, while at the same time introducing fewer hinges."
+
-
*Help on server: YES with snapshots; some context-sensitive help.
+
-
*Does alignment involve sequence comparison? UNCLEAR.
+
-
*Rigid alignment: YES (optional)
+
-
*Flexible alignment: YES (optional)
+
-
*Structure neighbors (pre-calculated): YES
+
-
*Pairwise alignment including uploaded models: YES
+
-
*Visualization: Jmol or Chime. See ''Special features''.
+
-
*Color by deviation: '''NO'''. (Colors identify twist/hinge boundaries.)
+
-
*Offered by RCSB? YES
+
-
*Special features: Produces a morph between the two aligned chains (at the link "Interpolating between ..."). Produces a sequence alignment. Offers a RasMol script to color each rigid segment distinctly (separated by twists/hinges).
+
-
 
+
-
Notes from the publication: With 10 "difficult examples"<ref>Fischer,D., Elofsson,A., Rice,D. and Eisenberg,D. (1996) Assessing the performance of fold recognition methods by means of a comprehensive benchmark. In Pacific Symposium on Biocomputing. pp. 300–318.</ref> FATCAT produced results comparable (length, RMSD) to the rigid alignment servers DALI, VAST, CE with no twists in 8 cases. This shows that FATCAT is not biased to introduce twists (hinges). Hinges were introduced in two of the difficult cases, producing arguably better alignments. In a comparison with ''FlexProt''<ref>PMID: 12112693</ref>, FATCAT obtained similar RMSD's and aligned lengths with fewer twists (hinges).
+
-
 
+
-
===FlexProt===
+
-
 
+
-
*Server: [http://bioinfo3d.cs.tau.ac.il/FlexProt/ FlexProt].
+
-
*Publication (2002)<ref>PMID: 12112693</ref>
+
-
*Rigid alignment: YES (Results include alignment for 0 hinges, but only a well-aligning subset of residues are aligned.)
+
-
*Flexible alignment: YES (Results are given for various numbers of hinges.)
+
-
*Visualization: '''NONE''' (You can download PDB files.)
+
-
*Ligands: '''Discarded'''.
+
-
*Special features: Assigns a distinct chain name to each rigid segment separated by a hinge, facilitating informative coloring.
+
-
Note: [[#FATCAT|FATCAT]] provides evidence that it out-performs FlexProt.
+
-
 
+
-
===MAMMOTH===
+
-
*Server: [http://ub.cbm.uam.es/servers/mammoth/mammoth.php mammoth MAMMOTH] (MAtching Molecular Models Obtained from THeory)
+
-
*Publication (2002)<ref name="mammoth">PMID: 12381844</ref>
+
-
*Help on server: Little or none.
+
-
*Does alignment involve sequence comparison? NO: They state that this is a "sequence-independent structural alignment".
+
-
*Rigid alignment: YES.
+
-
*Flexible alignment: NO.
+
-
*Multiple alignment: YES.
+
-
*Structure neighbors (pre-calculated): NO.
+
-
*Pairwise alignment including uploaded models: YES
+
-
*Visualization: None (you can download a PDB file and a RasMol script. '''PDB file lacks MODEL/ENDMDL delimiters. PDB file has no chain names. There is a PDB file with chains A and B in the downloadable file ''rasmol.tcl'' but this is not a Jmol-ready file.''')
+
-
*Color by deviation: '''NO'''.
+
-
*Offered by RCSB? YES
+
-
 
+
-
===TM-Align===
+
-
 
+
-
*Server: [http://zhanglab.ccmb.med.umich.edu/TM-align/ TM-align]
+
-
*Publication (2005)<ref name="mammoth">PMID: 15849316</ref>
+
-
*Help on server: A little.
+
-
*Does alignment involve sequence comparison? UNCLEAR.
+
-
*Rigid alignment: YES.
+
-
*Flexible alignment: NO.
+
-
*Multiple alignment: You can download the software to run on linux.
+
-
*Structure neighbors (pre-calculated): NO.
+
-
*Pairwise alignment including uploaded models: YES
+
-
*Visualization: None (you can download a script for RasMol that contains PDB coordinates. '''PDB file lacks MODEL/ENDMDL delimiters. PDB file has no chain names. File does not run as a script in Jmol due to REMARK lines that are not legal Jmol commands.''')
+
-
*Color by deviation: '''NO'''.
+
-
*Offered by RCSB? YES
+
-
 
+
-
===VAST===
+
-
Rigid alignment: ONLY (according to FATCAT<ref name="fatcat" />)
+
-
 
+
-
==Structural Alignment Software==
+
-
 
+
-
==Examples==
+
-
 
+
-
===Example Requiring Flexibility===
+
-
 
+
-
'''This example requires flexibility for a good alignment''': [[2bbm]]:A vs. [[1cfc]]:A. Length: 148. 97% sequence identity (145/148), 99% similar. These files contain calmodulin. In [[2bbm]] (''Drosophila''), the two calcium-binding domains are wrapped around a peptide. In [[1cfc]] (''Xenopus''), there is no calcium and no peptide, and the linker between the two domains is flexible.
+
-
*CE:
+
-
**4.8 &Aring; RMSD.
+
-
**38.5% sequence identity in structure-based sequence alignment. Aligned/gap positions = 109/47.
+
-
**Uses old, unremediated PDB files (1cfc has no chain A).
+
-
*FATCAT:
+
-
**5 hinges(twists): 140 residues aligned, RMSD 2.08 &Aring;.
+
-
*FlexProt:
+
-
**0 hinges: 49 residues aligned, RMSD 2.94 &Aring;.
+
-
**1 hinge: 84 residues aligned, RMSD 2.97 &Aring;.
+
-
**2 hinges: 102 residues aligned, RMSD 2.82 &Aring;.
+
-
**3 hinges: 118 residues aligned, RMSD 2.60 &Aring;.
+
-
**4 hinges: 134 residues aligned, RMSD 2.62 &Aring;.
+
-
 
+
-
===Example for Rigid Alignment===
+
-
[[1fsz]] is the bacterial cell division protein FtsZ, length 334 residues with coordinates (372 in crystallized protein). It has structural similarity to mammalian tubulin found in [[1tub]] chain A, length 440.
+
-
*CE
+
-
**3.2 &Aring; RMSD for 305 residues. The structural alignment has 96 unaligned "gap" residues: one large gap of ~30 residues, and ten smaller gaps of 8 residues or less.
+
-
**12.5% sequence identity.
+
-
*FATCAT
+
-
**3.02 &Aring; RMSD RIGID alignment includes 298 residues.
+
-
**FLEXIBLE alignment introduced ZERO twists (hinges), so gave the same result as the rigid alignment.
+
-
**The structure-based sequence alignment has many gaps, looking similar to that generated by CE.
+
-
*MAMMOTH
+
-
**4.0 &Aring; (?) with 298 aligned residues (?) (Labeling in results is unclear.)
+
-
**Structure-based sequence alignment is displayed.
+
-
 
+
-
==References==
+
-
<references />
+

Current revision

Proposed title: Knots in Proteins: Methods

This page describes technical details of the methods used for the morph at Knots in proteins.

The PDB file that displays the animated morph of the protein knot being "shrunk" is Image:1yve i knot morph.pdb. This file contains 51 models. Model 1 is a slightly smoothed backbone trace of the knotted segment of chain I in 1yve, namely 234 residues, sequence numbers 312-545. Model 51 is a simplified knot, which is the result of shrinking the backbone while keeping the ends fixed in space. These models were kindly provided by William R. Taylor (Division of Mathematical Biology, National Institute for Medical Research, the Ridgeway, Mill Hill, London UK). These models represent the residues with alpha carbon atoms only, and all residues are labeled GLY.

For the website developed in 2000, the coordinates of these models were reduced (using the MS-DOS PDB Tools program shrink.exe) so the visualization plug-in MDL Chime would automatically connect the alpha carbon atoms (some being unrealistically distant from each other) with a backbone trace. However, Jmol did not connect consecutive alpha carbon atoms due to the mixture of unrealistically small and large distances between them. Therefore, the following script was employed to generate bonds (wireframe) in Jmol between consecutive alpha carbon atoms: 

select all # all 51 models
wireframe off
backbone off
spacefill off
connect delete # deletes all the automatically generated inter-atomic bonds
# now connect consecutive alpha carbons. model 1 contains 234 alpha carbons. model 51 contains 14.
for (i=1; i<=234; i++)
{
	ii = i + 1;
	connect (atomno=i) (atomno=ii);
}
wireframe 0.05 # display the new connections
color group

The morph for the comparison example that has no knot, 1tph chain "1", is in Image:1tph 1 noknot morph.pdb. Similar methods were employed, adjusting the maximum residue count in the above script from 234 to 242.

Proteopedia Page Contributors and Editors (what is this?)

Eric Martz

Retrieved from "http://52.214.119.220/wiki/index.php/User:Eric_Martz/Sandbox_11"
Personal tools