410 likes | 751 Views
Applications of Homology Modeling. Hanka Venselaar. This seminar…. Homology Modeling… What? Why? When? How? And a few real world examples…. EEC syndrome. No structure:. ?. Sequence:.
E N D
Applications of Homology Modeling Hanka Venselaar
This seminar…. Homology Modeling… • What? • Why? • When? • How? • And a few real world examples….
EEC syndrome No structure: ? Sequence: MSQSTQTNEFLSPEVFQHIWDFLEQPICSVQPIDLNFVDEPSEDGATNKIEISMDCIRMQDSDLSDMWPQYTNLGLLNSMDQQIQNGSSSTSPYNTDHAQNSVTAPSPYAQPSSTFDALSPSPAIPSNTDYPGPHSFDVSFQQSSTAKSATWTYSTELKKLYCQIAKTCPIQIKVMTPPPQGAVIRAMPVYKKAEHVTEVVKRCPNHELSREFNEGQIAPPSHLIRVEGNSHAQYVEDPITGRQSVLVPYEPPQVGTEFTTVLYNFMCNSSCVGGMNRRPILIIVTLETRDGQVLGRRCFEARICACPGRDRKADEDSIRKQQVSDSTKNGDGTKRPFRQNTHGIQMTSIKKRRSPDDELLYLPVRGRETYEMLLKIKESLELMQYLPQHTIETYRQQQQQQHQHLLQKQTSIQSPSSYGNSSPPLNKMNSMNKLPSVSQLINPQQRNALTPTTIPDGMGANIPMMGTHMPMAGDMNGLSPTQALPPPLSMPSTSHCTPPPPYPTDCSIVSFLARLGCSSCLDYFTTQGLTTIYQIEHYSMDDLASLKIPEQFRHAIWKGILDHRQLHEFSSPSHLLRTPSSASTVSVGSSETRGERVIDAVRFTLRQTISFPPRDEWNDFNFDMDARRNKQQRIKEEGE EEC syndrome
Homology modeling in short… Prediction of structure based upon a highly similar structure • 2 basic assumptions: • Structure defines function • During evolution structures are more conserved than sequence Use one structure to predict another
Homology modeling % identity * O # residues * Actually, modelling is possible, but we cannot get an alignment… Example: by 80 residues 30% identity sufficient
NSDSECPLSHDG || || | || NSYPGCPSSYDG NSDSECPLSHDG ? Model sequence Unknown structure Known structure Back bone copied Homology modeling in short… Prediction of structure based upon a highly similar structure Model! Copy backbone and conserved residues Add sidechains, Molecular Dynamics simulation on model Known structure
NSDSECPLSHDGYCLHDGVC || || | ||||| ||| NSYPGCPSSYDGYCLNGGVC 1: Template recognition and initial alignment • BLAST your sequence against PDB • Best hit normally template • Initial alignment
2: Alignment correction 1: Template recognition and initial alignment
CPISRTAAS-FRCW CPISRTG-SMFRCW CPISRTA--TFRCW CPISRTAASHFRCW CPISRTGASIFRCW CPISRTA---FRCW CPISRTGASIFRCW CPISRTGASIFRCW CPISRTA---FRCW CPISRT---AFRCW Correct alignment 2: Alignment correction • Functional residues conserved • Use multiple sequence alignments • Deletions shift gaps Multipe sequence alignment Sequence with known structure Your sequence Both are possible
E I E E V V A P C C C S R R M R G L M P P 2: Alignment correction -A-V F-D- • Core residues conserved • Use multiple sequence alignments • Deletions in your sequence shift gaps Known structure FDICRLPGSAEAV Model FNVCRMP---EAI Model FNVCR---MPEAI Correct alignment
3: Backbone generation 2: Alignment correction 1: Template recognition and initial alignment
3: Backbone generation • Making the model…. • Copy backbone of template to model • Make deletions as discussed • (Keep conserved residues)
4: Loop modeling 2: Alignment correction 1: Template recognition and initial alignment 3: Backbone generation
4: Loop modeling Known structure GVCMYIEA---LDKYACNC Your sequence GECFMVKDLSNPSRYLCKC Loop library, try different options
5: Sidechain modeling 2: Alignment correction 1: Template recognition and initial alignment 3: Backbone generation 4: Loop modeling
5: Side-chain modeling • Several options • Libraries of preferred rotamers based upon backbone conformation
6: Model optimization 2: Alignment correction 1: Template recognition and initial alignment 3: Backbone generation 4: Loop modeling 5: Sidechain modeling
6: Model optimization • Molecular dynamics simulation • Remove big errors • Structure moves to lowest energy conformation
2: Alignment correction 1: Template recognition and initial alignment 3: Backbone generation 4: Loop modeling 5: Sidechain modeling 7: Model validation 6: Model optimization
7: Model Validation • Second opinion by PDBreport /WHATIF • Errors in active site? new alignment/ template • No errors? Model!
2: Alignment correction 1: Template recognition and initial alignment 3: Backbone generation 4: Loop modeling 8: Iteration 8: Iteration 5: Sidechain modeling 8: Iteration 8: Iteration 7: Model validation 6: Model optimization
2: Alignment correction 1: Template recognition and initial alignment 3: Backbone generation 4: Loop modeling 8: Iteration 8: Iteration 5: Sidechain modeling Model! 8: Iteration 8: Iteration 7: Model validation 6: Model optimization
Alignment Modeling Correction 8 steps of homology modeling 1: Template recognition and initial alignment 2: Alignment correction 3: Backbone generation 4: Loop modeling 5: Side-chain modeling 6: Model optimization 7: Model validation 8: Iteration
EEC syndrome Structure! P63 MSQSTQTNEFLSPEVFQHIWDFLEQPICSVQPIDLNFVDEPSEDGATNKIEISMDCIRMQDSDLSDMWPQYTNLGLLNSMDQQIQNGSSSTSPYNTDHAQNSVTAPSPYAQPSSTFDALSPSPAIPSNTDYPGPHSFDVSFQQSSTAKSATWTYSTELKKLYCQIAKTCPIQIKVMTPPPQGAVIRAMPVYKKAEHVTEVVKRCPNHELSREFNEGQIAPPSHLIRVEGNSHAQYVEDPITGRQSVLVPYEPPQVGTEFTTVLYNFMCNSSCVGGMNRRPILIIVTLETRDGQVLGRRCFEARICACPGRDRKADEDSIRKQQVSDSTKNGDGTKRPFRQNTHGIQMTSIKKRRSPDDELLYLPVRGRETYEMLLKIKESLELMQYLPQHTIETYRQQQQQQHQHLLQKQTSIQSPSSYGNSSPPLNKMNSMNKLPSVSQLINPQQRNALTPTTIPDGMGANIPMMGTHMPMAGDMNGLSPTQALPPPLSMPSTSHCTPPPPYPTDCSIVSFLARLGCSSCLDYFTTQGLTTIYQIEHYSMDDLASLKIPEQFRHAIWKGILDHRQLHEFSSPSHLLRTPSSASTVSVGSSETRGERVIDAVRFTLRQTISFPPRDEWNDFNFDMDARRNKQQRIKEEGE EEC syndrome
Arginine • Loss of negative charge • Loss of interaction with the DNA Mutation RS Serine
Another real world example: Mutation analysis HFE
Transferrin receptor (dimer) binds iron/transferrin complex HFE – complex: -Signaling and regulation of iron in bloodstream. -Expressed in liver and colon. -Mutations cause iron deposition disease “Hereditary Hemachromatosis“ HFE β2-microglobulin Facilitates trafficking of HFE to the cellmembrane
Hereditary Hemachromatosis 3 occuring mutations • C280Y • D41H • L161P L161P D41H C280Y
Loss of cystein bridge • Disturbing of β2-microglobulin binding domain • No trafficking to membrane Mutation C260Y
Introduction of additional negative charge • Disturbing of hydrogen bridges • Loss of stability in this area Mutation H41D
Loss hydrophobic interactions • Major disturbance of the helix • Less interaction of the helix with the transferrin receptor Mutation L161P
Seriousness of mutation Seriousness of the disease D41H L161P C260Y D41H L161P C260Y Conclusion: the seriousness of the mutation is related to the seriousness of the disease and can be explained by analyzing the mutations with the 3D structure.
To conclude…. Homology Modeling… • What? Prediction of an unknown structure based on an homologous and known structure • Why? To answer biological and medical questions when the “real” structure is unknown • When? A template with enough identity must be available • How? 8 Steps • Real world examples: mutations in EEC syndrome and HFE can be explained