- By
**Anita** - Follow User

- 318 Views
- Uploaded on

Download Presentation
## PowerPoint Slideshow about 'Protein Tertiary Structure Comparison' - Anita

Download Now**An Image/Link below is provided (as is) to download presentation**

Download Now

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript

### Protein Tertiary Structure Comparison

Dong Xu

Computer Science Department

271C Life Sciences Center

1201 East Rollins Road

University of Missouri-Columbia

Columbia, MO 65211-2060

E-mail: xudong@missouri.edu

573-882-7064 (O)

http://digbio.missouri.edu

Lecture Outline

- Why structural alignment
- Technical definition
- SSAP
- DALI
- Fast search
- Protein families

Structure Is Better Conserved during Evolution

Structure can adopt a

wide range of mutations.

Physical forces favor

certain structures.

Concept of fold.

Number of fold is limited.

Currently ~1000

Total: 1,000s ~10,000s TIM barrel

Alignment of Protein Structure

- Three-dimensional structure of one protein compared against three-dimensional structure of second protein
- Atoms (protein backbones) fit together as closely as possible to minimize the average deviation

Why Align Structures? (1)

- Additional measure of protein similarity
- Structure generally preserved better than sequence over the course of evolution
- Provide more information on the relationship between proteins than what sequence alignment can offer
- Allows classification of proteins based on structural similarities

Why Align Structures? (2)

- Basis for protein fold identification (prediction)
- Sometimes sequence similarity between two proteins exists, but is not strong enough to produce an unambiguous alignment (gold standard for sequence comparison).
- Pinpoint the active sites more accurately.
- Allows identification of common sub-structures of interest

Why Align Structures? (4)

Illustrate interesting evolutionary/functional relationship between proteins:

Two ferredoxins, 1DOI and

1AWD, are aligned structurally,

showing an insertion in 1DOI

that contains potassium-ion

binding sites. This may be the

result of adaptations to the high

salt environment of the Dead Sea.

Lecture Outline

- Why structural alignment
- Technical definition
- SSAP
- DALI
- Fast search
- Protein families

T

Find a transformation

to achieve the best

superposition

Simple case – two closely related proteins with the same number of amino acids.

Types ofStructure Comparison

- Sequence-dependent vs. sequence-independent structural alignment
- Global vs. local structural alignment
- Pairwise vs. multiple structural alignment

Sequence-dependent Structure Comparison (1)

Given two sets of 3-D points :

P={pi}, Q={qi} , i=1,…,n;

rmsd(P,Q) = √ S i|pi - qi |2 /n

(root mean square deviation)

Find a 3-D rigid transformation T* such that:

rmsd( T*(P), Q ) = minT√ S i|T(pi) - qi |2 /n

Sequence-dependent Structure Comparison (2)

2

2

6

6

5

5

7

7

1

1

3

3

4

4

2

2

1

1

4

4

5

5

7

7

3

3

6

6

1234567

ASCRKLE

¦¦¦¦¦¦¦

ASCRKLE

Minimize rmsd

of distances 1-1,...,7-7

Sequence-dependent Structure Comparison (3)

- Can be solved in O(n) time.
- Useful in comparing structures of the same protein solved in different methods, under different conformation, through dynamics.
- Evaluation protein structure prediction.

Sequence-independent Structure Comparison

T

Correspondence is Unknown!Given two configurations of points in the three dimensional space,

find T which produces “largest” superimpositions of corresponding 3-D points.

Order-Dependent vs. Order-Independent Comparison

Alignment (order dependent): a correspondence between elements of two sequences with order (topology) kept (typical structural alignment)

bipartite matching (order-independent): one-to-one matching

residues

of protein

sequence

FESYTTHRGHR

:::::::: ::

FESYTTHRPHR

FSEYTTHRGHR

: ::::: ::

FESYTTHRPHR

Evaluating Structural Alignments

1. Number of amino acid correspondences created.

2. RMSD of corresponding amino acids

3. Percent identity in aligned residues

4. Number of gaps introduced

5. Size of the two proteins

6. Conservation of known active site environments …

No universally agreed upon criteria. It depends on what you are using the alignment for.

Output

1ABR:B - ABRIN-A

1BAS:_ - BASIC FIBROBLAST GROWTH FACTOR (BFGF)

Seq. identity = 10% RMSD = 1.9Å

Lecture Outline

- Why structural alignment
- Technical definition
- SSAP
- DALI
- Fast search
- Protein families

How to recognize structural similarities

- By eye (SCOP)
- 2. Algorithmically
- point-based methods use properties of points (distances) to establish correspondence
- Dynamic programming (SSAP)
- Distance matrix (DALI)
- secondary structure-based methods use vectors representing secondary structures to establish correspondences (LOCK).
- Image processing based method.

Structural Comparison Algorithms

- Due to the high compute complexity, practical algorithms rely on heuristics
- Fully automated structure analysis has not been as successful as analyses with human intervention in taking in to account the biological implications

SSAP

- SSAP: Secondary Structure Alignment Program
- Incorporates double dynamic programming to produce a structural alignment between two proteins

Basic Ideas of SSAP

The similarity between residue i in molecule A and residue k in molecule B is characterised in terms of their structural surroundings

This similarity can be quantified into a score, Sik

Based on this similarity score and some specified gap penalty, dynamic programming is used to find the optimal structural alignment

Scoring Function of SSAP (1)

Distance between residue i & j in molecule A ; dAi,j

Similarity for two pairs of residues, ij in A & kl in B ;

a,b constants

l

j

i

k

Scoring Function of SSAP (2)

Similarity between residue i in A and residue k in B ;

Si,k is big if the distances from residue i in A to the 2n nearest neighbours are similar to the corresponding distances around k in B

HSERAHVFIM..

GQ-VMAC-NW..

A :

B :

k=4

Alignment Gaps in SSAPThis works well for small structures and local structural alignments - however, insertions and deletions cause problems unrelated distances

The actual SSAP algorithm uses Dynamic programming on two levels, first to find which distances to compare Sik, then to align the structures using these scores

Steps in SSAP (1)

- 1)Calculate vectors from C of one amino acid to set of nearby amino acids
- Vectors from two separate proteins compared
- Difference (expressed as an angle) calculated, and converted to score
- 2)Matrix for scores of vector differences from one protein to the next is computed.

Steps in SSAP (2)

- 3) Optimal alignment found using global dynamic programming, with a constant gap penalty
- 4)Next amino acid residue considered, optimal path to align this amino acid to the second sequence computed

Steps in SSAP (3)

- 5) Alignments transferred to summary matrix
- If paths cross same matrix position, scores are summed
- If part of alignment path found in both matrices, evidence of similarity

Steps in SSAP (4)

- 6) Dynamic programming alignment is performed for the summary matrix
- Final alignment represents optimal alignment between the protein structures
- Resulting score converted so it can be compared to see how closely related two structures are

Lecture Outline

- Why structural alignment
- Technical definition
- SSAP
- DALI
- Fast search
- Protein families

Distance Matrix Approach

- Uses graphical procedure similar to dot plots
- Identifies residues that lie most closely together in three-dimensional structure
- Two sequences with similar structure can have dot plots superimposed

Distance Matrix

- Similar 3D structures have similar inter-residue distances

DALI

- Distance Alignment Tool (DALI)
- Uses distance matrix method to align protein structures
- Assembly step uses Monte Carlo simulation to find submatrices that can be aligned

Structural Analysis Algorithms – DALI (1)

- DALI is based on distance matrices – 2D matrices containing all pairwise distances between points of a molecule
- Distance matrices of two molecules are compared to find regions of similar patterns of distances, which indicate similarities in their 3D structure
- Key algorithm steps:
- Divide distance matrices into overlapping sub-matrices of fixed size
- Search through two matrices (of two molecules) to find similar patterns
- Assemble matching pairs of sub-matrices in to larger sets to maximize their similarity score

Structural Analysis Algorithms – DALI (2)

- Assembly of aligned sub-matrices is done using a Monte Carlooptimization
- Monte Carlo optimization is an iterative improvement by a random walk exploration of the search space, with occasional excursions in to non-optimal territory (i.e. occasionally, a move that reduces the overall score is carried out)
- The occasional non-optimal moves help avoid getting “trapped” in local optima of the score function, improving the chance of finding the global optimum

Lecture Outline

- Why structural alignment
- Technical definition
- SSAP
- DALI
- Fast search
- Protein families

Fast Structural Similarity Search

- Compare types and arrangements of secondary structures within two proteins
- If elements similarly arranged, three-dimensional structures are similar
- LOCK, VAST and SARF are programs that use these fast methods

Structural Analysis Algorithms – LOCK

- Both SSAP and DALI deal only with points (atoms) of the molecules
- LOCK uses a hierarchical approach
- Larger secondary structures such as helixes and strands are represented using vectors and dealt with first
- Individual residues are dealt with afterwards
- Assumes large secondary structures provide most stability and function to a protein, and are most likely to be preserved during evolution

LOCK Algorithm

- Key algorithm steps:
- Represent secondary structures as vectors
- Obtain initial superposition by computing local alignment of the secondary structure vectors (using dynamic programming)
- Compute residue superposition by performing a greedy search to try to minimize root mean square deviation (a RMS distance measure) between pairs of nearest backbone atoms from the two proteins
- Identify “core” (well aligned) atoms and try to improve their superposition (possibly at the cost of degrading superposition of non-core atoms)
- Steps 2, 3, and 4 require iteration at each step

ProteinDBS

Shyu, Chi, Scott, Xu. Nucleic Acid Research. 32, W572 - CW575, 2004

Comparison between different methods

- CATH
- Fully automated
- SSAP
- SCOP
- Based on subjective interpretation of evolutionary history of proteins
- FSSP
- DALI
- Agreement between CATH and SCOP may be at most 60%.
- FSSP vs CATH 40%
- FSSP vs SCOP 60%

Lecture Outline

- Why structural alignment
- Technical definition
- SSAP
- DALI
- Fast search
- Protein families

Structure Families (1)

Homologous family: evolutionarily related with a significant sequence identity;

Superfamily: different families whose structural and functional features suggest common evolutionary origin;

Fold: different superfamilies having same major secondary structures in same arrangement and with same topological connections (energetics favoring certain packing arrangements);

Class: secondary structure composition.

6 Classes of Protein Structures (1)

1) Class : bundles of helices connected by loops on surface of proteins

2) Class : antiparallel sheets, usually two sheets in close contact forming sandwich

3) Class /: mainly parallel sheets with intervening helices; may also have mixed sheets (metabolic enzymes)

6 Classes of Protein Structures (2)

4) Class + : mainly segregated helices and anti-parallel sheets

5) Multi-domain ( and ) proteins more than one of the above four domains

6) Membrane and cell-surface proteins and peptides excluding proteins of the immune system

Reading Assignments

- Suggested reading:
- Contemporary approaches to protein structure classification. Mark B. Swindells, et al. BioEssay. Volume 20, Issue 11, 1998, Pages: 884-891
- Optional reading:
- The structural alignment between two proteins: Is there a unique answer? Adam Godzik, Protein Science (1996), 5 1325-1338
- Protein Structure Similarities. Patrice Koehl, Current Opinions in Structural Biology (2001), 11 348-353

Develop a program that can perform protein structural alignment using SSAP:

The Ca coordinates of two proteins (A and B) of will be sent to the mailing list

Calculate the similarity matrix between residue i in A and residue k in B (let n = 4, a = b = 1):

Perform dynamic programming on Si,k, and retrieve the alignment to print out.

Project AssignmentProject Phase III Report

- Due on 11/17, send me through email
- Write on top of Phase II report.
- 7-30 Pages
- As a draft of the final report
- Free style in writing (use 11pt font or larger)
- Present key results
- Software implementation
- Benchmark (computing time)
- Computational data
- Interpret the meaning of the data

Download Presentation

Connecting to Server..