390 likes | 521 Vues
This document delves into the PLAPA method for modeling phosphate localization and energy minimization in proteins. It outlines the step-by-step process including the surface generation tool usage, phosphate ball placement, energy calculations, and ranking of positions. The results showcase the evaluation of five test cases across different protein domains, highlighting challenges faced and modifications made to optimize outcomes. The future directions aim to explore unresolved examples and improve structural alignment methods for more accurate prediction and modeling of protein structures.
E N D
PLAPA Phosphate Localization And Prediction Algorithm
How does it work? • Generates a surface using Matt’s surface generation tool in PLOP • Places a phosphate “ball” of -2 charge at each point • Minimizes the position of the ball • Calculates energy • Ranks all positions and chooses the lowest energy point – no averages among points
What Have I Done? • Used five test cases to evaluate different strategies • 1gxc – FHA domain • 1g6g – FHA domain • 1fpr – SH2 domain • 1bm2 – SH2 domain • 1een – SH2 domain
What Have I Done? • Started by using a phosphate ball • This did not work at either of the trial radii of 1.4, 2.0, and 3.0 A, or the OPLSAA phosphate radius of 4.83 A • Why?? (All following examples with 1BM2 only, Native is always BLUE, Decoy RED)
Radius vs Distance from Native Post Minimization Distance from Native Radius of Sphere in A
Energy of Native vs Decoy at Different Radii Energy in Kcals Radius of Sphere in A
Energy of Native vs Decoy After Subtraction of LJ Energy in Kcals Radius of Sphere in A
Energy of Native vs Decoy Using only Electrostatics Energy in Kcals Radius of Sphere in A
1G6G Fix • Added capping group (ACE) • Lowest energy moved to within 5 A of native phosphate position in native vs decoy trials
Final Output for 1BM2 Energy in Kcals Distance from Native in A
Future Directions • Investigate the two examples that are not working, haven’t had time due to CASP…
CASP Target T0196 Possible Elongation Factor from Pyrococcus furiosus
Initial Alignment * Poly His at the beginning of a sequence is usually a His tag…
Ok, so it doesn’t quite match • Genomic search – is it only a piece? • Bacterial gene prediction • Operon search (functional proteins are grouped onto operons in prokaryotes) • Similar operons = more sequences for the alignment • It’s a real protein…
Structural Alignment • Schrodinger Structural Alignment Tool • ~jacobson/bin/dinger/utilities/structalign pdb1 pdb2 pdbN • Must use PDB files with only one chain • ~jacobson/bin/dinger/utilities/getpdb pdbID:chain (CAPS) • BUGS • Rotation matrix a bit screwy for some pdbs • Some of your pdbs will appear in alignment but will not appear as files at the end of the run • This can be fixed by creating a .fasta file with the structural alignment, loading into chimera, and aligning the structures
Structural Alignment • Whole protein structural alignment performed with CE and Schrodinger • CE provided a good alignment, but output pdb file useless since the other domains influenced position • Schrodinger script could not handle this • Domain only alignment based on sequence alignment • CE could not do this • Schrodinger could…with bugs
Interlude • Here is where I bitch about how chimera crashed my system about 5 times so I downloaded the linux version, in which I could not figure out how to change the background to white…so its black.
Different Partial Templates Different Full Templates
Further Directions • Build longer models • Look over surface of final model to see whether there is a hydrophobic patch for helix packing • Get your suggestions…
Brought to you by… People like you!