- Chimera-users - RBVI Mailing Lists

Re: [Chimera-users] Scripting to write / save multiple models into single pdf file
by Elaine Meng 08 Mar '10

08 Mar '10

On Mar 7, 2010, at 5:12 AM, Sumitro Harjanto wrote: > Hi Meng, > I would like to use chimera to help me with a process of combining > receptors and ligands that were modelled separately (in separate > files) > > So here is my intended framework: > > For each iteration: > • open “template.pdb” template file (which contains the receptor/ > ligand complex) I used to model the receptor and ligand > • open the receptor_model.pdb > • open the ligand_model.pdb > • match receptor to template > • match ligand to template > • close template.pdb > • save the models in their current position into a single new pdb > file. > > I was exploring the feasibility of using the command “write” but I > am not too sure if it will allow me to save the complex as a single > pdb file..? > > On a side note, if it’s possible, any suggestion on how to optimise > the run time since I might be running this for a large number of > iterations ..? I read somewhere about switching off the GUI to save > time, does it help..? > > Thanks! > Cheers, > Sumitro Hi Sumitro, You can't save multiple models in one file with the "write" command, but you can use "write" multiple times, one for each model you want to save. If you make sure to save each model "relative" to the same model, they will be in the correct positions with respect to each other when you reopen the files. See the "relative" keyword: <http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/write.html> Or, if you really need to have them in a single file, you can use the "combine" command to create a new single model from multiple models. You would do that after the matching, but before saving the PDB file. <http://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/combine.html> I can't answer about how to script for multiple models, or how much faster it might be to use noGUI, so I'm sending this to the chimera- users list. In the future, please send Chimera questions to the chimera-users(a)cgl.ucsf.edu list, because I might be away from work or (like now) others may be better able to answer the question. I hope this helps, ----- Elaine C. Meng, Ph.D. meng(a)cgl.ucsf.edu UCSF Computer Graphics Lab and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco http://www.cgl.ucsf.edu/home/meng/index.html

2 1

Problem with energy minimization
by Bill McIntire 05 Mar '10

05 Mar '10

I used Dock Prep for my modeled structure, and get no error comments. However, when I try to energy-minimize, I get the following message: No MMTK name for atom "HN1" in standard residue "PRO" PRO is at the N-terminus. I cannot find any help in manual. Please help with this problem. Thanks Bill McIntire

2 2

cap color from command line
by Jean-Paul Armache 05 Mar '10

05 Mar '10

Hello, Thank you for all your previous answers. I would like to ask you whether there is a way to change the cap color not only from the Surface capping tool, but the command line as well and declare it maybe in a script (like for example coloring of the maps and models). Thanks, Jean-Paul

2 1

make one chain from two chains ID in PDB
by Dong, Jianbo 04 Mar '10

04 Mar '10

Hello, I'm working with a PDB file, in this file, one protein peptide was recognized as two chains (B and C), I want to make it as one chain, so it will be easier to do some interface analysis and modeling. I guess I can open the file as WORD file and change all the C to B one by one, but it's too much work to edit it that way. May I get some help for this? Is there a easier way to do this? Many thanks, Jianbo

4 5

positive and negative overlaps
by Kshatresh Dutta Dubey 03 Mar '10

03 Mar '10

Dear all, Can anyone tell me the significance of positive and negative overlaps results in structure analysis module of chimera, as i got 140 contacts that contains both positive and negative overlap values? Any help is appreciated. Thanks in advance Kshatresh Dutta Dubey

2 2

show axis ?
by Stefan Bäckström 03 Mar '10

03 Mar '10

Hi. It is possible to show the box of a volume, and I wonder if it is possible to show the X, Y and Z axis as lines ? Best, Stefan Bäckström Medical Biochemistry Göteborg University Sweden

2 1

groups
by Jean-Paul Armache 03 Mar '10

03 Mar '10

Hello, I am often juggling between different sets that are to be turned on/off together. Is there a way of grouping maps/pdbs so that instead of turning 20 pdbs on, i could just click on a group to which they belong ? If there is something like that, could you point me to it ? And if not, would it be possible to have it in the Model Panel? Best regards, Jean-Paul Armache

3 2

chimera-users@cgl.ucsf.edu
by Jozef Lewandowski 03 Mar '10

03 Mar '10

Hi everybody, I am trying to get a list of all hydrogen atoms within certain distance from a specific atom. I found an old thread suggesting using findclash for that. When I use: findclash #0.265:4@c overlap -8 hbond 0 ignoreIntraRes false log true bondSeparation 0 everything works fine and I get the list of all the atoms within the desired distance limits. However, I get all types of atoms. To select only for contacts to hydrogens, I have tried to add the "test" option but when I type: findclash #0.265:4@c overlap -8 hbond 0 ignoreIntraRes false log true bondSeparation 0 test @H the procedure does not report on any contacts. Am I misunderstanding the syntax? What is the way to get only the distances between :4@c and all hydrogens within ~10A? Moreover, I am ultimately interested in the sum of all reported distances (or more specifically sum of (1/r)6 where r is the distance). Is there simple way to get this value calculated in Chimera (e.g. simple script)? Unfortunately, I have no idea about Python coding, so right now I have to import the list to Excel and calculate what I need in there (which is fine for a few atoms but slightly tedious when you have to check several hundred of them...). Any help would be appreciated. Thanks, Józef

3 5

Sum of distances
by Jozef Lewandowski 03 Mar '10

03 Mar '10

Hi everybody, I am trying to get a list of all hydrogen atoms within certain distance from a specific atom. I found an old thread suggesting using findclash for that. When I use: findclash #0.265:4@c overlap -8 hbond 0 ignoreIntraRes false log true bondSeparation 0 everything works fine and I get the list of all the atoms within the desired distance limits. However, I get all types of atoms. To select only for contacts to hydrogens, I have tried to add the "test" option but when I type: findclash #0.265:4@c overlap -8 hbond 0 ignoreIntraRes false log true bondSeparation 0 test @H the procedure does not report on any contacts. Am I misunderstanding the syntax? What is the way to get only the distances between :4@c and all hydrogens within ~10A? Moreover, I am ultimately interested in the sum of all reported distances (or more specifically sum of (1/r)^6 where r is the distance). Is there simple way to get this value calculated in Chimera (e.g. simple script)? Unfortunately, I have no idea about Python coding, so right now I have to import the list to Excel and calculate what I need in there (which is fine for a few atoms but slightly tedious when you have to check several hundred of them...). Any help would be appreciated. Thanks, Józef

1 0

constraints and restraints
by Elaine Meng 02 Mar '10

02 Mar '10

Hello, The Chimera release notes describe what is new or changed in a release: <http://www.cgl.ucsf.edu/chimera/docs/relnotes.html> ... clicking the "1.4" link: <http://www.cgl.ucsf.edu/chimera/docs/relnotes/1.4.html> The meaning of "constraint" "restraint" and "tether" may change based on which program you are using. Each program should explain it in its own documentation. In programs I've used, "constraints" meant atoms are completely fixed in place (frozen, cannot move) whereas "restraints" use a force constant like a spring, and allowed some movement. I don't have any opinion on the technical definition of "tethering." Again, don't rely on my descriptions as final because different programs may use the words to mean different things. Chimera minimization only allows "fixing" some of the atoms so they can't move at all. It does not allow holding atoms in place with force constants. http://www.cgl.ucsf.edu/chimera/docs/ContributedSoftware/minimize/minimize.… > I hope this helps, Elaine ----- Elaine C. Meng, Ph.D. meng(a)cgl.ucsf.edu UCSF Computer Graphics Lab and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco http://www.cgl.ucsf.edu/home/meng/index.html > nice to give you a letter! > is there a new one ,really? > what is the new feature in the least software? > > I want constraint ,restraint and tether some atoms and backbone. > how can I do constraiting ,restrainting and tethering in > software.how can i set value for constraiting ,restrainting and > tethering. > what is the difference between them? > > thank you! > >

1 0