Before you begin
Input Data Format
Step by step
Submitting your query
- Go to the box entitled "Input Target-Template Alignment". Pick a way to provide your input. You can either upload a file or copy and paste your data into a text area. Either way, make sure your input sequences are in FASTA or PIR format. Your target and template sequences should be labelled as ">target" and ">template", respectively (without the quotes).
- Go to the box entitled "Input Template Structure". Pick a way to provide your input. You can either enter a 5-character PDB code (the last character is the chain identifier), upload a file or copy and paste your data into a text area. Make sure that any data you provide here is in PDB format.
Go to the box entitled "Options". If you are submitting a query for the first time, do not do anything here,
simply leave the default settings unchanged and move on to the next step.
If you are running MEDELER, you have the option of choosing the type of model you will generate. A more complete model takes longer to generate than just the core model. In case you are submitting a protein that is not an integral membrane protein, you may wish to activate the option to "Start modelling from a bigger minimal core".
You have a choice of different iMembrane search methods. Your choice influences the speed of the search and the number of hits found. The different choices are briefly explained below the radio buttons in the "Options" box. In general, it is fine to leave the search in "auto" mode. You may also choose different cut-off values for your search. The default values are quite loose, but still manage to filter out absolutely terrible hits. In any case, your search results are ranked. If you are running MEDELLER, only the top hit will be used to annotate your template protein, which is then in turn used to annotate the models of your target protein.
- The last box: "Submit". Press the "Submit" button once, and only once please, to submit your query. Your browser may stay inert from a few seconds up to several minutes (in some cases), while it uploads your input to our web site. If you are unsure whether you have successfully clicked the button, have a look at your browser's status bar. If this starts with something like "Sending", "Waiting" or "Uploading" then everything is fine and you simply need to wait. The screen should automatically refresh and forward you to your results.
Viewing your results
The results screen consists of 4 main frames:
- Top-left: The Query and Match structure viewers.
- Top-right: The Sequence Alignment viewer.
- Bottom half: In MEDELLER mode, this area lists all your generated models and some links where you can view or download these. In iMembrane mode, the Results list shows all the hits for your query in a tabular format, one query/hit pair per row. The links in the hits table control what is being shown in the above 3 frames and allow you to download individual hits to your computer.
The structure viewers
These frames show three-dimensional representations of your models (MEDELLER) or your annotated input structure and its database match (iMembrane). Colours represent the annotation produced by iMembrane.
In MEDELLER mode, the left structure is coloured by modelling confidence, i.e. the model core has high confidence (coloured blue) whereas
loop coordinates have lower confidence and ab-initio loops have the lowest conficence (coloured red).
You can turn the molecule by clicking and dragging it with the left mouse button. Right-clicking gives various view options. Consult the Jmol manual for details.
The sequence alignment viewer
This shows the sequence of your input protein, aligned to its corresponding database hit (iMembrane) or your model aligned to your template (MEDELLER).
The lines entitled "membrane layer" and "membrane contact" show the annotation for the database hit, as produced by iMembrane.
The colour key below the sequences explains the meaning of the colouring in both the sequence alignment as well as the structure viewers. Refer to the section "Understanding iMembrane annotation" below for details.
The results list
The text at the top of the Results list mentions your result's Result ID. Make a note of this or bookmark the link to your result page to come back to it later. Note, however, that your result will only be saved for 24 hours before being deleted from the server.
Clicking on the links named "Layer" or "Contact" in each row will change the structures and sequences being displayed in the structure viewer and sequence viewer frames. For the exact meaning of "Layer" and "Contact", refer to the section "Understanding iMembrane annotation" below.
In MEDELLER mode, the list of generated models is shows here: core (a reliable subset of the template co-ordinates), hiacc (high accuracy model, with high accuracy database loops), hicov (high coverage model, with high and low accuracy database loops) and complete (the complete model with high and low accuracy database loops and all remaining gaps filled with ab initio methods). The complete model is the only model containing sidechains, terminal residues not present in the input model as well as very long gap regions (>30 residues). The largest generated model (usually the complete model) is shown by default.
The result table (iMembrane only) has several columns. Their meaning is explained below.
- dbid: The name given to the particular database molecule that matched the above query molecule. This is simply a PDB code, simulation number, chain identifier and a duplicate chain number (some chains appear more than once in a PDB structure, e.g. via rotational operators). The simulation number identifies the molecular dynamics simulation that underlies iMembrane's annotation for this molecule. Clicking on the dbid will bring up the particular CGDB entry associated with this molecule.
- E-value: The BLAST E-value corresponding to this match. Values below 1.0E-20 (ten to the power of -20) usually assure a good quality structure-based annotation, while values below 1.0E-40 (ten to the power of -40) tend to assure a good sequence-based annotation. Higher values do not guarantee a bad annotation.
- TM-score: The TM-score corresponding to the final alignment of the query and match molecules. Values above 0.5 tend to indicate an acceptable match.
- pannot: The percentage of your query protein that could be annotated using this match (measured on the sequence-based "membrane contact" annotation). Currently this is always equal to the coverage. See below.
- pcov: The percentage coverage, i.e. fraction of query residues aligned
[100 * count(aligned residues) / length(query)].
pid <= pcov.
- pid: The percentage sequence identity between the query and match molecules, computed from the final alignment between the two. This number is relative to the length of the query molecule
[100 * count(identical residues) / length(query)].
pid <= pcov.
- RMSD: The Root Mean Square Deviation of the structural alignment between the query and match proteins. This is not an indicator of annotation quality.
- View Results: This column contains links to change the molecules shown in the "Query", "Match" and "Sequence Alignment" frames. Clicking on the "Contact" or "Layer" link in any of the rows will display the corresponding query/match pair, coloured by either "membrane layer" or "membrane contact". See the section "Understanding iMembrane annotation" below for detail on the meaning of these annotation types.
- Download Results: This column contains links to download individual query/match pairs and all associated files to your computer. Clicking on the link will download a g-zipped tar archive (file extension "tar.gz"). Make sure to give the archive an appropriate name if you intend to keep it. Most modern archiving programmes support tar.gz archives. The author recommends 7-zip to Windows users. Linux and Mac users can unzip the archive with the "tar zxf <filename>" command.
Understanding iMembrane annotation
The annotation produced by iMembrane comes in two variants:
- Membrane contact
- Membrane layer
This is the most basic type of annotation. The data is merely a descretised version of the annotation obtained from the original molecular dynamics experiments performed on the database match protein. The annotation distinguishes between residues that are in direct contact with the inner (lipid tails) or outer (polar head groups) membrane layers and those residues that are not at all in contact with the membrane. By aligning the query sequence to the match sequence, this annotation can be transferred over to the query protein, in regions where the two proteins are similar.
This is a more abstract and simplified model of the membrane. By using the above membrane contact data, one can fit parallel plains onto the database protein. These planes represent the limits between the different layers of the lipid bilayer. The protein's residues are essentially classified into three categories: Those in the middle layer of the membrane, where the hydrophobic tails of the membrane lipids reside; the residues in the outer membrane layers, where the polar head groups of the membrane lipids reside; the residues outside the membrane. Note that, unlike the contact annotation, being in a certain layer does not mean that the residue in question is actually in contact with this environment. For example, residues in the core of the protein can be classified as being in the middle of the membrane, if the are positioned on the same level as the lipid tails along the bilayer normal.