 Untitled.png (105.95 KiB) Viewed 5821 times
Help with GNOM/Gasbor
Help with GNOM/Gasbor
Can someone walk me through the steps of getting data into and out of GNOM and into Gasbor. The samples are of HSA 5mg/ml taken in a xeuss 2.0 beamline at 1209 sd, high resolution for 1200 seconds. The files were circle gathered using foxtrot and then normalized and background subtracted using mantid plot. Now I have .dat files and would like to do some abinitio modeling, the probelm is that some of the samples are at higher temperatures, and as I understand Dammif is not very good with unfolded systems and Gasbor is more appropriate.
Use PRIMUSqt and GASBOR online
1. Make sure you have the latest ATSAS version installed.
2. Open your backgroundsubtracted data in PRIMUS (on Windows make sure you use primusqt; it is called "SAS Data Analysis" in the Start menu).
3. In the bottom right panel go to the 'Analysis' tab and click Distance Distribution'. This will start the GNOM GUI.
3.1. If PRIMUS was unable to automatically estimate R_{g} from the Guinier approximation there will be an additional step where you have to choose a Guinier region manually.
4. Increase the last point of 'Range' such that the whole experimental data are displayed (you can look up the number of data points in the parent window under 'To').
5. If necessary  adjust 'D_{max}'; make sure your p(r) is smooth, positive and fits the experimental data.
6. 'Save' your p(r) and the fit to the experimental data in a file with a *.out extension.
7. Go to GASBOR online.
Run GASBOR 510 times to make sure that you always get about the same model and all your models fit your input data. To do so go to My Projects and click 'rerun with new id' several times.
9. If the fit is good and the model makes sense  deposit it to SASBDB.
Check out the EMBO practical course on solution scattering from biological macromolecules 19  26 November 2018  Hamburg, Germany
2. Open your backgroundsubtracted data in PRIMUS (on Windows make sure you use primusqt; it is called "SAS Data Analysis" in the Start menu).
3. In the bottom right panel go to the 'Analysis' tab and click Distance Distribution'. This will start the GNOM GUI.
3.1. If PRIMUS was unable to automatically estimate R_{g} from the Guinier approximation there will be an additional step where you have to choose a Guinier region manually.
4. Increase the last point of 'Range' such that the whole experimental data are displayed (you can look up the number of data points in the parent window under 'To').
5. If necessary  adjust 'D_{max}'; make sure your p(r) is smooth, positive and fits the experimental data.
6. 'Save' your p(r) and the fit to the experimental data in a file with a *.out extension.
7. Go to GASBOR online.
 'GNOM file': choose your *.out file;
 'Angular units': choose inverse nanometres (1/nm) or inverse angstroms (1/A).
 'Symmetry': if the molecular weight estimated from the SAXS data points to a monomer  choose 'P1'.
 'Anisometry': if you think your protein is rather elongated than flat  choose 'Prolate', otherwise leave 'Unknown'.
 'Number of dummy residues': number of residues in the protein you've measured (including Histag etc.).
 'Mode': 'Reciprocal space'. (If you have evidence your protein is a monomer:dimer mixture in solution  you may try 'Symmetry: P2' and choose 'Mode: oligomermonomer mixture'.)
Run GASBOR 510 times to make sure that you always get about the same model and all your models fit your input data. To do so go to My Projects and click 'rerun with new id' several times.
9. If the fit is good and the model makes sense  deposit it to SASBDB.
Check out the EMBO practical course on solution scattering from biological macromolecules 19  26 November 2018  Hamburg, Germany
Re: Hel p with GNOM/Gasbor
Thank you this really helps, but i still have questions. In the distance distribution analysis tab I have total quality estimate which I guess is a percentage value. Guinier Rg/I(0) , p(r) Rc/I(0) and porod volume, what are these values and what should they be. For system I'm using a Xeuss 2.0 by xenocs in High resolution. How do I know the range and what system im dealing with, I know I'm using HSA. What does Rmin=0 and Rmax=0 mean and points and Alpha. I know it's a lot of questions but i would greatly appreciate it.
 Attachments

 Untitled.png (40.33 KiB) Viewed 5802 times
GNOM parameters
 Left plot log(I) vs. s: fit of the Fourier transform of p(r) (blue line) to the experimental SAS data (black dots).
The fit is computed like this: I(s) = 4Π ∫_{o}^{Dmax} p(r)sin(sr)/sr dr  Right plot p(r) vs. r: pair distance distribution function.
 Total quality estimate ranges between 1.0 (best) and 0.0 (worst).
 Guinier R_{g}/I(0): Radius of gyration and Forward scattering estimated from the fit in the left plot using Guinier approximation.
 p(r) R_{g}/I(0): Radius of gyration and Forward scattering estimated from the pair distance distribution function (right plot).
If the experimental data are in inverse angstroms (Å^{1}), R_{g} is provided in angstroms (Å). I(0) is provided in units of I(s).  Porod Volume: Excluded volume of the hydrated particle computed by DATPOROD.
If the experimental data are in inverse angstroms (Å^{1}), the volume is provided in cubic angstroms (Å^{3}).
For globular proteins Porod volume in cubic angstroms divided by 1.6 roughly corresponds to the molecular weight in daltons (Da).  System: Arbitrary Monodisperse  monodisperse systems of globular particles; GNOM can treat experimental data from other types of systems.
 Range: first and last points of the experimental data used by GNOM. By default the first point is the first point of the Guinier interval, the last point is s = 8/R_{g}.
 Setup: experimental setup, either Point Collimation (no smearing) or Slit Collimation.
 D_{max}: maximum intraparticle distance.
 Rmin=0: force r(0) = 0 (typically you want this).
 Rmax=0: force r(D_{max}) = 0 (typically you want this).
 Points: number of p(r) data points.
 Alpha: regularization parameter.
Larger alpha values yield a smoother p(r) but worse fit to the experimental data.
Smaller alpha values improve the fit but may result in a wobbly p(r).  Autognom: run DATGNOM and return to the initial state of this Wizard (if you've changed any values).
Re: Help with GNOM/Gasbor
Thank you AL,
You have no idea how much this has helped.
But I have a question, is it absolutely necessary to "Increase the last point of 'Range' such that the whole experimental data are displayed". My data is pretty noisy at high q. What do you recommend?
You have no idea how much this has helped.
But I have a question, is it absolutely necessary to "Increase the last point of 'Range' such that the whole experimental data are displayed". My data is pretty noisy at high q. What do you recommend?
Re: My data is pretty noisy at high q
If your experimental SAXS data at q > 8/R_{g} is just noise around zero (check by plotting your data on an absolute scale) there is not much point in using GASBOR. Just use DAMMIF.
Re: Help with GNOM/Gasbor
What i meant to say is that if its necessary to increase the whole range of data when calculating the pair distance distribution. Autognom gives the best quality estimate but does truncate the data.