
Hi, I just installed chimerax RC-1.8 and tried to use "modeller" (tools->sequence->modeller comparative). The resulting error message is attached to this email. If I use chimerax 1.7 with the very same data modeller works. Am I doing something wrong or did I find a bug? cheers Roberto Follows the complete log: === modeller comparative aligned_1.fasta:2 numModels 5 fast false multichain false hetPreserve true hydrogens false waterPreserve falseWebservices job id: GN1G33GHX39IPCAE Modeller job (ID GN1G33GHX39IPCAE) finished Modeller error output Traceback (most recent call last): File "ModellerModelling.py", line 74, in <module> a.make() File "/usr/lib64/python3.8/site-packages/modeller/automodel/automodel.py", line 141, in make self.homcsr(exit_stage) File "/usr/lib64/python3.8/site-packages/modeller/automodel/automodel.py", line 624, in homcsr self.check_alignment(aln) File "/usr/lib64/python3.8/site-packages/modeller/automodel/automodel.py", line 577, in check_alignment aln.check() File "/usr/lib64/python3.8/site-packages/modeller/alignment.py", line 213, in check self.check_structure_structure(io=io) File "/usr/lib64/python3.8/site-packages/modeller/alignment.py", line 222, in check_structure_structure return f(self.modpt, io.modpt, self.env.libs.modpt, eqvdst) _modeller.SequenceMismatchError: read_te_291E> Sequence difference between alignment and pdb : Modeller run output MODELLER 10.5, 2024/01/19, r12711 PROTEIN STRUCTURE MODELLING BY SATISFACTION OF SPATIAL RESTRAINTS Copyright(c) 1989-2024 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 wilkins.cgl.ucsf.edu 4.18.0-513.24.1.el8_9.x86_64 x86_64 Date and time of compilation : 2024/01/20 00:40:16 MODELLER executable type : x86_64-intel8 Job starting time (YY/MM/DD HH:MM:SS): 2024/06/06 03:33:20 environ____W> The class 'environ' is deprecated; use 'Environ' instead openf___224_> Open $(LIB)/restyp.lib openf___224_> Open ${MODINSTALL10v5}/modlib/resgrp.lib rdresgr_266_> Number of residue groups: 2 openf___224_> Open ${MODINSTALL10v5}/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 ${MODINSTALL10v5}/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 ${MODINSTALL10v5}/modlib/radii.lib Dynamically allocated memory at amaxlibraries [B,KiB,MiB]: 253994 248.041 0.242 openf___224_> Open ${MODINSTALL10v5}/modlib/af_mnchdef.lib rdwilmo_274_> Mainchain residue conformation classes: APBLE openf___224_> Open ${MODINSTALL10v5}/modlib/mnch.lib rdclass_257_> Number of classes: 5 openf___224_> Open ${MODINSTALL10v5}/modlib/mnch1.lib openf___224_> Open ${MODINSTALL10v5}/modlib/mnch2.lib openf___224_> Open ${MODINSTALL10v5}/modlib/mnch3.lib openf___224_> Open ${MODINSTALL10v5}/modlib/xs4.mat rdrrwgh_268_> Number of residue types: 21 automodel__W> The class 'automodel' is deprecated; use 'AutoModel' instead openf___224_> Open alignment.ali Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 282425 275.806 0.269 Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 297517 290.544 0.284 Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 312917 305.583 0.298 Dynamically allocated memory at amaxalignment [B,KiB,MiB]: 343717 335.661 0.328 Dynamically allocated memory at amaxsequence [B,KiB,MiB]: 344285 336.216 0.328 read_te_291E> Sequence difference between alignment and pdb : x (mismatch at alignment position 1) Alignment XSLSDKDKAAVRALWSKIGKSADAIGNDALSRMIVVYPQTKTYFSHWPDVTPGSP PDB SLSDKDKAAVRALWSKIGKSADAIGNDALSRMIVVYPQTKTYFSHWPDVTPGSPH Match * * Alignment residue type 24 (X, UNK) does not match pdb residue type 16 (S, SER), for align code 1pbx.cif_2/A (atom file 1pbx.cif_2), pdb residue number " 1", chain "A" 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 1pbx.cif_2/A No output models from Modeller; see log for Modeller text output.