Script mode - spins from a sequence file

Image sequence_600x600

Alternatively to setting up the molecule, residue, and spin containers via 3D structural data, a plain text columnar formatted file can be used. This is useful for when no 3D structure exists for the molecule. It also has the advantage that the residue and atom names need not conform to the PDB standard. An example for reading sequence data is:

# Create a data pipe.
pipe.create(pipe_name='R1 1200', pipe_type='relax_fit')

# Set up the 15N spins.
sequence.read(file='noe.500.out', mol_name_col=1, res_num_col=2, res_name_col=3, spin_num_col=4, spin_name_col=5)
spin.element(element='N', spin_id='@N*')
spin.isotope('15N', spin_id='@N')

Here the molecule, residue, and spin information is extracted from the “noe.500.out” file which could look like:

# mol_name          res_num  res_name  spin_num  spin_name  value             error               
Ap4Aase_new_3_mol1  1        GLY       1         N          None              None                
Ap4Aase_new_3_mol1  2        PRO       11        N          None              None                
Ap4Aase_new_3_mol1  3        LEU       28        N          None              None                
Ap4Aase_new_3_mol1  4        GLY       51        N          0.03892194698453  0.01903177024613    
Ap4Aase_new_3_mol1  5        SER       59        N          0.31240422567912  0.01859693729836    
Ap4Aase_new_3_mol1  6        MET       71        N          0.42850831873249  0.0252585632304     
Ap4Aase_new_3_mol1  7        ASP       91        N          0.53054928103134  0.02799062314416    
Ap4Aase_new_3_mol1  8        SER       104       N          0.56528429775819  0.02170612146773    
Ap4Aase_new_3_mol1  9        PRO       116       N          None              None                
Ap4Aase_new_3_mol1  40       TRP       685       N          0.65394813490548  0.03830061886537    
Ap4Aase_new_3_mol1  40       TRP       698       NE1        0.67073879732046  0.01426066343831

The file can contain columns for the molecule name, the residue name and number, and the spin name and number in any order though not all are needed. For example for a single protein system, the molecule name, residue name and spin number are nonessential. Or for an organic molecule, the molecule name, residue name and number and spin number could be nonessential. The subsequent user functions in the above example are used to set up the spin containers appropriately for a model-free analysis. These are not required in the automatic analysis of GUI as these user functions will be presented to you when adding relaxation data, or when clicking on the heteronucleus and proton buttons (“X isotope” and “H isotope”).

In the GUI, the creation of molecule, residue, and spin containers from a sequence file is also available via the “Load spins” wizard within the spin viewer window (vide supra).

The relax user manual (PDF), created 2020-08-26.