Jma at 21:57, 6 January 2010

2010-01-06T21:57:18Z

← Older revision		Revision as of 21:57, 6 January 2010
Line 135:		Line 135:
	<br>		<br>

	~~-- PaoloRossi - 14 Dec 2009~~		<br>

Prossi at 21:22, 15 December 2009

2009-12-15T21:22:11Z

← Older revision		Revision as of 21:22, 15 December 2009
Line 40:		Line 40:
	rm -f t.aco		rm -f t.aco
	</pre>		</pre>
	==== '''~~Sequence~~ file''' ====		==== '''The sequence file''' ====

	The example sequence file, shown below in abbreviated form, now includes the doubled sequence separated by a 48-element (LL2) linker: <br>		The example sequence file, shown below in abbreviated form, now includes the doubled sequence separated by a 48-element (LL2) linker: <br>
Line 77:		Line 77:

	</pre>		</pre>
	==== CALC.cya file ====		==== The CALC.cya file ====

	The working directory that contains all the files necessary to start the calcualtion is ready and the CALC.cya modified for the presence of RDC constraints is used:<br>		The working directory that contains all the files necessary to start the calcualtion is ready and the CALC.cya modified for the presence of RDC constraints is used:<br>
Line 108:		Line 108:

	</pre>		</pre>
	==== Init.cya file <br> ====		==== The Init.cya file <br> ====

	The new intstructions for symmetric dimer treatment are included in the init.cya file, in the example below molecule indentity and molecule symmetry are parametrized and adjustable (probably the parameters have been optimized so results tweaking the weight_ide and weight_sym should be carefully).		The new intstructions for symmetric dimer treatment are included in the init.cya file, in the example below molecule indentity and molecule symmetry are parametrized and adjustable (probably the parameters have been optimized so results tweaking the weight_ide and weight_sym should be carefully).

Prossi: Created page with '== '''Introduction''' == '''IMPORTANT DISCLAIMER:''' A number of NESG NMR groups are currently beta-testing sites for CYANA-3.0 and the information in this page is in…'

2009-12-15T21:06:07Z

Created page with '== '''Introduction''' == '''IMPORTANT DISCLAIMER:''' A number of NESG NMR groups are currently beta-testing sites for CYANA-3.0 and the information in this page is in…'

New page

== '''Introduction''' ==

'''IMPORTANT DISCLAIMER:''' A number of NESG NMR groups are currently beta-testing sites for CYANA-3.0 and the information in this page is intended for use by the licensed members of the NESG consortium, other beta testers, and is to be used in accordance to the program licensing agreement.

 

The following page describes the setup and analysis of an automated structure determination run for symmertic homodimers starting from NOE peaklists in the framework of CYANA-3.0. 

 Results from CYANA should always be accompanied by energy refinement in CNS or NIH-XPLOR.  Further information about the program and publication references can be found in the [http://www.cyana.org/wiki/index.php/Main_Page CYANA WIKI] page.

 

== Automated Homodimer Structure Calculation Setup ==

The structure calculation with automated noesy assignments and intermolecular restraints follows the canonical CYANA recipe.  Simple annealing calculations starting from a set of constraints that include RDC are easily derived by simplifying the scripts below and following the demo scripts. The program requires a sequence file (name.seq), a proton assignment list (name.prot), a noesy peaklist set (name.peaks), a CALC.cya script and an init.cya script.   

An helpful script called '''DIMERIZE''', to be run from within CYANA, can be used to symmetrize various files such as the proton (.prot), angle and NOE constraint files: 
<pre>!/bin/sh

# Chemical shift list
cp RpR324.prot RpR324d.prot
awk 'NF>=5 { printf("%6d%8.3f%8.3f %-5s%4d\n",$1+1000,$2,$3,$4,$5+200) }' $* RpR324.prot >> RpR324d.prot

# Additional distance constraints
echo "
read upl D5.upl
read lol D5.lol
molecules symmetrize
write upl D5dimer.upl
write lol D5dimer.lol
quit" | cyana32

# angle constraints
echo "
read dihe.aco
shift +200
writeaco t.aco
quit" | ancoma
cat dihe.aco t.aco > dihed.aco
rm -f t.aco
</pre>
==== '''Sequence file''' ====

The example sequence file, shown below in abbreviated form, now includes the doubled sequence separated by a 48-element (LL2) linker: 
<pre>MET 1
THR 2
SER 3
THR 4
PHE 5
ASP 6
.
.
.
HIS 100
HIS 101

PL 102
LL2 103
LL2 104
.
.
LL2 148
LL2 149
LP 150

MET 201
THR 202
SER 203
THR 204
PHE 205
ASP 206
.
.
.
HIS 300
HIS 301

</pre>
==== CALC.cya file ====

The working directory that contains all the files necessary to start the calcualtion is ready and the CALC.cya modified for the presence of RDC constraints is used: 
<pre>peaks  := ali5.peaks,aro5.peaks,n3.peaks,xfilt.peaks # names of NOESY peak lists
prot  := RpR324d.prot # names of chemical shift lists
restraints  := ssa.cya,hbonds.upl,hbonds.lol,manual.upl,dihed.aco # additional (non-NOE) constraints
tolerance  := 0.04,0.025,0.3 # chemical shift tolerances: HX2-HX1-X1
calibration := # NOE calibration parameters
structures  := 100,20 # number of initial, final structures
steps  := 10000 # number of torsion angle dynamics steps
rmsdrange  := 10..80 # residue range for RMSD calculation
randomseed  := 56231 # random number generator seed

ssa
noeassign peaks=$peaks prot=$prot autoaco
</pre>
==== Defining inter-chain only peaks ====

Notice the peak contain a list of intermolecular NOEs called xfilt.peaks.  It is possible to instruct the program to make only interchain NOE assignments for the peaklist by using the 9 U xeasy color-code as shown below, the remaining peaklists can be treated as usual as a mixture of inter and intra chain NOEs. 
<pre># Number of dimensions 3
#FORMAT xeasy3D
#INAME 1 h
#INAME 2 C
#INAME 3 H
#CYANAFORMAT hCH
1 8.072 35.125 4.110 9 U 1.000E+06 0.000E+00 e 0 0 0 0
2 2.584 33.247 4.446 9 U 3.610E+05 0.000E+00 e 0 0 0 0
3 6.801 42.287 3.199 9 U 4.280E+05 0.000E+00 e 0 0 0 0
4 6.799 42.276 2.940 9 U 3.440E+05 0.000E+00 e 0 0 0 0

</pre>
==== Init.cya file ====

The new intstructions for symmetric dimer treatment are included in the init.cya file, in the example below molecule indentity and molecule symmetry are parametrized and adjustable (probably the parameters have been optimized so results tweaking the weight_ide and weight_sym should be carefully).
<pre>name:=RpR324
rmsdrange:=10-80
cyanalib
nproc:=8
read seq $name.seq
molecules define 1..101 201..301
molecule identity
weight_ide=0.03
molecule symdist "CA 1..101" "CA 201..301"
weight_sym=0.0025

</pre>
The weight of the symmetry and indentity constraints does not appear to be parametrized in the target function as a penalty.  Further explanation will hopefully become available with the official release.

The use of MPI routine is highly recommended due to the length of the calculation.

 

== Output analysis ==

The output analysis is carried out in the usual manner as for CYANA -2.1. 

 

-- PaoloRossi - 14 Dec 2009

Homodimer Structure Calculation with Symmetry Constraints - Revision history

Jma at 21:57, 6 January 2010

Prossi at 21:22, 15 December 2009

Prossi: Created page with '== '''Introduction''' == '''IMPORTANT DISCLAIMER:''' A number of NESG NMR groups are currently beta-testing sites for CYANA-3.0 and the information in this page is in…'