<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
<title></title>
</head>
<body bgcolor="#ffffff" text="#000000">
Hi Joe,<br>
<br>
yes, I agree on this. It's on my list already and hopefully I do it
soon.<br>
<br>
Thanks for the feedback!<br>
<br>
Pavel.<br>
<br>
<br>
On 9/30/09 4:41 PM, Joe Krahn wrote:
<blockquote cite="mid:4AC3ECC1.2010805@niehs.nih.gov" type="cite">
<pre wrap="">A few more comments on writing peak coordinates, after going through
peaks based on the current PHENIX log file:
Having actual coordinates present as atoms makes it easy to measure
distances to other atoms. Less importantly, if there are two peaks close
to each other, it is not always obvious which peak matches the log.
It may be sufficient to write the peak coordinates in the log file. That
would avoid the need for adding another output file. It would not be
that hard to parse into either a PDB file, or some sort of
program-specific list or script.
Thanks,
Joe Krahn
Krahn, Joe (NIH/NIEHS) [C] wrote:
</pre>
<blockquote type="cite">
<pre wrap="">Pavel Afonine wrote:
</pre>
<blockquote type="cite">
<pre wrap="">Hi Joe,
</pre>
<blockquote type="cite">
<pre wrap="">PHENIX has an option to find positive and negative peaks, but only lists
them in the log file along with the residue they are closest to.
</pre>
</blockquote>
<pre wrap="">I have done it this way having the following in mind:
- you load you refined model and maps in your favorite graphics program,
Coot for example;
- you take that +/- peaks list;
- and you go residue-by-residue and check the flagged peaks.
</pre>
<blockquote type="cite">
<pre wrap="">Is it
possible to write out the actual peak coordinates, such as to a PDB file?
</pre>
</blockquote>
<pre wrap="">It's possible of course, but I have a couple of questions:
- why to do this?
- what do I put as atom name, residue name, occupancy, B-factor, atom
chemical type ?
Thanks again for your feedback!
Pavel.
</pre>
</blockquote>
<pre wrap="">Most of the time, I optimize a model by stepping through peaks, rather
than walking through all residues. In most software, it is easy to step
through residues in a PDB file, whereas converting the PHENIX log output
into something that can automate this is likely program specific.
My standard approach is to get a peak list in the form of a PDB file,
sorted from the highest peak, and work through those peaks until they
become too weak. I can fix model errors, or add a water at the peak
position, or decide that it is something other than water. Having an
atom at the peak makes it easy to add waters, which is good for
lower-resolution or noisy maps, where I don't trust automated waters.
Many programs write peaks as PDB files, usually with B-factor set from
the difference-map sigma level. Some programs write them out as waters,
which is probably safest to avoid problems with non-standard names. The
occupancy is usually 1.0, but it could be set to zero to flag that they
are not actual model coordinates.
This approach used to be very common. Some new programs have built-in
peak search utilities, so it may be getting less common for that reason.
Thanks,
Joe Krahn
_______________________________________________
phenixbb mailing list
<a class="moz-txt-link-abbreviated" href="mailto:phenixbb@phenix-online.org">phenixbb@phenix-online.org</a>
<a class="moz-txt-link-freetext" href="http://www.phenix-online.org/mailman/listinfo/phenixbb">http://www.phenix-online.org/mailman/listinfo/phenixbb</a>
</pre>
</blockquote>
<pre wrap=""><!---->
_______________________________________________
phenixbb mailing list
<a class="moz-txt-link-abbreviated" href="mailto:phenixbb@phenix-online.org">phenixbb@phenix-online.org</a>
<a class="moz-txt-link-freetext" href="http://www.phenix-online.org/mailman/listinfo/phenixbb">http://www.phenix-online.org/mailman/listinfo/phenixbb</a>
</pre>
</blockquote>
</body>
</html>