- Chimera-users - RBVI Mailing Lists

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

color zone on command line
by Kyle Morris 24 Dec '16

24 Dec '16

Hi Chimera dev, Following on from this post back from 2009: http://plato.cgl.ucsf.edu/pipermail/chimera-users/2009-January/003514.html <http://plato.cgl.ucsf.edu/pipermail/chimera-users/2009-January/003514.html> How would one also split map in this python script? Thanks! Kyle

2 1

Re: [Chimera-users] Chimera-users Digest, Vol 164, Issue 26
by Healey, Joe 23 Dec '16

23 Dec '16

I might be able to offer a bit of help here, the Chimera team and Jaime helped me with basically this exact problem. I have a script here https://github.com/jrjhealey/bioinfo-tools/blob/master/strucfit.py that you could take from. You'll probably want to look specifically at lines 216-240. (Note - you'll also need to look at the relevant module imports for match and pychimera etc. The script will need Jaime's implementation of pychimera too (he's super helpful!) https://github.com/insilichem/pychimera The script gives output as below, as the script automatically finds the best matching PDB file (that's what the first 150 lines or so is looking for), there is additional information there. Best Model I'm Hit Comparing RMSD Assorted homology metrics 2k4r PAK_01796_model5.pdb 2.81843275676 37.5 13 0.00035 21.5 2k4r PAK_01796_model4.pdb 12.6869844869 37.5 13 0.00035 21.5 2k4r PAK_01796_model3.pdb 14.4133021316 37.5 13 0.00035 21.5 2k4r PAK_01796_model1.pdb 0.459177533295 37.5 13 0.00035 21.5 2k4r PAK_01796_model2.pdb 15.1861411912 37.5 13 0.00035 21.5 Hope that's useful! Happy to offer more help (though I'm sure the chimera team have it covered!) Joe M.Sc. B.Sc. (Hons) MSRB PhD Student MOAC CDT, Senate House University of Warwick Coventry CV47AL Mob: +44 (0) 7536 042620 | Email: J.R.J.Healey(a)warwick.ac.uk Jointly working in: Waterfield Lab<http://www2.warwick.ac.uk/fac/med/research/tsm/microinfect/staff/waterfield…> (WMS Microbiology and Infection Unit) and the Gibson Lab<http://www2.warwick.ac.uk/fac/sci/chemistry/research/gibson/gibsongroup/> (Warwick Chemistry) Twitter: @JRJHealey<https://twitter.com/JRJHealey> | Website: MOAC Page<http://www2.warwick.ac.uk/fac/sci/moac/people/students/2013/joseph_healey> | ORCID: orcid.org/0000-0002-9569-6738 ________________________________ From: chimera-users-bounces(a)cgl.ucsf.edu <chimera-users-bounces(a)cgl.ucsf.edu> on behalf of chimera-users-request(a)cgl.ucsf.edu <chimera-users-request(a)cgl.ucsf.edu> Sent: 22 December 2016 20:00 To: chimera-users(a)cgl.ucsf.edu Subject: Chimera-users Digest, Vol 164, Issue 26 Send Chimera-users mailing list submissions to chimera-users(a)cgl.ucsf.edu To subscribe or unsubscribe via the World Wide Web, visit http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users or, via email, send a message with subject or body 'help' to chimera-users-request(a)cgl.ucsf.edu You can reach the person managing the list at chimera-users-owner(a)cgl.ucsf.edu When replying, please edit your Subject line so it is more specific than "Re: Contents of Chimera-users digest..." Today's Topics: 1. Saving Output From Match Align (Muk, Sanychen) ---------------------------------------------------------------------- Message: 1 Date: Thu, 22 Dec 2016 16:11:00 +0000 From: "Muk, Sanychen" <smuk(a)coh.org> To: "chimera-users(a)cgl.ucsf.edu" <chimera-users(a)cgl.ucsf.edu> Subject: [Chimera-users] Saving Output From Match Align Message-ID: <A69D81CB04F0734E96302AA9E36F9B233E521BEE(a)ppwexch2kx03.coh.org> Content-Type: text/plain; charset="iso-8859-1" Hello, I'm working on a script to superimpose pdb structures onto one another and then run the match -> align tool to align them to get the RMSD between residues of the pdb structures. I know that from the multialign interviewer window, I can then go to structure acess match and then save the attributes as a text file that will contain the residue number and RMSD. Is there a command that I can use to output this attributes file automatically? Thank you, and happy holidays. Best, Sanychen Muk --------------------------------------------------------------------- *SECURITY/CONFIDENTIALITY WARNING: This message and any attachments are intended solely for the individual or entity to which they are addressed. This communication may contain information that is privileged, confidential, or exempt from disclosure under applicable law (e.g., personal health information, research data, financial information). Because this e-mail has been sent without encryption, individuals other than the intended recipient may be able to view the information, forward it to others or tamper with the information without the knowledge or consent of the sender. If you are not the intended recipient, or the employee or person responsible for delivering the message to the intended recipient, any dissemination, distribution or copying of the communication is strictly prohibited. If you received the communication in error, please notify the sender immediately by replying to this message and deleting the message and any accompanying files from your system. If, due to the security risks, you do not wi! sh to receive further communications via e-mail, please reply to this message and inform the sender that you do not wish to receive further e-mail from the sender. (fpc5p) ---------------------------------------------------------------------

1 0