0
EMBL Hamburg Biological
Small Angle Scattering
BioSAXS
  Home  >  Software  >  OLIGOMER  >  Manual

OLIGOMER manual

OLIGOMER

Written by D. I. Svergun and P. V. Konarev.
Post all your questions about OLIGOMER to the ATSAS Forum.

© ATSAS Team, 1994-2004

Manual

Introduction

The program OLIGOMER fits an experimental scattering curve from a multicompoment mixture of proteins to find the volume fractions of each component in the mixture.

The input data for oligomer is experimental SAXS or SANS scattering curve from a mixture of several dictinct components and a set of form-factors for each component of this mixture.

The experimental scattering intensity I(s) from a mixture of K different particles (components) is written as:

I(s) = ∑(wi×Ii(s))

where wi and Ii(s) are the volume fraction and the scattering intensity from the i-th component, respectively.

Given the intensities from the components (form-factors), OLIGOMER finds the volume fractions by solving a system of linear equations using the algorithm of nonnegative or unconstrained least-squares to minimize the discrepancy between the experimental and calculated scattering curves.

X-ray or neutron form-factors of proteins with known structure or of dummy atoms models can be computed using the programs CRYSOL and CRYSON respectively, which both are part of the ATSAS package.

Running OLIGOMER

Prior to running the application a file containing the Scattering intensity from the form-factors must be prepared in the format:

1st line: comment

from 2nd line: scattering vector (s), iForm1, iForm2 ...

here, s = 4 * π * sin (θ)/λ, [1/Angstroms]

Iform1, Iform2 etc are the scattering intensities of the components. They should be normalized in such a way that iForm(s=0) is equal (or proportional to) the square of the molecular weight MW (or total scattering length).

OLIGOMER estimates molecular weights for calculation of volume fractions and radii of gyration from experimental or model scattering curves of components of the mixture. The latter values are used only for determination of the average value of radius of giration for all mixtures.

Therefore, one must enter the real values of molecular weight for each component (for example, based on information about primary sequence) in following cases:

  • scattering curves from components (i.e. sets of form-factors) are standardized to Iexp(0) = 1
  • scattering curves from components calculated by CRYDAM or CRYSOL from dummy atom models built by programs DAMMIN or GASBOR with beads not equal size
  • sets of form-factors are standardized to equally values of Iexp(0)

If the iForm curves and not properly scaled, the fit not be affected but the volume fractions will be incorrect.

If one just uses predicted CRYSOL curves from different oligomers, no additional normalization of the form factors is required.

Please note: To quickly create form-factor file from pdb files and/or from scattering data files (either from ASCII *.dat files where column 1 will be taken as S-axis, column 2 as intensities or from GNOM output files where desmeared curve - "Ireg" column will be taken for intensity) one can use program "FFMAKER" (runs on Windows) that automatically creates form-factor file (one can use interactive or batch mode to run the program). To see usage information type: "FFMAKER /help". There is no internal scaling in FFMAKER, the curves will correspond to CRYSOL output files (if calculated from PDB models) or to the input scattering curves (as they are).

Command-Line Arguments and Options

The following program options can be given as keys with parameters:

Short OptionDescription
/ff Name of the form factor input file.
/dat Name of the experimental data input file.
/un Angular units (1: A-1=default; 2: nm-1; 3: 2pi*s A-1; 4: 2pi*s nm-1)
/smin minimum value of angular axis s
/smax maximum value of angular axis s
/out the name of the output file, default is oligomer.log
/cst adds constant as additional component. Default value is 10 times greater than maximum I(0) value from the form-factor file. This choice of constant will keep the volume fraction of the "constant" component equal to 0 providing correction of the scattering curve from parasitic constant background. In batch mode one needs to provide only the option "/cst" (without value for the constant). The constant value can be changed only in interactive mode.
/pl plots the fitted curves using sasplot
/br brief output for .log file

Interactive Configuration

OLIGOMER reads experimental data and form-factors (both should be in ASCII format). After starting OLIGOMER you may specify:

Prompt Possible value(s) Default value Description
Program option 0 or 1 0 option 0 - permits to run OLIGOMER for multiple data sets using one form-factor file, option 1 - permits to use multiple form-factor files for fitting one data file.
Input data, form-factor file name filename ff.dat Input file should contain intensity values from at least one component. It can be prepared by FFMAKER program from *.pdb, *.dat and *.out (GNOM) files or manually by composing the following columns: 1st column - s-axis, 2nd column - intensity from component 1, 3rd column - intensity from component 2 etc.
Use constant as additional component Y or N N if 'Yes' - additional constant component will be added, that permits to improve the fit at higher angles, if required
Calculated MW and Rg X1,X2 var default values are estimated from form-factor intensity curves
Experimental data file name filename lys.dat The name of the file containing the experimental SAXS profile from multicomponent system
Angular units in the input file :
4*pi*sin(theta)/lambda [1/angstrom] (1)
4*pi*sin(theta)/lambda [1/nm ] (2)
2* sin(theta)/lambda [1/angstrom] (3)
2* sin(theta)/lambda [1/nm ] (4)
1-4 1 Formula for the scattering vector in the data file and its units.
Range for evaluation of scattering X1,X2 var default values are estimated from the full-range of experimental data file
Use of non-negativity condition Y or N Y if 'Yes' - the volume fractions values will be required to be positive
Output fit file filename *.fit by default: experimental file name + .fit extenstion.
Plot the result Y or N Y if 'Yes' - the fit curve and exp. curve will be plotted
Process another data set Y or N Y if 'Yes' - the fitting procedure can be repeated for another data set.

Runtime Output

On runtime, two lines of output will be generated:

 Chisquare  MW   Rg   Volume fractions+-errors
 1.5     22.5    1.7     0.600+-0.001 0.400+-0.001

The fields can be interpreted as follows, top-left to bottom-right:

Field Description
Chisquare The quality of the fit (discrepancy value)
MW The averaged molecular mass of the system
Rg The averaged radius of gyration of the system
Volume fractions+-errors The volume fractions of the components in the mixture and errors of volume fractions. The sequence of the components is the same as specified in the form-factor file

OLIGOMER Input Files

The experimental data file is composed as follows:

SAXS experimental data
experimental description
0.177473E+00  0.737508E+02  0.169651E+01
0.180129E+00  0.728071E+02  0.155912E+01
.....            ....          ....

The first two lines are reserved for sample description and further details. The first column in the third row is the s-value, the second column are the scattering intensities and the third column are the experimental errors.

The form factor files consists of the individual form factor of the components. If they are derived by the program CRYSOL the second column of the calculated form factors has to be used.

0.00E+00  4.30E+08  3.31E+09  ...
1.00E-02  4.23E+08  3.19E+09  ...
2.00E-02  4.03E+08  2.87E+09  ...
 ....        ...       ...    ...

The first column is the s-value, the second is the first formfactor and the third column the second form factor.

The two input files can have different s-values and data length.

OLIGOMER Output Files

OLIGOMER creates a output file named oligomer.log as default or as specified using the /out option. If the file already exists, the results are appended to the file.

         ***   O L I G O M E R ***
================================================================================
                                              01-Sep-2009   17:26:17
 Option 0: a set of form-factors and several sets of experimental data
 Form-factor file              test.dat
 Real oligomer weights          20849   58193
 Oligomer radii of gyration    22.57  34.37
 Experimental data file        Merge00.dat    Range of Scattering angle:    0.02    0.27
 Using non-negativity condition
 Output file                   Merge00.fit
 ChiSquare     <MW>      <Rg>        Volume fractions +- errors
--------------------------------------------------------------------------------
Merge00.da    2.17      52869     33.82      0.143+-0.004  0.857+-0.002
================================================================================

If the brief /br option is used only the last line with the results is plotted.

If the batch mode is used, olig_summary.txt file is created, where the information about the time (when the program was launched) and the parameters of batch mode are stored.

OLIGOMER creates also *.fit file that contains 5 columns: S-axis (1st column), experimental intensity Iexp(2nd column), fit from OLIGOMER Ifit (3rd column), errors of intensity Ierr(4th column), the difference curve (Iexp-Ifit) (5th column).

Examples

Example 1

For fitting the experimental data set m01.dat using the formfactor file ff.dat with the angular axis in nm-1

$> oligomer  /ff  ff.dat  /dat  m01.dat  /un  2

Example 2

For processing all experimental data files with the names beginning with "m" and with "dat" extension using the formfactor file ff2.dat. The fitting will be applied from smin=0.5 for the experimental curves. The fits will be plotted and output oligomer.log file will have the brief format.

$> oligomer  /ff  ff2.dat  m*.dat  /smin 0.5  /pl  

  Last modified: July 29, 2011

© BioSAXS group 2011