Re: [ccp4bb] Fd-3m with pdbset - test pdb file

2010-09-06 Thread Herman . Schreuder
Dear bulletin board,

Tim asked me for the pdb file I used for testing (coot version 0.6),
which I have attached here. It is a totally artifical pdb file with one
random small ligand in it. Using a symmetry radius of 50 Angstrom, one
can get an idea of the symmetry in 'F d -3 m'. If in Tim's case with a
real protein coot comes to a grinding halt, it looks to me more like a
question of memory management than a problem with the symmetry
definitions.

Best,
Herman

-Original Message-
From: CCP4 bulletin board [mailto:ccp...@jiscmail.ac.uk] On Behalf Of
herman.schreu...@sanofi-aventis.com
Sent: Monday, September 06, 2010 2:44 PM
To: CCP4BB@JISCMAIL.AC.UK
Subject: Re: [ccp4bb] Fd-3m with pdbset

Hi Tim,

The string 'Fd-3m:1' is recognized by pdbset, but still produces 'F d -3
m :'. I did a test with a single small ligand with coot and if you
delete the ':' from the CRYST1 card to make 'F d -3 m' coot happily
accepts it and seems to produce the correct symmetry mates. I hope that
this way you are able to display your molecule. I agree that there must
be a bug somewhere.

Best,
Herman

-Original Message-
From: CCP4 bulletin board [mailto:ccp...@jiscmail.ac.uk] On Behalf Of
Tim Gruene
Sent: Monday, September 06, 2010 1:38 PM
To: CCP4BB@JISCMAIL.AC.UK
Subject: [ccp4bb] Fd-3m with pdbset

Hello,
I tried to assign a cell with space group F d -3 m to a PDB file using
pdbset.
Since it did not recognize the string 'F d -3 m' (nor 'F d -3 m :1') , I
used the space group number 227. This resulted in the space group string
'F d -3 m :'
in the PDB-file instead of 'F d -3 m :1' as it is listed in syminfo.lib,
and hence coot refused to display the symmetry mates for the molecule.

When I add the missing '1' in the resulting PDB-file, coot does display
the cell, but when I ask to display the symmetry mates, coot chokes and
consumed all my cpu power (for 8 molecules present in the unit cell...)
until I killed it.

I was wondering if this is a known bug in syminfo.lib (and/or the
dependent
programs) and if there is a way around it (without using
phenix.get_me_out_of_here, I could of course use xp instead of coot for
such a small molecule).

Kind regards, Tim

P.S.: The bug in pdbset persits also in pdbset from ccp4-6.1.3

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

GPG Key ID = A46BEE1A


set.pdb
Description: set.pdb


Re: [ccp4bb] Reverse Translatase

2010-09-06 Thread Greg Alushin

Hi Jacob-

What an intriguing proposition.  I can think of multiple reasons why  
such a system would not exist, but there is a mechanistic one which is  
most fundamental, having to do with the nature of the genetic code.


Say that there is a cellular machine which  would unfold a protein and  
expose it to some sort of reading system (already a hard problem).   
There is now the issue of transforming the amino acid information into  
nucleic acid information.  For simplicity let's assume that this  
system only uses one codon per amino acid, bypassing the degeneracy  
problem.


How would the cell then read off the amino acid sequence?  It seems  
the simplest solution would be analogous to translation, i.e. having  
adaptor molecules analogues to tRNAs which would guide an enzyme that  
synthesized the nucleic acid.  Otherwise, one would have to invoke the  
idea of a single enzyme recognizing every amino acid, which seems  
impossible to me.


As we know, the problem of protein-protein recognition is relatively  
complex.  At a minimum, one would need 20 adaptor proteins to  
recognize the 20 canonical amino acids: however, it seems unlikely  
that recognition of a single amino acid would be robust enough to  
select for the correct adaptor molecule.


So, let's say instead these adaptors recognize 2 amino acids at a time  
(still probably not robust enough).  Then, one would need 2^20  
adaptors, already a far greater number of gene products than that  
present in any genome than I know of...


It might be tempting to draw an analogy between this system and the  
immune system, where an incredible diversity is generated from a small  
number of genes.  However, diverse immune proteins all take the same  
input sequence (say antigen recognition) and lead to a single  
response, whereas this system has a 1 to 1 correspondence between  
inputs (protein sequence) and outputs (nucleic acid sequences),  and  
thus there is no way that a randomization system could generate the  
required diversity.


Cheers,
-Greg Alushin
Nogales lab
UC Berkeley

On Sep 6, 2010, at 7:12 PM, Michael Thompson wrote:


Jacob,

The idea is enticing, but don't forget that there are multiple  
degenerate codons for a given amino acid. Once the protein is  
synthesized, the specific codon information is lost.


I think that's a fundamental problem.

Keep the ideas coming,

Mike Thompson




- Original Message -
From: "Jacob Keller" 
To: CCP4BB@JISCMAIL.AC.UK
Sent: Monday, September 6, 2010 6:36:14 PM GMT -08:00 US/Canada  
Pacific

Subject: [ccp4bb] Reverse Translatase

Dear Crystallographers,

does anyone know of any conceptual reason why a reverse translatase  
enzyme
(protein-->nucleic acid) could not exist? I can think of so many  
things for
which such an enzyme would be helpful, both to cells and to  
scientists...!
Unless there is something I am missing, it would seem to me  
conceptually

almost impossible that it *not* exist.

Best Regards,

Jacob Keller


***
Jacob Pearson Keller
Northwestern University
Medical Scientist Training Program
Dallos Laboratory
F. Searle 1-240
2240 Campus Drive
Evanston IL 60208
lab: 847.491.2438
cel: 773.608.9185
email: j-kell...@northwestern.edu
***

--
Michael C. Thompson

Graduate Student

Biochemistry & Molecular Biology Division

Department of Chemistry & Biochemistry

University of California, Los Angeles

mi...@chem.ucla.edu


Re: [ccp4bb] Reverse Translatase

2010-09-06 Thread Michael Thompson
Jacob,

The idea is enticing, but don't forget that there are multiple degenerate 
codons for a given amino acid. Once the protein is synthesized, the specific 
codon information is lost.

I think that's a fundamental problem.

Keep the ideas coming,

Mike Thompson




- Original Message -
From: "Jacob Keller" 
To: CCP4BB@JISCMAIL.AC.UK
Sent: Monday, September 6, 2010 6:36:14 PM GMT -08:00 US/Canada Pacific
Subject: [ccp4bb] Reverse Translatase

Dear Crystallographers,

does anyone know of any conceptual reason why a reverse translatase enzyme 
(protein-->nucleic acid) could not exist? I can think of so many things for 
which such an enzyme would be helpful, both to cells and to scientists...! 
Unless there is something I am missing, it would seem to me conceptually 
almost impossible that it *not* exist.

Best Regards,

Jacob Keller


***
Jacob Pearson Keller
Northwestern University
Medical Scientist Training Program
Dallos Laboratory
F. Searle 1-240
2240 Campus Drive
Evanston IL 60208
lab: 847.491.2438
cel: 773.608.9185
email: j-kell...@northwestern.edu
***

-- 
Michael C. Thompson

Graduate Student

Biochemistry & Molecular Biology Division

Department of Chemistry & Biochemistry

University of California, Los Angeles

mi...@chem.ucla.edu


[ccp4bb] Reverse Translatase

2010-09-06 Thread Jacob Keller

Dear Crystallographers,

does anyone know of any conceptual reason why a reverse translatase enzyme 
(protein-->nucleic acid) could not exist? I can think of so many things for 
which such an enzyme would be helpful, both to cells and to scientists...! 
Unless there is something I am missing, it would seem to me conceptually 
almost impossible that it *not* exist.


Best Regards,

Jacob Keller


***
Jacob Pearson Keller
Northwestern University
Medical Scientist Training Program
Dallos Laboratory
F. Searle 1-240
2240 Campus Drive
Evanston IL 60208
lab: 847.491.2438
cel: 773.608.9185
email: j-kell...@northwestern.edu
***


Re: [ccp4bb] Can CCP4 calculate the R factor based on the given Fo, Fc, and sigFo by a mtz file?

2010-09-06 Thread Hailiang Zhang
Thanks a lot!

> Hailiang Zhang wrote:
>> Hi,
>>
>> I want to calculate the R factor based on the given Fo, Fc, and sigFo by
>> a
>> mtz file. Can some CCP4 tools do this? Thanks!
>>
>> Best Regards, Hailiang
>>
>
> Yes- rstats
>
> something like
>
> rstats hklin sfall.mtz hklout rescaled.mtz  < LABIN FP=FP SIGFP=SIGFP FC=FC PHIC=PHIC
> TITLE  FP column scaled to FC
> #RESOLUTION 999.0 2.6! default is 1 to 100 Ang
> PRINT LAST! default is LAST
> CYCLES 10   ! default is 6
> !LIST 0
> OUTPUT FOFC! default is FOFC
> SCALE 1.0   ! 52.5  ! default is 1.0
> !REJECT DELTA 4000  ! default is no rejections
> WEIGHTING_SCHEME NONE ! default is NONE
> WIDTH_OF_BINS RTHETA=0.005 FBINR=100! defaults are .01 and 1000
> PROCESS FCAL   ! default is FCAL
> eof-rstats
>
>
> This will rescale Fc to Fo or vice versa, depending on "process".
> But if they are already optimally scaled that shouldn't make any
> difference.
> it prints out R and weighted R in bins and overall.
>
> OR, scaleit, using
> labinFP=FP SIGFP=SIGFP   FPH1=Fc
> And take R-merge on F as R??
>
> OR, use fft to calculate an Fc,Phic map and use sfall to calculate
> structure factors
> from the map (should precisely replicate the input Fc) giving it the Fobs
> for reference
>
>


Re: [ccp4bb] Can CCP4 calculate the R factor based on the given Fo, Fc, and sigFo by a mtz file?

2010-09-06 Thread Edward A. Berry

Hailiang Zhang wrote:

Hi,

I want to calculate the R factor based on the given Fo, Fc, and sigFo by a
mtz file. Can some CCP4 tools do this? Thanks!

Best Regards, Hailiang



Yes- rstats

something like

rstats hklin sfall.mtz hklout rescaled.mtz  <

[ccp4bb] Can CCP4 calculate the R factor based on the given Fo, Fc, and sigFo by a mtz file?

2010-09-06 Thread Hailiang Zhang
Hi,

I want to calculate the R factor based on the given Fo, Fc, and sigFo by a
mtz file. Can some CCP4 tools do this? Thanks!

Best Regards, Hailiang


Re: [ccp4bb] Fd-3m with pdbset

2010-09-06 Thread George M. Sheldrick

I was going to say that we could adopt the convention that the 
default setting of a centrosymmetric space group is the one with the
origin on an inversion centre, since this is assumed by all programs
that I am aware of, thus making ':2' for 'second setting' unnecessary
and keeping the name to 11 characters. At first glance this would
solve the problem for Ian's two 'long' space groups ('P 42/n n m :2' 
and 'I 41/a m d :2') since only the second of the two settings given
in volume A of International Tables has the origin on an inversion
centre. However as chance would have it, both these space groups have
another alternative non-equivalent setting that also has the origin 
on an inversion centre but this third setting is not given in volume 
A of the IT!
   
Perhaps I should mention that with the SHELX method of specifying
the space group symmetry using the symmetry operators, alternative 
settings, specified in IT or not, cause no problems. But then we  
would not have had this thread (or the H3/R3 and P21221 threads).

George 

Prof. George M. Sheldrick FRS
Dept. Structural Chemistry,
University of Goettingen,
Tammannstr. 4,
D37077 Goettingen, Germany
Tel. +49-551-39-3021 or -3068
Fax. +49-551-39-22582


On Mon, 6 Sep 2010, Ian Tickle wrote:

> It seems that it's the new addition of ':n' which has pushed it over
> the limit of the 11 char field.  Has this been thought through?   I
> can see that the ':n' qualifier is necessary to disambiguate alternate
> origins.  Without the redundant formatting spaces 9 chars would be
> sufficient even for 'P 42/n n m :2' which is indeed the longest one at
> 13 chars (also 'I 41/a m d :2').
> 
> -- Ian
> 
> On Mon, Sep 6, 2010 at 2:57 PM, Martyn Winn  wrote:
> >> > There seem to be a few CCP4 Fortran files that have character spgnam*10
> >> > from the good old days ...
> >>
> >> What good old days were those? The Cryst1-format for PDB-files in the ccp4
> >> documentation says
> >> 56-66    Space group symbol, left justified (not used)
> >> which is eleven characters. Were there days when this was different?
> >>
> >> Tim
> >
> > I remember when we didn't bother with spacegroup on the
> > CRYST1 line at all, let alone formatting it correctly.
> >
> > And what about 'P 42/n n m :2' which is more than 11 characters?
> >
> > I realise my answer doesn't help with Coot problems, but I will clean up
> > CCP4 Fortran programs at least.
> >
> > m
> >
> > --
> > ***
> > *                                                                     *
> > *               Dr. Martyn Winn                                       *
> > *                                                                     *
> > *   STFC Daresbury Laboratory, Daresbury, Warrington, WA4 4AD, U.K.   *
> > *   Tel: +44 1925 603455    E-mail: martyn.w...@stfc.ac.uk            *
> > *   Fax: +44 1925 603634    Skype name: martyn.winn                   *
> > *             URL: http://www.ccp4.ac.uk/martyn/                      *
> > ***
> >
> 
> 

Re: [ccp4bb] Refmac bref MIXED failed

2010-09-06 Thread Lionel Costenaro
Hi Garib,

I tried with refmac 5.6 and the same inputs; it worked fine.

Thanks,
Lionel

2010/9/6 Garib N Murshudov 

> Could you please try refmac 5.6. I think we have fixed this problem. It is
> available from:
>
> www.ysbl.york.ac.uk/refmac/latest_refmac.html
>
> You need to take experimental version (it is now fairly stable)
>
> Regards
> Garib
>
> On 6 Sep 2010, at 15:01, Lionel Costenaro wrote:
>
> > Hello,
> >
> > When I run Refmac with the mixed iso/anisotripic B-value option (refi
> bref MIXED), refmac failed with the error message below.
> > My pdb comes from previous refmac runs with anisotropic B values for all
> atoms (H generated but not written).
> > I just removed from the pdb the ANISOU lines for all water molecules.
> >
> > Refmac with "refi bref ANIS" (all anisotropic B) works fine (normal
> termination) with the same input data and pdb - with all ANISOU lines.
> > CCP4Interface 2.0.5
> > Refmac_5.5.0102
> >
> > Any comments, help, fix are welcome.
> >
> > Best regards,
> > Lionel
> >
> > **
> >  CGMAT cycle number =  1
> >
> > .
> >  Trying gamma equal0.00
> >  Not converging with gamma equal0.00
> >  Trying gamma equal   5.503E-02
> >  Not converging with gamma equal   5.503E-02
> >  Trying gamma equal   0.1155000
> > .
> >  Not converging with gamma equal2655.967
> >  Trying gamma equal2921.619
> >  Not converging with gamma equal2921.619
> >  Trying gamma equal3213.836
> >  Gamma decreased to3155.393
> >
> ***
> > * Information from CCP4Interface script
> >
> ***
> > The program run with command: /all/programas/ccp4-6.1.2/bin/refmac5 XYZIN
> "/data/refmac8-coot-0Wiso.pdb" XYZOUT "/tmp/lcocri/21_2_pdb_1.tmp" HKLIN
> "/data/4_unique1.mtz" HKLOUT "/tmp/lcocri/21_3_mtz_1.tmp" LIBIN
> "/data/ligand.cif" LIBOUT "/data/21_lib.cif"
> > has failed with error message
> > child killed: SIGABRT
> >
> ***
> >
> >
> > #CCP4I TERMINATION STATUS 0 child killed: SIGABRT
> > #CCP4I TERMINATION TIME 06 sep 2010  15:28:12
> > #CCP4I MESSAGE Task failed
> >
> >
> >
> >
> >
>
>


-- 
___

Interactions faibles entre protéines en solution
Costenaro L. (2010)
Editions Universitaires Européennes, ISBN 9786131506048

Livre disponible en francais sur amazon:
http://www.amazon.com/s/ref=nb_sb_noss?url=search-alias%3Daps&field-keywords=ISBN+9786131506048&x=0&y=0

___

Dr. Lionel Costenaro

Institut de Biologia Molecular de Barcelona - CSIC
and Institute for Research in Biomedicine
Parc Cientific de Barcelona
C/ Baldiri i Reixac 10-12
E-08028 Barcelona, Spain

Phone  +34 93 403 49 57
Fax  +34 93 403 49 79
Emaillionel.costen...@ibmb.csic.es ,  lionel.costen...@irbbarcelona.org
URL  www.ibmb.csic.es ,  www.irbbarcelona.org


Re: [ccp4bb] Fd-3m with pdbset

2010-09-06 Thread Ian Tickle
It seems that it's the new addition of ':n' which has pushed it over
the limit of the 11 char field.  Has this been thought through?   I
can see that the ':n' qualifier is necessary to disambiguate alternate
origins.  Without the redundant formatting spaces 9 chars would be
sufficient even for 'P 42/n n m :2' which is indeed the longest one at
13 chars (also 'I 41/a m d :2').

-- Ian

On Mon, Sep 6, 2010 at 2:57 PM, Martyn Winn  wrote:
>> > There seem to be a few CCP4 Fortran files that have character spgnam*10
>> > from the good old days ...
>>
>> What good old days were those? The Cryst1-format for PDB-files in the ccp4
>> documentation says
>> 56-66    Space group symbol, left justified (not used)
>> which is eleven characters. Were there days when this was different?
>>
>> Tim
>
> I remember when we didn't bother with spacegroup on the
> CRYST1 line at all, let alone formatting it correctly.
>
> And what about 'P 42/n n m :2' which is more than 11 characters?
>
> I realise my answer doesn't help with Coot problems, but I will clean up
> CCP4 Fortran programs at least.
>
> m
>
> --
> ***
> *                                                                     *
> *               Dr. Martyn Winn                                       *
> *                                                                     *
> *   STFC Daresbury Laboratory, Daresbury, Warrington, WA4 4AD, U.K.   *
> *   Tel: +44 1925 603455    E-mail: martyn.w...@stfc.ac.uk            *
> *   Fax: +44 1925 603634    Skype name: martyn.winn                   *
> *             URL: http://www.ccp4.ac.uk/martyn/                      *
> ***
>


Re: [ccp4bb] Problem with NCS detection in Parrot

2010-09-06 Thread LEGRAND Pierre
Hi Florian,

I have found that, in some cases, playing with -ncs-mask-filter-radius can help 
parrot to successfully detect and apply ncs averaging.
In one recent case, with 8 molecules where NCS was not found without this 
option, values from 14 to 30 angstroms worked (but gave maps of variable 
quality).
I don't know what's the default value for this.

Best regards,

Pierre Legrand, PROXIMA-1
Phone-office   +33 1 69 35 96 60
Phone-beamline +33 1 69 35 97 38
http://www.synchrotron-soleil.fr/Recherche/LignesLumiere/PROXIMA1

-Original Message-
From: CCP4 bulletin board on behalf of Florian Schmitzberger
Sent: Fri 9/3/2010 9:34 PM
To: CCP4BB@JISCMAIL.AC.UK
Subject: [ccp4bb] Problem with NCS detection in Parrot
 
Dear All,

I am encountering a problem when using Parrot (for combined density  
modification and non crystallographic symmetry (NCS) averaging) in  
ccp4 6.1.13, run via ccp4i.

Parrot does not detect the (2-fold) NCS present among my heavy atom  
substructure with 20 seleniums (the pdb was output by Phaser-EP,  
single chain ID, and is read by Parrot from what I can tell). I have  
tried a to split the NCS related heavy atoms into separate chains, but  
Parrot does still not appear to detect any NCS (error message:  
WARNING: No NCS found from heavy atoms).

The Professs program seems to detect the NCS readily. Unfortunately, I  
don't think it is possible to input externally determined NCS  
operators into Parrot.

Regards,

Florian

---
Florian Schmitzberger
Biological Chemistry and Molecular Pharmacology
Harvard Medical School
250 Longwood Avenue, SGM 130
Boston, MA 02115, US
Tel: 001 617 432 5602


Re: [ccp4bb] Refmac bref MIXED failed

2010-09-06 Thread Boaz Shaanan
Hi,

I'm not sure whether it's related to the error you're seeing but as far as I 
know, in order for refmac to do mixed B-refinement (that is Biso and Baniso) 
you have to set yourself by editing (!!) the incoming pdb such that whichever 
atoms you'd like to have Biso refined should have only Biso in the input file. 
I suggested a while ago that this be fixed to become more elegant (say by 
specifying range of residues in the gui) but apparently it still works only the 
old way.

   Cheers,

  Boaz

- Original Message -
From: Lionel Costenaro 
Date: Monday, September 6, 2010 17:09
Subject: [ccp4bb] Refmac bref MIXED failed
To: CCP4BB@JISCMAIL.AC.UK

> Hello,
> 
> When I run Refmac with the mixed iso/anisotripic B-value option 
> (refi bref
> MIXED), refmac failed with the error message below.
> My pdb comes from previous refmac runs with anisotropic B values 
> for all
> atoms (H generated but not written).
> I just removed from the pdb the ANISOU lines for all water molecules.
> 
> Refmac with "refi bref ANIS" (all anisotropic B) works fine (normal
> termination) with the same input data and pdb - with all ANISOU lines.
> CCP4Interface 2.0.5
> Refmac_5.5.0102
> 
> Any comments, help, fix are welcome.
> 
> Best regards,
> Lionel
> 
> **
>  CGMAT cycle number 
> =  1
> 
> .
>  Trying gamma equal    0.00
>  Not converging with gamma equal    0.00
>  Trying gamma equal   5.503E-02
>  Not converging with gamma equal   5.503E-02
>  Trying gamma equal   0.1155000
> .
>  Not converging with gamma equal    2655.967
>  Trying gamma equal    2921.619
>  Not converging with gamma equal    2921.619
>  Trying gamma equal    3213.836
>  Gamma decreased to    3155.393
> ***
> * Information from CCP4Interface script
> ***
> The program run with command: /all/programas/ccp4-
> 6.1.2/bin/refmac5 XYZIN
> "/data/refmac8-coot-0Wiso.pdb" XYZOUT 
> "/tmp/lcocri/21_2_pdb_1.tmp" HKLIN
> "/data/4_unique1.mtz" HKLOUT "/tmp/lcocri/21_3_mtz_1.tmp" LIBIN
> "/data/ligand.cif" LIBOUT "/data/21_lib.cif"
> has failed with error message
> child killed: SIGABRT
> ***
> 
> 
> #CCP4I TERMINATION STATUS 0 child killed: SIGABRT
> #CCP4I TERMINATION TIME 06 sep 2010  15:28:12
> #CCP4I MESSAGE Task failed
> 

Boaz Shaanan, Ph.D.
Dept. of Life Sciences
Ben-Gurion University of the Negev
Beer-Sheva 84105
Israel
Phone: 972-8-647-2220 ; Fax: 646-1710
Skype: boaz.shaanan‎


Re: [ccp4bb] Refmac bref MIXED failed

2010-09-06 Thread Garib N Murshudov
Could you please try refmac 5.6. I think we have fixed this problem. It is 
available from:

www.ysbl.york.ac.uk/refmac/latest_refmac.html

You need to take experimental version (it is now fairly stable)

Regards
Garib

On 6 Sep 2010, at 15:01, Lionel Costenaro wrote:

> Hello, 
> 
> When I run Refmac with the mixed iso/anisotripic B-value option (refi bref 
> MIXED), refmac failed with the error message below.
> My pdb comes from previous refmac runs with anisotropic B values for all 
> atoms (H generated but not written).
> I just removed from the pdb the ANISOU lines for all water molecules.
> 
> Refmac with "refi bref ANIS" (all anisotropic B) works fine (normal 
> termination) with the same input data and pdb - with all ANISOU lines.
> CCP4Interface 2.0.5
> Refmac_5.5.0102
> 
> Any comments, help, fix are welcome.
> 
> Best regards, 
> Lionel
> 
> **
>  CGMAT cycle number =  1
>  
> .
>  Trying gamma equal0.00
>  Not converging with gamma equal0.00
>  Trying gamma equal   5.503E-02
>  Not converging with gamma equal   5.503E-02
>  Trying gamma equal   0.1155000
> .
>  Not converging with gamma equal2655.967
>  Trying gamma equal2921.619
>  Not converging with gamma equal2921.619
>  Trying gamma equal3213.836
>  Gamma decreased to3155.393
> ***
> * Information from CCP4Interface script
> ***
> The program run with command: /all/programas/ccp4-6.1.2/bin/refmac5 XYZIN 
> "/data/refmac8-coot-0Wiso.pdb" XYZOUT "/tmp/lcocri/21_2_pdb_1.tmp" HKLIN 
> "/data/4_unique1.mtz" HKLOUT "/tmp/lcocri/21_3_mtz_1.tmp" LIBIN 
> "/data/ligand.cif" LIBOUT "/data/21_lib.cif" 
> has failed with error message
> child killed: SIGABRT
> ***
> 
> 
> #CCP4I TERMINATION STATUS 0 child killed: SIGABRT
> #CCP4I TERMINATION TIME 06 sep 2010  15:28:12
> #CCP4I MESSAGE Task failed
> 
> 
> 
> 
> 


Re: [ccp4bb] fail to install ccp4i task for arp warp

2010-09-06 Thread Huw Jenkins
On 31 Aug 2010, at 14:56, Jochen Kuper wrote:

> sorry to bother you with this but I just can't figure out what goes wrong ...
> 
> I am running the latest version of snow leopard on a MBP unibody with the 
> 64bit kernel.
> 
> I have an up to date fink installation for native 64bit. 

I've seen this too - it seems to be due to problems with 64bit bltwish. There 
is a workaround here: 

http://proteincrystallography.org/ccp4bb/message13191.html

but that unfortunately requires a 32 bit fink installation of CCP4 alongside 
the 64bit one.

I can confirm that this works and once installed ARP/wARP and 64bit CCP4 work 
together fine. 



Huw
--
Dr Huw Jenkins
Astbury Centre for Structural Molecular Biology 
University of Leeds


[ccp4bb] Refmac bref MIXED failed

2010-09-06 Thread Lionel Costenaro
Hello,

When I run Refmac with the mixed iso/anisotripic B-value option (refi bref
MIXED), refmac failed with the error message below.
My pdb comes from previous refmac runs with anisotropic B values for all
atoms (H generated but not written).
I just removed from the pdb the ANISOU lines for all water molecules.

Refmac with "refi bref ANIS" (all anisotropic B) works fine (normal
termination) with the same input data and pdb - with all ANISOU lines.
CCP4Interface 2.0.5
Refmac_5.5.0102

Any comments, help, fix are welcome.

Best regards,
Lionel

**
 CGMAT cycle number =  1

.
 Trying gamma equal0.00
 Not converging with gamma equal0.00
 Trying gamma equal   5.503E-02
 Not converging with gamma equal   5.503E-02
 Trying gamma equal   0.1155000
.
 Not converging with gamma equal2655.967
 Trying gamma equal2921.619
 Not converging with gamma equal2921.619
 Trying gamma equal3213.836
 Gamma decreased to3155.393
***
* Information from CCP4Interface script
***
The program run with command: /all/programas/ccp4-6.1.2/bin/refmac5 XYZIN
"/data/refmac8-coot-0Wiso.pdb" XYZOUT "/tmp/lcocri/21_2_pdb_1.tmp" HKLIN
"/data/4_unique1.mtz" HKLOUT "/tmp/lcocri/21_3_mtz_1.tmp" LIBIN
"/data/ligand.cif" LIBOUT "/data/21_lib.cif"
has failed with error message
child killed: SIGABRT
***


#CCP4I TERMINATION STATUS 0 child killed: SIGABRT
#CCP4I TERMINATION TIME 06 sep 2010  15:28:12
#CCP4I MESSAGE Task failed


Re: [ccp4bb] Fd-3m with pdbset

2010-09-06 Thread Martyn Winn
> > There seem to be a few CCP4 Fortran files that have character spgnam*10
> > from the good old days ...
> 
> What good old days were those? The Cryst1-format for PDB-files in the ccp4
> documentation says
> 56-66Space group symbol, left justified (not used)
> which is eleven characters. Were there days when this was different?
> 
> Tim

I remember when we didn't bother with spacegroup on the
CRYST1 line at all, let alone formatting it correctly.

And what about 'P 42/n n m :2' which is more than 11 characters?

I realise my answer doesn't help with Coot problems, but I will clean up
CCP4 Fortran programs at least.

m

-- 
***
* *
*   Dr. Martyn Winn   *
* *
*   STFC Daresbury Laboratory, Daresbury, Warrington, WA4 4AD, U.K.   *
*   Tel: +44 1925 603455E-mail: martyn.w...@stfc.ac.uk*
*   Fax: +44 1925 603634Skype name: martyn.winn   * 
* URL: http://www.ccp4.ac.uk/martyn/  *
***


Re: [ccp4bb] Fd-3m with pdbset

2010-09-06 Thread Tim Gruene
On Mon, Sep 06, 2010 at 02:45:49PM +0100, Martyn Winn wrote:
> 
> Yes, I can reproduce this, and I think I have found the bug.
> It is simply that PDBSET dimensions the spacegroup name to 10
> characters. Increasing that to 20 allows you to use either
> 
> spac 'F d -3 m :1'
> or
> spac 227
> 
> However, it still won't recognise 'F d -3 m' because this is not in the
> syminfo file (at least, not in a field that is read).
> 
> There seem to be a few CCP4 Fortran files that have character spgnam*10
> from the good old days ...

What good old days were those? The Cryst1-format for PDB-files in the ccp4
documentation says
56-66Space group symbol, left justified (not used)
which is eleven characters. Were there days when this was different?

Tim

> 
> Cheers
> Martyn
> 
> On Mon, 2010-09-06 at 13:38 +0200, Tim Gruene wrote:
> > Hello,
> > I tried to assign a cell with space group F d -3 m to a PDB file using 
> > pdbset.
> > Since it did not recognize the string 'F d -3 m' (nor 'F d -3 m :1') , I 
> > used
> > the space group number 227. This resulted in the space group string 'F d -3 
> > m :'
> > in the PDB-file instead of 'F d -3 m :1' as it is listed in syminfo.lib, and
> > hence coot refused to display the symmetry mates for the molecule.
> > 
> > When I add the missing '1' in the resulting PDB-file, coot does display the
> > cell, but when I ask to display the symmetry mates, coot chokes and 
> > consumed all
> > my cpu power (for 8 molecules present in the unit cell...) until I killed 
> > it.
> > 
> > I was wondering if this is a known bug in syminfo.lib (and/or the dependent
> > programs) and if there is a way around it (without using
> > phenix.get_me_out_of_here, I could of course use xp instead of coot for 
> > such a
> > small molecule).
> > 
> > Kind regards, Tim
> > 
> > P.S.: The bug in pdbset persits also in pdbset from ccp4-6.1.3
> > 
> 
> -- 
> ***
> * *
> *   Dr. Martyn Winn   *
> * *
> *   STFC Daresbury Laboratory, Daresbury, Warrington, WA4 4AD, U.K.   *
> *   Tel: +44 1925 603455E-mail: martyn.w...@stfc.ac.uk*
> *   Fax: +44 1925 603634Skype name: martyn.winn   * 
> * URL: http://www.ccp4.ac.uk/martyn/  *
> ***

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

GPG Key ID = A46BEE1A



signature.asc
Description: Digital signature


Re: [ccp4bb] Fd-3m with pdbset

2010-09-06 Thread Tim Gruene
Dear Herman,

thanks for the tip. Unfortunately it does not help with my version of coot
(linux, 0.6.1). Coot eats up all my memory (8GB RAM + 4GB swap) until it is
being terminated by the OS.

Which version of coot are you using? Could you send me your PDB-file that works
for you, please (you could even send it plain text to the bb to avoid an
attachment...) ?

The same happens when I read in an .ins-file which does not contain the space
group name or symbol, but provides the symmetry operators themselves.

Cheers, Tim

P.S.: I'll skip sending a copy to the coot-list since Paul is also reading this
one, I am sure.

On Mon, Sep 06, 2010 at 02:43:34PM +0200, herman.schreu...@sanofi-aventis.com 
wrote:
> Hi Tim,
> 
> The string 'Fd-3m:1' is recognized by pdbset, but still produces 'F d -3
> m :'. I did a test with a single small ligand with coot and if you
> delete the ':' from the CRYST1 card to make 'F d -3 m' coot happily
> accepts it and seems to produce the correct symmetry mates. I hope that
> this way you are able to display your molecule. I agree that there must
> be a bug somewhere.
> 
> Best,
> Herman
> 
> -Original Message-
> From: CCP4 bulletin board [mailto:ccp...@jiscmail.ac.uk] On Behalf Of
> Tim Gruene
> Sent: Monday, September 06, 2010 1:38 PM
> To: CCP4BB@JISCMAIL.AC.UK
> Subject: [ccp4bb] Fd-3m with pdbset
> 
> Hello,
> I tried to assign a cell with space group F d -3 m to a PDB file using
> pdbset.
> Since it did not recognize the string 'F d -3 m' (nor 'F d -3 m :1') , I
> used the space group number 227. This resulted in the space group string
> 'F d -3 m :'
> in the PDB-file instead of 'F d -3 m :1' as it is listed in syminfo.lib,
> and hence coot refused to display the symmetry mates for the molecule.
> 
> When I add the missing '1' in the resulting PDB-file, coot does display
> the cell, but when I ask to display the symmetry mates, coot chokes and
> consumed all my cpu power (for 8 molecules present in the unit cell...)
> until I killed it.
> 
> I was wondering if this is a known bug in syminfo.lib (and/or the
> dependent
> programs) and if there is a way around it (without using
> phenix.get_me_out_of_here, I could of course use xp instead of coot for
> such a small molecule).
> 
> Kind regards, Tim
> 
> P.S.: The bug in pdbset persits also in pdbset from ccp4-6.1.3
> 
> --
> --
> Tim Gruene
> Institut fuer anorganische Chemie
> Tammannstr. 4
> D-37077 Goettingen
> 
> GPG Key ID = A46BEE1A

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

GPG Key ID = A46BEE1A



signature.asc
Description: Digital signature


Re: [ccp4bb] Fd-3m with pdbset

2010-09-06 Thread Martyn Winn
Yes, I can reproduce this, and I think I have found the bug.
It is simply that PDBSET dimensions the spacegroup name to 10
characters. Increasing that to 20 allows you to use either

spac 'F d -3 m :1'
or
spac 227

However, it still won't recognise 'F d -3 m' because this is not in the
syminfo file (at least, not in a field that is read).

There seem to be a few CCP4 Fortran files that have character spgnam*10
from the good old days ...

Cheers
Martyn

On Mon, 2010-09-06 at 13:38 +0200, Tim Gruene wrote:
> Hello,
> I tried to assign a cell with space group F d -3 m to a PDB file using pdbset.
> Since it did not recognize the string 'F d -3 m' (nor 'F d -3 m :1') , I used
> the space group number 227. This resulted in the space group string 'F d -3 m 
> :'
> in the PDB-file instead of 'F d -3 m :1' as it is listed in syminfo.lib, and
> hence coot refused to display the symmetry mates for the molecule.
> 
> When I add the missing '1' in the resulting PDB-file, coot does display the
> cell, but when I ask to display the symmetry mates, coot chokes and consumed 
> all
> my cpu power (for 8 molecules present in the unit cell...) until I killed it.
> 
> I was wondering if this is a known bug in syminfo.lib (and/or the dependent
> programs) and if there is a way around it (without using
> phenix.get_me_out_of_here, I could of course use xp instead of coot for such a
> small molecule).
> 
> Kind regards, Tim
> 
> P.S.: The bug in pdbset persits also in pdbset from ccp4-6.1.3
> 

-- 
***
* *
*   Dr. Martyn Winn   *
* *
*   STFC Daresbury Laboratory, Daresbury, Warrington, WA4 4AD, U.K.   *
*   Tel: +44 1925 603455E-mail: martyn.w...@stfc.ac.uk*
*   Fax: +44 1925 603634Skype name: martyn.winn   * 
* URL: http://www.ccp4.ac.uk/martyn/  *
***


[ccp4bb] Post-doc vacancy in the Integral Membrane Protein group, SGC, Oxford

2010-09-06 Thread Liz Carpenter
Dear All,

The Integral Membrane Protein group at the SGC in Oxford is looking for a 
post-doctoral research scientist with 
experience in purification and crystallisation of membrane proteins. In the 
past year we have developed a pipeline
for the production of human membrane proteins in the baculovirus/insect cell 
expression system.  We are now optimizing  
purification and crystallisation protocols for a series of well-expressed and 
purified membrane proteins in a broad spectrum 
of families, including GPCRs, channels, solute carriers, ABC transporters and 
enzymes. 

We are looking for an enthusiastic and highly motivated scientist to join our 
team of 12 membrane protein scientists to help with this work.   

For further information on this post, please follow this link 
http://www.sgc.ox.ac.uk/jobs/H910012/H910012.html.  

For informal enquiries please  contact Dr Liz Carpenter, the Principal 
Investigator for the IMP group:  liz.carpen...@sgc.ox.ac.uk.

Liz


Re: [ccp4bb] Problem with NCS detection in Parrot

2010-09-06 Thread Kevin Cowtan

Florian Schmitzberger wrote:

Dear All,

I am encountering a problem when using Parrot (for combined density 
modification and non crystallographic symmetry (NCS) averaging) in ccp4 
6.1.13, run via ccp4i.


Parrot does not detect the (2-fold) NCS present among my heavy atom 
substructure with 20 seleniums (the pdb was output by Phaser-EP, single 
chain ID, and is read by Parrot from what I can tell). I have tried a to 
split the NCS related heavy atoms into separate chains, but Parrot does 
still not appear to detect any NCS (error message: WARNING: No NCS found 
from heavy atoms).


Can I try it?

The Professs program seems to detect the NCS readily. Unfortunately, I 
don't think it is possible to input externally determined NCS operators 
into Parrot.


Yup, you can, but you have to 'run and edit com file'.
 ncs-operator alpha,beta,gamma,x,y,z,x,y,z
First set of x,y,z are the centre of mass of the source molecule, second 
are the target molecule. This info is provided by recent versions of 
professs.


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


Re: [ccp4bb] Fd-3m with pdbset

2010-09-06 Thread Herman . Schreuder
Hi Tim,

The string 'Fd-3m:1' is recognized by pdbset, but still produces 'F d -3
m :'. I did a test with a single small ligand with coot and if you
delete the ':' from the CRYST1 card to make 'F d -3 m' coot happily
accepts it and seems to produce the correct symmetry mates. I hope that
this way you are able to display your molecule. I agree that there must
be a bug somewhere.

Best,
Herman

-Original Message-
From: CCP4 bulletin board [mailto:ccp...@jiscmail.ac.uk] On Behalf Of
Tim Gruene
Sent: Monday, September 06, 2010 1:38 PM
To: CCP4BB@JISCMAIL.AC.UK
Subject: [ccp4bb] Fd-3m with pdbset

Hello,
I tried to assign a cell with space group F d -3 m to a PDB file using
pdbset.
Since it did not recognize the string 'F d -3 m' (nor 'F d -3 m :1') , I
used the space group number 227. This resulted in the space group string
'F d -3 m :'
in the PDB-file instead of 'F d -3 m :1' as it is listed in syminfo.lib,
and hence coot refused to display the symmetry mates for the molecule.

When I add the missing '1' in the resulting PDB-file, coot does display
the cell, but when I ask to display the symmetry mates, coot chokes and
consumed all my cpu power (for 8 molecules present in the unit cell...)
until I killed it.

I was wondering if this is a known bug in syminfo.lib (and/or the
dependent
programs) and if there is a way around it (without using
phenix.get_me_out_of_here, I could of course use xp instead of coot for
such a small molecule).

Kind regards, Tim

P.S.: The bug in pdbset persits also in pdbset from ccp4-6.1.3

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

GPG Key ID = A46BEE1A


[ccp4bb] Fd-3m with pdbset

2010-09-06 Thread Tim Gruene
Hello,
I tried to assign a cell with space group F d -3 m to a PDB file using pdbset.
Since it did not recognize the string 'F d -3 m' (nor 'F d -3 m :1') , I used
the space group number 227. This resulted in the space group string 'F d -3 m :'
in the PDB-file instead of 'F d -3 m :1' as it is listed in syminfo.lib, and
hence coot refused to display the symmetry mates for the molecule.

When I add the missing '1' in the resulting PDB-file, coot does display the
cell, but when I ask to display the symmetry mates, coot chokes and consumed all
my cpu power (for 8 molecules present in the unit cell...) until I killed it.

I was wondering if this is a known bug in syminfo.lib (and/or the dependent
programs) and if there is a way around it (without using
phenix.get_me_out_of_here, I could of course use xp instead of coot for such a
small molecule).

Kind regards, Tim

P.S.: The bug in pdbset persits also in pdbset from ccp4-6.1.3

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

GPG Key ID = A46BEE1A



signature.asc
Description: Digital signature