User:Boris/Temp/APB

Note! This assignment is currently inactive. Major and minor unannounced changes may be made at any time.

 
 

   

Assignment 4 - Homology modeling

How could the search for ultimate truth have revealed so hideous and visceral-looking an object?

Max Perutz (on his first glimpse of the Hemoglobin structure)

   

Where is the hidden beauty in structure, and where, the "ultimate truth"? In the previous assignments we have studied sequence conservation in APSES family domains and we have seen homologues in all fungal species. This is an ancient protein family that had already duplicated to several paralogues at the time the cenancestor of all fungi lived, more than 600,000,000 years ago, in the Vendian period of the Proterozoic era of Precambrian times.

In order to understand how specific residues in the sequence contribute to the putative function of the protein, and why and how they are conserved throughout evolution, we would need to study an explicit molecular model of an APSES domain protein, bound to its cognate DNA sequence. Explanations of a protein's observed properties and functions can't rely on the general fact that it binds DNA, we need to consider details in terms of specific residues and their spatial arrangement. In particular, it would be interesting to correlate the conservation patterns of key residues with their potential to make specific DNA binding interactions. Unfortunately, no APSES domain structures in complex with bound DNA has been solved up to now, and the experimental evidence we have considered in Assignment 2 (Taylor et al., 2000) is not sufficient to unambiguously define the details of how a DNA double helix might be bound. Moreover, at least two distinct modes of DNA binding are known for proteins of the winged-helix superfamily, of which the APSES domain is a member.

In this assignment you will (1) construct a molecular model of the Mbp1 orthologue in your assigned organism, (2) identify similar structures of distantly related domains for which protein-DNA complexes are known, (3) define whether the available evidence allows you to distinguish between different modes of ligand binding, and (4) assemble a hypothetical complex structure.

For the following, please remember the following terminology:

Target

The protein that you are planning to model.

Template

The protein whose structure you are using as a guide to build the model.

Model

The structure that results from the modeling process. It has the Target sequence and is similar to the Template structure.

 

A brief overview article on the construction and use of homology models is linked to the resource section at the bottom of this page. That section also contains links to other sites and resources you might require.

Preparation, submission and due date

Read carefully.

Be sure you have understood all parts of the assignment and cover all questions in your answers! Sadly, we see too many assignments which, arduously effected, nevertheless intimate nescience of elementary tenets of molecular biology. If the sentence above did not trigger an urge to open a dictionary, you are trying to guess, rather than confirm possibly important information.

Review the guidelines for preparation and submission of BCH441 assignments.

The due date for the assignment is Monday, November 12 at 10:00 in the morning.

   

(1) Preparation

Template choice and sequence (1 marks)

 
Often more than one related structure can be found in the PDB. We have touched on principles of selecting template structures in the lecture and there is a short summary of template choice principles on this Wiki. One can either search the PDB itself through its Advanced Search interface; for example one can search for sequence similarity with a BLAST search, or search for structural similarity by accessing structures according to their CATH or SCOP classification. But one can always also use the BLAST interface at the NCBI, since the sequences contained in PDB files are accessible as a database subsection on the BLAST menu.

• Use the NCBI BLAST interface to identify all PDB files that are clearly homologous to your target APSES domain, if you haven't already done so in Assignment 2. Document that you have searched in the correct subsection of the database by selecting "pdb" on the database options menu. For the hits you find, consider how these coordinate sets differ and which features would make each more or less suitable for your task by commenting briefly on

• sequence similarity to your target
• size of expected model (length of alignment)
• presence or absence of ligands
• experimental method and quality of the data set

Then choose the template you consider the most suitable and note why you have decided to use this template.

• Retrieve the most suitable template structure coordinate file from the PDB.

It is not straightforward at all how to number sequence in such a project. The "natural" numbering is to start a sequential numbering from the start-codon of the full length protein and go sequentially from there. However, this does not map well with other numbering schemes we have encountered. As you know the first residue of the APSES domain as the CDD defines it is not Residue 1 of the Mbp1 protein. The first residue of the e.g. 1MB1 FASTA file is the first residue of Mbp1 protein, but the last five residues are an artifical His tag. Is H125 of 1MB1 thus equivalent to R125 in MBP1_SACCE? The N-terminus of the Mbp1 crystal structure is disordered. The first residue in the structure is ASN 3, therefore N is the first residue in a FASTA sequence derived from the cordinate section of the PDB file (the ATOM records; whereas the SEQRES records start with MET ... and so on. You need to remember: a sequence number is not absolute, but derived from a particular context.

The homology model will be based on an alignment of target and template. Thus we have to define the target sequence. As discussed in class, PDB files have an explicit and an implied sequence and these do not necessarily have to be the same. To compare the implied and the explicit sequence for the template, you need to extract sequence information from coordinates. One way to do this is via the Web interface for WhatIf, a crystallography and molecular modeling package that offers many useful tools for coordinate manipulation tasks.

• Navigate to the Administration sub-menu of the WhatIf Web server. Follow the link to Make sequence file from PDB file. Enter the PDB-ID of your template into the form field and Send the request to the server. The server accesses the PDB file and extracts sequence information directly from the ATOM   records of the file. The results will be returned in PIR format. Copy the results, edit them to FASTA format and save them in a text-only file. Make sure you create a valid FASTA formatted file! Use this implied sequence to check if and how it differs from the sequence ...

• ... listed in the SEQRES records of the coordinate file;
• ... given in the FASTA sequence for the template, which is provided by the PDB;
• ... stored in the protein database of the NCBI.

and record your results.

• In a table, establish how the sequence numbers in the coordinate section of your template(*) correspond to your target sequence numbering.

(*) These residue numbers are important, since they are referenced e.g. by VMD when you visualize the structure. The easiest way to list them is via the Sequence Viewer extension of VMD..

Don't do this for every residue individually but define ranges. Look at the correspondence of the first and last residue of target and template sequence and take indels into account. Establishing sequence correspondence precisely is crucially important! For example, when a publication refers to a residue by its sequence number, you have to be able to relate that number to the residue numbers of the model as well as your target sequence..

   

The input alignment (1 mark)

 

The sequence alignment between target and template is the single most important factor that determines the quality of your model.

No comparative modeling process will repair an incorrect alignment; it is useful to consider a homology model rather like a three-dimensional map of a sequence alignment rather than a structure in its own right. In a homology modeling project, typically the largest amount of time should be spent on preparing the best possible alignment. Even though automated servers like the SwissModel server will align sequences and select template structures for you, it would be unwise to use these only because they are convenient. You should take advantage of the much more sophisticated alignment methods available. Analysis of wrong models can't be expected to produce right results.

The best possible alignment is usually constructed from a multiple sequence alignment that includes at least the target and template sequence and other related sequences as well. The additional sequences are an important aid in identifying the correct placement of insertions and deletions. Your alignment should have been carefully reviewed by you and wherever required, manually adjusted to move insertions or deletions between target and template out of the secondary structure elements of the template structure.

In the case of Mbp1 genes however, all orthologues we have considered have no indels in the APSES domain regions. Evolutionary pressure on the APSES domains has selected against indels in the more than 600 million years these sequences have evolved independently in their respective species.

Accordingly, all we need to do is to write the APSES domain sequences one under the other.

• Copy the FASTA formatted sequence for the APSES domain of your organism's Mbp1 orthologue from the sequences defined in Assignment 3 and save it as FASTA formatted text file. This is your target sequence. Compare this with the FASTA formatted file you have extracted from the PDB coordinate set. This is your template sequence. Then generate a multi-FASTA formatted file that contains both sequences, and pad the sequence(s) where required with hyphens as gap characters, so that target and template sequences have exactly the same length are aligned. Refer to the Fallback data if you are not sure about the format. (1 mark)

 
 

(2) Homology model

   

(2.1) SwissModel (1 mark)

 

Access the Swissmodel server at http://swissmodel.expasy.org . Navigate to the Alignment Interface.

 

• Paste your alignment for target and model into the form field. Refer to the Fallback Data file if you are not sure about the format. Make sure to select the correct option for the alignment input format on the form.

(You have to choose the correct format, and, if e.g. you choose a CLUSTAL format, you have to include a header line and a blank line. In the past we have seen problems with uploading alignments that have not been saved as "text only" and including periods i.e. "." in sequence names of CLUSTAL formatted alignments. Underscores appear to be safe.

• Click submit alignment and on the returned page define your target and template sequence. For the template sequence define the PDB ID of the coordinate file. Enter the correct Chain-ID.

Recently the PDB has undergone a "remediation" process in which archived coordinate files were altered by the database to conform to new format standards. One of the changes was to assign a chain identifier of "A" to all chains that did not previously have a chain identifier. SwissModel uses a derivative of coordinate sets from the PDB (a dataset they call ExPDB). Apparently the PDB proper and ExPDB have now gone out of synchrony; when I entered the (correct, according to PDB) chain designation "A" for 1MB1, SwissModel rejected the alignment with a nondescript error message. When I entered an underscore "_" instead, which would be the designation for a chain without explicit chain identifier, such as the pre-remidation versio of the coordinates, the alignment was accepted and processed. I have e-mailed SwissModel about the problem which may or may not be corrected while you are working on your assignments. If your template chain has the chain identifier "A" and your alignment gets rejected, enter an underscore instead.

Enter the correct chain ID into the form-field even if you think it already appears there, don't simply accept the preloaded default. There is a bug in SwissModel's parser code that may cause incorrect strings to be sent to the server from that field. I have e-mailed SwissModel about the problem which may or may not be corrected while you are working on your assignments.

• Click submit alignment and review the alignment on the returned page. Make sure it has been interpreted correctly by the server. The conserved residues have to be lined up and matching. Then click submit alignment again, to start the modeling process.

• The resulting page returns information about the resulting model. Save the model coordinates on your computer. Read the information on what is being returned by the server (click on the red questionmark icon). Paste the Anolea profile into your assignment.

Do not paste a screenshot of the result, but copy and paste the image from the Web-page! You do not need to submit the actual coordinate files with your assignment.

(1 mark)

 
In case you do not wish to submit the modelling job yourself, or have insurmountable problems when using the SwissModel interface, you may access the result files from the Fallback Data file. Document the problems and note this in your assignment.

(3) Model analysis

   

(3.1) The PDB file (1 mark)

 

Open your model coordinates in a text-editor (make sure you view the PDB file in a fixed-width font) and consider the following questions: (Alternatively, view the coordinates linked to the Fallback Data file.)

 

• What is the residue number of the first residue in the model? What should it be, based on the alignment? If the putative DNA binding region was reported to be residues 50-74 in the Mbp1 protein, which residues of the model correspond to that?

(1 mark)

[...]

   

(3.2) first visualization (2 marks)

 

In assignment 2 you have already studied the 1MB1 coordinate file and compared it to your organism's Mbp1 APSES domain, Since a homology model inherits its structural details from the template, the model should look very similar to the original structure but contain the sequence of the target.

 

• Save your model coordinates to your harddisk and visualize the structure in VMD. (Alternatively, copy and save the coordinates linked to the Fallback Data file to your harddisk.) Make an informative stereo view (see below), and paste it into your assignment.

• Discuss briefly which parts of the model may be unreliable and color these (if any) distinctly in your submitted image.

(2 marks)

 

Stereo-view of a subdomain within the 1MB1 structure that includes residues 36 to 76. The color gradient ramps from blue (36) to green (76).

   

(3.3) modeling a DNA ligand (4 marks)

 

One of the really interesting questions we can discuss with reference to our model is how sequence variation might be converted into changing DNA recognition sites, and then lead to changed cognate DNA binding sequences. But in order to address this, we would need to add a plausible model for how DNA is bound to APSES domains.

Since there is currently no software available that would accurately model such a complex from first principles, we will base a model of on homology modeling as well. This means we need to find a similar structure for which the complex structure is known. However, you may remember from the third assignment that the APSES domains in fungi seem to be a relatively small family. And there is no structure available of a protein-DNA complex. Now what?

Remember that homologous sequences can have diverged to the point where their sequence similarity is no longer recognizable, however their structure may be quite well conserved. Thus if we could find similar structures in the PDB, these might provide us with some plausible hypotheses for how DNA is bound by APSES domains. We thus need a tool similar to BLAST, but not for the purpose of sequence alignment, but for structure alignment. A kind of BLAST for structures. Very similar to using BLAST, we might not want to search with the entire protein, if all we are interested in is a subdomain that binds to DNA. Attempting to match all structural elements in addition to the ones we are actually interested in is likely to make the search less specific - we would find false positives that are similar to some irrelevant part of our structure. However, defining too small of a subdomain would also lead to a loss of specificity: in the extreme it is easy to imagine that the search for e.g. a single helix would retrieve very many hits that would be quite meaningless. The arrangement of the residues from 50 to 74 that we have already discussed in Assignment 2 suggests that the compact subdomain from 36 to 76 (see the image above) might be a useful structure to search with: it contains the residues we are interested in and enough of connected secondary structure elements to be structurally meaningful.

At the NCBI, VAST is provided a search tool for structural similarity search.

At the EBI there are a number of very well designed structure analysis tools linked off the Structural Analysis page. As part of its MSD Services, MSDfold provides a convenient interface for structure searches.

However we have also read previously that the APSES domains are members of a much larger superfamily, the "winged helix" DNA binding domains (CATH 1.10.10.10), of which hundreds of structures have been solved. They represent one branch of the tree of helix-turn-helix (HTH) DNA binding modules. (A recent review on HTH proteins is linked from the resources section at the bottom of this page). Winged Helix domains typically bind their cognate DNA with a "recognition helix" which precedes the beta hairpin and binds into the major groove; additional stabilizing interactions are provided by the edge of the beta strand binding into the minor groove. This is good news: once we have determined that the APSES domain is actually an example of a larger group of transcription factors, we can pick one of these for which a DNA complex structure is known. CATH does not provide information on complexes, but we can serach the PDB with CATH codes in the following way:

• Navigate to the PDB home page and follow the link to Advanced Search
• In the options menu for Choose a Query Type select Structure Features → CATH classification. A window will open that allows you to navigate down through the CATH tree. The interface is awkward because it does not display the actual CATH codes along with the class names, but you can view the class names on the CATH page linked above. Click on tthe arrows before "Mainly Alpha"→"Orthogonal Bundle"→"ARC repressor mutant, subunit A" then click on the link to "winged helix repressor DNA binding domain". As of this writing, this subquery matches 295 structures.

Click on the (+) button behind the subquery to add an additional query. Select the option "Structure Summary"→"Moelcule / Chain type". In the option menus that pop up, select "Contains Protein → Yes", "Contains DNA → Yes""Contains RNA → Ignore". This selects files that contain Protein-DNA complexes.

• Check the box below this subquery to "Remove Similar Sequences at 90% identity" and click on "Evaluate Query". As of this writing, seventy complexes were returned.
• In the left-hand menu, under the Tabulate section, click on the "Collage" function to display icons of the structure files. This is a fast way to obtain an overview of the structures that have been returned. First of all you may notice that in fact not all of the structures are really different, despite [...]

This is good news: once we have determined that the APSES domain is actually an example of a larger group of transcription factors, we can pick one of these for which a DNA complex structure is known. I have picked one such structure from the list of hits that were returned by SSM: it is the Elk-1 transcription factor.

1DUX). Note how the "recognition helix" inserts into the major groove of the DNA molecule. The color gradient ramps from blue (34) to green (84). Note how the first helix of the "helix-turn-helix" architecture serves only to position the recognition helix and makes few interactions by itself.

Now all that is left to do is to bring the DNA molecule into the correct orientation for our model and then to combine the two files. We need to superimpose the Elk-1 protein/DNA complex onto our model.

Structure superposition

There are quite a number of superposition servers available on the Web, a remarkably comprehensive overview can be found in Wikipedia. However, overengineering and black-box mentality makes our task more difficult than it need be: most tools do not allow users to specify particular alignment zones but attempt to automatically define the zones of residues to be supoerimposed according to some geometric target function. Almost none return the actual rotation matrix and translation vector that is used for the superposition. And almost none transform the coordinates of heteratoms such as solvent, ligands or DNA molecules along with the protein coordinates. An exception that I have found to be very useable is the Local-Global Alignment server (LGA), written by Adam Zemla. The procedure is quite straightforward:

• Define the structure to be rotated (1DUX in this case). This is a dimer, so download the file from the PDB and manually edit to contain only DNA chains A and B and protein chain C.
• Define the structure to be held constant (1MB1 in this case). Download from PDB.
• Use the "browse" option to define both files as input on the LGA inpput form
• Use the option to have both coordinate sets included in your output: -o2
• Submit

The results arrive per e-mail. I have linked the resulting PDB file to the Fallback Data page. If you run this analysis on your own, you may want to review the types of edits the edits I made to the PDB file to get it displayed correctly in Rasmol.

 

• Save the superimposed coordinates in a file, open and view in Rasmol and note how well the "recognition helix" and adjacent beta strands superimpose! (Alternatively, copy and save the coordinates from the c to your harddisk.) Make an informative view, divergent stereo and paste it into your assignment. (4 marks)

 
 

(4) Summary of Resources

 

Links

• Review (PDF, restricted) Manuel Peitsch on Homology Modeling
• Review (PDF, restricted) Aravind et al. Helix-turn-helix domains (background reading, not required reading)
• Assigned Organisms
• PDB file format
• Wikipedia on Structural Superposition (although the article is called "Structural Alignment")

• Fallback Data page

Alignments

• Mbp1 proteins T-Coffee aligned (text version)

   

[End of assignment]

If you have any questions at all, don't hesitate to mail me at boris.steipe@utoronto.ca or post your question to the Course Mailing List

<Tasks: review location of fallback files; rewrite SwissModel interface section ...>