<div dir="ltr">
<div><span lang="en">
</span><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;line-height:normal"><font size="2"><span style="font-family:arial,sans-serif"><span lang="EN">Dear Elaine.</span></span></font></p><p class="MsoNormal" style="margin:0cm 0cm 0.0001pt;line-height:normal;font-size:11pt"><span style="font-family:arial,sans-serif"><br><span style="font-size:12pt" lang="EN"><span></span></span></span></p><span style="font-family:arial,sans-serif">
</span></div><div><span style="font-family:arial,sans-serif"><span lang="en">Thank you very much, I really appreciate it.</span></span><span style="font-family:arial,sans-serif"><span lang="en"></span></span>
</div><div><span style="font-family:arial,sans-serif"><span lang="en"></span></span></div><span style="font-family:arial,sans-serif"><span lang="en">I am confused with some data.<br>A few months ago <span style="color:rgb(0,0,255)">Tom Goddard </span>sent me a phyton code to get the vertex coordinates and the values in Reply log:<br><br><span style="color:rgb(0,0,255)">In Windows:<br> chimera /share/SurfaceColor/__init__.py<br>or on Mac in:<br>Chimera.app/Contents/Resources/share/SurfaceColor/__init__.py<br><br>in the volume_values () routine change<br><br> def volume_values (self, surface_piece):<br><br> p = surface_piece<br> xf = p.model.openState.xform<br> v = p.geometry [0]<br> n = p.normals<br> return self.offset_values (v, n, xf)<br>to<br> def volume_values (self, surface_piece):<br><br> p = surface_piece<br> xf = p.model.openState.xform<br> v = p.geometry [0]<br> n = p.normals<br> values, outside = self.offset_values (v, n, xf)<br>print '\ n'.join ('%. 5g% .5g% .5g% .5g '% (x, y, z, val) for (x, y, z), val in zip (v, values))<br> return values, outside<br></span><br><span style="color:rgb(56,118,29)">To obtain the output vertex coordinates and potential values:<br>Commands:<br>open 1avx<br>display<br>~ longbond<br>split<br>delete solvent<br>surf # 0.1<br><br>Tools> Surface/Binding analysis>Coulombic surface coloring<br> </span><span style="color:rgb(56,118,29)"> > Check the box compute grid<br> </span><span style="color:rgb(56,118,29)"> > Apply<br>>Surface color>Click color<br></span></span></span><div><span style="font-family:arial,sans-serif"><span style="color:rgb(56,118,29)"><span lang="en">= All values Appear in reply log<span lang="en"> (</span><span lang="en">attached file</span><span lang="en"></span>: 001_xyz-val_ESP.txt)<br></span></span></span></div><div><div><span style="font-family:arial,sans-serif"><span style="color:rgb(56,118,29)"><img src="cid:ii_k4xnvgi32" alt="imagen.png" style="margin-right:0px" width="110" height="284"><br></span></span></div><span style="font-family:arial,sans-serif"><span lang="en"><span style="color:rgb(56,118,29)">I suppose that the first three numbers are the vertex coordinates and the last value is the ESP Coulombic (Kcal / mol.e). </span><br></span></span></div><div><span style="font-family:arial,sans-serif"><span lang="en"></span><br><span lang="en"></span></span></div><span style="font-family:arial,sans-serif"><span lang="en"><span style="color:rgb(255,153,0)">Value at mouse position:<br>> Surface color> Options<br> > Check the box Report value at mouse position<br></span></span></span><div><span style="font-family:arial,sans-serif"><span lang="en"><span style="color:rgb(255,153,0)"> > Click color:<br></span></span></span></div><div><div><span style="font-family:arial,sans-serif"><img src="cid:ii_k4xnpgdy0" alt="imagen.png" width="377" height="22"></span></div><span style="font-family:arial,sans-serif"><span lang="en"><span style="color:rgb(255,153,0)"></span></span></span></div><div><span style="font-family:arial,sans-serif"><span lang="en"><span style="color:rgb(255,153,0)">All the electrostatic potential values appear by Coulombic [Kcal / (mol.e)] at the mouse position, I suppose they are the same values shown in replylog = file 001_xyz-val_ENG.txt.</span><br></span></span></div><div><span style="font-family:arial,sans-serif"><span lang="en"><br></span></span></div><span style="font-family:arial,sans-serif"><span lang="en"><span style="color:rgb(153,0,255)">However, when I open the surfval.py script<br>Command:<br>surfvalues # 0.1 # 1 ~/Desktop /002_xyz-val_script-val_ENG.txt<br><br></span></span></span><div><span style="font-family:arial,sans-serif"><span lang="en"><span style="color:rgb(153,0,255)">A file is saved (</span></span><span style="color:rgb(153,0,255)"><span lang="en">attached file</span></span><span lang="en"><span style="color:rgb(153,0,255)">: 002_xyz-val_script-val_ESP.txt) with the same vertex coordinates. The ESP values in Kcal/(mol.e) are different from 001_xyz-val_ESP.txt.</span></span><br><span lang="en"><span style="color:rgb(153,0,255)"></span></span></span><div><div id="m_-3947095932993323525gmail-tw-target"><div id="m_-3947095932993323525gmail-kAz1tf"><div id="m_-3947095932993323525gmail-tw-target-text-container"><pre id="m_-3947095932993323525gmail-tw-target-text" style="text-align:left" dir="ltr"><span style="font-family:arial,sans-serif"><span lang="en"><span style="color:rgb(153,0,255)"></span></span></span></pre></div></div></div></div></div><div><div><span style="font-family:arial,sans-serif"><img src="cid:ii_k4xnsa271" alt="imagen.png" style="margin-right:0px" width="272" height="119"><br><br></span></div><span style="font-family:arial,sans-serif"><span lang="en"><span style="color:rgb(153,0,255)"></span></span></span></div><div><span style="font-family:arial,sans-serif"><span style="color:rgb(120,63,4)"><span lang="en">Besides, if I save the values from Volume viewer:</span><br><span lang="en"></span></span></span></div><span style="font-family:arial,sans-serif"><span style="color:rgb(120,63,4)"><span lang="en">Volume viewer> Tools> Values at atom position <br>> Click Histogram> Render / Select by attribute<br></span></span></span><div><span style="font-family:arial,sans-serif"><span style="color:rgb(120,63,4)"><span lang="en"> > File> Save attribute> Click box “Restrict save to current selection, if any”> File name: 003_@val_ENG.txt> Save</span></span></span><span style="color:rgb(120,63,4)"><span style="font-family:arial,sans-serif"><span lang="en"> (attached file:
003_@val_ENG.txt)</span></span></span></div><div><span style="color:rgb(120,63,4)">
The ESP values --> Kcal/(mol.e)
<br><span style="font-family:arial,sans-serif"><span lang="en"></span></span><span style="font-family:arial,sans-serif"><span lang="en"></span></span></span></div><div><span style="font-family:arial,sans-serif"><span style="color:rgb(120,63,4)"><span lang="en"><br></span></span></span></div><span style="font-family:arial,sans-serif"><span lang="en">My question is:<br></span></span><div><span style="font-family:arial,sans-serif"><span lang="en">The files 001_xyz-val_ESP.txt, 002_xyz-val_script-val_ENP.txt, 003_@val_ESP.txt, and the values at mouse position have units in Kcal/(mol.e)?</span></span></div><div><span style="font-family:arial,sans-serif"><span lang="en"><br></span></span></div><div><span style="font-family:arial,sans-serif"><span lang="en">
If I'm wrong, let me know please. <br></span></span></div><div><span style="font-family:arial,sans-serif"><span lang="en"><br></span></span></div><div><span style="font-family:arial,sans-serif"><span lang="en">Regards.</span></span></div><div><span style="font-family:arial,sans-serif"><span lang="en"><br></span></span></div><div><span style="font-family:arial,sans-serif"><span lang="en">Fer.<br></span></span></div><span style="font-family:arial,sans-serif">
</span><div><div dir="ltr" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="text-align:center"><font face="monospace, monospace"><br></font></div></div></div></div></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El lun., 30 de dic. de 2019 a la(s) 11:56, Elaine Meng (<a href="mailto:meng@cgl.ucsf.edu" target="_blank">meng@cgl.ucsf.edu</a>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Fernando,<br>
You already calculated the values of Coulombic ESP on the interface: it is displayed. However, I guess you mean that you want to write the values that you have already shown with coloring to some file. There is a “surfval.py” python script in the Chimera scripts page that is supposed to do that, although it may be the whole surface (not just the interface), see:<br>
<<a href="http://www.cgl.ucsf.edu/trac/chimera/wiki/Scripts" rel="noreferrer" target="_blank">http://www.cgl.ucsf.edu/trac/chimera/wiki/Scripts</a>><br>
<br>
The Coulombic Surface Coloring tool (and command) also has an option to open the grid of Coulombic ESP values as as separate model. This grid is what is interpolated to get the values for surface coloring. This separate grid model can be saved using the Volume Viewer tool’s File menu.<br>
<<a href="http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/coulombic/coulombic.html" rel="noreferrer" target="_blank">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/coulombic/coulombic.html</a>><br>
<<a href="http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/coulombic/coulombic.html#volume" rel="noreferrer" target="_blank">http://www.rbvi.ucsf.edu/chimera/docs/ContributedSoftware/coulombic/coulombic.html#volume</a>><br>
<<a href="http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/coulombic.html" rel="noreferrer" target="_blank">http://www.rbvi.ucsf.edu/chimera/docs/UsersGuide/midas/coulombic.html</a>><br>
<br>
Getting polar/nonpolar surface area is the same question you asked in September. Here were the responses:<br>
<<a href="http://www.cgl.ucsf.edu/pipermail/chimera-users/2019-September/016139.html" rel="noreferrer" target="_blank">http://www.cgl.ucsf.edu/pipermail/chimera-users/2019-September/016139.html</a>><br>
<<a href="http://www.cgl.ucsf.edu/pipermail/chimera-users/2019-September/016159.html" rel="noreferrer" target="_blank">http://www.cgl.ucsf.edu/pipermail/chimera-users/2019-September/016159.html</a>><br>
<br>
I hope this helps,<br>
Elaine<br>
-----<br>
Elaine C. Meng, Ph.D.<br>
UCSF Chimera(X) team<br>
Department of Pharmaceutical Chemistry<br>
University of California, San Francisco<br>
<br>
> On Dec 30, 2019, at 9:35 AM, Fernando Villa <<a href="mailto:fer.vdl1928@gmail.com" target="_blank">fer.vdl1928@gmail.com</a>> wrote:<br>
> <br>
> Hi all Chimera users,<br>
> <br>
> <br>
> I would like to know how to calculate the values of Coulombic ESP on the interface between two proteins:<br>
> <br>
> <br>
> At the moment I have done the following:<br>
> <br>
> *UCSF Chimera version 1.14<br>
> <br>
> Commands:<br>
> <br>
> open 1avx<br>
> display<br>
> ~longbond <br>
> split<br>
> delete solvent <br>
> surf #0.1<br>
> <br>
> Tools>Surface/Binding analysis>Coulombic surface coloring<br>
> >Check the box compute grid<br>
> >Apply<br>
> >Surface color>Options<br>
> >Check the box Report value at mouse position<br>
> >Click color:<br>
> <br>
> All the electrostatic potential values appear by Coulombic Kcal/(Mol.e) at the mouse position:<br>
> <br>
> <imagen.png><br>
> <br>
> <br>
> Commands:<br>
> ~surface #0.1<br>
> measure buriedArea #0.1 #0.2<br>
> select #0.1@/buriedSASArea> 1<br>
> surf sel<br>
> <br>
> <br>
> Volume viewer>Tools>Values at atom position>Click Histogram<br>
> <br>
> <br>
> <imagen.png><br>
> <br>
> <br>
> <br>
> Is there a way to save all the ESP Coulombic values in the selected interface surface in a file?<br>
> <br>
> Is it possible to calculate the polar and non-polar surface area in the interface?<br>
> <br>
> Regards,<br>
> <br>
> Fer.<br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> El lun., 16 de sep. de 2019 a la(s) 17:39, Tom Goddard (<a href="mailto:goddard@sonic.net" target="_blank">goddard@sonic.net</a>) escribió:<br>
> Hi Fernando,<br>
> <br>
> 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<br>
> <br>
> chimera/share/SurfaceColor/__init__.py<br>
> <br>
> or on Mac in <br>
> <br>
> Chimera.app/Contents/Resources/share/SurfaceColor/__init__.py<br>
> <br>
> in the volume_values() routine change<br>
> <br>
> def volume_values(self, surface_piece):<br>
> <br>
> p = surface_piece<br>
> xf = p.model.openState.xform<br>
> v = p.geometry[0]<br>
> n = p.normals<br>
> return self.offset_values(v, n, xf)<br>
> <br>
> to<br>
> <br>
> def volume_values(self, surface_piece):<br>
> <br>
> p = surface_piece<br>
> xf = p.model.openState.xform<br>
> v = p.geometry[0]<br>
> n = p.normals<br>
> values, outside = self.offset_values(v, n, xf)<br>
> print '\n'.join('%.5g %.5g %.5g %.5g' % (x, y, z, val) for (x,y,z), val in zip(v, values))<br>
> return values, outside<br>
> <br>
> Use the Surface Color tool and the vertices and values will be printed to the reply log.<br>
> <br>
> Tom<br>
> <br>
> <br>
>> On Sep 16, 2019, at 1:20 PM, Fernando Villa <<a href="mailto:fer.vdl1928@gmail.com" target="_blank">fer.vdl1928@gmail.com</a>> wrote:<br>
>> <br>
>> Dear Tom and Chimera users,<br>
>> <br>
>> In ChimeraX or Chimera, is it possible to calculate the number of vertices and their respective total charges?<br>
>> <br>
>> 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.<br>
>> <br>
>> 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)<br>
>> <br>
>> e.g.<br>
>> <image.png><br>
>> <br>
>> Now my question is:<br>
>> 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?<br>
>> <br>
>> I would greatly appreciate your help<br>
>> <br>
>> best regards,<br>
>> <br>
>> Fer<br>
>> <br>
>> <br>
>> <br>
>> El lun., 16 de sep. de 2019 a la(s) 13:03, Tom Goddard (<a href="mailto:goddard@sonic.net" target="_blank">goddard@sonic.net</a>) escribió:<br>
>> Hi Fernando,<br>
>> <br>
>> 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.<br>
>> <br>
>> Tom<br>
>> <br>
>> <br>
>>> On Sep 11, 2019, at 7:27 PM, Fernando Villa <<a href="mailto:fer.vdl1928@gmail.com" target="_blank">fer.vdl1928@gmail.com</a>> wrote:<br>
>>> <br>
>>> Dear all Chimera users<br>
>>> <br>
>>> In UCSF ChimeraX version 0.91 (2019-08-30)<br>
>>> 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)<br>
>>> <br>
>>> Then I input the commands:<br>
>>> <br>
>>> >surface<br>
>>> >color electrostatic #1 map #2 palette -10,red:0,white:10,blue<br>
>>> <br>
>>> <image.png><br>
>>> Is it possible to calculate the red and blue area in Å2 of the molecule?<br>
>>> <br>
>>> Best regards,<br>
>>> <br>
>>> Fernando<br>
>>> <br>
>>> <br>
>>> ATTE<br>
>>> Fernando Villa Díaz<br>
>>> <br>
<br>
</blockquote></div>