Praktikum Biophysikalische Chemie: Difference between revisions
No edit summary |
|||
Line 59: | Line 59: | ||
# First residue number | # First residue number | ||
# First residue name | # First residue name | ||
=== Write the CYANA script to execute the structure calculation === | |||
# First atom name | # First atom name | ||
# Second residue number | # Second residue number |
Revision as of 18:25, 9 January 2009
Download and install the course material
Store the file PBCPraktikumCyana.tgz in your home directory and execute the setup command:
tar zxf PBCPraktikumCyana.tgz setup_cyana
Make a new directory for the structure calculation, and change into it:
mkdir vaso cd vaso
Write the sequence file
Use a text editor to write a new file called vaso.seq that contains the peptide sequence, one upper-case residue name per line, given in the standard three-letter code for amino acids (except for cysteine residues that are involved in a disulfide bond, which are denoted by "CYSS"), e.g.
CYSS TYR PHE GLN ASN CYSS PRO ARG GLY
Write the initialization script
Use a text editor to write a new initialization script, init.cya, for the program CYANA with the following content:
cyanalib read seq vaso.seq
These two commands will be executed automatically whenever the program CYANA is started. The cyanalib command reads the standard residue library of CYANA, and the command read seq vaso.seq reads the polypeptide sequence.
Write the NOE distance restraint file
Use a text editor to write a new file, vaso.upl, that contains the upper distance bounds derived from NOESY cross peaks, using the standard CYANA nomenclature for atoms in proteins and the same format as in the following example:
91 THR HB 93 GLN QB 5.50 80 SER HB2 81 ILE H 4.22 80 SER HB3 81 ILE H 4.22 81 ILE HA 84 LEU H 4.01 81 ILE HA 84 LEU HB2 4.47 81 ILE HA 81 ILE QG2 3.46 81 ILE HA 81 ILE HG12 3.77 28 VAL HA 39 LEU HG 3.97 52 SER H 52 SER HB2 3.96 52 SER H 52 SER HB3 3.96 99 SER QB 101 VAL H 5.50 43 SER H 43 SER QB 3.12 43 SER QB 48 GLU H 4.07 42 GLU HA 43 SER QB 5.50 43 SER QB 48 GLU HB2 3.95
Each line specifies an upper bound on the distance between two hydrogen atoms. The data in the 7 columns are:
- First residue number
- First residue name
Write the CYANA script to execute the structure calculation
- First atom name
- Second residue number
- Second residue name
- Second atom name
- Upper distance bound in Å
Residue and atom names are given in upper case letters. The exact number of spaces between different items is irrelevant, but the "TAB" key should not be used.
Note that the backbone amide hydrogen is called "H", not "HN".
Degenerate groups of atoms, e.g. methyl groups, and diastereotopic pairs of hydrogen atoms, e.g. HB2/HB3 in serine, are referred to by "pseudoatoms" whose names are derived from the names of the hydrogen atoms that they represent by changing the first letter from "H" to "Q" and omitting the last digit. For instance, "HB2" and "HB3" are represented by a pseudoatom called "QB".