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.