Re: [ccp4bb] Merohedral twining for P212121. -

2010-06-29 Thread Loes Kroon-Batenburg

Dear Frank,

DIRAX is very good at finding the twin lattices in case of non-merohedral 
twinning.  In the reticular case you might want to use the LEPAGE-TWIN routine 
in PLATON to find the correct subcell.
My collegue, Martin Lutz, suggests to change the measurement temperature in 
order to change the cell parameters and thus the amount of overlap.
For integration of the data you may want to use our EVAL15 (J. Appl. Cryst. 43 
(2010) 70-82).

Have a look at  http://www.crystal.chem.uu.nl/distr/
EVAL15 can deconvolute reflections with significant overlap, or otherwise 
outputs the summed intensties of the two reflections. This method is 
particularly suitable for your reticular case. Scaling can be done with TWINABS 
and refinement can be done with SHELX (hklf5 format).


Loes.

P.S. a collection of twinning literature can be found at 
http://www.cryst.chem.uu.nl/lutz/twin/gen_twin.html


On 06/23/10 11:46, Frank von Delft wrote:

My experience with pseudo-merohedral twinning (it was actually the
reticular case with half the spots overlapped and the other
non-overlapped half on a pseudo C-centred lattice) is that the degree
of splitting varies widely over the diffraction pattern. In some
places there was complete overlap, in others you see elongation of the
spots, in others partial separation, and in others complete separation
(and of course all shades in-between), with around 50-50 intensity
split. In this situation the mosaicity becomes meaningless! I'm not
aware of any software that can handle this kind of thing successfully
(and certainly the data we did manage to get turned out to be
garbage!).


Both DIRAX or SAINT should be able to handle it, you'll need SADABS to
scale it. (The latter two are in the Bruker software.)

phx.



--
__

Dr. Loes Kroon-Batenburg
Dept. of Crystal and Structural Chemistry
Bijvoet Center for Biomolecular Research
Utrecht University
Padualaan 8, 3584 CH Utrecht
The Netherlands

E-mail : l.m.j.kroon-batenb...@uu.nl
phone  : +31-30-2532865
fax: +31-30-2533940
__


Re: [ccp4bb] Lodovico...

2010-06-29 Thread Vellieux Frederic

Dear all,

A tribute to Lodovico (pictures + letters + a song) can be found on 
http://www.crystalerice.org/Erice2010/Lodovico/index.htm


Fred.


Re: [ccp4bb] f2mtz error: problems reading reflection

2010-06-29 Thread Eleanor Dodson

Can you attach the first 50 lines of your hkl file?
Eleanor

Yogesh Gupta wrote:

Dear Experts,

After a new installation on Mac OS 10.6, i am getting this error (related to
f2mtz) during the Find SITES step by SHELX in Autosharp.

Data line--- LABO H K L FA SIGFA ALPHA
 Number of columns to be read in:   6

 Data line--- CTYP H H H F  Q P
 Data line--- FORM '(3F4.0,2F8.2,F4.0)'

 $TEXT:Warning: $$ comment $$
 WARNING:  Format I replaced with F
 $$

FORMATTED  OLD file opened on unit   1

Logical name: HKLIN, Filename: SHELX/1_fa.hkl


 $TEXT:Warning: $$ comment $$
 WARNING:  PROJECTNAME not assigned
 $$

 $TEXT:Warning: $$ comment $$
 WARNING:  CRYSTALNAME not assigned
 $$


 $TEXT:Warning: $$ comment $$
 WARNING:  DATASETNAME not assigned
 $$
*** Read error
Check FORMAT -- especially the need for REALs
 f2mtz:  problems reading reflection   9

 f2mtz:  problems reading reflection   9


I am running ccp4-6.1.13. It looks like a FORTRAN error but do not know how
to fix this or which file to edit particularly. SHARP otherwise running fine
if i exclude SHELX step by feeding known sites.

Thanks for your help.

Yogesh



Re: [ccp4bb] Fast,medium,slow axis do not match

2010-06-29 Thread Kevin Cowtan
You could feed *both* maps through mapmask with AXIS X Y Z to convert 
them to the same axis order.


You may also have a problem with the maps having different XYZLIM 
ranges. In that case, using XYZLIM MATCH on the second run of mapmask to 
match the second map to the first should fix it.


Kevin

Hailiang Zhang wrote:

Hi there,

I was generating the atomic mask using ccpr-sfall, and generating real
maps using ccp4-fft, and then ccp4-overlap these maps to calculate the cc
values. However, ccp4-overlapmap frequently complaints that the
Fast,medium,slow axis do not match for these maps. Following are the
script I was using, and thanks for any advice (I can't manually change the
keywords values one by one since everything needs automated):


sfall xyzin ${OUTTMPDIR}${PDB}-Bsmall-0B.pdb \
  mapout ${OUTTMPDIR}${PDB}-Bsmall-mask.map\
 END-sfall
TITL Toxd Atom map from final coordinates
MODE ATMMAP RESMOD
GRID ${GRID}
SYMM ${spcgrp}
FORM NGAU 5
END
END-sfall

...



fft ${OUTLOGDIR}${PDB}-0B-fft.log\
hklin ${OUTTMPDIR}${PDB}-0B.mtz \
mapout ${OUTTMPDIR}${PDB}-0B.map\
 END-fft
TITL data_toxd
GRID ${GRID}
XYZLIM  ASU
RESO 100.0 ${resolution}
LABI F1=FOBS PHI=PHIFMODEL
END
END-fft
*



--
EMAIL DISCLAIMER http://www.york.ac.uk/docs/disclaimer/email.htm


[ccp4bb] measure of anamolous signal

2010-06-29 Thread Vandu Murugan
Dear all,
   I have collected a 2.7 angstrom home source data with Cu-Kalpha source
for a protein with 6 cysteines, with a multiplicity of around 23.  I need to
know, is there any significant anamolous signal present in the data set,
since there is no good model for my protein.  Can any one tell, which
program to run, and what parameter to see?  Thanks in advance.

cheers,
Murugan


Re: [ccp4bb] measure of anamolous signal

2010-06-29 Thread Vellieux Frederic

Vandu Murugan wrote:

Dear all,
   I have collected a 2.7 angstrom home source data with Cu-Kalpha 
source for a protein with 6 cysteines, with a multiplicity of around 
23.  I need to know, is there any significant anamolous signal present 
in the data set, since there is no good model for my protein.  Can any 
one tell, which program to run, and what parameter to see?  Thanks in 
advance.


cheers,
Murugan
When processing the data, ensure that Bijvoet pairs are not merged. The 
data processing software should provide you with an R-ano value and that 
is already a start. The values provided should tell you if you have an 
anomalous signal or not. You may also have to play with the number of 
frames to integrate in order to obtain an optimal anomalous signal in 
the resulting data set.


There are several publications describing Sulphur SAD, but one of them 
which is freely available on the Web can be found here:
http://www.stoe.com/pages/brochure/labnote_genix_cu.pdf , Schiltz M (pp 
4-6).


Good luck!

Fred.


Re: [ccp4bb] measure of anamolous signal

2010-06-29 Thread Tim Gruene
Dear Murugan,

you can use the program hkl2map from Thomas Schneider, available at 
http://webapps.embl-hamburg.de/hkl2map/
It's a graphical interface to the programs shelx c/d/e which are available from
http://shelx.uni-ac.gwdg.de/SHELX/index.html

With SAD data you want to look at the d/sigma line at the end of the shelxc
output. Where that drops below about 1.3 is approximately where your anomalous
signal ends. You might get slightly improved statistics with xprep instead of
shelxc, but xprep is not free and you have to get a copy from Bruker-AXS.

Tim

On Tue, Jun 29, 2010 at 02:35:29PM +0530, Vandu Murugan wrote:
 Dear all,
I have collected a 2.7 angstrom home source data with Cu-Kalpha source
 for a protein with 6 cysteines, with a multiplicity of around 23.  I need to
 know, is there any significant anamolous signal present in the data set,
 since there is no good model for my protein.  Can any one tell, which
 program to run, and what parameter to see?  Thanks in advance.
 
 cheers,
 Murugan

-- 
--
Tim Gruene
Institut fuer anorganische Chemie
Tammannstr. 4
D-37077 Goettingen

GPG Key ID = A46BEE1A



signature.asc
Description: Digital signature


Re: [ccp4bb] measure of anamolous signal

2010-06-29 Thread Graeme Winter
Hi Murugan,

One useful indicator of raw anomalous signal is the ANOMPLOT graph
from Scala - this shows the differences between reflections compared
with the expected differences. If the gradient of the plot is 1
there's no more differences that you would expect. If the gradient is
more than one there is (or may be.) - also check out the merging
statistics as a function of batch, if there's significant radiation
damage this may mess things up.

Scala writes out the gradient (assuming you told it anomalous on) in the summary

Another rule-of-thumb is the resolution limit where cc-anom is  0.5.

The most practical indicator of the anomalous signal is of course the
success or failure of the subsequent phasing :o)

Best wishes,

Graeme

On 29 June 2010 10:05, Vandu Murugan wandumuru...@gmail.com wrote:
 Dear all,
    I have collected a 2.7 angstrom home source data with Cu-Kalpha source
 for a protein with 6 cysteines, with a multiplicity of around 23.  I need to
 know, is there any significant anamolous signal present in the data set,
 since there is no good model for my protein.  Can any one tell, which
 program to run, and what parameter to see?  Thanks in advance.

 cheers,
 Murugan



Re: [ccp4bb] measure of anamolous signal

2010-06-29 Thread Dirk Kostrewa

 Dear Murugan,

Am 29.06.10 11:05, schrieb Vandu Murugan:

Dear all,
   I have collected a 2.7 angstrom home source data with Cu-Kalpha 
source for a protein with 6 cysteines, with a multiplicity of around 
23.  I need to know, is there any significant anamolous signal present 
in the data set, since there is no good model for my protein.  Can any 
one tell, which program to run, and what parameter to see?  Thanks in 
advance.


cheers,
Murugan
estimating the quality of the anomalous signal is not trivial, and 
several quality indicators have been discussed (see for example Fu, Rose 
 Wang, Acta Cryst D60, 499-506 (2004), or Zwart, Acta Cryst D61, 
1437-1444 (2005)).


If you process your data with XDS, there are two quality indicators for 
the anomalous signal, given both in the CORRECT.LP and in the XSCALE.LP 
file. One is the correlation of anomalous differences between two 
randomly chosen subsets that should have the same anomalous difference 
due to crystallographic symmetry, called RANOM. The other describes the 
absolute anomalous difference divided by their standard deviation, 
called SIGANO. A typical rule of thumb (and the one that I use) is, that 
RANOM should be ~ 30%, and SIGANO should be ~ 1.2. However, these 
indicators might not be realistic (see references above) and therefore 
should be taken with a grain of salt.


Good luck,

Dirk.

--

***
Dirk Kostrewa
Gene Center Munich, A5.07
Department of Biochemistry
Ludwig-Maximilians-Universität München
Feodor-Lynen-Str. 25
D-81377 Munich
Germany
Phone:  +49-89-2180-76845
Fax:+49-89-2180-76999
E-mail: kostr...@genzentrum.lmu.de
WWW:www.genzentrum.lmu.de
***


[ccp4bb] P-CUBE - Last chance to sign up for the Workshop and 1st Annual Meeting in Grenoble- Registration is closing soon!

2010-06-29 Thread Peer Mittl

Dear All,

Don't miss out on the P-CUBE meeting in Grenoble in September! Sign up today 
underwww.p-cube.eu.  The registration deadline is June 30. 2010.

See you in Grenoble
P-CUBE Management Team


***
Dr. Jutta Tatzel
Program Manager P-CUBE

Department of Biochemistry
University of Zurich
Winterthurerstr. 190
CH-8057 Zurich

Tel +41 44 635 5593
Email:j.tat...@bioc.uzh.ch


Re: [ccp4bb] measure of anamolous signal

2010-06-29 Thread Frank von Delft
I've found the scala CC-anom significantly underestimates the anomalous 
signal, relative to e.g. xprep.  I don't know why that is, but the 
latter seems to agree with what shelxd is happy with.


Cheers
phx




On 29/06/2010 10:35, Graeme Winter wrote:

Hi Murugan,

One useful indicator of raw anomalous signal is the ANOMPLOT graph
from Scala - this shows the differences between reflections compared
with the expected differences. If the gradient of the plot is 1
there's no more differences that you would expect. If the gradient is
more than one there is (or may be.) - also check out the merging
statistics as a function of batch, if there's significant radiation
damage this may mess things up.

Scala writes out the gradient (assuming you told it anomalous on) in the summary

Another rule-of-thumb is the resolution limit where cc-anom is  0.5.

The most practical indicator of the anomalous signal is of course the
success or failure of the subsequent phasing :o)

Best wishes,

Graeme

On 29 June 2010 10:05, Vandu Muruganwandumuru...@gmail.com  wrote:
   

Dear all,
I have collected a 2.7 angstrom home source data with Cu-Kalpha source
for a protein with 6 cysteines, with a multiplicity of around 23.  I need to
know, is there any significant anamolous signal present in the data set,
since there is no good model for my protein.  Can any one tell, which
program to run, and what parameter to see?  Thanks in advance.

cheers,
Murugan

 


Re: [ccp4bb] measure of anamolous signal

2010-06-29 Thread Phil Evans
I would say CC-anom  0.3 or even 0.2 (note that the scala CC-anom is defined 
on I not F)

Phil

On 29 Jun 2010, at 14:37, Frank von Delft wrote:

 I've found the scala CC-anom significantly underestimates the anomalous 
 signal, relative to e.g. xprep.  I don't know why that is, but the latter 
 seems to agree with what shelxd is happy with.
 
 Cheers
 phx
 
 
 
 
 On 29/06/2010 10:35, Graeme Winter wrote:
 Hi Murugan,
 
 One useful indicator of raw anomalous signal is the ANOMPLOT graph
 from Scala - this shows the differences between reflections compared
 with the expected differences. If the gradient of the plot is 1
 there's no more differences that you would expect. If the gradient is
 more than one there is (or may be.) - also check out the merging
 statistics as a function of batch, if there's significant radiation
 damage this may mess things up.
 
 Scala writes out the gradient (assuming you told it anomalous on) in the 
 summary
 
 Another rule-of-thumb is the resolution limit where cc-anom is  0.5.
 
 The most practical indicator of the anomalous signal is of course the
 success or failure of the subsequent phasing :o)
 
 Best wishes,
 
 Graeme
 
 On 29 June 2010 10:05, Vandu Muruganwandumuru...@gmail.com  wrote:
   
 Dear all,
I have collected a 2.7 angstrom home source data with Cu-Kalpha source
 for a protein with 6 cysteines, with a multiplicity of around 23.  I need to
 know, is there any significant anamolous signal present in the data set,
 since there is no good model for my protein.  Can any one tell, which
 program to run, and what parameter to see?  Thanks in advance.
 
 cheers,
 Murugan
 
 


Re: [ccp4bb] Fast,medium,slow axis do not match

2010-06-29 Thread Eleanor Dodson
Alternatively you can force the fft to generate a map with the sfall 
axiis order (sfall has a fixed axis order governed by the spacegroup)

From the fft documentation
You can set

AXIS fast medium slow

Eleanor

Kevin Cowtan wrote:
You could feed *both* maps through mapmask with AXIS X Y Z to convert 
them to the same axis order.


You may also have a problem with the maps having different XYZLIM 
ranges. In that case, using XYZLIM MATCH on the second run of mapmask to 
match the second map to the first should fix it.


Kevin

Hailiang Zhang wrote:

Hi there,

I was generating the atomic mask using ccpr-sfall, and generating real
maps using ccp4-fft, and then ccp4-overlap these maps to calculate the cc
values. However, ccp4-overlapmap frequently complaints that the
Fast,medium,slow axis do not match for these maps. Following are the
script I was using, and thanks for any advice (I can't manually change 
the

keywords values one by one since everything needs automated):


sfall xyzin ${OUTTMPDIR}${PDB}-Bsmall-0B.pdb \
  mapout ${OUTTMPDIR}${PDB}-Bsmall-mask.map\
 END-sfall
TITL Toxd Atom map from final coordinates
MODE ATMMAP RESMOD
GRID ${GRID}
SYMM ${spcgrp}
FORM NGAU 5
END
END-sfall

...



fft ${OUTLOGDIR}${PDB}-0B-fft.log\
hklin ${OUTTMPDIR}${PDB}-0B.mtz \
mapout ${OUTTMPDIR}${PDB}-0B.map\
 END-fft
TITL data_toxd
GRID ${GRID}
XYZLIM  ASU
RESO 100.0 ${resolution}
LABI F1=FOBS PHI=PHIFMODEL
END
END-fft
*





Re: [ccp4bb] measure of anamolous signal

2010-06-29 Thread James Holton
Answering the question should I even bother trying? can be 
complicated, but I get asked that a lot at the beamline!  I recently 
incorporated a number of data quality formulas into a little interactive 
web page here:

http://bl831.als.lbl.gov/xtalsize.html

which focuses on calculating how many crystals of a given size you will 
need to average together to attain a specified goal (either weak spots 
(high res) or subtle differences (anomalous)) given radiation damage 
limits. 


For those interested in explanation:

A quick-and-dirty way to estimate your Bijvoet ratio is with the formula:
0.75*fpp*sqrt(sites/MW)
where sites is the number of sites you expect per MW, where MW is 
in Daltons (amu).  Note that MW can represent the protein monomer, 
asymmetric unit, or whatever, as long as sites refers to the same 
thing. fpp is the expected number of anomalous electrons, which you 
can get from the CCP4 program crossec or several websites.


Now, since you didn't mention your molecular weight, I will guess that 
it is about 30 kDa, which means your Bijvoet ratio at CuKa (where the 
fpp of sulfur is ~0.56 electron) will be about 0.6%.


Now, detecting a 0.6% difference requires that the two things you are 
subtracting (I+ and I-) be measured to at least a precision of 0.42% 
(because the error in the difference will then be: sqrt(0.42%^2+0.42%^2) 
= 0.6%).  This is the BARE minimum (where the signal is equal to the 
noise), but even to reach this goal, your signal-to-noise ratio must be 
1/0.42% = 240.  This is challenging!  Typical data sets have I/sig(I) ~ 
30 in their best bins (see Deiderichs 2010 
http://dx.doi.org/10.1107/S0907444910014836).  A multiplicity of 23 can 
push I/sig(I) to 30*sqrt(23) = 143, which is good, but still not close 
to 240.  You will probably need a multiplicity of ~65 to measure Bijvoet 
differences to an accuracy of 0.6%


Then again, if your protein is only 10 kDa, then your Bijvoet ratio is 
~1% and you will only need I/sig(I)  140.


-James Holton
MAD Scientist

Vandu Murugan wrote:

Dear all,
   I have collected a 2.7 angstrom home source data with Cu-Kalpha 
source for a protein with 6 cysteines, with a multiplicity of around 
23.  I need to know, is there any significant anamolous signal present 
in the data set, since there is no good model for my protein.  Can any 
one tell, which program to run, and what parameter to see?  Thanks in 
advance.


cheers,
Murugan


Re: [ccp4bb] measure of anamolous signal

2010-06-29 Thread Tim Gruene
On Tue, Jun 29, 2010 at 09:30:30AM -0700, James Holton wrote:
 Answering the question should I even bother trying? can be  
 complicated, [...]

 -James Holton
 MAD Scientist


Since 'trying' may only take a semi-experienced user about 30min until they
could have a poly-Alanine trace of their protein, I would say the answer is
definitely Yes irrespective of what statistics may tell you - unless you don't
have a protein in your crystal, but e.g.  DNA/ RNA...

Tim


 Vandu Murugan wrote:
 Dear all,
I have collected a 2.7 angstrom home source data with Cu-Kalpha  
 source for a protein with 6 cysteines, with a multiplicity of around  
 23.  I need to know, is there any significant anamolous signal present  
 in the data set, since there is no good model for my protein.  Can any  
 one tell, which program to run, and what parameter to see?  Thanks in  
 advance.

 cheers,
 Murugan

-- 
--
Tim Gruene
Institut fuer anorganische Chemie
Tammannstr. 4
D-37077 Goettingen

GPG Key ID = A46BEE1A



signature.asc
Description: Digital signature


[ccp4bb] Postdoc position

2010-06-29 Thread Howard S. Young
Applications are invited for a postdoctoral position to work on the  
structure determination of membrane transport proteins critical in  
cardiovascular function. The transporters we study have been linked to  
human disorders such as cardiomyopathy, ischemia-reperfusion injury,  
and various cancers. Our aim is to understand the molecular regulation  
of the transporters with the goal of rational design of novel  
inhibitors and activators. We seek an enthusiastic candidate  
experienced in macromolecular x-ray or electron crystallography.  
Several projects are currently available, including 2D and 3D crystals  
of membrane transport proteins and 3D crystals of soluble regulatory  
proteins. Experience with membrane protein biochemistry would be an  
asset. In addition, the successful candidate should be well-versed in  
one or more of the following areas: molecular biology, protein  
purification, crystallization, and membrane protein biochemistry.


This is an opportunity to join a lab experienced in the expression and  
crystallization of membrane proteins [See Biophysical Journal (2006)  
90:4213-23, JBC (2007) 282:9748-57, JBC (2008) 283:4145-54, and JBC  
(2009) 284:13513-8]. Our department provides interaction with numerous  
structural biologists. In addition, the University of Alberta Membrane  
Protein Research Group provides a dedicated forum for membrane protein  
research that includes 14 research laboratories  
(http://www.mprg.med.ualberta.ca/index.php ).


The University of Alberta, located in Edmonton, Alberta (Canada), is  
home to a large and interactive community of biomedical scientists  
http://www.med.ualberta.ca/; support and facilities for structural  
biology are excellent. A state-of-the-art facility is available for  
electron microscopy at the National Institute for Nanotechnology,  
University of Alberta  
(http://nint-innt.nrc-cnrc.gc.ca/research/micro/index_e.html), and a  
newly commissioned synchrotron facility is located a short drive away.  
Edmonton, having a population of approximately 1 million, offers a  
cosmopolitan environment with world class performing arts, sports,  
culinary and recreational opportunities.  In addition, the city’s  
proximity to the Rocky Mountains including the towns of Jasper and  
Banff is an additional bonus.


Salary is commensurate with training and experience at competitive  
rates. A comprehensive medical and dental benefit package is included  
with salary. Postdoctoral fellowship awards are also available through  
several competitive funding agencies, including the Canadian  
Institutes for Health Research and Alberta Innovates – Health  
Solutions. The University of Alberta hires on the basis of merit. We  
are committed to the principle of equity in employment. We welcome  
diversity and encourage applications from all qualified women and men,  
including persons with disabilities, members of visible minorities,  
and Aboriginal persons. The records arising from this competition will  
be managed in accordance with provisions of the Alberta Freedom of  
Information and Protection of Privacy Act (FOIPP).


Please direct CV's or inquiries to Dr. Howard S. Young at the  
following email address: hyo...@ualberta.ca.


Howard S. Young, PhD
Associate Professor  AHFMR Senior Scholar
Biochemistry, University of Alberta
Edmonton, Alberta T6G 2H7
(780)492-3931

National Institute for Nanotechnology,
National Research Council of Canada

Membrane Protein Research Group
(http://www.mprg.med.ualberta.ca/group-members.php)

CONFIDENTIALITY WARNING
This communication is intended for the use of the recipient to which  
it is addressed, and may contain confidential, personal, and/or  
privileged information. Please contact me immediately if you are not  
the intended recipient of this communication, and do not copy,  
distribute, or take action relying on it. Any communication received  
in error, or subsequent reply, should be deleted or destroyed.


Re: [ccp4bb] measure of anamolous signal

2010-06-29 Thread Bernhard Rupp
I second the hkl2map/SHELXCDE approach. Two complete examples 
explaining how to do this for MAD and S-SAD cases are in my book.
I wish to emphasize the importance of 
a) running enough trials 
b) careful selection of resolution cutoffs
c) look at the solution distribution
d) play with SHELXE parameters. 
The hkl2map graphs are enormously helpful for this purpose.

BR
-Original Message-
From: CCP4 bulletin board [mailto:ccp...@jiscmail.ac.uk] On Behalf Of Tim
Gruene
Sent: Tuesday, June 29, 2010 2:29 AM
To: CCP4BB@JISCMAIL.AC.UK
Subject: Re: [ccp4bb] measure of anamolous signal

Dear Murugan,

you can use the program hkl2map from Thomas Schneider, available at
http://webapps.embl-hamburg.de/hkl2map/
It's a graphical interface to the programs shelx c/d/e which are available
from http://shelx.uni-ac.gwdg.de/SHELX/index.html

With SAD data you want to look at the d/sigma line at the end of the shelxc
output. Where that drops below about 1.3 is approximately where your
anomalous signal ends. You might get slightly improved statistics with xprep
instead of shelxc, but xprep is not free and you have to get a copy from
Bruker-AXS.

Tim

On Tue, Jun 29, 2010 at 02:35:29PM +0530, Vandu Murugan wrote:
 Dear all,
I have collected a 2.7 angstrom home source data with Cu-Kalpha 
 source for a protein with 6 cysteines, with a multiplicity of around 
 23.  I need to know, is there any significant anamolous signal present 
 in the data set, since there is no good model for my protein.  Can any 
 one tell, which program to run, and what parameter to see?  Thanks in
advance.
 
 cheers,
 Murugan

--
--
Tim Gruene
Institut fuer anorganische Chemie
Tammannstr. 4
D-37077 Goettingen

GPG Key ID = A46BEE1A


Re: [ccp4bb] detection of anomalous signal using HKL2000

2010-06-29 Thread Phil Evans
Can anyone explain what Zbyszek Otwinowski means by Chi squared? I can't find a 
definition in any of his papers (though I may have missed it). Is there a 
reference?

It doesn't seem obviously related to the chi squared distribution (In 
probability theory and statistics, the chi-square distribution (also 
chi-squared orχ²-distribution) with k degrees of freedom is the distribution of 
a sum of the squares of k independent standard normal random variables.  
http://en.wikipedia.org/wiki/Chi-square_distribution)

Phil


On 29 Jun 2010, at 21:14, Felix Frolow wrote:

 Graphical information from Scalepack as used in HKL2000 is of unprecedented
 help to detect anomalous signal. Anomalous detection for S anomalous data 
 using
 CHI**2 and Rfactor statistics for reflections with averaged and separated 
 Bijvoet pairs is attached. It is very well described in HKL2000 manual. There 
 is nothing special about data collection
 (strategy was used) and measurement was relatively fast (4 h on MicroMax007 
 and RaxisIV++).
 BTW Rfactor is 1.9% with separated Bijvoet reflections and 2.6% with averaged
 
 
 BTW James Holton website calculate for this case 0.078 crystal
 
 Dr  Felix Frolow
 Professor of Structural Biology and Biotechnology
 Department of Molecular Microbiology
 and Biotechnology
 Tel Aviv University 69978, Israel
 
 Acta Crystallographica D, co-editor
 
 e-mail: mbfro...@post.tau.ac.il
 Tel:  ++972-3640-8723
 Fax: ++972-3640-9407
 Cellular: 0547 459 608
 
 
 
 
 
 
 
 
 Quoting Bernhard Rupp b...@ruppweb.org:
 
 I second the hkl2map/SHELXCDE approach. Two complete examples
 explaining how to do this for MAD and S-SAD cases are in my book.
 I wish to emphasize the importance of
 a) running enough trials
 b) careful selection of resolution cutoffs
 c) look at the solution distribution
 d) play with SHELXE parameters.
 The hkl2map graphs are enormously helpful for this purpose.
 
 BR
 -Original Message-
 From: CCP4 bulletin board [mailto:ccp...@jiscmail.ac.uk] On Behalf Of Tim
 Gruene
 Sent: Tuesday, June 29, 2010 2:29 AM
 To: CCP4BB@JISCMAIL.AC.UK
 Subject: Re: [ccp4bb] measure of anamolous signal
 
 Dear Murugan,
 
 you can use the program hkl2map from Thomas Schneider, available at
 http://webapps.embl-hamburg.de/hkl2map/
 It's a graphical interface to the programs shelx c/d/e which are available
 from http://shelx.uni-ac.gwdg.de/SHELX/index.html
 
 With SAD data you want to look at the d/sigma line at the end of the shelxc
 output. Where that drops below about 1.3 is approximately where your
 anomalous signal ends. You might get slightly improved statistics with xprep
 instead of shelxc, but xprep is not free and you have to get a copy from
 Bruker-AXS.
 
 Tim
 
 On Tue, Jun 29, 2010 at 02:35:29PM +0530, Vandu Murugan wrote:
 Dear all,
   I have collected a 2.7 angstrom home source data with Cu-Kalpha
 source for a protein with 6 cysteines, with a multiplicity of around
 23.  I need to know, is there any significant anamolous signal present
 in the data set, since there is no good model for my protein.  Can any
 one tell, which program to run, and what parameter to see?  Thanks in
 advance.
 
 cheers,
 Murugan
 
 --
 --
 Tim Gruene
 Institut fuer anorganische Chemie
 Tammannstr. 4
 D-37077 Goettingen
 
 GPG Key ID = A46BEE1A
 
 
 
 
 Dr  Felix Frolow
 Professor of Structural Biology and Biotechnology
 Department of Molecular Microbiology
 and Biotechnology
 Tel Aviv University 69978, Israel
 
 Acta Crystallographica D, co-editor
 
 e-mail: mbfro...@post.tau.ac.il
 Tel:  ++972-(0)-3640-8723
 Cel:  ++972-0547-459-608
 Fax:  ++972-(0)-3640-9407detectionOfAnomalousSignalBy Scalepack-HKL2000.pdf


Re: [ccp4bb] density quiz

2010-06-29 Thread Bernhard Rupp
Looks atypical for 3-me-Lys, M3L, (usually in natively methylated products) 
http://www.ruppweb.org/garland/gallery/Ch2/pages/Biomolecular_Crystallograph
y_Fig_2-27.htm

nor fits dimethyllysine, MYL, (via chem. red)

http://www.ruppweb.org/garland/gallery/Ch4/pages/Biomolecular_Crystallograph
y_Fig_4-19.htm

BR

PS: ...and don’t forget to remove those waters later ;-)

-Original Message-
From: CCP4 bulletin board [mailto:ccp...@jiscmail.ac.uk] On Behalf Of J.
Preben Morth
Sent: Tuesday, June 29, 2010 1:44 PM
To: CCP4BB@JISCMAIL.AC.UK
Subject: Re: [ccp4bb] density quiz

Hi Christian
One guess would be one of the possible post translational  
modifications seen for Lysine. acetylation or a triple methylation.  
However propionylations and butyrylations have also been observed on  
lysine.
It looks like you could build in waters, you may want to leave this  
lysine till later and see if more clear density show up.
good luck

Preben



On 29/06/2010, at 21.52, Christian Biertuempfel wrote:



 Dear all,
 Does anyone have an idea what this density at a lysine residue could  
 be?
  The structure is a protein-DNA complex and I suspect that something
 left-over from the DNA synthesis reacted with the amino group even
 though the DNA substrate was purified.

 christian

 p.s. The protein has seen the following substances during purification
 and crystallization: Leupeptin, Pepstatin, AEBSF, MES, Tris, HEPES,  
 DTT,
 EDTA, imidazole, glycerol, PEG 2K-MME, Mg2+, Na+, K+, Ca2+, Cl-.


 ___

 Dr. Christian Biertümpfel
 Laboratory of Molecular Biology

 NIDDK/National Institutes of Health  phone: +1 301 402  
 4647
 9000 Rockville Pike, Bldg. 5, Rm. B1-03  fax:   +1 301 496  
 0201
 Bethesda, MD 20892-0580
 USA
 ___
  
 density_quiz1 
 .png 
  
  
 density_quiz2 
 .pngdensity_quiz3.pngdensity_quiz4.pngdensity_quiz5.png

Jens Preben Morth, Ph.D
Aarhus University
Department of Molecular Biology
Gustav Wieds Vej 10 C
DK - 8000 Aarhus C
Tel. +45 8942 5257, Fax. +45 8612 3178
j...@mb.au.dk
website: http://person.au.dk/en/j...@mb.au.dk


Re: [ccp4bb] detection of anomalous signal using HKL2000

2010-06-29 Thread Felix Frolow

Quoting Phil Evans p...@mrc-lmb.cam.ac.uk:



Can anyone explain what Zbyszek Otwinowski means by Chi squared?


If I understand properly, CHI**2 value as used in Scalepack is:
SUM(I-Ij)**2/SUMsigma(I)**2 (I have to use formula editor to write  
it properly, but the idea is clear) and is useful in exhibiting and  
detection of systematic errors.
Intuitively this value will be close to 1.0 if only counting  
statistics contribute to errors in measurements of I.
Drawing CHI**2 distribution as a function of various values such as  
frame number, shell of resolution, intensity
of reflection etc. may show very interesting things related to the  
status of systematic errors during data collection.
In the case of absence of systematic errors and radiation decay and in  
presence of a stable X-ray beam these distributions will be  
featureless and their value will be dependent on a counting  
statistics, quality

of a detector and absorption of a crystal mainly.
Systematic errors of various sources change these distributions in a  
sensible way. Despite the fact that it is impossible to correct  
systematic errors, understanding of CHI**2 distributions lead to
better understanding of the experiment and in most cases these  
systematic errors can be eliminated leading to a perfect data shaped  
mostly by a counting statistics.
CHI**2  distributions a la Otwinowsky - Minor (or else, I also do  
not know if Zbyszek Otwinovsky developed it by himself or adopted from  
earlier sources),is as used in Scalepack an instrument of a great power.

FF

Dr  Felix Frolow
Professor of Structural Biology and Biotechnology
Department of Molecular Microbiology
and Biotechnology
Tel Aviv University 69978, Israel

Acta Crystallographica D, co-editor

e-mail: mbfro...@post.tau.ac.il
Tel:  ++972-3640-8723
Fax: ++972-3640-9407
Cellular: 0547 459 608




I can't find a definition in any of his papers (though I may have  
missed it). Is there a reference?


It doesn't seem obviously related to the chi squared distribution  
(In probability theory and statistics, the chi-square distribution  
(also chi-squared orχ²-distribution) with k degrees of freedom is  
the distribution of a sum of the squares of k independent standard  
normal random variables.   
http://en.wikipedia.org/wiki/Chi-square_distribution)


Phil


On 29 Jun 2010, at 21:14, Felix Frolow wrote:


Graphical information from Scalepack as used in HKL2000 is of unprecedented
help to detect anomalous signal. Anomalous detection for S  
anomalous data using
CHI**2 and Rfactor statistics for reflections with averaged and  
separated Bijvoet pairs is attached. It is very well described in  
HKL2000 manual. There is nothing special about data collection
(strategy was used) and measurement was relatively fast (4 h on  
MicroMax007 and RaxisIV++).
BTW Rfactor is 1.9% with separated Bijvoet reflections and 2.6%  
with averaged



BTW James Holton website calculate for this case 0.078 crystal

Dr  Felix Frolow
Professor of Structural Biology and Biotechnology
Department of Molecular Microbiology
and Biotechnology
Tel Aviv University 69978, Israel

Acta Crystallographica D, co-editor

e-mail: mbfro...@post.tau.ac.il
Tel:  ++972-3640-8723
Fax: ++972-3640-9407
Cellular: 0547 459 608








Quoting Bernhard Rupp b...@ruppweb.org:


I second the hkl2map/SHELXCDE approach. Two complete examples
explaining how to do this for MAD and S-SAD cases are in my book.
I wish to emphasize the importance of
a) running enough trials
b) careful selection of resolution cutoffs
c) look at the solution distribution
d) play with SHELXE parameters.
The hkl2map graphs are enormously helpful for this purpose.

BR
-Original Message-
From: CCP4 bulletin board [mailto:ccp...@jiscmail.ac.uk] On Behalf Of Tim
Gruene
Sent: Tuesday, June 29, 2010 2:29 AM
To: CCP4BB@JISCMAIL.AC.UK
Subject: Re: [ccp4bb] measure of anamolous signal

Dear Murugan,

you can use the program hkl2map from Thomas Schneider, available at
http://webapps.embl-hamburg.de/hkl2map/
It's a graphical interface to the programs shelx c/d/e which are available
from http://shelx.uni-ac.gwdg.de/SHELX/index.html

With SAD data you want to look at the d/sigma line at the end of  
the shelxc

output. Where that drops below about 1.3 is approximately where your
anomalous signal ends. You might get slightly improved statistics  
with xprep

instead of shelxc, but xprep is not free and you have to get a copy from
Bruker-AXS.

Tim

On Tue, Jun 29, 2010 at 02:35:29PM +0530, Vandu Murugan wrote:

Dear all,
  I have collected a 2.7 angstrom home source data with Cu-Kalpha
source for a protein with 6 cysteines, with a multiplicity of around
23.  I need to know, is there any significant anamolous signal present
in the data set, since there is no good model for my protein.  Can any
one tell, which program to run, and what parameter to see?  Thanks in

advance.


cheers,
Murugan


--
--
Tim Gruene
Institut fuer anorganische Chemie
Tammannstr. 4
D-37077 

Re: [ccp4bb] vacuum adapter for Dry shipper Cryo Dewar

2010-06-29 Thread James Holton


The part you are looking for is known as a valve operator, at least in 
the parlance of our supplier.   Here's a cut  paste of the item we 
purchased for pumping out our shipping dewars with a nominal 30mm O.D.:


V1000 Series Operator for 1 tube.
Part Number:  V1085-3-25 (KF-25 flange) @ $356.00 each

From:
DLH INDUSTRIES INC.
P.O. Box 562
San Luis Obispo, California 93406
(805) 781-3565   FAX (805) 781-3566
Email: sa...@cryocomp.com
Web Site: http://www.cryocomp.com
**

So far, we have tried to rescue 3-4 different people's dewars by pumping 
them down with this thing, and only one of them has worked.  Most of the 
time the vacuum fails because of a crack or leak on the inside of the 
dewar somewhere, and that cannot easily be fixed.  Still, the one 
success pays for the tool.


-James Holton
MAD Scientist


On 6/28/2010 9:34 PM, Frank von Delft wrote:
They seem to have a finite life-time:  we had two of these, bought at 
the same time, self-immolate the way you describe within months of 
another, after ~5 years of good service.


We asked TW about it and they said it would cost more to get it fixed 
than just get a new one -- then again, they would have said that, 
wouldn't they.  If you figure out how to re-suck the vacuum, do tell 
the list!


Cheers
phx.

P.S. Sorry, we don't have an adaptor -- you probably figured :)



On 29/06/2010 00:07, Gunnar Olovsson wrote:

Dear ccp4bb,

We have a *Taylor-Wharton *CX100 Dry shipper Dewar that has lost it's 
vacuum insulation. The vacuum sealing plug popped out of the valve 
(standard 30mm diameter).


Has anyone found a vacuum adapter that fits so that one can 
re-establish vacuum insulation at home?


I would ideally need to find a vacuum adapter that fits the valve of 
the Dewar and that ends in a KF16 Flange (preferably), so that we can 
re-establish vacuum insulation ourselves.


 Thank you very much in advance for your help!

 Sincerely - Gunnar


Dr. Gunnar Olovsson
Life Sciences Centre (4th floor)
University of British Columbia
2350 Health Sciences Mall
Vancouver,BC  V6T1Z3  Canada
gun...@byron.biochem.ubc.ca mailto:gun...@byron.biochem.ubc.ca




[ccp4bb] prediction of anomalous signal

2010-06-29 Thread James Holton

On 6/29/2010 12:55 PM, Felix Frolow wrote:

BTW James Holton website calculate for this case 0.078 crystal

   


So far, every single one of these reports has come down to a 
misinterpretation of the web interface.  I am trying to make it clearer, 
so I am interested in what values were entered into the various fields 
when someone gets a result that they think is not consistent with a real 
experiment that they did!


http://bl831.als.lbl.gov/xtalsize.html

BTW, this is a JavaScript program, so it runs on your web browser, not 
my server.  This means I am not spying on you, but it also means that 
there is no way for me to tell what people are typing into it.


-James Holton
MAD Scientist


Re: [ccp4bb] prediction of anomalous signal

2010-06-29 Thread James Holton
Woops!  Sorry in case I just confused a lot of people.  The number Felix 
reports: n_xtals = 0.078 means that the structure solution should have 
been easy (I.E. it could have been done with a crystal having only 8% 
of the volume used), and indeed it sounds like it was a slam dunk.  So, 
it sounds like Felix got a good prediction from it.


However, I am still very interested if anyone can solve a structure from 
a case where my predictor says that you can't.


-James Holton
MAD Scientist

On 6/29/2010 4:18 PM, James Holton wrote:

On 6/29/2010 12:55 PM, Felix Frolow wrote:

BTW James Holton website calculate for this case 0.078 crystal



So far, every single one of these reports has come down to a 
misinterpretation of the web interface.  I am trying to make it 
clearer, so I am interested in what values were entered into the 
various fields when someone gets a result that they think is not 
consistent with a real experiment that they did!


http://bl831.als.lbl.gov/xtalsize.html

BTW, this is a JavaScript program, so it runs on your web browser, not 
my server.  This means I am not spying on you, but it also means that 
there is no way for me to tell what people are typing into it.


-James Holton
MAD Scientist


Re: [ccp4bb] Density changes from Positive to negative after ligand addition and refinement

2010-06-29 Thread Frank von Delft
Negative (or positive) density tells you what you that what you've 
modelled is wrong.


If you have 1.15A X-ray data and it's telling you what you think is Ca2+ 
is not Ca2+, then I would think your X-ray data wins and your previous 
structures lose (a lot can happen from one structure to the next...).  I 
must say, your 2fo-fc looks like a perfectly good Mg2+, since it is the 
same size (ish) as the surrounding O and C atoms;  a Ca2+ would have 
given a much larger green blob.


phx.



On 30/06/2010 02:35, xaravich ivan wrote:

Dear CCP4BB,

I have come across something that might be pretty obvious to 
experienced people but is making me crazy.
I have this great 1.15 angs data and I know that I have a Calcium ion 
(pics attached) from previous structures of the same protein, that I 
have solved. Rightly when I add waters with Arp solvent it does not 
put water at that positive density.
Now whenever I have tried to put the calcium, and refine the structure 
it is giving me a negative density at the metal site. I csn see that 
the 2fc-fo is clear there, but why negative density. This is just the 
start of my refinement and I have to refine multiple ligands in the 
structure and I would like to get past this issue before that.


I tried putting atom at the pointer, adding water and renaming it 
according to the naming convention in the PDB for Calcium, but nothing.


Your suggestions would be invaluable, as always.

Ivan



Re: [ccp4bb] Density changes from Positive to negative after ligand addition and refinement

2010-06-29 Thread Pavel Afonine

Hi Ivan,

these Dale Tronrud's slides might help you with understanding the maps:

http://www.ccp4.ac.uk/courses/stwk10/talk_files/Dale_The_Wonderful_World_of_Maps.pdf

Good luck!
Pavel.


On 6/29/10 6:35 PM, xaravich ivan wrote:

Dear CCP4BB,

I have come across something that might be pretty obvious to 
experienced people but is making me crazy.
I have this great 1.15 angs data and I know that I have a Calcium ion 
(pics attached) from previous structures of the same protein, that I 
have solved. Rightly when I add waters with Arp solvent it does not 
put water at that positive density.
Now whenever I have tried to put the calcium, and refine the structure 
it is giving me a negative density at the metal site. I csn see that 
the 2fc-fo is clear there, but why negative density. This is just the 
start of my refinement and I have to refine multiple ligands in the 
structure and I would like to get past this issue before that.


I tried putting atom at the pointer, adding water and renaming it 
according to the naming convention in the PDB for Calcium, but nothing.


Your suggestions would be invaluable, as always.

Ivan



Re: [ccp4bb] Density changes from Positive to negative after ligand addition and  refinement

2010-06-29 Thread Frederic VELLIEUX
Could it be that for some reason (like the components of your solutions) you 
have Mg2+ in that site? Also, Magnesium 
is a common contaminant, trace amounts are usually present in one chemical or 
another.

Subtle differences in the site might make it a better site for another 
lighter ion (we still do not fully 
understand the subtility of protein structures, and of their modifications).

Fred.

Fred.

 Message du 30/06/10 03:36
 De : xaravich ivan 
 A : CCP4BB@JISCMAIL.AC.UK
 Copie à : 
 Objet : [ccp4bb] Density changes from Positive to negative after ligand 
 addition and  refinement
 
 Dear CCP4BB,
 
 I have come across something that might be pretty obvious to experienced
 people but is making me crazy.
 I have this great 1.15 angs data and I know that I have a Calcium ion (pics
 attached) from previous structures of the same protein, that I have solved.
 Rightly when I add waters with Arp solvent it does not put water at that
 positive density.
 Now whenever I have tried to put the calcium, and refine the structure it is
 giving me a negative density at the metal site. I csn see that the 2fc-fo is
 clear there, but why negative density. This is just the start of my
 refinement and I have to refine multiple ligands in the structure and I
 would like to get past this issue before that.
 
 I tried putting atom at the pointer, adding water and renaming it according
 to the naming convention in the PDB for Calcium, but nothing.
 
 Your suggestions would be invaluable, as always.
 
 Ivan
 

 [ Screen shot 2010-06-29 at 6.15.29 PM.png (375.3 Ko) ]
 [ Screen shot 2010-06-29 at 6.21.22 PM.png (383.5 Ko) ]