Structure Refinement Using XPLOR-NIH: Difference between revisions

From NESG Wiki
Jump to navigation Jump to search
No edit summary
No edit summary
 
(22 intermediate revisions by the same user not shown)
Line 37: Line 37:




==== I.  The first four files needed are NMR restraints in Xplor/CNS format.  ====
==== The NMR restraints should be in Xplor/CNS format.  ====


If you have restraint files in another format, one way to convert is using PDBStat (reference)
If you have restraint files in another format, one way to convert is using [[PdbStat]]


==== II.  Creating protein structure file, (psf), using mkpsf_prot.inp script  ====
==== The input PDB&nbsp;files should be in Xplor-NIH format<br>  ====
 
One way to convert them from another format, such as Cyana or CNS is using [[PdbStat]]
 
=== Creating protein structure file, (psf) ===
 
==== Editing and running the mkpsf_prot.inp script  ====


The X-Plor psf file specifies the atom parameters and connectivity for all amino acids.&nbsp; The xplor script mkpsf_prot.inp is edited with the protein sequence in three letter code and used to generate the psf files by running:  
The X-Plor psf file specifies the atom parameters and connectivity for all amino acids.&nbsp; The xplor script mkpsf_prot.inp is edited with the protein sequence in three letter code and used to generate the psf files by running:  
Line 48: Line 54:
Download the script [[Media:Mkpsf_prot.inp|mkpsf_prot.inp]].<br> Edit the script with your sequence and to match your input xplor PDB files:  
Download the script [[Media:Mkpsf_prot.inp|mkpsf_prot.inp]].<br> Edit the script with your sequence and to match your input xplor PDB files:  


#Have chain A or no chain ID to match your input PDB file  
#Have chain A or no chain ID to match your input PDB file.
#Add the TOPPAR:toph19.pep line if your input structures have NH3+ and COO- at the N- and C-term with the atoms HT1, HT2, HT3 and OT1, OT2, respectively.  
#Add the TOPPAR:toph19.pep line if your input structures have NH<sup>3+</sup> and COO- at the N- and C-termini with the atoms HT1, HT2, HT3 and OT1, OT2, respectively.
#Fix the Histidine patch to match the protonation state / tautomer  
#Fix the Histidine patch to match the protonation state / tautomer.
#If you protein sequence starts at residues 43, for example, add extra residues to the sequence and uncomment the "dele sele (resid 1:42)" line by removing the&nbsp;!.  
#If you protein sequence starts at residues 43, for example, add extra residues to the sequence and uncomment the "dele sele (resid 1:42)" line by removing the "!".  
#Add cis Proline patch if you have any
#Add cis Proline patch if you have any cis Prolines.<br>


==== III.&nbsp; Edit simulated annealling script  ====
=== Edit and run the simulated annealing script<br> ===


Download the script prot_sa_refine_nogyr.inp (for refinement without the radius of gyration, such as proteins that aren't globular, or you haven't yet determined the residues that should be included in the "core")<br>
==== Edit simulated annealing script ====


or prot_sa_refine.inp  
Download the script [[Media:Prot_sa_refine_nogyr.inp|prot_sa_refine_nogyr.inp]] (for refinement without the radius of gyration, R<sub>gyr</sub>, such as proteins that aren't globular, or you haven't yet determined the residues that should be included in the "core")<br>


or [[Media:Prot_sa_refine.inp|prot_sa_refine.inp]] to include the R<sub>gyr</sub> term.


<br>


Edit the script:
Edit the script:  


#To match your input .tbl file names.&nbsp; Remove the "!" from prot_hbond.tbl if you have hbond restraints.
#To match your input PDB and .tbl file names.&nbsp; Remove the "!" from prot_hbond.tbl if you have hbond restraints.  
#Edit the radius of gyration line for your "core" residues.
#Edit the radius of gyration line for your "core" residues.


For example, if your core residues are 13-92 , then the number of ordered residues, N<sup>residues</sup> is 80.&nbsp; Use the equation below to determine the radiue of gyration.&nbsp; In this case R<sub>gyr</sub> = 11.6.
For example, if your core residues are 13-92 , then the number of ordered residues, N<sub><sup>residues</sup></sub> is 80.&nbsp; Use the equation below to determine the radiue of gyration.&nbsp; In this case R<sub>gyr</sub> = 11.6.  
<pre>    Rgyr = 2.2 * (N<sub>residues</sub>^ 0.38)  
<pre>    Rgyr = 2.2 * (Nresidues^ 0.38)  
</pre>
</pre>  
The line in the prot_sa_refine.inp script should look like this:
The line in the prot_sa_refine.inp script should look like this:  
<pre>    assign ( resid 13:92 ) 50.0 11.6
<pre>    assign ( resid 13:92 ) 50.0 11.6
</pre>
</pre>  
The number 50 is a scaling constant, so don't change it.
The number 50 is a scaling constant, so don't change it.


==== Run the Xplor-NIH refinement<br> ====


 
Edit the number of structures that your are&nbsp; refining and run like this:  
3.&nbsp; Edit the number of structures that your are&nbsp; refining and run like this:
<pre>xplor &lt; prot_sa_refine.inp &gt; sa.out &amp;
<pre>xplor &lt; prot_sa_refine.inp &gt; sa.out &amp;
</pre>
</pre>
Check the sa.out for the line "ERR" to find errors.&nbsp; Typically there will be errors if the .psf file atom names do not match the input pdb file names.&nbsp; This can help you track them down.<br>


== References  ==
== References  ==


<references/>
<references/>

Latest revision as of 19:32, 30 November 2009

Introduction

Here we describe the protocol for Xplor-NIH refinement of NMR structures.  The Xplor-NIH molecular structure determination package is available at http://nmr.cit.nih.gov/xplor-nih/[1][2]


The restrained simulated annealing protocol uses many of the updated features of the Xplor-NIH software, including:

  • the IVM module for torsion angle and rigid body dynamics
  • a radius of gyration term to represent the weak packing potential
  • database potentials of mean force to refine against:
  1. Cα/Cβ chemical shifts
  2. multidimensional torsion angles,
  3. backbone hydrogen bonding database 
  • RDCs can be included but aren't described here

The topology and parameter files protein.top and protein.par, which were designed to agree with bond lengths and angles from the CSDX force field, are used.


Protocol

Setup

The files needed to run the Xplor-NIH refinement are:

  • prot_noe.tbl
  • prot_dihe.tbl
  • prot_hbond.tbl (optional)
  • prot_shiftsCACB.tbl
  • prot.psf (described below), can be obtained from mkpsf_prot.inp
  • sa_refine.inp (edited to reflect your protein parameters)
  • xplor format PDB files to be refined (described below)


The NMR restraints should be in Xplor/CNS format.

If you have restraint files in another format, one way to convert is using PdbStat

The input PDB files should be in Xplor-NIH format

One way to convert them from another format, such as Cyana or CNS is using PdbStat

Creating protein structure file, (psf)

Editing and running the mkpsf_prot.inp script

The X-Plor psf file specifies the atom parameters and connectivity for all amino acids.  The xplor script mkpsf_prot.inp is edited with the protein sequence in three letter code and used to generate the psf files by running:

xplor < mkpsf_prot.inp > temp.out

Download the script mkpsf_prot.inp.
Edit the script with your sequence and to match your input xplor PDB files:

  1. Have chain A or no chain ID to match your input PDB file.
  2. Add the TOPPAR:toph19.pep line if your input structures have NH3+ and COO- at the N- and C-termini with the atoms HT1, HT2, HT3 and OT1, OT2, respectively.
  3. Fix the Histidine patch to match the protonation state / tautomer.
  4. If you protein sequence starts at residues 43, for example, add extra residues to the sequence and uncomment the "dele sele (resid 1:42)" line by removing the "!".
  5. Add cis Proline patch if you have any cis Prolines.

Edit and run the simulated annealing script

Edit simulated annealing script

Download the script prot_sa_refine_nogyr.inp (for refinement without the radius of gyration, Rgyr, such as proteins that aren't globular, or you haven't yet determined the residues that should be included in the "core")

or prot_sa_refine.inp to include the Rgyr term.


Edit the script:

  1. To match your input PDB and .tbl file names.  Remove the "!" from prot_hbond.tbl if you have hbond restraints.
  2. Edit the radius of gyration line for your "core" residues.

For example, if your core residues are 13-92 , then the number of ordered residues, Nresidues is 80.  Use the equation below to determine the radiue of gyration.  In this case Rgyr = 11.6.

    Rgyr = 2.2 * (Nresidues^ 0.38) 

The line in the prot_sa_refine.inp script should look like this:

    assign ( resid 13:92 ) 50.0 11.6

The number 50 is a scaling constant, so don't change it.

Run the Xplor-NIH refinement

Edit the number of structures that your are  refining and run like this:

xplor < prot_sa_refine.inp > sa.out &

Check the sa.out for the line "ERR" to find errors.  Typically there will be errors if the .psf file atom names do not match the input pdb file names.  This can help you track them down.

References

  1. Schweiters, C.D.; Kuszewski, J.J.; Tjandra, N.; Clore, G.M. “The Xplor-NIH NMR Molecular Structure Determination Package,” J. Magn. Reson. 2003, 160,66-74
  2. Schweiters, C.D.; Kuszewski, J.J.; Clore, G.M. “Using Xplor-NIH for NMR molecular structure determination,” Progr. NMR Spectroscopy. 2006, 48, 47-62