[Chimera-users] Electrostatic Potential

Tom Goddard goddard at sonic.net
Mon Sep 16 15:39:04 PDT 2019


Hi Fernando,

  I don't understand what you need vertex coordinates for.  But you could put a print statement into the Chimera Python code that does surface coloring to output vertex coordinates and potential values.  This would be in file

	 chimera/share/SurfaceColor/__init__.py

or on Mac in 

	Chimera.app/Contents/Resources/share/SurfaceColor/__init__.py

in the volume_values() routine change

    def volume_values(self, surface_piece):

        p = surface_piece
        xf = p.model.openState.xform
        v = p.geometry[0]
        n = p.normals
        return self.offset_values(v, n, xf)

to

    def volume_values(self, surface_piece):

        p = surface_piece
        xf = p.model.openState.xform
        v = p.geometry[0]
        n = p.normals
        values, outside = self.offset_values(v, n, xf)
	print '\n'.join('%.5g %.5g %.5g %.5g'  % (x, y, z, val) for (x,y,z), val in zip(v, values))
        return values, outside

Use the Surface Color tool and the vertices and values will be printed to the reply log.

  Tom


> On Sep 16, 2019, at 1:20 PM, Fernando Villa <fer.vdl1928 at gmail.com> wrote:
> 
> Dear Tom and Chimera users,
> 
> In ChimeraX or Chimera, is it possible to calculate the number of vertices and their respective total charges?
> 
> I calculated the electrostatic potential surface (Coulombic) with another program, it gave me the following results (attached in this email) and I calculated the area in red and blue, above and below a threshold = 0.
> 
> Unfortunately, these results do not help me, because I need to obtain an electrostatics calculations from Adaptive Poisson-Boltzmann Solver (Calculations that I have already made for study proteins with PDB2PQR and APBS, in interface with Chimera: .pqr and .dx files and emulated in Chimera X)
> 
> e.g.
> <image.png>
> 
> Now my question is:
> Can I obtain these results of vertices and charges in a file from Chimera or Chimera X (possibly .cvs, .txt, etc.), so that I can make the corresponding calculations?
> 
> I would greatly appreciate your help
> 
> best regards,
> 
> Fer
> 
> 
> 
> El lun., 16 de sep. de 2019 a la(s) 13:03, Tom Goddard (goddard at sonic.net <mailto:goddard at sonic.net>) escribió:
> Hi Fernando,
> 
>   Instead of thinking about the coloring, it would be clearer to ask how much surface area has electrostatic potential value >= 5.  Chimera does not have any code to compute that, although it is something we could add to ChimeraX.
> 
> 	Tom
> 
> 
>> On Sep 11, 2019, at 7:27 PM, Fernando Villa <fer.vdl1928 at gmail.com <mailto:fer.vdl1928 at gmail.com>> wrote:
>> 
>> Dear all Chimera users
>> 
>> In UCSF ChimeraX version 0.91 (2019-08-30)
>> I opened a .pqr file (generated with pdb2pqr in Chimera 1.14)  and then I opened a .dx file (generated with APBS in Chimera 1.14)
>> 
>> Then I input the commands:
>> 
>> >surface
>> >color electrostatic #1 map #2 palette -10,red:0,white:10,blue
>> 
>> <image.png>
>> Is it possible to calculate the red and blue area in Å2 of the molecule?
>> 
>> Best regards,
>> 
>> Fernando
>> 
>> 
>> ATTE
>> Fernando Villa Díaz
>> 
>> 
>> _______________________________________________
>> Chimera-users mailing list: Chimera-users at cgl.ucsf.edu <mailto:Chimera-users at cgl.ucsf.edu>
>> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users <http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users>
> 
> <Test_MEP_Protein.xlsx>_______________________________________________
> Chimera-users mailing list: Chimera-users at cgl.ucsf.edu
> Manage subscription: http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://plato.cgl.ucsf.edu/pipermail/chimera-users/attachments/20190916/02b021d1/attachment.html>


More information about the Chimera-users mailing list