bugrelax - Bugs: bug #21598, Considering molecule numbers when...

 
 
Show feedback again

bug #21598: Considering molecule numbers when writing pyMol (or Molmol) macros

Submitted by:  Martin Ballaschk <mab>
Submitted on:  Wed 05 Feb 2014 10:16:04 AM UTC  
 
Category: relax's source codeSpecific analysis category: Model-free analysis
Priority: 5 - NormalSeverity: 1 - Wish
Status: PostponedAssigned to: None
Originator Name: Open/Closed: Open
Release: 3.1.5Operating System: All systems

Add a New Comment (Rich MarkupRich Markup):
   

You are not logged in

Please log in, so followups can be emailed to you.

 

Tue 11 Nov 2014 04:36:09 PM UTC, comment #2:

This is being postponed as no reasonable solution can currently be found.

Edward d Auvergne <bugman>
Project Administrator
Mon 17 Feb 2014 05:05:43 PM UTC, comment #1:

This issue is discussed here:
http://thread.gmane.org/gmane.science.nmr.relax.devel/4956

Troels E. Linnet <tlinnet>
Project Member
Wed 05 Feb 2014 10:16:04 AM UTC, original submission:

Hi Edward.

Relax automatically creates mappings of the various model-free parameters onto the loaded structure by generatiing .pml (or Molmol .mac) scripts/macros. To use them, the PDB structure has to be opened in PyMol, and then the mapping script has to be run.

The problem: when loading the mapping script, all chains of the current molecule are treated the same, i.e. the values are not only mapped to chain A of my multimer, but also onto chain B, C, etc.

The reason is that these mappings are based on residue numbers only. To make one's life easier, all present molecules should be treated individually.

E.g., Instead of:
select pept_bond, (name ca,n AND resi 2) or (name ca,c AND resi 1)

it should read:
select pept_bond, (name ca,n AND resi 2 AND chain A) or (name ca,c AND resi 1 AND chain A)

in the pymol script.

To fix this, one just would have take into account the molecule number that is read by relax:

structure.read_pdb(file='XYZ.pdb', read_mol=1)

"Molecule 1" would translate into "chain A", "molecule 2" to "chain B" , etc in the PyMol script, by looping over the molecules present, assuming all present molecules have been loaded from the same pdb. If the different molecules are loaded from different files, the molecule-chain mapping would make little sense. One way to circumvent this problem could be something like a "multimer" flag to tell relax to specify molecule/chain numbers.

I don't know what the scripts would look like if there are several molecules loaded into relax at the same time. If there is no seperate treatment for them, a fix like this would probably help.

Cheers,
Martin

Martin Ballaschk <mab>

 

(Note: upload size limit is set to 1024 kB, after insertion of the required escape characters.)

Attach File(s):
   
   
Comment:
   

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -unavailable- added by bugman (Posted a comment)
  • -unavailable- added by tlinnet (Posted a comment)
  • -unavailable- added by mab (Submitted the item)
  •  

    Do you think this task is very important?
    If so, you can click here to add your encouragement to it.
    This task has 0 encouragements so far.

    Only logged-in users can vote.

     

    Please enter the title of George Orwell's famous dystopian book (it's a date):

     

     

    Follows 1 latest change.

    Date Changed By Updated Field Previous Value => Replaced By
    Tue 11 Nov 2014 04:36:09 PM UTCbugmanStatusNone=>Postponed
    Show feedback again

    Back to the top


    Powered by Savane 3.1-cleanup