Flex-EM

A method for protein structure fitting and refinement guided by cryoEM density.

The Monte Carlo (MC) rigid fitting protocol (Mod-EM) is documented in the MODELLER manual: the density class. The scripts below run the Conjugate Gradients (CG) and Molecular Dynamics (MD) refinement protocols of Flex-EM. Note that you will need MODELLER 9v2 or later to use these scripts.

2 files have to be edited by the user:

1. flex-em.py: the user has to set the input parameters (the EM map parameters and the name of the input structure) and define the type of optimization (CG or MD).
2. rigid.txt: the user has to define the rigid bodies. The format of this file is the following (see example-rigid.txt for an example):
   • A comment line describing what kind of a rigid body : "domain", "helix", "beta", "loop", followed by lines describing the rigid bodies : each line describes one rigid body by specifying the initial and final residue of each of the segments in that rigid body.
   • For example: 36 50 55 100 means that residues 36-50 and 55-100 will be included in the rigid body. (Note that before assigning the rigid bodies the structure has to be renumbered from 1 because MODELLER is going to renumber it during the optimization).
   • If you add "nr" after the rigid body the position of this rigid body will not be randomized prior to running of the CG minimization.

Scripts available for download:

• CG.py
• MD.py
• rigid.py
• flex-em.py
• example-rigid.txt

Reference: M. Topf, K. Lasker, B. Webb, H. Wolfson, W. Chiu, A. Sali. Structure 16, 295-307, 2008

Benchmark:

download benchmark*

Notes:

• To minimize bias, the benchmark density maps were not produced with Mod-EM, but with pdb2vol in SITUS (Wriggers, et al. 1999), which uses a different Gaussian smoothing technique.
• Cα RMSDs were calculated for the benchmark using the superpose method in MODELLER 9v0.

For further information, please contact Maya Topf at