Working With Metal Ions

From NESG Wiki
Revision as of 15:46, 4 January 2010 by Gliu (talk | contribs)
Jump to navigation Jump to search

How to Identify a Protein with metal ion?

  1. Bioinformatic:
  2. Experimental:

'Running CYANA/AUTOSTRUCTURE for Proteins with Metal Ions'

CYANA/AUTOSTRUCTURE can handle structure calculation for protein with metal ions.

  1. Check CYANA.lib for special residues: (ME, PL, LL5, LP, CYSZ,HISZ,...,CYS-, CYSS,...).
  2. Modify the protein's sequence file to have the required special residues included
  3. Metal-ligand geometry parameters. Use reliable bond length and bond angle parameters to derive upper/lower distance limits and angle constraints, and treated as addtional manual constraints during CYANA calculation.

Running CNS for Proteins with Metal Ions

The CNS refinement for proteins with metal ions can also be performed with the new WaterrefCNS script, using the appropriate metal ion nomenclature in both the pdb coordinates and distance constraints.

For manual CNS refinement please do the refinement as following:

1.  Set environment for CNS1.1 by run source /farm/users/gliu/alias.cns

      alias cns1 /farm/software/cns/cns_solve_1.1/intel-i686-linux_g77/bin/cns
      setenv CNS_TOPPAR /farm/data/gliu/cns1/

In addition to the topology and parameter files, the metal ion parameter file ion.top is required. An example can be found in "/farm/users/gliu/projects/cns_cuttha_cis" with all required input files. 2.  Prepared required files as described above (final.tbl, final cns format PDB files and put in xplorPDB dir with name as sa_#.pdb) except the PDB file should include the metal ions with format according to CNS library ion.top. cp sa_1.pdb as  template.pdb, input files for creating mtf file Note that alignment is important. eg:

ATOM   1249  OT2 ALA    83      69.296  13.232   5.744  1.00  0.00
ATOM   1250 ZN+2 ZN2   150      63.086  13.789 -10.407  1.00  0.00      zinc

3.  Run generate_h2o.inp once to create temp_h2o.pdb and temp_h2o.mtf. The extra proton atom in the ligand residues, eg. HIS HD1 or CYS S, are removed by editing the generate_h2o.inp; cis proline is also defined here in the generate_h2o.inp (resid is the residue number prior the proline).

{* any special prosthetic group patches can be applied here *}
{===>}
delete select (name hg and resname cys and resid 61) end
delete select (name hg and resname cys and resid 85) end
delete select (name hd1 and resname his and resid 46) end
delete select (name he2 and resname his and resid 83) end
 patch cisp
           reference=1=( resid 13 )
 end

4.  Edit generate_1.inp to remove the extra proton as did above.

5.  Run generate_20.com, this will run generate_#.inp 20 times, updating each pdb number and this creates cnsPDB/sa_cns_#.pdb

6.  Edit and run re_h2o_cu.inp, the refined pdb is kept in refinedPDB, or

7.  Use subcns to submit cns refinement by using PBS: eg, type "sh subcns". Before run subcns , make a folder " com" contains the following file. Type getpdb to get refined pdb files in refinedPDB after it finished.

  • cns.sc: PBS submission
  • cutc_h2o.mtf: mtf file created as descrive above
  • topology and parameter files: parallhdg5.3C.pro, parallhdg5.3.pro, topallhdg5.3.pro
  • re_h2o_cu.inp: input file for cns refinement