· Delarue web servers|
· Nomad Flow-chart
· Normal Mode calculation
· Examples (Movies)
· Submit a job (from PDB file)
· Split trajectory (for MR)
· Generate decoys (MixMod)
· Elastic Energy (Perturb. Anal.)
· Overlap coefficients
· Submit a job
· Include Profit
· X-Ray refinement
· Standard refinement
· Screening Mol. Repl. Solns.
· EM refinement
· Get Structure Factors
· Submit a job (no NCS)
· Submit job with NCS
· Docking refinement
· Submit a job
· Force field methods
· Energy minimization
· Gromacs NMA
· Job queue status
Flexible refinement in EM maps
Performs refinement of a model against phased phased structure factors
along Normal Mode directions of low frequency.
In addition to the PDB and the formatted structure factors data files, you must choose:
The engine of refinement is a Conjugate-Gradient algorithm that will optimize the amplitudes of the modes to fit experimental data.
- The crystallographic symmetry operators and cell parameters.
- The resolution limits of the data.
- The number of modes to use (increasing frequency order).
- The cutoff for interactions in the Elastic Network Model.
This can also make use of (cryo) electron microscopy data, expressed in reciprocal space.
For this, you need first to (Fourier) transform the electron density map of electron microscopy
data into structure factors (reciprocal space) in your space group, e.g. P1.
See the Documents section on the left for more details on how to invert
your map into structure factors (Get Struct. Fact.).
Input data formats
- The job title is just for your own identification, but note that
it will show up in the public job queue (but your results will not be public).
- The coordinate files should be in PDB format, with only a single
structure (no multiple models). Atoms marked with alternate residue
flags will be removed. By default the refinement only uses C-alpha atoms,
and the refined output will only consist of C-alpha coordinates.
- The symmetry data is a text file that describes the box and symmetry operators.
On the first line there should be six numbers, representing unit cell dimensions (a,b,c)
in Ångströms and the cell angles (alpha,beta,gamma) in degrees.
The format is (6f8.3).
This is followed by symmetry operators, using the AMoRe program format. Each operator
is followed by a '*', and the list is terminated with an 'end' record. A couple of examples
(you should NOT just copy these):
- The reflection data is a free-format text file with one
reflection per line, including the phase.
Columns are separated by at least one space, and
there are 5 columns per line. The first
three are integers corresponding to the h,k,l indices. Column 4 is
the amplitude of the reflection and column 5 the phase in degrees.
Reflections outside the low/high limits specified by you will be discarded.
An example: Fobs_env.data.
- The refinement is carried out in reciprocal space using normal mode amplitudes.
Low-frequency modes offer the advantage to allow for collective and large-
In the order of 10 modes should work well for normal proteins, but you
can experiment with higher values, although the execution will be slower.
If you choose to include modes 1-6, the program will also attempt rigid-body
In any case, you should watch closely the output file which displays some information about possible CA-CA distances violations. If there are too many of these violations, the model should be discarded and another job (with a smaller number of modes) submitted, until this kind of behaviour disappears.
- The number of modes is usually 7-16 or 7-21. If you want to refine the position and orientation
of the model in the map as well, you can use the first 6 modes (i.e. modes 1 to 16 for instance). Their
frequency is damped to soften the effect of this refinement using (lambda(k),k=1,6) = lambda(7)/5., as suggested in
Hinsen et al., Biophys. J. (2005) 88:818-827.
- The STRONG option for the force constant puts a weight of 1000 to the energy of CA(i)-CA(i+1) distances,
compared to 1 for all other (non-bonded) distances (NORMAL option).
Only CAs from the same subunit are taken into account. For P(i)-P(i+1) distances (if nucleic acids are present), a weight of 30 is taken for distances between bonded nucleotides.