<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<span style="color: rgb(0, 0, 0); font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;">Elaine: I think the approach with the fake atoms and measure mapvalues might be what I need. </span><br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
Tom: I had anticipated I would probably need to use a Python script to obtain the values however I was unable to find the ChimeraX command that would return the value. </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
I might have enough information to give this a test. </div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255);" class="elementToProof">
Thank you!</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Elaine Meng <meng@cgl.ucsf.edu><br>
<b>Sent:</b> Thursday, February 2, 2023 1:30 PM<br>
<b>To:</b> Lohry, David <dlohry@vols.utk.edu><br>
<b>Cc:</b> chimerax-users@cgl.ucsf.edu <chimerax-users@cgl.ucsf.edu><br>
<b>Subject:</b> Re: [chimerax-users] Obtain electrostatic potential from PQR/DX via script</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:11pt;">
<div class="PlainText">[You don't often get email from meng@cgl.ucsf.edu. Learn why this is important at
<a href="https://aka.ms/LearnAboutSenderIdentification">https://aka.ms/LearnAboutSenderIdentification</a> ]<br>
<br>
Hi David,<br>
There is a "measure mapvalues" command that will get values of a map (e.g. density or electrostatic potential) at atomic center positions and assign them as an atom attribute. Then the resulting attribute can be written to a file with the "save" command.<br>
<br>
<<a href="https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Fcommands%2Fmeasure.html%23mapvalues&data=05%7C01%7Cdlohry%40vols.utk.edu%7Cd209a42f93ca469bea7608db054b92d6%7C515813d9717d45dd9eca9aa19c09d6f9%7C0%7C0%7C638109594368478433%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=jhgNkhPlNnDEfmKqaeVVEkFXKyvQt3h54SY03P38J5c%3D&reserved=0">https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Fcommands%2Fmeasure.html%23mapvalues&data=05%7C01%7Cdlohry%40vols.utk.edu%7Cd209a42f93ca469bea7608db054b92d6%7C515813d9717d45dd9eca9aa19c09d6f9%7C0%7C0%7C638109594368478433%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=jhgNkhPlNnDEfmKqaeVVEkFXKyvQt3h54SY03P38J5c%3D&reserved=0</a>><br>
<<a href="https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Fcommands%2Fsave.html%23attributes&data=05%7C01%7Cdlohry%40vols.utk.edu%7Cd209a42f93ca469bea7608db054b92d6%7C515813d9717d45dd9eca9aa19c09d6f9%7C0%7C0%7C638109594368478433%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=imoOH4eeaQJ3kc%2Fh%2F4AM1XDfzXa83tnq1ICFxnGnvww%3D&reserved=0">https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Fcommands%2Fsave.html%23attributes&data=05%7C01%7Cdlohry%40vols.utk.edu%7Cd209a42f93ca469bea7608db054b92d6%7C515813d9717d45dd9eca9aa19c09d6f9%7C0%7C0%7C638109594368478433%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=imoOH4eeaQJ3kc%2Fh%2F4AM1XDfzXa83tnq1ICFxnGnvww%3D&reserved=0</a>><br>
<br>
However, that is different than what you are viewing, the map values at surface vertices (or probably even 1.4 A outward from surface vertices since the default electrostatic coloring uses that offset to approximate the ESP at the solvent-accessible surface
even though the solvent-excluded surface is what is actually displayed). These different surface definitions are shown in a diagram here:<br>
<<a href="https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Fcommands%2Fsurface.html&data=05%7C01%7Cdlohry%40vols.utk.edu%7Cd209a42f93ca469bea7608db054b92d6%7C515813d9717d45dd9eca9aa19c09d6f9%7C0%7C0%7C638109594368478433%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=wWqwUxxmQiQ046gc0TBMjKd%2FwElgWAtQ65r9tc1JeBw%3D&reserved=0">https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Fcommands%2Fsurface.html&data=05%7C01%7Cdlohry%40vols.utk.edu%7Cd209a42f93ca469bea7608db054b92d6%7C515813d9717d45dd9eca9aa19c09d6f9%7C0%7C0%7C638109594368478433%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=wWqwUxxmQiQ046gc0TBMjKd%2FwElgWAtQ65r9tc1JeBw%3D&reserved=0</a>><br>
<br>
These surface values are probably more relevant to understanding molecular interactions than the potential at the atomic centers, which would be singularities if those were the same atoms with charges used to calculate the potential. So an additional step
is to make the surface vertices into "fake atoms" (markers) before using measure mapvalues etc. as mentioned at the beginning of this message.<br>
<br>
The command to make markers out of a surface is "marker fromMesh":<br>
<<a href="https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Fcommands%2Fmarker.html%23create&data=05%7C01%7Cdlohry%40vols.utk.edu%7Cd209a42f93ca469bea7608db054b92d6%7C515813d9717d45dd9eca9aa19c09d6f9%7C0%7C0%7C638109594368478433%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=kofUdQwsNq43TvB%2F0z58ZxQibAmeP9dKIV9Bzup2b3E%3D&reserved=0">https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Fcommands%2Fmarker.html%23create&data=05%7C01%7Cdlohry%40vols.utk.edu%7Cd209a42f93ca469bea7608db054b92d6%7C515813d9717d45dd9eca9aa19c09d6f9%7C0%7C0%7C638109594368478433%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=kofUdQwsNq43TvB%2F0z58ZxQibAmeP9dKIV9Bzup2b3E%3D&reserved=0</a>><br>
<br>
So I was able to do this on a small example, first using "coulombic" to calculate an ESP map since I didn't have a dx handy, but your procedure would be analogous, simply opening your dx instead of using "coulombic". In the following example, the molecular
surface is #1.1 and the fake atoms made from it are #2:<br>
<br>
open 1gcn<br>
coulombic protein map true<br>
surf<br>
marker fromMesh #1.1 edge 0.05<br>
measure mapvalues #1.2 atoms #2 attribute cesp<br>
save ~/Desktop/cesp.defattr attrName cesp models #2<br>
<br>
So now I have a plain text file cesp.defattr, format described here<br>
<<a href="https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Fformats%2Fdefattr.html&data=05%7C01%7Cdlohry%40vols.utk.edu%7Cd209a42f93ca469bea7608db054b92d6%7C515813d9717d45dd9eca9aa19c09d6f9%7C0%7C0%7C638109594368478433%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=fUhdJj%2FyT2XFO7ziNcYrqw6b3SfF822LqxrvYIHEGhw%3D&reserved=0">https://nam11.safelinks.protection.outlook.com/?url=https%3A%2F%2Frbvi.ucsf.edu%2Fchimerax%2Fdocs%2Fuser%2Fformats%2Fdefattr.html&data=05%7C01%7Cdlohry%40vols.utk.edu%7Cd209a42f93ca469bea7608db054b92d6%7C515813d9717d45dd9eca9aa19c09d6f9%7C0%7C0%7C638109594368478433%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=fUhdJj%2FyT2XFO7ziNcYrqw6b3SfF822LqxrvYIHEGhw%3D&reserved=0</a>><br>
<br>
...with contents:<br>
<br>
attribute: cesp<br>
recipient: atoms<br>
match mode: 1-to-1<br>
#2/M:1@M 1.6767187<br>
#2/M:2@M 1.9566636<br>
#2/M:3@M -0.15276606<br>
#2/M:4@M -0.13766149<br>
#2/M:5@M 3.8019018<br>
#2/M:6@M 3.8222048<br>
#2/M:7@M 1.3163166<br>
#2/M:8@M 1.6913735<br>
#2/M:9@M -12.664384<br>
#2/M:10@M -12.545011<br>
#2/M:11@M 0.17837371<br>
#2/M:12@M 0.13978194<br>
#2/M:13@M -13.340148<br>
#2/M:14@M -13.404252<br>
[.... etc ...]<br>
<br>
The main issue I can see is that you don't know where on the surface each fake atom is unless you do something to label or zoom in on it in the 3D display, e.g.<br>
<br>
size /M atomRadius 0.1<br>
view /M:1 :<5<br>
label /M:1 color yellow<br>
select /M:1<br>
<br>
...to see where on the surface that particular marker (fake atom) resides. Another issue is that the file will be quite large since there are lots of surface vertices.<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 Feb 1, 2023, at 7:10 PM, Lohry, David via ChimeraX-users <chimerax-users@cgl.ucsf.edu> wrote:<br>
><br>
> Hi everyone,<br>
><br>
> I'm performing an analysis where I want to select different regions on a molecule and obtain a list of all of the electrostatic potentials from that region. My input is a PQR file with a DX file. These were generated from a PDB and APBS.<br>
><br>
> The electrostatic potential can be viewed when the option "Report value at mouse position" is selected in the Surface Color tool (Tools...Volume Data...Surface Color...Options button) In the status bar at the bottom, it will show the value along with the
coordinates.<br>
><br>
> Is there a way to obtain this information via a ChimeraX command that I could script?<br>
><br>
> Thank you!<br>
<br>
</div>
</span></font></div>
</body>
</html>