Comqum-X

The following contains a short guideline to the software developed by Ulf Ryde for solving protein crystal structures guided by QM/MM energies. A more exhaustive guideline can be found at Ulf Ryde's personal homepage; http://signe.teokem.lu.se/~ulf/Methods/comqum_x.html
 * 1) First, we need a refined crystal structure (pdb file) including structure factors (CVS file). Both can be downloaded from the pdb data base http://www.rcsb.org
 * 2) Copy needed files cp /lunarc/nobackup/projects/bio/Data/Comqumx_lib/*.
 * 3) Modify the sedfile (see example below) and specify the source of coordinates (usually the pbd file - this can optionally be renamed e.g. to coord_pdb_id).
 * 4) pdb data (see http://signe.teokem.lu.se/~ulf/Methods/comqum_x.html)
 * 5) Alternate configurations. Delete (keep highest) or include; can be done by maestro or by changepdb, command occ. A good check if occupations are deleted is to take the major occupation and see if Rfree changes
 * 6) Insert END at the end of all files (we did not do this?)
 * 7) Do changepdb, command seg followed by w (formatting)
 * 8) get topology files for unusual residues (many can be found at the HIC-UP server)
 * 9) Continue from here..... (NOT DONE YET)
 * 10) Run cqx_sedscript; this will construct the CNS files "bindividual.inp  bindividual1.inp  generate.inp  make_cv.inp  mmrun2  mmrun3". The b-individual files contain information about the B-factors from the crystallographic data.
 * 11) run generate.inp
 *  cns generate.out 


 * grep ERR    generate.out


 * grep WRN    generate.out


 * grep INFO   generate.out


 * grep error  generate.out 


 * grep warning generate.out 
 * Note; if there are alternate configurations in the coordinate pdb file given in the CNS format (look for AAC1 or BAC1 in the coordinate pdb file), these have to be deleted by vi-editor command:
 * :%s/AAC1/B  /g


 * :%s/BAC1/B  /g</tt>
 * Note that the lowest points 10. - 12. in Ulf's guide are not required configurations
 * 1) run pdbtocomqum: Select system 2 as 1000.0 and system 3 as 0.0
 * Note: When running an X-ray with hydrogens, the control file contains $pointcharges. This should be deleted!
 * 1) Cont. list

sedfile example
s/cqx_space_group/p2(1)2(1)2(1)/ # Change! Look for SPACE GROUP </tt>

s/cqx_crystal_alpha/90/ # Default, check: Look for CRYST1 (unit cell dimension; angles) </tt>

s/cqx_crystal_beta/90/ </tt>

s/cqx_crystal_gamma/90/ </tt>

# Look for CRYST1 (unit cell dimension; lenghts) </tt>

s/cqx_crystal_a/32.040/ # Change! Look for CRYST1 (unit cell dimension; lenghts) </tt>

s/cqx_crystal_b/75.560/ </tt>

s/cqx_crystal_c/120.360/ </tt>

</tt>

s/cqx_coord_file/lpmo_5vg0_coord/ # rename pdb file to this name</tt>

# Find these under "Dictonaries" at HIC-UP server ("CNS/topology files and CNS parameter files") </tt>

s/cqx_top_file1/per.top/ </tt>

s/cqx_top_file2// </tt>

s/cqx_top_file3// </tt>

s/cqx_par_file1/per.par/ </tt>

s/cqx_par_file2// </tt>

<tt>s/cqx_par_file3// </tt>

<tt># Reflections file </tt>

<tt>s/cqx_reflex_file/lpmo_5vg0.cv/ # Change this one! </tt>

<tt> </tt>

<tt>s/cqx_test_percent/1.72/ # Change: Look for </tt>FREE R VALUE TEST SET SIZE under FIT TO DATA USED IN REFINEMENT.

<tt>s/cqx_resolution_up/32.0/ # Change: Look for RESOLUTION RANGE LOW under FIT TO DATA USED IN REFINEMENT. </tt>

<tt>s/cqx_resolution/1.10/ </tt>

<tt> </tt>

<tt>s/cqx_map_residues/(resid 1) or (resid 2)/ </tt>