Closed GoogleCodeExporter closed 9 years ago
One comment already: can't find the 134d joinedCoord file on sunfish, or has
that changed? Also no joy for 1gnc, 1msh, 2neo, 2ofg, ... didn't try the others
Original comment by wfvran...@gmail.com
on 27 Jan 2011 at 2:33
These entries have no experimental data in PDB. They're aren't processed at
BMRB but in Nijmegen at the above locations. Can you look at eg:
http://nmr.cmbi.ru.nl/NRG-CING/prep/C/34/134d/ ?
I attach in this comment the overall flow for NRG-CING.
Original comment by jurge...@gmail.com
on 31 Jan 2011 at 9:07
Attachments:
This is still an issue when I redid the preparation and validation steps for
this entry yesterday.
Original comment by jurge...@gmail.com
on 1 Feb 2011 at 10:07
Similar problem?
================================================================================
======================
| CING: Common Interface for NMR structure Generation version 0.95 (r918)
AW,JFD,GWV 2004-2011 |
================================================================================
======================
User: jd on: nmr.cmbi.umcn.nl (darwin/32bit/8cores/2.6.5) at: (9564)
Mon Jan 31 21:46:39 2011
DEBUG: Using:
DEBUG: inputDir:
file:///Library/WebServer/Documents/NRG-CING/input/ki
DEBUG: outputDir: /Library/WebServer/Documents/NRG-CING/data/ki/2kid
DEBUG: pdbConvention: .
DEBUG: restraintsConvention: .
DEBUG: archiveType: BY_CH23
DEBUG: projectType: CCPN
DEBUG: modelCount: None
DEBUG: storeCING2db: True
DEBUG: ranges: cv
Enforcing a redo
copying file: //Library/WebServer/Documents/NRG-CING/input/ki/2kid.tgz to:
2kid.tgz
Moving CCPN directory from [linkNmrStarData] to [2kid]
==> Importing data from Ccpn project "2kid"
Error loading file for: <ccp.molecule.ChemComp.NonStdChemComp ['protein',
'Phq']>
Reading: <open file
'/Volumes/tera4/NRG-CING/data/ki/2kid/2kid/ccp/molecule/ChemComp/protein+Phq+2ki
d_user_2010-12-08-13-35-37-016_00001.xml', mode 'r' at 0x7faf498>
Last xml tag read: _StorageUnit
Parser state was: checking validity
Object stack was empty
Error loading file for: <ccp.molecule.MolSystem.MolSystem ['2kid']>
Reading: <open file
'/Volumes/tera4/NRG-CING/data/ki/2kid/2kid/ccp/molecule/MolSystem/2kid+2kid_user
_2010-12-08-13-34-39-344_00002.xml', mode 'r' at 0x772bf98>
Last xml tag read: _StorageUnit
Parser state was: checking validity
Object stack was empty
Error loading file for: <ccp.nmr.NmrConstraint.NmrConstraintStore [1]>
Reading: <open file
'/Volumes/tera4/NRG-CING/data/ki/2kid/2kid/ccp/nmr/NmrConstraint/1+2kid_user_201
0-12-08-13-34-39-371_00001.xml', mode 'r' at 0x64c1f40>
Last xml tag read: NMRC.FixedAtomSet.atoms
Parser state was: reading
Current object was: <ccp.nmr.NmrConstraint.FixedAtomSet [1, 1]>
CING started at : Mon Jan 31 21:46:34 2011
CING stopped at : Mon Jan 31 21:46:45 2011
CING took : 11.009 s
Traceback (most recent call last):
File "/Users/jd/workspace35/cingStable/python/cing/Scripts/validateEntry.py", line 355, in <module>
status = main(*sys.argv[1:])
File "/Users/jd/workspace35/cingStable/python/cing/Scripts/validateEntry.py", line 183, in main
if not project.initCcpn(ccpnFolder=fileNameTgz, modelCount=modelCount):
File "/Users/jd/workspace35/cingStable/python/cing/PluginCode/Ccpn.py", line 2332, in initCcpn
if not ccpn.importFromCcpn(modelCount = modelCount):
File "/Users/jd/workspace35/cingStable/python/cing/PluginCode/Ccpn.py", line 317, in importFromCcpn
if not self.importFromCcpnMolecule(modelCount = modelCount):
File "/Users/jd/workspace35/cingStable/python/cing/PluginCode/Ccpn.py", line 412, in importFromCcpnMolecule
ccpnConstraintLists.extend(constraintStore.sortedConstraintLists())
File "/Users/jd/workspace35/ccpn/python/ccp/api/nmr/NmrConstraint.py", line 22236, in sortedConstraintLists
self.load()
File "/Users/jd/workspace35/ccpn/python/memops/api/Implementation.py", line 5030, in load
self.loadFrom(repositories[0])
File "/Users/jd/workspace35/ccpn/python/memops/api/Implementation.py", line 5064, in loadFrom
XmlIO.loadTopObject(repository.url.getDataLocation(), self)
File "/Users/jd/workspace35/ccpn/python/memops/format/xml/XmlIO.py", line 236, in loadTopObject
return loadFromStream(open(filePath), topObjId=topObjId, topObject=topObject)
File "/Users/jd/workspace35/ccpn/python/memops/format/xml/XmlIO.py", line 267, in loadFromStream
partialLoad=partialLoad)
File "/Users/jd/workspace35/ccpn/python/memops/xml/Implementation.py", line 4912, in loadFromStream
value = clazz.getByKey(oo,vv[1:-1])
File "/Users/jd/workspace35/ccpn/python/ccp/api/molecule/MolSystem.py", line 1583, in getByKey
obj1.load()
File "/Users/jd/workspace35/ccpn/python/memops/api/Implementation.py", line 5030, in load
self.loadFrom(repositories[0])
File "/Users/jd/workspace35/ccpn/python/memops/api/Implementation.py", line 5064, in loadFrom
XmlIO.loadTopObject(repository.url.getDataLocation(), self)
File "/Users/jd/workspace35/ccpn/python/memops/format/xml/XmlIO.py", line 236, in loadTopObject
return loadFromStream(open(filePath), topObjId=topObjId, topObject=topObject)
File "/Users/jd/workspace35/ccpn/python/memops/format/xml/XmlIO.py", line 273, in loadFromStream
partialLoad=partialLoad)
File "/Users/jd/workspace35/ccpn/python/memops/xml/Implementation.py", line 4993, in loadFromStream
obj.checkValid()
File "/Users/jd/workspace35/ccpn/python/ccp/api/molecule/MolSystem.py", line 20447, in checkValid
('chemCompVars', (self.linking, self.descriptor))
File "/Users/jd/workspace35/ccpn/python/memops/api/Implementation.py", line 3974, in getByNavigation
result.load()
File "/Users/jd/workspace35/ccpn/python/memops/api/Implementation.py", line 5030, in load
self.loadFrom(repositories[0])
File "/Users/jd/workspace35/ccpn/python/memops/api/Implementation.py", line 5064, in loadFrom
XmlIO.loadTopObject(repository.url.getDataLocation(), self)
File "/Users/jd/workspace35/ccpn/python/memops/format/xml/XmlIO.py", line 236, in loadTopObject
return loadFromStream(open(filePath), topObjId=topObjId, topObject=topObject)
File "/Users/jd/workspace35/ccpn/python/memops/format/xml/XmlIO.py", line 273, in loadFromStream
partialLoad=partialLoad)
File "/Users/jd/workspace35/ccpn/python/memops/xml/Implementation.py", line 4993, in loadFromStream
obj.checkValid()
File "/Users/jd/workspace35/ccpn/python/ccp/api/molecule/ChemComp.py", line 26753, in checkValid
+ ": %s" % (self,)
memops.general.Implementation.ApiError:
ccp.molecule.ChemComp.NonStdChemComp.checkValid:
constraint protein_DNA_RNA_are_linear_polymers violated: <ccp.molecule.ChemComp.NonStdChemComp ['protein', 'Phq']>
Original comment by jurge...@gmail.com
on 1 Feb 2011 at 11:15
For the 134d problem, you are sure that this file:
/Volumes/tera4/NRG-CING/data/34/134d/134d/ccp/molecule/MolSystem/134d+134d_user_
2010-11-10-09-52-07-265_00001.xml
exists? It could be CCPN project confusion.
Original comment by wfvran...@gmail.com
on 2 Feb 2011 at 3:28
I've rerun this entry. Now it's log is at:
http://nmr.cmbi.ru.nl/NRG-CING/data/34/134d/log_validateEntry/134d_2011-02-02_05
-29-06.log
Same problem though, right?
Original comment by jurge...@gmail.com
on 2 Feb 2011 at 4:35
Indeed, this is a fun one. The sequence is 31 residues long, but not all
residues have coordinates, and the gaps are flanked by residues with a sequence
insertion code. Because I set the PDB/author info in CCPN (for NRG reasons), it
got a bit confused with the sequence insertion code and was not using this for
the 'missing' residues, so there was duplication of seqCode/insertionCode pairs
in the sequence (which is not allowed).
Anyway fixed now, again have to check this into SF, not sure when that'll be up
and running.
Original comment by wfvran...@gmail.com
on 3 Feb 2011 at 4:23
Great, let me know when you had a chance to commit. The main sf.net site is up
now.
Original comment by jurge...@gmail.com
on 4 Feb 2011 at 8:59
But cvs access is still blocked I feel:
https://sourceforge.net/blog/update-sourceforgenet-attack/#comments
Time to move over to Google Code's svn ? ;-)
Original comment by jurge...@gmail.com
on 4 Feb 2011 at 3:15
Seriously... I'll keep an eye on it, will do it asap.
Original comment by wfvran...@gmail.com
on 4 Feb 2011 at 3:19
Is now in SF CVS
Original comment by wfvran...@gmail.com
on 22 Feb 2011 at 1:47
Analysis got the newly created CCPN project files for entry 134d in fine.
Thanks Wim!
Original comment by jurge...@gmail.com
on 23 Feb 2011 at 3:05
I'm having to reopen this issue for e.g. entry 1lcc.
Take a look at:
http://nmr.cmbi.ru.nl/NRG-CING/data/lc/1lcc/log_validateEntry/1lcc_2011-03-04_05
-19-17.log
I checked the waters in the input DOCR STAR and they are fine.
Traceback (most recent call last):
File "/Users/jd/workspace35/cingStable/python/cing/Scripts/validateEntry.py", line 357, in <module>
status = main(*sys.argv[1:])
File "/Users/jd/workspace35/cingStable/python/cing/Scripts/validateEntry.py", line 185, in main
if not project.initCcpn(ccpnFolder=fileNameTgz, modelCount=modelCount):
File "/Users/jd/workspace35/cingStable/python/cing/PluginCode/Ccpn.py", line 2334, in initCcpn
if not ccpn.importFromCcpn(modelCount = modelCount):
File "/Users/jd/workspace35/cingStable/python/cing/PluginCode/Ccpn.py", line 310, in importFromCcpn
if not self.importFromCcpnMolecule(modelCount = modelCount):
File "/Users/jd/workspace35/cingStable/python/cing/PluginCode/Ccpn.py", line 422, in importFromCcpnMolecule
if not len(ccpnMolSys.structureEnsembles):
File "/Users/jd/workspace35/ccpn/python/ccp/api/molecule/MolSystem.py", line 13467, in getStructureEnsembles
self.load()
File "/Users/jd/workspace35/ccpn/python/memops/api/Implementation.py", line 5030, in load
self.loadFrom(repositories[0])
File "/Users/jd/workspace35/ccpn/python/memops/api/Implementation.py", line 5064, in loadFrom
XmlIO.loadTopObject(repository.url.getDataLocation(), self)
File "/Users/jd/workspace35/ccpn/python/memops/format/xml/XmlIO.py", line 236, in loadTopObject
return loadFromStream(open(filePath), topObjId=topObjId, topObject=topObject)
File "/Users/jd/workspace35/ccpn/python/memops/format/xml/XmlIO.py", line 267, in loadFromStream
partialLoad=partialLoad)
File "/Users/jd/workspace35/ccpn/python/memops/xml/Implementation.py", line 4993, in loadFromStream
obj.checkValid()
File "/Users/jd/workspace35/ccpn/python/ccp/api/molecule/MolSystem.py", line 3933, in checkValid
+ ": %s" % (self,)
memops.general.Implementation.ApiError: ccp.molecule.MolSystem.Chain.checkValid:
constraint seqCode_and_seqInsertCode_form_secondary_key_for_residues violated: <ccp.molecule.MolSystem.Chain ['1lcc', 'E']>
I suggest we test read the CCPN project at BMRB so the problem can be caught
early on.
We've seen this issue at NRG too:
http://code.google.com/p/nmrrestrntsgrid/issues/detail?id=159
I've checked manually if I'm on the most recent versions of the files and code.
Indeed so with below output.
================================================================================
======================
| CING: Common Interface for NMR structure Generation version 0.95 (r941)
AW,JFD,GWV 2004-2011 |
================================================================================
======================
User: jd on: stella.cmbi.umcn.nl (darwin/64bit/2cores/2.6.6) at:
(53346) Fri Mar 4 19:36:59 2011
DEBUG: Using program arguments:
DEBUG: inputDir:
file:///Library/WebServer/Documents/NRG-CING/input/lc
DEBUG: outputDir: /Library/WebServer/Documents/NRG-CING/data/lc/1lcc
DEBUG: pdbConvention: .
DEBUG: restraintsConvention: .
DEBUG: archiveType: BY_CH23
DEBUG: projectType: CCPN
DEBUG: storeCING2db: True
DEBUG: ranges: cv
DEBUG: filterTopViolations: 1
DEBUG:
DEBUG: Using derived settings:
DEBUG: modelCount: None
DEBUG: isRemoteOutputDir: False
copying file: /Library/WebServer/Documents/NRG-CING/input/lc/1lcc.tgz to:
1lcc.tgz
Moving CCPN directory from [linkNmrStarData] to [1lcc]
==> Importing data from Ccpn project "1lcc"
Error loading file for: <ccp.molecule.MolSystem.MolSystem ['1lcc']>
Reading: <open file
'/Library/WebServer/Documents/NRG-CING/data/lc/1lcc/1lcc/ccp/molecule/MolSystem/
1lcc+1lcc_user_2011-02-25-08-47-48-991_00001.xml', mode 'r' at 0x10853fe40>
Last xml tag read: _StorageUnit
Parser state was: checking validity
Object stack was empty
CING started at : Fri Mar 4 19:36:55 2011
CING stopped at : Fri Mar 4 19:37:00 2011
CING took : 5.088 s
Traceback (most recent call last):
File "/Users/jd/workspace35/cing/python/cing/Scripts/validateEntry.py", line 393, in <module>
status = main(*sys.argv[1:])
File "/Users/jd/workspace35/cing/python/cing/Scripts/validateEntry.py", line 233, in main
if not project.initCcpn(ccpnFolder=fileNameTgz, modelCount=modelCount):
File "/Users/jd/workspace35/cing/python/cing/PluginCode/Ccpn.py", line 2343, in initCcpn
if not ccpn.importFromCcpn(modelCount = modelCount):
File "/Users/jd/workspace35/cing/python/cing/PluginCode/Ccpn.py", line 310, in importFromCcpn
if not self.importFromCcpnMolecule(modelCount = modelCount):
File "/Users/jd/workspace35/cing/python/cing/PluginCode/Ccpn.py", line 422, in importFromCcpnMolecule
if not len(ccpnMolSys.structureEnsembles):
File "/Users/jd/workspace35/ccpn/python/ccp/api/molecule/MolSystem.py", line 13467, in getStructureEnsembles
self.load()
File "/Users/jd/workspace35/ccpn/python/memops/api/Implementation.py", line 5030, in load
self.loadFrom(repositories[0])
File "/Users/jd/workspace35/ccpn/python/memops/api/Implementation.py", line 5064, in loadFrom
XmlIO.loadTopObject(repository.url.getDataLocation(), self)
File "/Users/jd/workspace35/ccpn/python/memops/format/xml/XmlIO.py", line 236, in loadTopObject
return loadFromStream(open(filePath), topObjId=topObjId, topObject=topObject)
File "/Users/jd/workspace35/ccpn/python/memops/format/xml/XmlIO.py", line 267, in loadFromStream
partialLoad=partialLoad)
File "/Users/jd/workspace35/ccpn/python/memops/xml/Implementation.py", line 4993, in loadFromStream
obj.checkValid()
File "/Users/jd/workspace35/ccpn/python/ccp/api/molecule/MolSystem.py", line 3933, in checkValid
+ ": %s" % (self,)
memops.general.Implementation.ApiError: ccp.molecule.MolSystem.Chain.checkValid:
constraint seqCode_and_seqInsertCode_form_secondary_key_for_residues violated: <ccp.molecule.MolSystem.Chain ['1lcc', 'E']>
Original comment by jurge...@gmail.com
on 4 Mar 2011 at 6:39
There is a problem in the NMR-STAR file for 1lcc... in the molecule
description, for chain E, it reads:
5 5 HOH 1 HOH 1 1
5 5 HOH 2 HOH 1 1
5 5 HOH 3 HOH 1 1
5 5 HOH 4 HOH 1 1
5 5 HOH 5 HOH 1 1
5 5 HOH 6 HOH 1 1
5 5 HOH 7 HOH 1 1
5 5 HOH 8 HOH 1 1
5 5 HOH 9 HOH 1 1
5 5 HOH 10 HOH 1 1
5 5 HOH 11 HOH 1 1
5 5 HOH 12 HOH 1 1
5 5 HOH 13 HOH 1 1
5 5 HOH 14 HOH 1 1
5 5 HOH 15 HOH 1 1
5 5 HOH 16 HOH 1 1
5 5 HOH 17 HOH 1 1
5 5 HOH 18 HOH 1 1
5 5 HOH 19 HOH 1 1
5 5 HOH 20 HOH 1 1
5 5 HOH 21 HOH 1 1
5 5 HOH 22 HOH 1 1
5 5 HOH 23 HOH 1 1
5 5 HOH 24 HOH 1 1
5 5 HOH 25 HOH 1 1
5 5 HOH 26 HOH 1 1
5 5 HOH 27 HOH 1 1
5 5 HOH 28 HOH 1 1
5 5 HOH 29 HOH 1 1
5 5 HOH 30 HOH 1 1
5 5 HOH 31 HOH 1 1
5 5 HOH 32 HOH 1 1
5 5 HOH 33 HOH 1 1
5 5 HOH 34 HOH 1 1
5 5 HOH 35 HOH 1 1
5 5 HOH 36 HOH 1 1
5 5 HOH 37 HOH 1 1
5 5 HOH 38 HOH 1 1
5 5 HOH 39 HOH 1 1
5 5 HOH 40 HOH 1 1
5 5 HOH 41 HOH 1 1
5 5 HOH 42 HOH 1 1
5 5 HOH 43 HOH 1 1
5 5 HOH 44 HOH 1 1
5 5 HOH 45 HOH 1 1
5 5 HOH 46 HOH 1 1
5 5 HOH 47 HOH 1 1
5 5 HOH 48 HOH 1 1
5 5 HOH 49 HOH 1 1
5 5 HOH 50 HOH 1 1
5 5 HOH 51 HOH 1 1
5 5 HOH 52 HOH 1 1
5 5 HOH 53 HOH 1 1
In the coordinates:
1 991 5 5 1 HOH O O 21.970 18.650 18.350 1.00 0.0 E A . 52 HOH O 1 1
1 992 5 5 1 HOH H1 H 22.400 19.540 18.450 1.00 0.0 E A . 52 HOH H1 1 1
1 993 5 5 1 HOH H2 H 22.170 18.090 19.160 1.00 0.0 E A . 52 HOH H2 1 1
1 994 5 5 2 HOH O O 9.910 20.180 29.550 1.00 0.0 E A . 53 HOH O 1 1
1 995 5 5 2 HOH H1 H 9.650 20.460 28.620 1.00 0.0 E A . 53 HOH H1 1 1
1 996 5 5 2 HOH H2 H 9.310 19.440 29.850 1.00 0.0 E A . 53 HOH H2 1 1
1 997 5 5 3 HOH O O 10.690 19.540 32.470 1.00 0.0 E A . 54 HOH O 1 1
1 998 5 5 3 HOH H1 H 10.760 18.700 31.940 1.00 0.0 E A . 54 HOH H1 1 1
1 999 5 5 3 HOH H2 H 10.340 19.330 33.390 1.00 0.0 E A . 54 HOH H2 1 1
1 1000 5 5 4 HOH O O 11.960 21.630 33.860 1.00 0.0 E A . 55 HOH O 1 1
1 1001 5 5 4 HOH H1 H 11.510 20.830 33.460 1.00 0.0 E A . 55 HOH H1 1 1
1 1002 5 5 4 HOH H2 H 11.500 21.880 34.710 1.00 0.0 E A . 55 HOH H2 1 1
1 1003 5 5 5 HOH O O 17.850 27.380 32.270 1.00 0.0 E A . 56 HOH O 1 1
1 1004 5 5 5 HOH H1 H 18.730 27.770 32.560 1.00 0.0 E A . 56 HOH H1 1 1
1 1005 5 5 5 HOH H2 H 17.230 28.110 32.000 1.00 0.0 E A . 56 HOH H2 1 1
1 1006 5 5 6 HOH O O 15.790 23.850 27.410 1.00 0.0 E A . 57 HOH O 1 1
1 1007 5 5 6 HOH H1 H 16.270 24.270 26.650 1.00 0.0 E A . 57 HOH H1 1 1
1 1008 5 5 6 HOH H2 H 14.880 23.560 27.110 1.00 0.0 E A . 57 HOH H2 1 1
1 1009 5 5 7 HOH O O 18.180 22.780 17.970 1.00 0.0 E A . 58 HOH O 1 1
1 1010 5 5 7 HOH H1 H 18.410 22.090 17.280 1.00 0.0 E A . 58 HOH H1 1 1
1 1011 5 5 7 HOH H2 H 18.400 22.430 18.880 1.00 0.0 E A . 58 HOH H2 1 1
1 1012 5 5 8 HOH O O 17.160 34.420 31.500 1.00 0.0 E A . 59 HOH O 1 1
1 1013 5 5 8 HOH H1 H 17.730 34.030 32.220 1.00 0.0 E A . 59 HOH H1 1 1
1 1014 5 5 8 HOH H2 H 16.880 33.690 30.880 1.00 0.0 E A . 59 HOH H2 1 1
1 1015 5 5 9 HOH O O 18.840 24.980 20.090 1.00 0.0 E A . 60 HOH O 1 1
1 1016 5 5 9 HOH H1 H 18.820 24.950 21.090 1.00 0.0 E A . 60 HOH H1 1 1
1 1017 5 5 9 HOH H2 H 19.650 25.470 19.790 1.00 0.0 E A . 60 HOH H2 1 1
1 1018 5 5 10 HOH O O 16.080 29.580 31.530 1.00 0.0 E A . 61 HOH O 1 1
1 1019 5 5 10 HOH H1 H 15.120 29.300 31.470 1.00 0.0 E A . 61 HOH H1 1 1
1 1020 5 5 10 HOH H2 H 16.220 30.090 32.380 1.00 0.0 E A . 61 HOH H2 1 1
1 1021 5 5 11 HOH O O 15.660 34.740 34.520 1.00 0.0 E A . 62 HOH O 1 1
1 1022 5 5 11 HOH H1 H 15.430 34.240 33.680 1.00 0.0 E A . 62 HOH H1 1 1
1 1023 5 5 11 HOH H2 H 15.980 34.090 35.210 1.00 0.0 E A . 62 HOH H2 1 1
1 1024 5 5 12 HOH O O 11.380 34.520 34.630 1.00 0.0 E A . 63 HOH O 1 1
1 1025 5 5 12 HOH H1 H 11.530 35.510 34.590 1.00 0.0 E A . 63 HOH H1 1 1
1 1026 5 5 12 HOH H2 H 11.640 34.190 35.540 1.00 0.0 E A . 63 HOH H2 1 1
1 1027 5 5 13 HOH O O 17.190 19.580 18.260 1.00 0.0 E A . 64 HOH O 1 1
1 1028 5 5 13 HOH H1 H 16.840 18.840 18.840 1.00 0.0 E A . 64 HOH H1 1 1
1 1029 5 5 13 HOH H2 H 16.530 20.340 18.260 1.00 0.0 E A . 64 HOH H2 1 1
1 1030 5 5 14 HOH O O 11.670 25.730 37.590 1.00 0.0 E A . 65 HOH O 1 1
1 1031 5 5 14 HOH H1 H 11.900 25.290 38.460 1.00 0.0 E A . 65 HOH H1 1 1
1 1032 5 5 14 HOH H2 H 10.990 25.160 37.110 1.00 0.0 E A . 65 HOH H2 1 1
1 1033 5 5 15 HOH O O 32.300 19.840 16.170 1.00 0.0 E A . 66 HOH O 1 1
1 1034 5 5 15 HOH H1 H 32.680 19.150 15.550 1.00 0.0 E A . 66 HOH H1 1 1
1 1035 5 5 15 HOH H2 H 32.870 19.910 16.990 1.00 0.0 E A . 66 HOH H2 1 1
1 1036 5 5 16 HOH O O 25.210 18.010 15.370 1.00 0.0 E A . 67 HOH O 1 1
1 1037 5 5 16 HOH H1 H 25.390 17.060 15.120 1.00 0.0 E A . 67 HOH H1 1 1
1 1038 5 5 16 HOH H2 H 24.260 18.090 15.680 1.00 0.0 E A . 67 HOH H2 1 1
1 1039 5 5 17 HOH O O 28.250 18.350 25.940 1.00 0.0 E A . 68 HOH O 1 1
1 1040 5 5 17 HOH H1 H 28.900 18.440 25.190 1.00 0.0 E A . 68 HOH H1 1 1
1 1041 5 5 17 HOH H2 H 28.200 17.390 26.240 1.00 0.0 E A . 68 HOH H2 1 1
1 1042 5 5 18 HOH O O 25.910 19.390 26.730 1.00 0.0 E A . 69 HOH O 1 1
1 1043 5 5 18 HOH H1 H 25.640 20.100 26.070 1.00 0.0 E A . 69 HOH H1 1 1
1 1044 5 5 18 HOH H2 H 26.720 18.910 26.390 1.00 0.0 E A . 69 HOH H2 1 1
1 1045 5 5 19 HOH O O 29.690 24.270 21.950 1.00 0.0 E A . 70 HOH O 1 1
1 1046 5 5 19 HOH H1 H 29.730 23.490 21.320 1.00 0.0 E A . 70 HOH H1 1 1
1 1047 5 5 19 HOH H2 H 29.130 24.990 21.550 1.00 0.0 E A . 70 HOH H2 1 1
1 1048 5 5 20 HOH O O 25.850 19.540 29.360 1.00 0.0 E A . 71 HOH O 1 1
1 1049 5 5 20 HOH H1 H 26.170 19.600 28.410 1.00 0.0 E A . 71 HOH H1 1 1
1 1050 5 5 20 HOH H2 H 25.730 18.580 29.610 1.00 0.0 E A . 71 HOH H2 1 1
1 1051 5 5 21 HOH O O 14.200 26.840 33.980 1.00 0.0 E A . 72 HOH O 1 1
1 1052 5 5 21 HOH H1 H 14.540 27.700 33.610 1.00 0.0 E A . 72 HOH H1 1 1
1 1053 5 5 21 HOH H2 H 13.230 26.740 33.780 1.00 0.0 E A . 72 HOH H2 1 1
1 1054 5 5 22 HOH O O 31.320 25.660 18.650 1.00 0.0 E A . 73 HOH O 1 1
1 1055 5 5 22 HOH H1 H 30.700 24.870 18.550 1.00 0.0 E A . 73 HOH H1 1 1
1 1056 5 5 22 HOH H2 H 30.770 26.480 18.820 1.00 0.0 E A . 73 HOH H2 1 1
1 1057 5 5 23 HOH O O 23.160 19.700 30.160 1.00 0.0 E A . 74 HOH O 1 1
1 1058 5 5 23 HOH H1 H 23.560 19.520 31.060 1.00 0.0 E A . 74 HOH H1 1 1
1 1059 5 5 23 HOH H2 H 23.880 19.660 29.470 1.00 0.0 E A . 74 HOH H2 1 1
1 1060 5 5 24 HOH O O 24.270 26.570 33.610 1.00 0.0 E A . 75 HOH O 1 1
1 1061 5 5 24 HOH H1 H 24.740 26.890 34.430 1.00 0.0 E A . 75 HOH H1 1 1
1 1062 5 5 24 HOH H2 H 24.810 25.840 33.190 1.00 0.0 E A . 75 HOH H2 1 1
1 1063 5 5 25 HOH O O 25.990 24.680 32.760 1.00 0.0 E A . 76 HOH O 1 1
1 1064 5 5 25 HOH H1 H 25.810 23.710 32.940 1.00 0.0 E A . 76 HOH H1 1 1
1 1065 5 5 25 HOH H2 H 26.940 24.800 32.480 1.00 0.0 E A . 76 HOH H2 1 1
1 1066 5 5 26 HOH O O 25.430 20.950 24.590 1.00 0.0 E A . 77 HOH O 1 1
1 1067 5 5 26 HOH H1 H 26.140 20.910 23.890 1.00 0.0 E A . 77 HOH H1 1 1
1 1068 5 5 26 HOH H2 H 24.550 21.170 24.170 1.00 0.0 E A . 77 HOH H2 1 1
1 1069 5 5 27 HOH O O 13.620 38.770 37.220 1.00 0.0 E A . 78 HOH O 1 1
1 1070 5 5 27 HOH H1 H 13.320 38.550 38.150 1.00 0.0 E A . 78 HOH H1 1 1
1 1071 5 5 27 HOH H2 H 14.610 38.650 37.170 1.00 0.0 E A . 78 HOH H2 1 1
So a lot less waters for entity_assembly_ID 5 in the coordinates than in the
molecule description. Because I reset the water sequence codes (on request!) to
what's in the coordinates, there is a problem here and some sequence codes end
up overlapping.
Original comment by wfvran...@gmail.com
on 8 Mar 2011 at 9:25
I looked some more at this... I introduced a bit of code a while ago to deal
with water molecules differently on reading the sequence from NMR-STAR, but
apparently it's changed again now in the joinedCoord files, is that correct?
Everything seems to work if I comment out that bit. I can check in this change,
but it is possible it will give problems with other entries (especially if the
water is handled differently in the coordinate NMR-STAR file of course).
Original comment by wfvran...@gmail.com
on 9 Mar 2011 at 10:04
This issue is fixed for entries:
134d
135d
136d
177d
1gnc
1msh
2neo
2ofg
I'm tagging this issue for 1lcc where it still does apply today:
http://nmr.cmbi.ru.nl/NRG-CING/data/lc/1lcc/log_validateEntry/1lcc_2011-03-28_20
-08-19.log
Let me look into 1lcc now.
Original comment by jurge...@gmail.com
on 30 Mar 2011 at 10:15
Note, I'll need to check this fix for all 35 water involved entries:
1ai0 1aiy 1cb9 1crp 1e08 1e0a 1ffj 1gxx 1hum 1lcc 1lcd 1mdi 1oa5 1oa6 1ocd 1ppx
1pun 1puq 1pus 1r2n 1r84 1rl5 1saf 1tum 1yn2 1z6c 1zrr 2frc 2hji 2hz8 2ke5 2kkh
2rmk 3usn 6i1b
They can be retrieved from NRG-CING simply by:
sort -n -t, -k3 nrgcing.cingresidue.csv > nrgcing.cingresidue.srt.csv
sort -n -t, -k1 nrgcing.cingentry.csv > nrgcing.cingentry.srt.csv
jd:nmr/~/tmp/ join -t, -2 3 -o 1.2,2.5 nrgcing.cingentry.srt.csv
nrgcing.cingresidue.srt.csv | grep HOH | sort -u
1ai0,HOH
1aiy,HOH
1cb9,HOH
1crp,HOH
1e08,HOH
1e0a,HOH
1ffj,HOH
1gxx,HOH
1hum,HOH
1lcc,HOH
1lcd,HOH
1mdi,HOH
1oa5,HOH
1oa6,HOH
1ocd,HOH
1ppx,HOH
1pun,HOH
1puq,HOH
1pus,HOH
1r2n,HOH
1r84,HOH
1rl5,HOH
1saf,HOH
1tum,HOH
1yn2,HOH
1z6c,HOH
1zrr,HOH
2frc,HOH
2hji,HOH
2hz8,HOH
2ke5,HOH
2kkh,HOH
2rmk,HOH
3usn,HOH
6i1b,HOH
Original comment by jurge...@gmail.com
on 30 Mar 2011 at 11:03
Today after 1lcc was updated by Chris and NRG-CING retried it still fails with:
constraint seqCode_and_seqInsertCode_form_secondary_key_for_residues violated:
My NRG-CING code can handle absent CCPN projects but not invalid projects. I
suggest we check the project before accepting it in NRG. Otherwise the
development cycles take too long and are too complex as well.
Wim, can you fix the actual problem or do you need input from us?
Original comment by jurge...@gmail.com
on 8 Apr 2011 at 2:22
To be honest. This is an NRG problem more than CING's.
Please see: NRG isssssue 272: Not accept invalid CCPN project anymore at:
http://code.google.com/p/nmrrestrntsgrid/issues/detail?id=272
Original comment by jurge...@gmail.com
on 11 Apr 2011 at 12:18
OK see my comment 15 - on my computer 1lcc works fine because I commented out a
bit of code that special-cases waters (the problem in 1lcc). That bit of code
was required before but not any more, I think, because we treat waters
alltogether differently now anyway... of course it is possible that other
entries will fail if I do take it out, I cannot guarantee that that won't
happen.
Let me know what you want me to do with this one.
Original comment by wfvran...@gmail.com
on 11 Apr 2011 at 12:25
If you have resources to fix, validate and check it, great. If not, let's not
do this halfway and break other parts. We'll have to accept the fact that FC
can't provide valid CCPN projects for all NRG entries. I do want NRG to realize
that though and not put it out in those cases but that's an NRG issue. The
issue here is closed I should think.
Original comment by jurge...@gmail.com
on 11 Apr 2011 at 12:39
Addressed in r968 but no real fix expected soon.
Original comment by jurge...@gmail.com
on 13 Apr 2011 at 10:32
Entries that have this problem: 177d 1gnc 1lcc 1lcd 1qch 1sae 1sah 1saj 1sak
1sal 3sak
They'll not make it into NRG-CING until fixed.
Original comment by jurge...@gmail.com
on 10 Oct 2011 at 8:15
Original issue reported on code.google.com by
jurge...@gmail.com
on 10 Nov 2010 at 1:08