The following sections shortly describe the method implemented in SREFLEX, how to run SREFLEX from the command-line, the required input and the produced output files.
If you use SREFLEX in your work, please cite:
Panjkovich A. and Svergun D.I. (2016) Deciphering conformational transitions of proteins by small angle X-ray scattering and normal mode analysis. Phys. Chem. Chem. Phys., 18, 5707-5719
The SREFLEX program uses normal mode analysis (NMA) in Cartesian space to estimate the flexibility of high-resolution models of biological macromolecules and improve their agreement with experimental small angle X-ray and neutron scattering (SAXS and SANS) data.
The method starts from a given structural conformation and a corresponding SAS profile in relative disagreement.
The structure is partitioned into pseudo-domains based on user input or automatically from the protein dynamics as predicted by NMA.
The algorithm proceeds hierarchically to first probe large rearrangements and progresses into smaller and more localized movements.
The output consists in a set of structural models representing possible conformational changes that improve the agreement with the experimental SAS profile. A mode to generate a pool of conformers from an initial structure file ("pool" mode) was also added.
Additionally, the following arguments can be forwared to CRYSOL/CRYSON to modify how theoretical intensities are calculated: --lm, --fb, --sm, --ns, --un, --dns, --dro, --cst and --eh. Please see the CRYSOL manual for more details.
SREFLEX reads SAS experimental data files (*.dat) in ascii format containing 3 columns: (1) experimental scattering vector, (2) experimental intensity and (3) experimental errors.
SREFLEX expects atomic coordinates in Protein Data Bank (PDB) format.
NMA calculations are based on backbone atoms (centroids): alpha carbons (CA) for proteins and sugar C1' for nucleotides.
Please note that SREFLEX parses centroid atoms from the ATOM record.
The program will not work properly if input coordinates lack backbone atoms.
Non-solvent HETATM entries are grouped by residue identifier and number.
These are then associated to the closest ATOM centroid in the structure for the application of rotations and translations.
Theoretical scattering calculations consider all non-H ATOM and non-solvent HETATM entries.
The restrained refinement stage that is performed initially by SREFLEX considers the input structure as a set of rigid bodies or pseudo-domains.
By default, SREFLEX will partition the structure based on its NMA and define (pseudo-) domains that can move with relative independence of each other.
The algorithm is explained in detail in the SREFLEX article.
This automatic partitioning procedure is the default behaviour if a single PDB file is provided as the <COORD FILE> argument.
Results may be improved if the user defines custom domains beforehand, these definitions can be based for example on structural or evolutionary information.
In such case, the coordinates for user-defined domains should be provided to SREFLEX as different PDB files in a comma separated list (no spaces) in place of the <COORD FILE> argument.
This will deactivate the automatic partitioning procedure and each PDB file provided by the user will be considered as a pseudo-domain or rigid-body during the initial restrained refinement stage.
Upon execution, SREFLEX will create a directory to store results.
The name of the output directory can be set through the prefix argument.
If the directory exists already, a number will be appended to create a new directory.
The examples are based on the known conformational change that adenylate kinase undergoes upon catalytic activity.
PDB entry 1ake (chain A) was used to simulate a SAXS profile for the closed conformation (available as example file closed1akeA.dat).
PDB entry 4ake (chain A) corresponds to the open conformational state (available as example file open4akeA.pdb).
SREFLEX will refine the 4ake open conformation until it matches the SAXS profile of the closed conformation.
The complete run should take a few minutes on a normal desktop computer.
To illustrate the definition of custom structural domains, open4akeA.pdb was separated into two files (open1.pdb and open2.pdb) corresponding to the SCOP domain classification of the 4ake PDB entry.
In this case, SREFLEX will refine the structure considering both PDB files as rigid-bodies for the initial restrained refinement stage, disabling the automatic structure partitioning approach.
The --prefix option is used to store the results in a new directory named customDomains.
Note that there are no spaces separating the PDB filenames at the <COORD FILE> argument: