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: | pett |
---|---|---|---|
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.