- Chimera-users - RBVI Mailing Lists

Web services (Modeller, Blast, etc.) running
by Elaine Meng 13 Jan '17

13 Jan '17

Hello users, I believe the RBVI web services called from Chimera are now running, thanks to the efforts of my coworkers. I just now verified running Modeller, Blast (requires Chimera 1.11.2 or later), and sequence realignment with Clustal Omega or Muscle. Best, Elaine ----- Elaine C. Meng, Ph.D. UCSF Computer Graphics Lab (Chimera team) and Babbitt Lab Department of Pharmaceutical Chemistry University of California, San Francisco

2 1

Deleting a surface on command line, not just hiding it
by Rebecca Swett 12 Jan '17

12 Jan '17

Hi all, I'm working on an animation, and I'd like to know if there's a way I can delete or close a calculated surface from the command line. For my purposes, I'd like to show the surface representation of the protein, remove it, then show a dynamics trajectory, and show the surface again at a separate point (while the trajectory is stopped). Currently, without being able to close the surface, it tries to recalculate the surface at each trajectory frame, which doesn't work for my purposes. Let me know if this is possible. Cheers, Rebecca This email message and any attachments are confidential and intended for use by the addressee(s) only. If you are not the intended recipient, please notify me immediately by replying to this message, and destroy all copies of this message and any attachments. Thank you.

2 2

Re: [Chimera-users] dipoles
by Sagar, Vatsala (NIH/NEI) [E] 12 Jan '17

12 Jan '17

Hi Tom, So if you align the pdb files of the rotated proteins, they have an exact match. But when you make the PQR files, it adds hydrogens differently to the rotated proteins, resulting in a charge difference in the said histidine. Should the program be adding hydrogens differently to the rotated proteins, and is there a way to remedy this> Thanks a lot, Vatsala Sagar From: "Sagar, Vatsala (NIH/NEI) [F]" <vatsala.sagar(a)nih.gov <mailto:vatsala.sagar@nih.gov>> Date: Friday, January 6, 2017 at 10:27 AM To: "chimera-bugs(a)cgl.ucsf.edu <mailto:chimera-bugs@cgl.ucsf.edu>" <chimera-bugs(a)cgl.ucsf.edu <mailto:chimera-bugs@cgl.ucsf.edu>> Subject: FW: [Chimera-users] dipoles Hi Tom, I am attaching the PQR files and raw PDB files of the crystal structure of human gamma D with two monomers rotated 90 degrees to one another. Please take a look at it and tell me what you think, Thanks, Vatsala From: Tom Goddard <goddard(a)sonic.net <mailto:goddard@sonic.net>> Date: Thursday, January 5, 2017 at 7:07 PM To: "chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu> BB" <chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu>>, "Sagar, Vatsala (NIH/NEI) [F]" <vatsala.sagar(a)nih.gov <mailto:vatsala.sagar@nih.gov>> Subject: Re: [Chimera-users] dipoles Hi Vatsala, I guess you are refering to the dipole.py script on the Chimera Python scripts web page: http://plato.cgl.ucsf.edu/trac/chimera/wiki/Scripts <http://plato.cgl.ucsf.edu/trac/chimera/wiki/Scripts> Yes it should give the same answer for a molecule and a copy of that molecule rotated by 90 degrees. Without have the structure and rotated PQR files I don’t know what is going wrong. Tom > On Jan 4, 2017, at 11:17 AM, Sagar, Vatsala (NIH/NEI) [F] wrote: > > Hi, > Do you know if there is a method to calculate dipoles in Chimera? I have generated a PQR structure in predicted protonation state and uploaded it onto Chimera. I then used a downloaded app called dimer.py to calculate the dipole moments of two of the same monomers rotated 90degrees. Technically they should both give the same dipole as the structures align 100% but according to the PDB coordinates are only rotated relative to one another. The dipole.py scipt is calculating them as two different dipoles. Can you help me with this? > Thanks, > Vatsala Sagar > _______________________________________________ > Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu <mailto:Chimera-users@cgl.ucsf.edu> > Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users <http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users>

2 1

Re: [Chimera-users] dipoles
by Tom Goddard 10 Jan '17

10 Jan '17

Hi Vatsala, If you run the Chimera pdb2pqr command or graphical interface on a molecule and a rotated copy of the same molecule and it gives different protonation states for histidine that seems like a bug. As described in the Chimera documentation this uses the PDB2PQR web service provided by the National Biomedical Computation Resource (NBCR). https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/pdb2pqr.html <https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/pdb2pqr.html> You might want to contact the developers of PDB2PQR to pursue this matter further http://www.poissonboltzmann.org/support/home/ <http://www.poissonboltzmann.org/support/home/> Chimera can also add hydrogens and compute charges using the Chimera DockPrep tool so maybe that is an option for you to avoid this problem in pdb2pqr. Tom > On Jan 10, 2017, at 7:35 AM, Sagar, Vatsala (NIH/NEI) [E] <vatsala.sagar(a)nih.gov> wrote: > > Hi Tom, > So if you align the pdb files of the rotated proteins, they have an exact match. But when you make the PQR files, it adds hydrogens differently to the rotated proteins, resulting in a charge difference in the said histidine. Should the program be adding hydrogens differently to the rotated proteins, and is there a way to remedy this> > Thanks a lot, > Vatsala Sagar > > From: "Sagar, Vatsala (NIH/NEI) [F]" <vatsala.sagar(a)nih.gov <mailto:vatsala.sagar@nih.gov>> > Date: Friday, January 6, 2017 at 10:09 AM > To: Tom Goddard <goddard(a)sonic.net <mailto:goddard@sonic.net>>, "chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu> BB" <chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu>>, "Sagar, Vatsala (NIH/NEI) [F]" <vatsala.sagar(a)nih.gov <mailto:vatsala.sagar@nih.gov>> > Subject: Re: [Chimera-users] dipoles > > Hi Tom, > I am attaching the PQR files and raw PDB files of the crystal structure of human gamma D with two monomers rotated 90 degrees to one another. Please take a look at it and tell me what you think, > Thanks, > Vatsala > From: Tom Goddard <goddard(a)sonic.net <mailto:goddard@sonic.net>> > Date: Thursday, January 5, 2017 at 7:07 PM > To: "chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu> BB" <chimera-users(a)cgl.ucsf.edu <mailto:chimera-users@cgl.ucsf.edu>>, "Sagar, Vatsala (NIH/NEI) [F]" <vatsala.sagar(a)nih.gov <mailto:vatsala.sagar@nih.gov>> > Subject: Re: [Chimera-users] dipoles > > Hi Vatsala, > > I guess you are refering to the dipole.py script on the Chimera Python scripts web page: > > http://plato.cgl.ucsf.edu/trac/chimera/wiki/Scripts <http://plato.cgl.ucsf.edu/trac/chimera/wiki/Scripts> > > Yes it should give the same answer for a molecule and a copy of that molecule rotated by 90 degrees. Without have the structure and rotated PQR files I don’t know what is going wrong. > > Tom > > >> On Jan 4, 2017, at 11:17 AM, Sagar, Vatsala (NIH/NEI) [F] wrote: >> >> Hi, >> Do you know if there is a method to calculate dipoles in Chimera? I have generated a PQR structure in predicted protonation state and uploaded it onto Chimera. I then used a downloaded app called dimer.py to calculate the dipole moments of two of the same monomers rotated 90degrees. Technically they should both give the same dipole as the structures align 100% but according to the PDB coordinates are only rotated relative to one another. The dipole.py scipt is calculating them as two different dipoles. Can you help me with this? >> Thanks, >> Vatsala Sagar >> _______________________________________________ >> Chimera-users mailing list: Chimera-users(a)cgl.ucsf.edu <mailto:Chimera-users@cgl.ucsf.edu> >> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users <http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users> >

1 0

Map conservation of my sequence compared to sequence alignment
by mfellner＠msu.edu 09 Jan '17

09 Jan '17

Good day I like the mapping sequence conservation tool in chimera, however I ran into a problem that IÂ could not overcome with anything IÂ found online. I have the consensus mapped onto my protein but that is not what IÂ want. IÂ want to show how the residues of my protein are conserved. For example at a particular residue the consensus is a C but in my structure its a W. at the moment the W is colored to be highly conserved based on the consenus of the C. The W is maybe 15% conserved, how can IÂ map the conservation of my sequence compared to all the other sequences. IÂ hope IÂ explained what my issue is Thank you for this great tool Matthias Fellner

3 2

dipoles
by Sagar, Vatsala (NIH/NEI) [F] 06 Jan '17

06 Jan '17

Hi, Do you know if there is a method to calculate dipoles in Chimera? I have generated a PQR structure in predicted protonation state and uploaded it onto Chimera. I then used a downloaded app called dimer.py to calculate the dipole moments of two of the same monomers rotated 90degrees. Technically they should both give the same dipole as the structures align 100% but according to the PDB coordinates are only rotated relative to one another. The dipole.py scipt is calculating them as two different dipoles. Can you help me with this? Thanks, Vatsala Sagar

2 3

Mutation in Chimera
by Islam, Tanzila 06 Jan '17

06 Jan '17

Hi I am learning Chimera and I have gone through soe online tutorials. The tutorials showed how to introduce point mutation by Rotamers command. But it only allows to replace 1 ammino acid residue at a time. I want to insert 4-5 ammino acid residues in a gap of a sequence. Is there any standard procedure, command in Chimera that will let me do this type of residue insertion? This is really crucial for my research work. Thanks Tanzila Islam Graduate Student Washington State University

2 6

Problem opening DOI from Protein Data Bank
by Tones,Jason R 05 Jan '17

05 Jan '17

Hello, I am a student and recently downloaded chimera to use for a class, when I clicked to open the DOI of a search on the protein database I picked the wrong file under chimera to use to open it and it wouldn't open. I tried uninstalling and reinstalling and it keeps going to the file that won't open what I'm trying to do. Thank you for any help or suggestions you may have for me. Jason

2 1

Measuring hexamer asymmetry
by Tom Goddard 05 Jan '17

05 Jan '17

Hi Elena, Ok. The “measure rotation” command shows the rotation axis defined by the rigid motion of one subunit onto another — it has nothing to do with the hexamer axis. I think you are asking for is that if you take subunit 1, rotate it 60 degrees about the hexamer axis, then compare the orientation and translation of subunit 2 to that moved copy of subunit 1. Then do the same for subunits 3,4,5,6. To do this take subunit 1 and make 5 rotated copies by 60,120,180,240,300 degrees and save that as a new hexamer PDB file. Now open the asymmetric hexamer and this symmetric hexamer. If the subunits are chain ids A,B,C,D,E,F then use “measure rotation #0:.B #1:.B”, “measure rotation #0:.C #1:.C", …. This will still show a rotation axis that is not the hexamer axis which indicates the rotation needed to make the exact symmetry copy of subunit 1 rotate onto the asymmetric subunit 2, and it will give the translation, rotation axis and angle in the reply log. Maybe you are not interested in how the symmetric copies of subunit 1 rotate when being aligned to asymmetric subunits 2,3,4…. Then you can just ignore this rotation information and get the translation from the reply log. Tom > On Jan 5, 2017, at 7:38 AM, Elena Zehr wrote: > > Hi Tom, > > Thank you very much for you response. Yes within the asymmetric hexamer subunits have different conformations. To test your advice on RMSD, I first generated hexamer by using six copies of the same subunit to align to the original hexametric structure, then I closed the original hexametric structure. The asymmetric hexamer now consisted of subunits with identical conformation, please note the hexamer is still asymmetric. Then I opened a monomer with the same conformation as the subunits within the newly-generated hexamer, and calculated rotation and translation matrices between the first protomer within the hexamer and the rest of the subunits (#2-6). The central axis for each of the calculations deviated from each other. > > My issue is the following: to make the calculations reliable I wish to have the same central axis running through the center of the hexamer and not the axes between #1 and #2 or #1 and #4, etc. which deviate from each other due to the asymmetry of the hexametric structure. Is it possible? For example I wish to calculate the central axis of the whole hexamer and then use it as a reference to calculate rotation and translation for each subunit (#2-5) around the ring with respect to the subunit #1 and the central axis of the whole hexamer. Is it possible to define parameters for the central axis in the command line together with “measure rotation #1 #2 or #1 and #4” ? > > > Many thanks, > Elena > > > On Jan 4, 2017, at 4:47 PM, Tom Goddard wrote: > >> Hi Elena, >> >> Elaine asked me to look at your problem measuring orientation of subunits of a hexamer. Does your hexamer have cyclic symmetry? If so I think the reason the displayed axes don’t appear to align with the hexamer symmetry axis is because your alignment of one subunit to another is not good. You can get the root mean square deviation of C-alpha positions for the alignment by looking in the Reply Log (menu Favorites) after running the match maker command, for example, >> >> RMSD between 308 pruned atom pairs is 0.361 angstroms; (across all 308 pairs: 0.361) >> >> From your picture with blue and tan complexes it does not appear they are aligned at all. If the sequences of the subunits are very different then match maker may not be able to align them. >> >> Tom >> >> >> >> Begin forwarded message: >> >>> From: Elena Zehr >>> Subject: Re: user-defined rotation axis >>> Date: January 4, 2017 at 1:09:10 PM PST >>> To: Elaine Meng >>> >>> Dea Elaine, >>> >>> Thank you for getting back to me so soon. Sorry for sending you my question again. >>> >>> I am getting the following (see attached image) Where the central axis between protomer A and B is different than the one between protomer A and F. I wish to have the same central axis that would be the central axis of the hexamer, otherwise I assume my measurements are not accurate. I get the same if I use the hexamers or separate monomers >>> >>> >> <axis.tiff> >>> On Jan 4, 2017, at 1:15 PM, Elaine Meng wrote: >>> >>> Begin forwarded message: >>> >>>> From: Elena Zehr >>>> Subject: Re: user-defined rotation axis >>>> Date: January 4, 2017 at 1:09:10 PM PST >>>> To: Elaine Meng >>>> >>>> Dea Elaine, >>>> >>>> Thank you for getting back to me so soon. Sorry for sending you my question again. >>>> >>>> I am getting the following (see attached image) Where the central axis between protomer A and B is different than the one between protomer A and F. I wish to have the same central axis that would be the central axis of the hexamer, otherwise I assume my measurements are not accurate. I get the same if I use the hexamers or separate monomers >>>> >>>> >>> <axis.tif> >> >

1 0

Showing color key programatically colouring from a volume...
by Pablo Conesa 31 Dec '16

31 Dec '16

Dear all We have this simple script to open a volume (#0) coloured using values from a second volume (#1): open betaGal.mrc open extra/MG_Chimera_resolution.vol volume #1 voxelSize 3.54 vol #1 hide scolor #0 volume #1 cmap rainbow reverseColors True For scolor command we are using a volume. https://www.cgl.ucsf.edu/chimera/docs/UsersGuide/midas/scolor.html#volume It would be nice to achieve something like this at start up (note the color key is based on the volume): I've did it dragging the mouse over the canvas while "Color key" utility is open. And, is there a posibility to round up the values. Some kind of: Psaudocode: Give the resulting color key for each round values display the color key I guess that if there is a way to get the color key, then rest can be done through commands. Cheers, Pablo.

2 3