ClusteringTool - 0.9

GibbsCluster takes as input a list of peptide sequences and attempts to cluster them into meaningful groups.
Sequence alignment and clustering are performed simultaneously by sampling the space of possible solutions using a Gibbs sampling strategy. Each cluster is represented by a position-specific scoring matrix (PSSM), and the algorithm aims at maximizing the information content of individual matrices while minimizing the overlap between distinct clusters. The server returns a detailed report on the optimal clustering solutions, including plots of the optimal number of clusters and graphical representations of the identified sequence motifs as sequence logos. For details on the algorithm refer to this paper

This page introduces the data formats, the parameters available to customize the analysis and some guidelines for the use of version 2.0. Users are welcome to contact the authors for any questions.

Paste a set of peptides (up to 50 amino acids in length), one sequence per line into the upper left window, or upload a file from your local disk.

The input file can be a plain list of peptides (Sample 1) or an annotated list (Sample 2). The annotation is carried over to the results, and may be useful for correlating a known classification with the clustering produced by the method. All input sequences must be in one-letter amino acid code. The allowed alphabet (not case sensitive) is as follows:

A C D E F G H I K L M N P Q R S T V W Y

The options are divided in three levels: Basic, Filtering and SOM. The only essential parameters to be set are in the Basic options; if you are using the server for the first time you may leave all the other options unchanged.

A brief explanation of each option can be visualized by hovering the mouse over the symbol next to each option in the submission form.

BASIC options

Job name:
This prefix is pre-pended to all files generated by the current run. If left empty, a system-generated number will be assigned as prefix.

Number of clusters:
You may provide a specific number of clusters (e.g. 3), or an interval of partitions (e.g. 2-5). In the second case, the method will suggest the optimal number of cluster it found in the data, given the parameter configuration of the job. Maximum number of clusters: 15.

FILTERING options

Make clustering moves at each iteration:
By default, simple shift moves are performed at each iteration, indel moves every 10 iterations, single peptide moves every 20 iterations, phase shift moves every 100 iterations.
You can alter this behavior by ticking this option; simple shift and phase shift moves become disabled, and single peptide moves are made at each iteration. This set-up is recommended for "nearly-aligned" data, where clustering and indels should be sampled more regularly than extensions at the termini. That is the case, for example, of sets of MHC class I ligands of different length, which would in most cases require central indels to model peptide bulging of long ligand.

SOM options

Number of iterations per sequence per temperature step:
This parameter ("I") specificies how long your clustering schedule should be. Note that total number of iterations is the results of "I" multiplied by the number of sequences times the number of temperature steps, and it will increase linearly the execution time.

Initial Monte Carlo temperature:
The temperature is a scalar, lowered by discreet steps as the iterations progress. The temperature influences the probability of accepting or rejecting the moves of the algorith. In the initial iterations (high temperature) the program is free to explore the landscape of solutions, and as the system cools off only moves that increase the energy will be accepted.

Number of temperature steps:
The number of steps in the cooling schedule (starting from the initial temperature specified above).

Interval between Indel moves:
Specifies how often to attempt introducing insertions and deletions (see glossary).

Interval between Single peptide moves:
Specifies how often to attempt moving a sequence between clusters (see glossary).

Interval between Phase shift moves:
Specifies how often to attempt shifting the alignment window of a single cluster.

Background amino acid frequencies:
Construction of PSSMs relies on calculating the frequency of a given residue at a given position, compared to the expected background frequency of that amino acid. You may use a flat background model identical for all amino acids (Flat), a pre-calculated distribution reflecting the relative frequency of each residue in naturally occurring proteins (Pre-calculated Uniprot), or determine the background model directly from the dataset you submitted (From data).

Preference for hydrophobic AAs at P1:
In the special case of MHC class II data, we have previously found helpful to guide the alignment by expressing a preference for hydrophobic residues at the P1 of the alignment.

Sequence weighting type:
Data redundancy may affect the quality of the clustering. You may use an explicit clustering of the sequences in a given group (Clustering), or use a faster heuristic that calculates the degree of variability at each column in the alignment (Heuristic, recommended); you may also disable sequence weighting for downweighting of redundant sequences (None).

Click on the "Submit query" button. The status of your job (either 'queued' or 'running') will be displayed and constantly updated until it terminates and the server output appears in the browser window.

At any time during the wait you may enter your e-mail address and simply leave the window. Your job will continue; you will be notified by e-mail when it has terminated. The e-mail message will contain the URL under which the results are stored; they will remain on the server for 24 hours for you to collect them.

• PSSM: Position-specific scoring matrix. A matrix of size L x A, where L is the length of the alignment window and A is the length of the alphabet; it stores the weight (or "preference") or any given amino acid at each position of the alignment.
• Indels: Insertions and deletions. An insertion adds one or more gaps in the alignment, a deletion removes one or more amino acids of a given sequence.
• Simple shift move: A move affecting the alignment core of a single peptide. The algorithm attempts to shift the alignment core of a peptide to the left or right, and accept/reject the move with a probability that depends on the temperature of the system.
• Single peptide move: A move that attempts to change the cluster of a single peptide sequence.
• Indel move: A move that introduces or remove an insertion or a deletion in a peptide sequence. In the GibbsCluster, indels can have any length up to the values specified as parameters, but there can only be one single indel stretch per peptide (i.e. a single sequence cannot contain indels at multiple positions)
• Phase shift move: A move affecting the alignment core of all peptide sequences in a cluster. In simple terms, the algorithm will attempt to shift en bloc the alignment window of a given cluster.