Opened 4 years ago
Last modified 4 years ago
#5274 assigned defect
Modeller: Number of residues in the alignment and pdb files are different
| Reported by: | Elaine Meng | Owned by: | Eric Pettersen |
|---|---|---|---|
| Priority: | normal | Milestone: | |
| Component: | Structure Prediction | Version: | |
| Keywords: | Cc: | ||
| Blocked By: | Blocking: | ||
| Notify when closed: | Platform: | all | |
| Project: | ChimeraX |
Description
The following bug report has been submitted:
Platform: macOS-10.15.7-x86_64-i386-64bit
ChimeraX Version: 1.3.dev202109212024 (2021-09-21 20:24:53 UTC)
Description
open 1maz
open uniprot:bclx_human
sequence associate /A bclx_human:1
modeller refine bclx_human:1:all-missing numModels 5 fast false adjacentFlexible 1 protocol standard
Log:
> open /Users/meng/Desktop/startup.cxc
> alias reset view orient; view initial
> alias start tool show $1
> alias whereprefs info path user unversioned config
> alias captut open help:user/tutorials/binding-sites.html#cap-example
> alias previewts toolshed url https://cxtoolshed-
> preview.rbvi.ucsf.edu;toolshed reload available
> alias normalts toolshed url https://cxtoolshed.rbvi.ucsf.edu;toolshed reload
> available
> alias btut open
> https://www.cgl.ucsf.edu/home/meng/chimerax/vdocs/user/tutorials/binding-
> sites.html
> alias stut open https://www.rbvi.ucsf.edu/chimerax/data/conservation-
> coloring/conservation-coloring.html
executed startup.cxc
UCSF ChimeraX version: 1.3.dev202109212024 (2021-09-21)
© 2016-2021 Regents of the University of California. All rights reserved.
How to cite UCSF ChimeraX
> open 1maz
1maz title:
X-ray structure of BCL-XL, an inhibitor of programmed cell death [more
info...]
Chain information for 1maz #1
---
Chain | Description | UniProt
A | Bcl-2-like protein 1 | BCLX_HUMAN
> open uniprot:bclx_human
Summary of feedback from opening bclx_human fetched from uniprot
---
notes | UniProt identifier bclx_human maps to entry Q07817
Alignment identifier is bclx_human
Opened UniProt bclx_human
> sequence associate /A bclx_human:1
Associated 1maz chain A to bclx_human with 12 mismatches and/or gaps
> ui tool show "Model Loops"
> modeller refine bclx_human:1:all-missing numModels 5 fast false
> adjacentFlexible 1 protocol standard
Web Service: Modeller9v8 is a Python wrapper that calls Modeller (v10.1) for
protein structure modeling
Opal service URL:
http://webservices.rbvi.ucsf.edu/opal2/services/Modeller9v8Service
Opal job id: appModeller9v8Service16323456707651186954005
Opal status URL prefix:
http://webservices.rbvi.ucsf.edu/appModeller9v8Service16323456707651186954005
stdout.txt = standard output
stderr.txt = standard error
Modeller job ID appModeller9v8Service16323456707651186954005 finished
Modeller error output
Traceback (most recent call last):
File "/usr/local/opal-local/bin/modeller9v8.py", line 332, in <module>
main()
File "/usr/local/opal-local/bin/modeller9v8.py", line 24, in main
VersionMap[cf["version"]](cf)
File "/usr/local/opal-local/bin/modeller9v8.py", line 34, in v2_run
execfile(fn)
File "ModellerModelling.py", line 88, in <module>
a.make()
File "/usr/lib64/python2.7/site-packages/modeller/automodel/loopmodel.py",
line 42, in make
AutoModel.make(self, exit_stage)
File "/usr/lib64/python2.7/site-packages/modeller/automodel/automodel.py",
line 141, in make
self.homcsr(exit_stage)
File "/usr/lib64/python2.7/site-packages/modeller/automodel/automodel.py",
line 625, in homcsr
self.check_alignment(aln)
File "/usr/lib64/python2.7/site-packages/modeller/automodel/automodel.py",
line 578, in check_alignment
aln.check()
File "/usr/lib64/python2.7/site-packages/modeller/alignment.py", line 213, in
check
self.check_structure_structure(io=io)
File "/usr/lib64/python2.7/site-packages/modeller/alignment.py", line 222, in
check_structure_structure
return f(self.modpt, io.modpt, self.env.libs.modpt, eqvdst)
_modeller.ModellerError: read_te_290E> Number of residues in the alignment and
pdb files are different: 143 291 For alignment entry: 1 1maz_1
Modeller run output
**** IGNORE FOLLOWING LICENSE KEY MESSAGE ****
check_lice_E> Invalid license key: xxx
Go to https://salilab.org/modeller/ to get a license key,
and then set the 'license' variable to it in the file
/usr/lib/modeller10.1/modlib/modeller/config.py
**** END IGNORED MESSAGE ****
MODELLER 10.1, 2021/03/12, r12156
PROTEIN STRUCTURE MODELLING BY SATISFACTION OF SPATIAL RESTRAINTS
Copyright(c) 1989-2021 Andrej Sali
All Rights Reserved
Written by A. Sali
with help from
B. Webb, M.S. Madhusudhan, M-Y. Shen, G.Q. Dong,
M.A. Marti-Renom, N. Eswar, F. Alber, M. Topf, B. Oliva,
A. Fiser, R. Sanchez, B. Yerkovich, A. Badretdinov,
F. Melo, J.P. Overington, E. Feyfant
University of California, San Francisco, USA
Rockefeller University, New York, USA
Harvard University, Cambridge, USA
Imperial Cancer Research Fund, London, UK
Birkbeck College, University of London, London, UK
Kind, OS, HostName, Kernel, Processor: 4, Linux franklin.cgl.ucsf.edu
3.10.0-1160.36.2.el7.x86_64 x86_64
Date and time of compilation : 2021/03/12 00:18:43
MODELLER executable type : x86_64-intel8
Job starting time (YY/MM/DD HH:MM:SS): 2021/09/22 14:21:10
environ____W> The class 'environ' is deprecated; use 'Environ' instead
openf___224_> Open $(LIB)/restyp.lib
openf___224_> Open ${MODINSTALL10v1}/modlib/resgrp.lib
rdresgr_266_> Number of residue groups: 2
openf___224_> Open ${MODINSTALL10v1}/modlib/sstruc.lib
Dynamically allocated memory at amaxlibraries [B,KiB,MiB]: 191566 187.076
0.183
Dynamically allocated memory at amaxlibraries [B,KiB,MiB]: 192094 187.592
0.183
openf___224_> Open ${MODINSTALL10v1}/modlib/resdih.lib
Dynamically allocated memory at amaxlibraries [B,KiB,MiB]: 240694 235.053
0.230
rdrdih__263_> Number of dihedral angle types : 9
Maximal number of dihedral angle optima: 3
Dihedral angle names : Alph Phi Psi Omeg chi1 chi2 chi3 chi4 chi5
openf___224_> Open ${MODINSTALL10v1}/modlib/radii.lib
Dynamically allocated memory at amaxlibraries [B,KiB,MiB]: 253994 248.041
0.242
openf___224_> Open ${MODINSTALL10v1}/modlib/af_mnchdef.lib
rdwilmo_274_> Mainchain residue conformation classes: APBLE
openf___224_> Open ${MODINSTALL10v1}/modlib/mnch.lib
rdclass_257_> Number of classes: 5
openf___224_> Open ${MODINSTALL10v1}/modlib/mnch1.lib
openf___224_> Open ${MODINSTALL10v1}/modlib/mnch2.lib
openf___224_> Open ${MODINSTALL10v1}/modlib/mnch3.lib
openf___224_> Open ${MODINSTALL10v1}/modlib/xs4.mat
rdrrwgh_268_> Number of residue types: 21
loopmodel__W> The class 'loopmodel' is deprecated; use 'LoopModel' instead
openf___224_> Open alignment.ali
Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 282617 275.993
0.270
Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 297709 290.731
0.284
Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 313109 305.771
0.299
Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 343909 335.849
0.328
Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 405509 396.005
0.387
Dynamically allocated memory at amaxsequence [B,KiB,MiB]: 406077 396.560 0.387
Dynamically allocated memory at amaxsequence [B,KiB,MiB]: 406957 397.419 0.388
Read the alignment from file : alignment.ali
Total number of alignment positions: 221
# Code #_Res #_Segm PDB_code Name
\-------------------------------------------------------------------------------
1 1maz_1 143 1 1maz_1
2 bclx_huma 221 1 bclx_human
check_a_343_> >> BEGINNING OF COMMAND
openf___224_> Open ./1maz_1.pdb
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 451582 440.998
0.431
Dynamically allocated memory at amaxsequence [B,KiB,MiB]: 451582 440.998 0.431
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 451753 441.165
0.431
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 452603 441.995
0.432
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 453861 443.224
0.433
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 455765 445.083
0.435
Dynamically allocated memory at amaxsequence [B,KiB,MiB]: 455765 445.083 0.435
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 455855 445.171
0.435
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 458711 447.960
0.437
Dynamically allocated memory at amaxsequence [B,KiB,MiB]: 458711 447.960 0.437
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 458846 448.092
0.438
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 463130 452.275
0.442
Dynamically allocated memory at amaxsequence [B,KiB,MiB]: 463130 452.275 0.442
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 463328 452.469
0.442
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 469754 458.744
0.448
Dynamically allocated memory at amaxsequence [B,KiB,MiB]: 469754 458.744 0.448
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 470051 459.034
0.448
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 479673 468.431
0.457
Dynamically allocated memory at amaxsequence [B,KiB,MiB]: 479673 468.431 0.457
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 480123 468.870
0.458
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 494573 482.981
0.472
Dynamically allocated memory at amaxsequence [B,KiB,MiB]: 494573 482.981 0.472
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 495248 483.641
0.472
Dynamically allocated memory at amaxsequence [B,KiB,MiB]: 495248 483.641 0.472
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 496256 484.625
0.473
Dynamically allocated memory at amaxcoordinates [B,KiB,MiB]: 517914 505.775
0.494
Dynamically allocated memory at amaxstructure [B,KiB,MiB]: 562492 549.309
0.536
read_te_290E> Number of residues in the alignment and pdb files are different:
143 291
For alignment entry: 1 1maz_1
x (mismatch at alignment position 144)
Alignment HLEPWIQENGGWDTFVELYG/
PDB HLEPWIQENGGWDTFVELYG-wwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwwww
Match ********************
Please check your alignment file header to be sure you correctly specified
the starting and ending residue numbers and chains. The alignment sequence
must match that from the atom file exactly.
Another possibility is that some residues in the atom file are missing,
perhaps because they could not be resolved experimentally. (Note that Modeller
reads only the ATOM and HETATM records in PDB, NOT the SEQRES records.)
In this case, simply replace the section of your alignment corresponding
to these missing residues with gaps.
read_te_288W> Protein not accepted: 1 1maz_1
No output models from Modeller; see log for Modeller text output.
OpenGL version: 4.1 INTEL-14.7.18
OpenGL renderer: Intel Iris Pro OpenGL Engine
OpenGL vendor: Intel Inc.Hardware:
Hardware Overview:
Model Name: MacBook Pro
Model Identifier: MacBookPro11,4
Processor Name: Quad-Core Intel Core i7
Processor Speed: 2.2 GHz
Number of Processors: 1
Total Number of Cores: 4
L2 Cache (per Core): 256 KB
L3 Cache: 6 MB
Hyper-Threading Technology: Enabled
Memory: 16 GB
Boot ROM Version: 427.140.8.0.0
SMC Version (system): 2.29f24
Software:
System Software Overview:
System Version: macOS 10.15.7 (19H1417)
Kernel Version: Darwin 19.6.0
Time since boot: 6:29
Graphics/Displays:
Intel Iris Pro:
Chipset Model: Intel Iris Pro
Type: GPU
Bus: Built-In
VRAM (Dynamic, Max): 1536 MB
Vendor: Intel
Device ID: 0x0d26
Revision ID: 0x0008
Metal: Supported, feature set macOS GPUFamily1 v4
Displays:
Color LCD:
Display Type: Built-In Retina LCD
Resolution: 2880 x 1800 Retina
Framebuffer Depth: 24-Bit Color (ARGB8888)
Main Display: Yes
Mirror: Off
Online: Yes
Automatically Adjust Brightness: No
Connection Type: Internal
Locale: (None, 'UTF-8')
PyQt5 5.15.2, Qt 5.15.2
Installed Packages:
alabaster: 0.7.12
appdirs: 1.4.4
appnope: 0.1.2
Babel: 2.9.1
backcall: 0.2.0
blockdiag: 2.0.1
certifi: 2021.5.30
cftime: 1.5.0
charset-normalizer: 2.0.6
ChimeraX-AddCharge: 1.1.4
ChimeraX-AddH: 2.1.10
ChimeraX-AlignmentAlgorithms: 2.0
ChimeraX-AlignmentHdrs: 3.2
ChimeraX-AlignmentMatrices: 2.0
ChimeraX-Alignments: 2.2.2
ChimeraX-AlphaFold: 1.0
ChimeraX-AltlocExplorer: 1.0.1
ChimeraX-AmberInfo: 1.0
ChimeraX-Arrays: 1.0
ChimeraX-Atomic: 1.29.1
ChimeraX-AtomicLibrary: 4.1.4
ChimeraX-AtomSearch: 2.0
ChimeraX-AtomSearchLibrary: 1.0
ChimeraX-AxesPlanes: 2.0
ChimeraX-BasicActions: 1.1
ChimeraX-BILD: 1.0
ChimeraX-BlastProtein: 2.0
ChimeraX-BondRot: 2.0
ChimeraX-BugReporter: 1.0
ChimeraX-BuildStructure: 2.6
ChimeraX-Bumps: 1.0
ChimeraX-BundleBuilder: 1.1
ChimeraX-ButtonPanel: 1.0
ChimeraX-CageBuilder: 1.0
ChimeraX-CellPack: 1.0
ChimeraX-Centroids: 1.2
ChimeraX-ChemGroup: 2.0
ChimeraX-Clashes: 2.1
ChimeraX-ColorActions: 1.0
ChimeraX-ColorGlobe: 1.0
ChimeraX-ColorKey: 1.5
ChimeraX-CommandLine: 1.1.5
ChimeraX-ConnectStructure: 2.0
ChimeraX-Contacts: 1.0
ChimeraX-Core: 1.3.dev202109212024
ChimeraX-CoreFormats: 1.1
ChimeraX-coulombic: 1.3.1
ChimeraX-Crosslinks: 1.0
ChimeraX-Crystal: 1.0
ChimeraX-CrystalContacts: 1.0
ChimeraX-DataFormats: 1.2.1
ChimeraX-Dicom: 1.0
ChimeraX-DistMonitor: 1.1.5
ChimeraX-DistUI: 1.0
ChimeraX-Dssp: 2.0
ChimeraX-EMDB-SFF: 1.0
ChimeraX-ExperimentalCommands: 1.0
ChimeraX-FileHistory: 1.0
ChimeraX-FunctionKey: 1.0
ChimeraX-Geometry: 1.1
ChimeraX-gltf: 1.0
ChimeraX-Graphics: 1.1
ChimeraX-Hbonds: 2.1.1
ChimeraX-Help: 1.2
ChimeraX-HKCage: 1.3
ChimeraX-IHM: 1.1
ChimeraX-ImageFormats: 1.2
ChimeraX-IMOD: 1.0
ChimeraX-IO: 1.0.1
ChimeraX-ItemsInspection: 1.0
ChimeraX-Label: 1.1
ChimeraX-ListInfo: 1.1.1
ChimeraX-Log: 1.1.4
ChimeraX-LookingGlass: 1.1
ChimeraX-Maestro: 1.8.1
ChimeraX-Map: 1.1
ChimeraX-MapData: 2.0
ChimeraX-MapEraser: 1.0
ChimeraX-MapFilter: 2.0
ChimeraX-MapFit: 2.0
ChimeraX-MapSeries: 2.1
ChimeraX-Markers: 1.0
ChimeraX-Mask: 1.0
ChimeraX-MatchMaker: 2.0.1
ChimeraX-MDcrds: 2.6
ChimeraX-MedicalToolbar: 1.0.1
ChimeraX-Meeting: 1.0
ChimeraX-MLP: 1.1
ChimeraX-mmCIF: 2.4
ChimeraX-MMTF: 2.1
ChimeraX-Modeller: 1.2
ChimeraX-ModelPanel: 1.2
ChimeraX-ModelSeries: 1.0
ChimeraX-Mol2: 2.0
ChimeraX-Morph: 1.0
ChimeraX-MouseModes: 1.1
ChimeraX-Movie: 1.0
ChimeraX-Neuron: 1.0
ChimeraX-Nucleotides: 2.0.2
ChimeraX-OpenCommand: 1.7
ChimeraX-PDB: 2.6.2
ChimeraX-PDBBio: 1.0
ChimeraX-PDBLibrary: 1.0.1
ChimeraX-PDBMatrices: 1.0
ChimeraX-PickBlobs: 1.0
ChimeraX-Positions: 1.0
ChimeraX-PresetMgr: 1.0.1
ChimeraX-PubChem: 2.1
ChimeraX-ReadPbonds: 1.0
ChimeraX-Registration: 1.1
ChimeraX-RemoteControl: 1.0
ChimeraX-ResidueFit: 1.0
ChimeraX-RestServer: 1.1
ChimeraX-RNALayout: 1.0
ChimeraX-RotamerLibMgr: 2.0
ChimeraX-RotamerLibsDunbrack: 2.0
ChimeraX-RotamerLibsDynameomics: 2.0
ChimeraX-RotamerLibsRichardson: 2.0
ChimeraX-SaveCommand: 1.5
ChimeraX-SchemeMgr: 1.0
ChimeraX-SDF: 2.0
ChimeraX-Segger: 1.0
ChimeraX-Segment: 1.0
ChimeraX-SelInspector: 1.0
ChimeraX-SeqView: 2.4.3
ChimeraX-Shape: 1.0.1
ChimeraX-Shell: 1.0
ChimeraX-Shortcuts: 1.1
ChimeraX-ShowAttr: 1.0
ChimeraX-ShowSequences: 1.0
ChimeraX-SideView: 1.0
ChimeraX-Smiles: 2.1
ChimeraX-SmoothLines: 1.0
ChimeraX-SpaceNavigator: 1.0
ChimeraX-StdCommands: 1.6
ChimeraX-STL: 1.0
ChimeraX-Storm: 1.0
ChimeraX-Struts: 1.0
ChimeraX-Surface: 1.0
ChimeraX-SwapAA: 2.0
ChimeraX-SwapRes: 2.1
ChimeraX-TapeMeasure: 1.0
ChimeraX-Test: 1.0
ChimeraX-Toolbar: 1.1
ChimeraX-ToolshedUtils: 1.2
ChimeraX-Tug: 1.0
ChimeraX-UI: 1.13
ChimeraX-uniprot: 2.2
ChimeraX-UnitCell: 1.0
ChimeraX-ViewDockX: 1.0.1
ChimeraX-VIPERdb: 1.0
ChimeraX-Vive: 1.1
ChimeraX-VolumeMenu: 1.0
ChimeraX-VTK: 1.0
ChimeraX-WavefrontOBJ: 1.0
ChimeraX-WebCam: 1.0
ChimeraX-WebServices: 1.0
ChimeraX-Zone: 1.0
colorama: 0.4.4
cxservices: 1.1
cycler: 0.10.0
Cython: 0.29.24
decorator: 5.1.0
distlib: 0.3.2
docutils: 0.17.1
filelock: 3.0.12
funcparserlib: 0.3.6
grako: 3.16.5
h5py: 3.4.0
html2text: 2020.1.16
idna: 3.2
ihm: 0.20
imagecodecs: 2021.4.28
imagesize: 1.2.0
ipykernel: 5.5.5
ipython: 7.23.1
ipython-genutils: 0.2.0
jedi: 0.18.0
Jinja2: 3.0.1
jupyter-client: 6.1.12
jupyter-core: 4.8.1
kiwisolver: 1.3.2
lxml: 4.6.3
lz4: 3.1.3
MarkupSafe: 2.0.1
matplotlib: 3.4.3
matplotlib-inline: 0.1.3
msgpack: 1.0.2
netCDF4: 1.5.7
networkx: 2.6.3
numexpr: 2.7.3
numpy: 1.21.2
openvr: 1.16.801
packaging: 21.0
ParmEd: 3.2.0
parso: 0.8.2
pexpect: 4.8.0
pickleshare: 0.7.5
Pillow: 8.3.2
pip: 21.2.4
pkginfo: 1.7.1
prompt-toolkit: 3.0.20
psutil: 5.8.0
ptyprocess: 0.7.0
pycollada: 0.7.1
pydicom: 2.1.2
Pygments: 2.10.0
PyOpenGL: 3.1.5
PyOpenGL-accelerate: 3.1.5
pyparsing: 2.4.7
PyQt5-commercial: 5.15.2
PyQt5-sip: 12.8.1
PyQtWebEngine-commercial: 5.15.2
python-dateutil: 2.8.2
pytz: 2021.1
pyzmq: 22.3.0
qtconsole: 5.1.1
QtPy: 1.11.1
RandomWords: 0.3.0
requests: 2.26.0
scipy: 1.7.1
setuptools: 57.5.0
sfftk-rw: 0.7.1
six: 1.16.0
snowballstemmer: 2.1.0
sortedcontainers: 2.4.0
Sphinx: 4.2.0
sphinx-autodoc-typehints: 1.12.0
sphinxcontrib-applehelp: 1.0.2
sphinxcontrib-blockdiag: 2.0.0
sphinxcontrib-devhelp: 1.0.2
sphinxcontrib-htmlhelp: 2.0.0
sphinxcontrib-jsmath: 1.0.1
sphinxcontrib-qthelp: 1.0.3
sphinxcontrib-serializinghtml: 1.1.5
suds-jurko: 0.6
tifffile: 2021.4.8
tinyarray: 1.2.3
tornado: 6.1
traitlets: 5.1.0
urllib3: 1.26.6
wcwidth: 0.2.5
webcolors: 1.11.1
wheel: 0.37.0
wheel-filename: 1.3.0
Change History (5)
comment:1 by , 4 years ago
| Component: | Unassigned → Structure Prediction |
|---|---|
| Owner: | set to |
| Platform: | → all |
| Project: | → ChimeraX |
| Status: | new → assigned |
| Summary: | ChimeraX bug report submission → Modeller: Number of residues in the alignment and pdb files are different |
comment:2 by , 4 years ago
follow-up: 4 comment:4 by , 4 years ago
I'm getting this error for everything I try to loop-model today: now for 6g2e, a totally different protein than those previously mentioned in this ticket.
comment:5 by , 4 years ago
Note:
See TracTickets
for help on using tickets.
Neither sequence is a subset of the other, which is undoubtedly where I'm getting tripped up in preparing the input for Modeller.