Week 6 and a bit of 7

Can't believe I'm already in week 7!

Well I finally have a program which can give some of the rotamer populations for the ILE side chains and some of the correlations look really nice! Though most are just noise. I'm now parameterising it to hopefully make it slightly more useful as a predictive tool ...

Currently having difficulty trying to find a way to combine a population which is defined by itself and then also defined as part of a group of rotamers ... this could be interesting

Week 5! (again a little late)

Well I have results! Hopefully I will be able to link these theoretical chemical shifts to ILE. It seemed like I can define the whole system but then I looked at the error (for now given as the standard deviation) and found out that some of the changed were not 'statistically significant'. oh well at least some of it. Now to find out if I can reasonably work out the populations of the ILE rotamers from the chemical shifts!

I seem to find coding python easier despite still writing disgusting code. It is becoming increasingly apparent that just because you can do maths doesn't mean that coding things will be a breeze. Still really good fun and quite different from wet lab stuff.

also have have i managed to finish 5 weeks already?

Week 4

 So the calculations for the magnetic shielding for each of the ILE rotamers are still running (3 days later) and I have finally managed to work out how to get the reference shielding to then calculate the the chemical shift.

A some what uneventful week just writing scripts to extract all the data from the log files.

week 3 but a bit late!

Well I am beginning to feel more confident in python and decided to implement the alignment of two structures as was done by W. Kabsch (1976). Nice bit of maths and also a fun introduction to eigenvectors! Though wasted two days on learning that np.eig wont work on symmetric matrices and one should use .eigh,.

Well at least I now have my alignments for the ILE rotamers  and calculated all my RSMDs to compare the geometry optimized structures against the ones I produced earlier.

So I found the alignment looked best (example shown below made using Pymol) when the alignment was done using only the fixed atoms to generate the rotation matrix.

Week 2 Many plots!

This week I have managed to make a Ramachandran plot for ILE using the top500 pdb files. It seems like the majority of the ILE phi psi angles fall in the beta-sheet and alpha-helix regions with very few in the left hand alpha-helix region. This was expected as the left hand alpha-helix is much rarer.

Following this I also made plots of the chi 1 and chi 2 angles according to their secondary structure. This gave some interesting results showing that in alpha helix the  x1 tended be in distinct regions. There was also some interesting results as when the back bone was in the beta sheet conformation there was a population in the chi 2 in +100.

This is just one of the 4 plots! also made the for ILE residers depending on the backbone conformation.

(disclaimer the figure is not finished yet!)

Based on these plots I have taken a selection of angle combination for phi, psi, chi 1 and chi2 and generated the structures. I hope to be able to optimise them soon but it seems I still have trouble with writing the input files for Gaussian 09. Oh well at least I understand for loops in python better!

First Week!

      Started putting together a series of structures for ILE using python ... I think I will never understand how for loops and dictionaries work! Oh well its good fun and sort of shows how useful computing can be since now its easy to produce any number of ILE rotamers.

Towards the end of this week I managed to measure the phi/psi angles for the 'top500' PDB files. Hopefully these will be able to give me an indication of the backbone angles for that are populated so that I can make chi 1 and ch 2 rotamers for each one.