Making Dihedral Restraints from Chemical Shifts
Introduction
Next we will generate restraints in a different manner; dihedral restraints from backbone chemical shifts. We will be using a program called DANGLE (Dihedral ANgles from Global Likelihood Estimates) which is embedded within Analysis. DANGLE estimates dihedral angles from chemical shifts in a similar manner to TALOS; i.e. it matches a chemical shift & sequence query to a structural database of known PHI/PSI angles and chemical shifts. However, DANGLE uses a different (Bayesian) method to produce an angle estimate and tolerance, compared to TALOS. The idea is to use Bayesian inference to infer what the range of likely PHI/PSI angles might be (using the chemical shifts) by checking all PHI/PSI combinations in 10 degree square bins to see how well such angles can be used to explain the data. Such an analysis allows for the user to see uncertainties in the angle predictions, including where the chemical shift to structure mapping is redundant and there are multiple regions in the Ramachandran plot which could explain the chemical shift data.
DANGLE
To run DANGLE select M: Structure: DANGLE: Predict Dihedrals. Note that at the top that the Chain should be set as "GI:A", the Shift List as "ShiftList 2:2" and "Max No. of Islands:" as 2. This simply specifies which data to use and how strict the analysis should be. Using two islands means that we will reject predictions that result in more than two discrete regions of the Ramachandran plot. To start the analysis press [Run Prediction] and accept "Run1" as the name for the job by pressing [OK] at the opportune moment. Please be aware that DANGLE will take several minutes to finish the calculation.
Once the calculation is over you will see the main table filled in with PHI and PSI backbone dihedral angle predictions and their associated error ranges:
Further, if you select a row in the main table you will see a plot in Ramachandran (PHI/PSI) space of where the likely angles are deemed to be. Click on the "7 Ser" row and note that there is a lot of red colour in the chart, indicating that DANGLE was not able to make a distinct choice of PHI/PSI: you should not use such a prediction in a structure calculation. Click on the [Next] button to get to "8 Lys". The prediction for this residue is somewhat better, and you could use this in a structure calculation (it has one discrete region) although the error bounds for such a dihedral restraint would be suitably large. Click on the "12 Glu" row. This residue has a very precise range of predicted PHI/PSI angles. Such a residue could be used in a structure calculation with a high degree of confidence and proportionately narrow error margins.
Note that DANGLE also predicts the secondary structure of the residues, but that this calculation is not made from the angles, but directly from the measured secondary structures in the shift-structure database. To make the restraints themselves set the Restraint Set to "1", which will place the PHI/PSI dihedral restraints with our existing distance restraints and press [Commit Constraints]. View the generated restraints by going to M: Structure: Restraints & Violations,{Restraints}. Note that if you have a structural model for your protein you can see how the model's angles match with the DANGLE prediction.
Previous: Making Distance Restraints