· 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 step
· X-Ray refinement
  · Standard refinement
  · Screening Mol. Repl. Solns.
· EM refinement
  · Documents/Examples
  · Get Structure Factors
  · Submit a job (no NCS)
  · Submit job with NCS
· Docking refinement
  · Submit a job (refinement)
  · Submit a job (scanning)
· Force field methods
  · Energy minimization
  · Gromacs NMA
· References
· Homes/Links
· Job queue status

Normal mode calculation and visualisation

Input
You need at least one pdb file to calculate Normal Modes using the Elastic Network Model, plus a couple of parameters described at the bottom of this form.
Options
Choose the number of modes you want, provided that there are less than 106 (=100+6); that is only to save some of our disk space: if you want more, please send us an e-mail.
The scale (amplitude) of the movement can be chosen by controlling the average rmsd of the trajectory (see below). The rmsd per residue is given on output for each chosen mode (as a plot).

Output
There will be one output PDB file per mode, named mode_#.pdb, directly loadable by PyMol or VMD: this is a concatenation of 30 PDB files separated by an ENDMDL card. This PDB file maybe split into 30 pieces for Molecular Replacement. There is also a tarred file with all modes.
Bonus
As a bonus, you will get a file loadable by Pymol to generate a picture of the Elastic Net associated to your protein 1GGG.pymol (thanks to Patrice Koehl).
Finally, there will be the modes.dat file that contains, for each mode, its frequency and eigenvector list, for your personal use.

Your email adress: (Recommended, for notification)



Job title: (Only alphanumerical characters, no space)



PDB file to calculate modes for
(Example: 1GGG.pdb). N.B. No Hydrogen atom(s)!



Processing the PDB (for proteins only)?
       
The Zacharias option keeps the C-alpha and models sidechains by up to 2 pseudoatoms as in ATTRACT program (Zacharias, Proteins, 2005, 60:252)


Number of modes to calculate (max 106 to save our disk):
(The first six ones are translation and rotation)


Distance weight parameter for elastic constant (Å):



ENM Cutoff to use for mode calculation (Ångström):



Average Rmsd (Angstrom) in output trajectories:



Method to use (advanced):
     



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 file should be in PDB format, with only a single structure (no multiple models). Atoms marked with alternate residue flags will be removed. Whatever atoms (ATOM card) are in the file will be used for the calculation (but not the HETATM ones). The length of each mode vector will be 3*natoms.
    It is not recommanded to include Hydrogen atoms in the PDB file.

  • All interactions are weighted by exp(-(d_ij(0)/d_0)^2), where d_0 is a distance-weight parameter. This effectively introduces a smoother "cutoff" value than in the original Tirion model. A value of 3.0 Angstrom works well for CA-only models, but you might prefer a larger value (5.0-10.0 Angstrom) for all-atom structures. See the Hinsen article in the References section.

  • A cutoff is used in the mode calculation. In the Tirion model (Elastic Network Model) only those pairs of atoms that are closer than the cutoff are linked by a spring of universal length. Ideally this should be choosen so that the weighting causes the interactions to be negligible outside the cutoff, but in practice a cutoff around 10 Angstrom works fine in almost all cases. In general you might be able to use a smaller cutoff value for all-atom calculations than CA-only calculation. See the Tirion article in the References section.

  • You can choose how many modes to calculate. The execution time for the sparse matrix solver is more or less proportional to the number of modes requested, while the full matrix solver execution is dominated by the reduction to tridiagonal form - the actual mode calculation is very fast in that case. Due to the memory requirements you cannot use the full solver for more than about 5000 atoms, and the sparse solver is much faster even for small systems (from about 500 atoms) if you only want the first (lowest-frequency) 10-100 modes. Remember that the first six modes correspond to rigid body motions (translation and rotation), so the structurally interesting ones begin at 7.
    You should check in the log file that the first 6 modes have an eigenvalue close to zero.

  • If you really want to you can override the automatic choice and specify either the sparse or full solver. Be aware that the full solver wont work for more than 5000 atoms, but it is faster for small systems or if you want a large number of mode vectors (or if you wanted a huge cutoff for some reason).

  • Some timing examples for 1ANF with 2860 atoms (Hessian dimension is 8580 by 8580): Calculating 20 modes takes 33 seconds and for 100 modes it increases to 1 minute and 40 seconds.



  Marc Delarue http://lorentz.dynstr.pasteur.fr
Page last modified 15:43 December 23, 2016.