Re: [ccp4bb] Matthews coefficient

2010-06-25 Thread Frederic VELLIEUX
To reply to the original question, can the calculated Matthews coefficient be 
wrong?:
have you read the paper by Matthews? There is a distribution of solvent content 
in macromolecular (protein) 
crystals. It has a peak and tails on both sides. Meaning that there are 
outliers in the distribution, i.e. 
crystals with very little solvent, and crystals with plenty of solvent. They 
are not really outliers since, as in 
any physical phenomenon, some data points fall outside 1 sigma or 2 sigmas from 
the mean. Exceptions do exist.

What you cannot have however, is an asymmetric unit that contains too many 
molecules/subunits than what is 
physically possible.

Molecular replacement calculations can be performed searching for 1, 2, 3,... N 
copies in the asymmetric unit. This 
is what I personally do because at first I do not know the content of the 
asymmetric unit in my crystal.

Fred.

 Message du 25/06/10 07:54
 De : xinghua qin 
 A : CCP4BB@JISCMAIL.AC.UK
 Copie à : 
 Objet : Re: [ccp4bb] Matthews coefficient
 
 Hi everyone:
 Thanks for all the responses.
 The Matthaws-cell content analysis programe in CCP4 package gives the
 results: 47% solvent content and 3 molecules in asu with 87% confidence.
 the space group is P3121. how to carry out self rotation function? can
 phaser do that work? If there areTwo moleculars in asu ,no clashes and the
 TF values are above 10.but if three, the clashes are too much and TF valur
 is about five.
 And how to calculate the solvent content? Is the calculated solvent content
 with the Matthews-cell content analysis programe not always right?
 Best regards
 
 Xinghua Qin
 On Fri, Jun 25, 2010 at 10:42 AM, Vineet Gaur wrote:
 
  Hi
  It would b good if u can mention the solvent content and space group along
  with the no. of molecules in asu/
 
  u can carry out self rotation function to check the no of mol in asu
  best
  vineet
  On Fri, Jun 25, 2010 at 6:42 AM, xinghua qin  wrote:
 
  hi CCPeers
  The Matthews coefficient of my protein is 3 calculated with
  matthews-cell content analysis CCP4 programe with 87% confidence , but when
  doing the refinement the third molecular couldn't get into the unit cell
  because of too many clashes.Deletion of the clashed AA did not work well,
  Then I used two molecules in the unit cell, After refining with Refmac, I
  found that the R factor is 0.29, R free is 0.40.I believe the value can be
  better with the real space refinement.But the question is that can the
  calculated Matthews coefficient be wrong?
 
  Best regards
 
  Xinghua Qin


Re: [ccp4bb] Matthews coefficient

2010-06-25 Thread Tim Gruene
Dear Xinghua,

the solvent content provided by yhe Matthews-program is certainly correct.
Whether or not it applies to you protein is a different issue. With only 2
molecules in the asymmetric unit the solvent content rises to 65%, and this is
stll perfectly fine for proteins. So carry out molecular replacement with only
two molecules and start model building and refinement. 

Your previous message sounds as though to started refinement directly after
molecular replacement. I would strongly discourage you from this practice (which
somehow seems to be stuck in the heads of many people). After molecular
replacement you should check the result with the model building program of your
choice and correct as many errors as possible before running a refinement
program.

Tim

On Fri, Jun 25, 2010 at 01:54:05PM +0800, xinghua qin wrote:
 Hi everyone:
  Thanks for all the responses.
  The Matthaws-cell content analysis programe in CCP4 package gives the
 results: 47% solvent content and 3 molecules in asu with  87% confidence.
 the space group is P3121. how to carry out self rotation function? can
 phaser do that work? If there areTwo moleculars in asu ,no clashes and the
 TF values are above 10.but if three, the clashes are too much and TF valur
 is about five.
 And how to calculate the solvent content? Is the  calculated solvent content
 with the Matthews-cell content analysis programe not always right?
 Best regards
 
  Xinghua Qin
 On Fri, Jun 25, 2010 at 10:42 AM, Vineet Gaur vineetgaur1...@gmail.comwrote:
 
  Hi
  It would b good if u can mention the solvent content and space group along
  with the no. of molecules in asu/
 
  u can carry out self rotation function to check the no of mol in asu
  best
  vineet
On Fri, Jun 25, 2010 at 6:42 AM, xinghua qin xtal...@gmail.com wrote:
 
  hi CCPeers
The Matthews coefficient  of my protein is 3 calculated with
  matthews-cell content analysis CCP4 programe with  87% confidence , but 
  when
  doing the refinement the third molecular couldn't get into the unit cell
  because of  too many clashes.Deletion of the clashed AA did not work well,
  Then I used two molecules in the unit cell, After refining with Refmac, I
  found that the R factor is 0.29, R free is 0.40.I believe the value can be
  better with the real space refinement.But the question is that can the
  calculated Matthews coefficient be wrong?
 
  Best regards
 
  Xinghua Qin
 
  --
  Xinghua Qin
  College of Biological Sciences
  No.2, Yuan Ming Yuan West Road
  Haidian District,Beijing,China,100094
  Tel: +86-10-62732672
 
 
 
 
 
 -- 
 Xinghua Qin
 College of Biological Sciences
 No.2, Yuan Ming Yuan West Road
 Haidian District,Beijing,China,100094
 Tel: +86-10-62732672

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

GPG Key ID = A46BEE1A



signature.asc
Description: Digital signature


Re: [ccp4bb] Matthews coefficient

2010-06-25 Thread Francois Berenger

Hello,

I am not a crystallographer, so I will ask
a maybe stupid question.

When you say After molecular replacement you should check the result
with the model building program of your choice and correct as many
errors as possible before running a refinement program.

Can this step be automated in some way (without any
human intervention)?

Are there programs to do this?
I am interested even if some program does only one part of the job.

Thanks a lot,
Francois.

Tim Gruene wrote:

Dear Xinghua,

the solvent content provided by yhe Matthews-program is certainly correct.
Whether or not it applies to you protein is a different issue. With only 2
molecules in the asymmetric unit the solvent content rises to 65%, and this is
stll perfectly fine for proteins. So carry out molecular replacement with only
two molecules and start model building and refinement. 


Your previous message sounds as though to started refinement directly after
molecular replacement. I would strongly discourage you from this practice (which
somehow seems to be stuck in the heads of many people). After molecular
replacement you should check the result with the model building program of your
choice and correct as many errors as possible before running a refinement
program.

Tim

On Fri, Jun 25, 2010 at 01:54:05PM +0800, xinghua qin wrote:

Hi everyone:
 Thanks for all the responses.
 The Matthaws-cell content analysis programe in CCP4 package gives the
results: 47% solvent content and 3 molecules in asu with  87% confidence.
the space group is P3121. how to carry out self rotation function? can
phaser do that work? If there areTwo moleculars in asu ,no clashes and the
TF values are above 10.but if three, the clashes are too much and TF valur
is about five.
And how to calculate the solvent content? Is the  calculated solvent content
with the Matthews-cell content analysis programe not always right?
Best regards

 Xinghua Qin
On Fri, Jun 25, 2010 at 10:42 AM, Vineet Gaur vineetgaur1...@gmail.comwrote:


Hi
It would b good if u can mention the solvent content and space group along
with the no. of molecules in asu/

u can carry out self rotation function to check the no of mol in asu
best
vineet
  On Fri, Jun 25, 2010 at 6:42 AM, xinghua qin xtal...@gmail.com wrote:


hi CCPeers
  The Matthews coefficient  of my protein is 3 calculated with
matthews-cell content analysis CCP4 programe with  87% confidence , but when
doing the refinement the third molecular couldn't get into the unit cell
because of  too many clashes.Deletion of the clashed AA did not work well,
Then I used two molecules in the unit cell, After refining with Refmac, I
found that the R factor is 0.29, R free is 0.40.I believe the value can be
better with the real space refinement.But the question is that can the
calculated Matthews coefficient be wrong?

Best regards

Xinghua Qin

--
Xinghua Qin
College of Biological Sciences
No.2, Yuan Ming Yuan West Road
Haidian District,Beijing,China,100094
Tel: +86-10-62732672





--
Xinghua Qin
College of Biological Sciences
No.2, Yuan Ming Yuan West Road
Haidian District,Beijing,China,100094
Tel: +86-10-62732672




Re: [ccp4bb] Matthews coefficient

2010-06-25 Thread Vellieux Frederic

Hi,

No such thing as a stupid question. If the resolution is sufficient, 
arp_warp does a good job. But it is a reconstruction and REFINEMENT 
program. But checking the proper packing (that there are indeed contacts 
in the 3 dimensions of space to form the crystal): takes only 10 seconds 
using a graphics program, so there is no need to automate that. And 
trying to do everything completely automatically (without any human 
intervention and checks) can lead to big problems (for example if the 
space group has been wrongly assigned). You may end up publishing a 
wrong structure if you rely entirely on totally automated 
crystallographic software without checking anything.


Fred.

Francois Berenger wrote:

Hello,

I am not a crystallographer, so I will ask
a maybe stupid question.

When you say After molecular replacement you should check the result
with the model building program of your choice and correct as many
errors as possible before running a refinement program.

Can this step be automated in some way (without any
human intervention)?

Are there programs to do this?
I am interested even if some program does only one part of the job.

Thanks a lot,
Francois.


Re: [ccp4bb] Matthews coefficient

2010-06-25 Thread Francois Berenger

Vellieux Frederic wrote:

Hi,

No such thing as a stupid question. If the resolution is sufficient, 
arp_warp does a good job. But it is a reconstruction and REFINEMENT 
program. But checking the proper packing (that there are indeed contacts 
in the 3 dimensions of space to form the crystal): takes only 10 seconds 
using a graphics program, so there is no need to automate that.


Well, if my colleague sends you his 300,000 computer-generated models
to use in MR, I am sure you will like to automate the job. ;)

 And
trying to do everything completely automatically (without any human 
intervention and checks) can lead to big problems (for example if the 
space group has been wrongly assigned). You may end up publishing a 
wrong structure if you rely entirely on totally automated 
crystallographic software without checking anything.


Sure, an expert is always useful to check things at some moment.

And, I also agree that things made only in the computer have
to be checked using some experimental data in order to be relied on.

Francois.


Fred.

Francois Berenger wrote:

Hello,

I am not a crystallographer, so I will ask
a maybe stupid question.

When you say After molecular replacement you should check the result
with the model building program of your choice and correct as many
errors as possible before running a refinement program.

Can this step be automated in some way (without any
human intervention)?

Are there programs to do this?
I am interested even if some program does only one part of the job.

Thanks a lot,
Francois.


Re: [ccp4bb] Matthews coefficient

2010-06-25 Thread Nat Echols
On Fri, Jun 25, 2010 at 12:41 AM, Francois Berenger beren...@riken.jpwrote:

 When you say After molecular replacement you should check the result
 with the model building program of your choice and correct as many
 errors as possible before running a refinement program.

 Can this step be automated in some way (without any
 human intervention)?

 Are there programs to do this?
 I am interested even if some program does only one part of the job.


Both ARP/wARP and phenix.autobuild will do this (again, given data of
sufficiently high quality).  But I disagree with the original statement.
 Molecular replacement is typically done at reduced resolution, and
particularly where the search model is different in conformation or sequence
from the protein in the crystal, there will be many errors, even if the
solution is generally excellent.  Attempting to correct each of these
manually is a waste of time, since the refinement or rebuilding programs
will do it faster and better.  If there are severe problems like overlapping
molecules, this will become obvious very quickly whether or not automation
is used.  Is it even possible to solve, rebuild, and refine a structure from
MR in the wrong space group?  I would expect the R-factors to be
unreasonably high and the maps uninterpretable.

I do agree that inspecting the lattice visually after MR is a good idea (I
like PyMOL for this).  I don't know whether there are programs that will do
this automatically - it would be a useful tool to have.

As always, there will still be corner cases, where the search model is too
remote for the automated steps to work and manual fixing is required, or you
get packing clashes before all copies are found and need to prune the search
model, or the crystal resolution is too poor, and so on.  It is also
possible to waste time trying to improve a nonsensical solution - until your
R-free is well below 0.5, you can't be certain that you've solved the
structure, and if the automated methods don't push it below 0.4, it's
probably worth checking your MR solution to make sure it's already as good
as possible.

The more general point, that no one should publish or deposit a structure
without ever thoroughly checking the maps, validation, and lattice contacts,
I of course agree with - but I haven't heard anyone propose doing away with
human intervention altogether, nor are any of the automated pipelines
anywhere close to making this feasible.

-Nat


Re: [ccp4bb] Matthews coefficient

2010-06-25 Thread Christian Biertuempfel
Dear Xinghua,
Of course, you can have 2 molecules in the asu. However, if things don't
it is worth to double-check your space group or check for twinning.

Cheers
christian


xinghua qin wrote:
 hi CCPeers
   The Matthews coefficient  of my protein is 3 calculated with
 matthews-cell content analysis CCP4 programe with  87% confidence , but
 when doing the refinement the third molecular couldn't get into the unit
 cell because of  too many clashes.Deletion of the clashed AA did not
 work well, Then I used two molecules in the unit cell, After refining
 with Refmac, I found that the R factor is 0.29, R free is 0.40.I believe
 the value can be better with the real space refinement.But the question
 is that can the calculated Matthews coefficient be wrong?
  
 Best regards
  
 Xinghua Qin
 
 -- 
 Xinghua Qin
 College of Biological Sciences
 No.2, Yuan Ming Yuan West Road
 Haidian District,Beijing,China,100094
 Tel: +86-10-62732672

-- 
___

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
___


Re: [ccp4bb] Matthews coefficient

2010-06-25 Thread Bernhard Rupp
It is out of principle misleading to think of the Matthews coefficient as one 
hard number. Xinghua mentions it actually: 'with 87 % confidence' (although 
'confidence as in confidence interval' is a statistically too hard term here 
imho). The Matthews probability calculator gives a more complete picture of the 
probable values and their distribution.   
 
http://www.ruppweb.org/mattprob/default.html

BR

-Original Message-
From: CCP4 bulletin board [mailto:ccp...@jiscmail.ac.uk] On Behalf Of Christian 
Biertuempfel
Sent: Friday, June 25, 2010 9:34 AM
To: CCP4BB@JISCMAIL.AC.UK
Subject: Re: [ccp4bb] Matthews coefficient

Dear Xinghua,
Of course, you can have 2 molecules in the asu. However, if things don't
it is worth to double-check your space group or check for twinning.

Cheers
christian


xinghua qin wrote:
 hi CCPeers
   The Matthews coefficient  of my protein is 3 calculated with
 matthews-cell content analysis CCP4 programe with  87% confidence , but
 when doing the refinement the third molecular couldn't get into the unit
 cell because of  too many clashes.Deletion of the clashed AA did not
 work well, Then I used two molecules in the unit cell, After refining
 with Refmac, I found that the R factor is 0.29, R free is 0.40.I believe
 the value can be better with the real space refinement.But the question
 is that can the calculated Matthews coefficient be wrong?
  
 Best regards
  
 Xinghua Qin
 
 -- 
 Xinghua Qin
 College of Biological Sciences
 No.2, Yuan Ming Yuan West Road
 Haidian District,Beijing,China,100094
 Tel: +86-10-62732672

-- 
___

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
___


[ccp4bb] Matthews coefficient

2010-06-24 Thread xinghua qin
hi CCPeers
  The Matthews coefficient  of my protein is 3 calculated with matthews-cell
content analysis CCP4 programe with  87% confidence , but when doing the
refinement the third molecular couldn't get into the unit cell because of
too many clashes.Deletion of the clashed AA did not work well, Then I used
two molecules in the unit cell, After refining with Refmac, I found that the
R factor is 0.29, R free is 0.40.I believe the value can be better with the
real space refinement.But the question is that can the calculated Matthews
coefficient be wrong?

Best regards

Xinghua Qin

-- 
Xinghua Qin
College of Biological Sciences
No.2, Yuan Ming Yuan West Road
Haidian District,Beijing,China,100094
Tel: +86-10-62732672


Re: [ccp4bb] Matthews coefficient

2010-06-24 Thread xinghua qin
Hi everyone:
 Thanks for all the responses.
 The Matthaws-cell content analysis programe in CCP4 package gives the
results: 47% solvent content and 3 molecules in asu with  87% confidence.
the space group is P3121. how to carry out self rotation function? can
phaser do that work? If there areTwo moleculars in asu ,no clashes and the
TF values are above 10.but if three, the clashes are too much and TF valur
is about five.
And how to calculate the solvent content? Is the  calculated solvent content
with the Matthews-cell content analysis programe not always right?
Best regards

 Xinghua Qin
On Fri, Jun 25, 2010 at 10:42 AM, Vineet Gaur vineetgaur1...@gmail.comwrote:

 Hi
 It would b good if u can mention the solvent content and space group along
 with the no. of molecules in asu/

 u can carry out self rotation function to check the no of mol in asu
 best
 vineet
   On Fri, Jun 25, 2010 at 6:42 AM, xinghua qin xtal...@gmail.com wrote:

 hi CCPeers
   The Matthews coefficient  of my protein is 3 calculated with
 matthews-cell content analysis CCP4 programe with  87% confidence , but when
 doing the refinement the third molecular couldn't get into the unit cell
 because of  too many clashes.Deletion of the clashed AA did not work well,
 Then I used two molecules in the unit cell, After refining with Refmac, I
 found that the R factor is 0.29, R free is 0.40.I believe the value can be
 better with the real space refinement.But the question is that can the
 calculated Matthews coefficient be wrong?

 Best regards

 Xinghua Qin

 --
 Xinghua Qin
 College of Biological Sciences
 No.2, Yuan Ming Yuan West Road
 Haidian District,Beijing,China,100094
 Tel: +86-10-62732672





-- 
Xinghua Qin
College of Biological Sciences
No.2, Yuan Ming Yuan West Road
Haidian District,Beijing,China,100094
Tel: +86-10-62732672