[phenixbb] Creating files with positive and negative map peaks

Jennifer Ekstrom ekstrom at msu.edu
Wed Jul 23 14:26:47 PDT 2008


Hi Pavel,
I'm actually looking for a way to get that peak data into the the format 
given by cns.  I've pasted an example below... 

Thanks,
Jennifer

......................

REMARK DATE:16-Jan-2008  04:53:47       created by user: ekstrom
REMARK VERSION:1.1
ATOM   6706 PEAK PEAK    1       2.936 -26.732   8.585  1.00 12.12      PEAK
ATOM   6707 PEAK PEAK    2     -29.892 -17.631  49.107  1.00 11.30      PEAK
ATOM   6708 PEAK PEAK    3     -22.298 -14.560  50.705  1.00  9.38      PEAK
ATOM   6709 PEAK PEAK    4      -0.284 -37.687  22.843  1.00  8.50      PEAK
ATOM   6710 PEAK PEAK    5     -26.184 -10.173  38.207  1.00  8.43      PEAK
ATOM   6711 PEAK PEAK    6     -27.455 -11.757  37.651  1.00  7.88      PEAK
ATOM   6712 PEAK PEAK    7       0.251 -36.207  17.111  1.00  7.54      PEAK
ATOM   6713 PEAK PEAK    8     -36.967   7.877  42.542  1.00  7.44      PEAK
ATOM   6714 PEAK PEAK    9       1.491 -38.169  16.878  1.00  7.44      PEAK
ATOM   6715 PEAK PEAK   10     -27.501 -10.922  37.959  1.00  7.35      PEAK
ATOM   6716 PEAK PEAK   11       4.477 -45.133  14.678  1.00  7.21      PEAK
ATOM   6717 PEAK PEAK   12      -0.963 -25.803  18.167  1.00  7.04      PEAK
ATOM   6718 PEAK PEAK   13     -26.482 -13.615  36.693  1.00  6.93      PEAK
ATOM   6719 PEAK PEAK   14     -26.015 -10.239  30.585  1.00  6.86      PEAK
ATOM   6720 PEAK PEAK   15     -30.002 -26.868  20.550  1.00  6.69      PEAK
ATOM   6721 PEAK PEAK   16     -27.396 -11.906  38.698  1.00  6.65      PEAK
ATOM   6722 PEAK PEAK   17     -31.373  -6.638  28.568  1.00  6.64      PEAK
ATOM   6723 PEAK PEAK   18      -8.221 -55.948   9.682  1.00  6.53      PEAK
ATOM   6724 PEAK PEAK   19     -28.291 -23.724  37.848  1.00  6.28      PEAK
ATOM   6725 PEAK PEAK   20     -30.468 -30.466  41.446  1.00  6.18      PEAK
ATOM   6726 PEAK PEAK   21     -31.063 -29.443  40.172  1.00  6.12      PEAK
ATOM   6727 PEAK PEAK   22     -26.755 -11.332  32.107  1.00  6.07      PEAK
ATOM   6728 PEAK PEAK   23       1.773 -36.238  17.792  1.00  6.05      PEAK
ATOM   6729 PEAK PEAK   24     -29.413 -25.920  38.032  1.00  6.04      PEAK
ATOM   6730 PEAK PEAK   25     -28.233 -22.681  37.519  1.00  6.03      PEAK
ATOM   6731 PEAK PEAK   26     -16.125 -21.401  26.461  1.00  5.99      PEAK
ATOM   6732 PEAK PEAK   27       0.732 -34.913  18.222  1.00  5.98      PEAK
ATOM   6733 PEAK PEAK   28     -27.046  -9.637  32.604  1.00  5.90      PEAK
ATOM   6734 PEAK PEAK   29     -29.019 -24.193  39.115  1.00  5.87      PEAK
ATOM   6735 PEAK PEAK   30       2.977 -35.032  16.644  1.00  5.76      PEAK
END


Pavel Afonine wrote:
> Hi Jennifer,
>
> phenix.refine reports the list of highest peaks and deepest holes 
> found at difference map. Look towards the end of .log file for lines 
> like this:
>
> ========== residual map mFobs-DFmodel: highest peaks and deepst holes 
> =========
>
>
>                            
> ----------peaks----------                          
>
> Number of peaks found at mFobs-DFmodel map (map cutoff=3.50 sigma)= 
> 21         
> Filter by distance & map next to the model:
>    mapped sites are within: 1.002 - 3.227
>    number of sites selected in [dist_min= 0.70, dist_max= 6.00]: 20 
> from: 21
>    mapped sites are within: 1.002 - 3.227
>
> peak=    4.154 closest distance to pdb=" N   VAL     2 " =    2.881
> peak=    3.572 closest distance to pdb=" O   VAL     2 " =    2.663
> peak=    3.895 closest distance to pdb=" CB  ARG     5 " =    2.743
> peak=    4.186 closest distance to pdb=" N   GLU     7 " =    2.237
> peak=    4.958 closest distance to pdb=" OE1 GLU     7 " =    2.920
>
> Same for negative peaks.
>
> Please let us know if you have any questions!
>
> Pavel.
>  
>
>
> On 7/23/2008 1:22 PM, Jennifer Ekstrom wrote:
>> I would like to obtain a file with a listing of positive and negative 
>> difference map peaks, similar in formatting and content to the files 
>> that would come out of cns: 
>>
>>          positive peaks in: <output_root>_positive.peaks
>>          negative peaks in: <output_root>_negative.peaks
>>
>>
>> We need this for use in DDQ at the jcsg validation server 
>> (http://www.jcsg.org/prod/scripts/validation/sv2.cgi).  Is there an easy 
>> way to create these peak files in phenix.refine?  We'd prefer to not 
>> re-refine in CNS as we have *many* alternate conformations... 
>>
>> Thanks, Jennifer
>> _______________________________________________
>> phenixbb mailing list
>> phenixbb at phenix-online.org
>> http://www.phenix-online.org/mailman/listinfo/phenixbb
>>   




More information about the phenixbb mailing list