Opened 5 months ago

Last modified 5 months ago

#17899 accepted defect

UniProt sequences should retain numbering

Reported by: Tristan Croll Owned by: pett
Priority: normal Milestone:
Component: Sequence Version:
Keywords: Cc:
Blocked By: Blocking:
Notify when closed: Platform: all
Project: ChimeraX

Description

The following bug report has been submitted:
Platform:        Linux-6.8.0-59-generic-x86_64-with-glibc2.35
ChimeraX Version: 1.9 (2024-12-11 19:11:19 UTC)
Description
If I:
- have a model (e.g. an AlphaFold structure) open containing a fragment of a larger protein
- open the associated UniProt entry
- renumber the protein residues to match the UniProt numbering:

"renumber /A start 647"

... then the protein model is renumbered correctly, but the complete associated sequence is renumbered to start from 647 (so residue 647 of the protein is now 1284 of the sequence).

Log:
Startup Messages  
---  
note | available bundle cache has not been initialized yet  
  
Running command "isolde shorthand"  

> isolde shorthand
    
    
    Initialising ISOLDE-specific command aliases:
    Alias	Equivalent full command
    -------------------------------------------------
    st	isolde step {arguments}
    aw	isolde add water {arguments}
    awsf	isolde add water {arguments} sim false
    al	isolde add ligand {arguments}
    aa	isolde add aa $1 sel {arguments}
    ht	isolde mod his sel {arguments}
    so	setattr sel atoms occupancy {arguments}
    ab	isolde adjust bfactors {arguments}
    ss	isolde sim start sel
    rt	isolde release torsions sel {arguments}
    rd	isolde release distances sel {arguments}
    ra	rd; rt
    pf	isolde pepflip sel
    cf	isolde cisflip sel
    cbb	color bfactor {arguments}
    cbo	color byattr occupancy {arguments}
    cbc	color {arguments} bychain; color {arguments} byhet
    cs	clipper set contourSensitivity {arguments}
    

  
UCSF ChimeraX version: 1.9 (2024-12-11)  
© 2016-2024 Regents of the University of California. All rights reserved.  
How to cite UCSF ChimeraX  
Running command "open
perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_1_seed_000.pdb"  

> open perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_1_seed_000.pdb

Chain information for
perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_1_seed_000.pdb #1  
---  
Chain | Description  
A | No description available  
B | No description available  
  
Computing secondary structure  
Running command "open
perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_2_seed_000.pdb"  

> open perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_2_seed_000.pdb

Chain information for
perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_2_seed_000.pdb #2  
---  
Chain | Description  
A | No description available  
B | No description available  
  
Running command "open
perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_3_seed_000.pdb"  

> open perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_3_seed_000.pdb

Chain information for
perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_3_seed_000.pdb #3  
---  
Chain | Description  
A | No description available  
B | No description available  
  
Computing secondary structure  
[Repeated 1 time(s)]

> matchmaker #2,3 to #1

Computing secondary structure  
[Repeated 2 time(s)]  Parameters  
---  
Chain pairing | bb  
Alignment algorithm | Needleman-Wunsch  
Similarity matrix | BLOSUM-62  
SS fraction | 0.3  
Gap open (HH/SS/other) | 18/18/6  
Gap extend | 1  
SS matrix |  |  | H | S | O  
---|---|---|---  
H | 6 | -9 | -6  
S |  | 6 | -6  
O |  |  | 4  
Iteration cutoff | 2  
  
Matchmaker
perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_1_seed_000.pdb, chain
B (#1) with
perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_2_seed_000.pdb, chain
B (#2), sequence alignment score = 1288.8  
RMSD between 241 pruned atom pairs is 0.278 angstroms; (across all 247 pairs:
0.761)  
  
Matchmaker
perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_1_seed_000.pdb, chain
B (#1) with
perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_3_seed_000.pdb, chain
B (#3), sequence alignment score = 1288.8  
RMSD between 243 pruned atom pairs is 0.455 angstroms; (across all 247 pairs:
0.593)  
  

> cbc

Running command "ui tool show "AlphaFold Error Plot""  

> ui tool show "AlphaFold Error Plot"

Running command "alphafold pae #3 file
/home/tcroll/colabfold/predictions/perk_short_loop_tmem33/output/perk_loop_tmem33_scores_alphafold2_multimer_v3_model_3_seed_000.json"  

> alphafold pae #3 file
> /home/tcroll/colabfold/predictions/perk_short_loop_tmem33/output/perk_loop_tmem33_scores_alphafold2_multimer_v3_model_3_seed_000.json

> alphafold contacts #3/A toAtoms #3/B

Found 23 residue or atom pairs within distance 3  
Running command "hide #2 models"  

> hide #2 models

Running command "hide #1 models"  

> hide #1 models

> show #3

Running command "show #2 models"  

> show #2 models

Running command "show #1 models"  

> show #1 models

Running command "hide #!3 models"  

> hide #!3 models

Running command "hide #2 models"  

> hide #2 models

Running command "ui tool show "AlphaFold Error Plot""  

> ui tool show "AlphaFold Error Plot"

Running command "alphafold pae #1 file
/home/tcroll/colabfold/predictions/perk_short_loop_tmem33/output/perk_loop_tmem33_scores_alphafold2_multimer_v3_model_1_seed_000.json"  

> alphafold pae #1 file
> /home/tcroll/colabfold/predictions/perk_short_loop_tmem33/output/perk_loop_tmem33_scores_alphafold2_multimer_v3_model_1_seed_000.json

> alphafold contacts #1/B toAtoms #1/A

Found 15 residue or atom pairs within distance 3  

> show #1

> clipper init

Running command "cofr center showPivot t"  

> cofr centerOfView showPivot true

Running command "camera ortho"  

> camera ortho

Running command "lighting simple"  

> lighting simple

> renumber /A start 647

150 residues renumbered  

> open Q9NZJ5 fromDatabase uniprot

Summary of feedback from opening Q9NZJ5 fetched from uniprot  
---  
notes | Alignment identifier is Q9NZJ5  
Associated
perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_1_seed_000.pdb chain A
to Q9NZJ5 with 0 mismatches  
Associated
perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_2_seed_000.pdb chain A
to Q9NZJ5 with 0 mismatches  
Associated
perk_loop_tmem33_unrelaxed_alphafold2_multimer_v3_model_3_seed_000.pdb chain A
to Q9NZJ5 with 0 mismatches  
  
Opened UniProt Q9NZJ5  

> select /A:THR

21 atoms, 18 bonds, 3 residues, 3 models selected  
Running command "~sel"  

> ~select

Nothing selected  
Running command "wait 1"  
Running command "~sel"  

> ~select

Nothing selected  
Running command "wait 1"  
Running command "~sel"  

> ~select

Nothing selected  
Running command "wait 1"  
Running command "~sel"  

> ~select

Nothing selected  
Running command "wait 1"  
Running command "~sel"  

> ~select

Nothing selected  
Running command "wait 1"  

> close

Running command "open
/home/tcroll/colabfold/predictions/perk_short_loop_tmem33/output/perk_loop_tmem33_unrelaxed_rank_001_alphafold2_multimer_v3_model_4_seed_000.pdb"  

> open
> /home/tcroll/colabfold/predictions/perk_short_loop_tmem33/output/perk_loop_tmem33_unrelaxed_rank_001_alphafold2_multimer_v3_model_4_seed_000.pdb

Chain information for
perk_loop_tmem33_unrelaxed_rank_001_alphafold2_multimer_v3_model_4_seed_000.pdb
#1  
---  
Chain | Description  
A | No description available  
B | No description available  
  
Associated
perk_loop_tmem33_unrelaxed_rank_001_alphafold2_multimer_v3_model_4_seed_000.pdb
chain A to Q9NZJ5 with 0 mismatches  
Computing secondary structure  

> cbc

> renumber /A&protein start 647

50 residues renumbered  

> select /A

414 atoms, 431 bonds, 50 residues, 1 model selected  

> renumber /A start 1

50 residues renumbered  

> renumber /A start 647

50 residues renumbered  




OpenGL version: 3.3.0 NVIDIA 535.230.02
OpenGL renderer: NVIDIA GeForce RTX 3070/PCIe/SSE2
OpenGL vendor: NVIDIA Corporation

Python: 3.11.4
Locale: en_GB.UTF-8
Qt version: PyQt6 6.7.1, Qt 6.7.1
Qt runtime version: 6.7.3
Qt platform: xcb

XDG_SESSION_TYPE=x11
DESKTOP_SESSION=ubuntu
XDG_SESSION_DESKTOP=ubuntu
XDG_CURRENT_DESKTOP=ubuntu:GNOME
DISPLAY=:1
Manufacturer: Dell Inc.
Model: XPS 8950
OS: Ubuntu 22.04
Architecture: 64bit ELF
Virtual Machine: none
CPU: 20 12th Gen Intel(R) Core(TM) i7-12700
Cache Size: 25600 KB
Memory:
	               total        used        free      shared  buff/cache   available
	Mem:            31Gi        12Gi       7.2Gi       788Mi        11Gi        17Gi
	Swap:          2.0Gi       1.5Gi       487Mi

Graphics:
	0000:01:00.0 VGA compatible controller [0300]: NVIDIA Corporation GA104 [GeForce RTX 3070 Lite Hash Rate] [10de:2488] (rev a1)	
	Subsystem: Dell GA104 [GeForce RTX 3070 Lite Hash Rate] [1028:c903]	
	Kernel driver in use: nvidia

Installed Packages:
    alabaster: 1.0.0
    anyio: 4.7.0
    appdirs: 1.4.4
    asttokens: 3.0.0
    auditwheel: 6.1.0
    autocommand: 2.2.2
    babel: 2.16.0
    backports.tarfile: 1.2.0
    beautifulsoup4: 4.12.3
    blockdiag: 3.0.0
    blosc2: 3.0.0
    build: 1.2.1
    certifi: 2024.8.30
    cftime: 1.6.4.post1
    charset-normalizer: 3.4.0
    ChimeraX-AddCharge: 1.5.18
    ChimeraX-AddH: 2.2.6
    ChimeraX-AlignmentAlgorithms: 2.0.2
    ChimeraX-AlignmentHdrs: 3.5
    ChimeraX-AlignmentMatrices: 2.1
    ChimeraX-Alignments: 2.16.1
    ChimeraX-AlphaFold: 1.0.1
    ChimeraX-AltlocExplorer: 1.1.2
    ChimeraX-AmberInfo: 1.0
    ChimeraX-Arrays: 1.1
    ChimeraX-Atomic: 1.58.8
    ChimeraX-AtomicLibrary: 14.1.11
    ChimeraX-AtomSearch: 2.0.1
    ChimeraX-AxesPlanes: 2.4
    ChimeraX-BasicActions: 1.1.2
    ChimeraX-BILD: 1.0
    ChimeraX-BlastProtein: 3.0.0
    ChimeraX-BondRot: 2.0.4
    ChimeraX-BugReporter: 1.0.1
    ChimeraX-BuildStructure: 2.13.1
    ChimeraX-Bumps: 1.0
    ChimeraX-BundleBuilder: 1.4.0
    ChimeraX-ButtonPanel: 1.0.1
    ChimeraX-CageBuilder: 1.0.1
    ChimeraX-CellPack: 1.0
    ChimeraX-Centroids: 1.4
    ChimeraX-ChangeChains: 1.1
    ChimeraX-CheckWaters: 1.4
    ChimeraX-ChemGroup: 2.0.1
    ChimeraX-Clashes: 2.3
    ChimeraX-Clipper: 0.24.0
    ChimeraX-ColorActions: 1.0.5
    ChimeraX-ColorGlobe: 1.0
    ChimeraX-ColorKey: 1.5.6
    ChimeraX-CommandLine: 1.2.5
    ChimeraX-ConnectStructure: 2.0.1
    ChimeraX-Contacts: 1.0.1
    ChimeraX-Core: 1.9
    ChimeraX-CoreFormats: 1.2
    ChimeraX-coulombic: 1.4.4
    ChimeraX-Crosslinks: 1.0
    ChimeraX-Crystal: 1.0
    ChimeraX-CrystalContacts: 1.0.1
    ChimeraX-DataFormats: 1.2.3
    ChimeraX-Dicom: 1.2.6
    ChimeraX-DistMonitor: 1.4.2
    ChimeraX-DockPrep: 1.1.3
    ChimeraX-Dssp: 2.0
    ChimeraX-EMDB-SFF: 1.0
    ChimeraX-ESMFold: 1.0
    ChimeraX-FileHistory: 1.0.1
    ChimeraX-FunctionKey: 1.0.1
    ChimeraX-Geometry: 1.3
    ChimeraX-gltf: 1.0
    ChimeraX-Graphics: 1.4.1
    ChimeraX-Hbonds: 2.5
    ChimeraX-Help: 1.3
    ChimeraX-HKCage: 1.3
    ChimeraX-IHM: 1.1
    ChimeraX-ImageFormats: 1.2
    ChimeraX-IMOD: 1.0
    ChimeraX-IO: 1.0.3
    ChimeraX-ISOLDE: 1.9
    ChimeraX-ItemsInspection: 1.0.1
    ChimeraX-IUPAC: 1.0
    ChimeraX-KVFinder: 1.2.1
    ChimeraX-Label: 1.1.14
    ChimeraX-LinuxSupport: 1.0.1
    ChimeraX-ListInfo: 1.2.2
    ChimeraX-Log: 1.2
    ChimeraX-LookingGlass: 1.1
    ChimeraX-Maestro: 1.9.1
    ChimeraX-Map: 1.3
    ChimeraX-MapData: 2.0
    ChimeraX-MapEraser: 1.0.1
    ChimeraX-MapFilter: 2.0.1
    ChimeraX-MapFit: 2.0
    ChimeraX-MapSeries: 2.1.1
    ChimeraX-Markers: 1.0.1
    ChimeraX-Mask: 1.0.2
    ChimeraX-MatchMaker: 2.1.6
    ChimeraX-MCopy: 1.0
    ChimeraX-MDcrds: 2.7.2
    ChimeraX-MedicalToolbar: 1.1
    ChimeraX-Meeting: 1.0.1
    ChimeraX-MLP: 1.1.1
    ChimeraX-mmCIF: 2.14.2
    ChimeraX-MMTF: 2.2
    ChimeraX-ModelArchive: 1.0
    ChimeraX-Modeller: 1.5.18
    ChimeraX-ModelPanel: 1.5
    ChimeraX-ModelSeries: 1.0.1
    ChimeraX-Mol2: 2.0.3
    ChimeraX-Mole: 1.0
    ChimeraX-Morph: 1.0.2
    ChimeraX-MouseModes: 1.2
    ChimeraX-Movie: 1.0
    ChimeraX-MutationScores: 1.0
    ChimeraX-Neuron: 1.0
    ChimeraX-Nifti: 1.2
    ChimeraX-NMRSTAR: 1.0.2
    ChimeraX-NRRD: 1.2
    ChimeraX-Nucleotides: 2.0.3
    ChimeraX-OpenCommand: 1.14
    ChimeraX-OrthoPick: 1.0.1
    ChimeraX-PDB: 2.7.6
    ChimeraX-PDBBio: 1.0.1
    ChimeraX-PDBLibrary: 1.0.4
    ChimeraX-PDBMatrices: 1.0
    ChimeraX-PickBlobs: 1.0.1
    ChimeraX-Positions: 1.0
    ChimeraX-PresetMgr: 1.1.2
    ChimeraX-PubChem: 2.2
    ChimeraX-QScore: 1.2
    ChimeraX-ReadPbonds: 1.0.1
    ChimeraX-Registration: 1.1.2
    ChimeraX-RemoteControl: 1.0
    ChimeraX-RenderByAttr: 1.6.2
    ChimeraX-RenumberResidues: 1.1
    ChimeraX-ResidueFit: 1.0.1
    ChimeraX-RestServer: 1.3.1
    ChimeraX-RNALayout: 1.0
    ChimeraX-RotamerLibMgr: 4.0
    ChimeraX-RotamerLibsDunbrack: 2.0
    ChimeraX-RotamerLibsDynameomics: 2.0
    ChimeraX-RotamerLibsRichardson: 2.0
    ChimeraX-SaveCommand: 1.5.1
    ChimeraX-SchemeMgr: 1.0
    ChimeraX-SDF: 2.0.2
    ChimeraX-Segger: 1.0
    ChimeraX-Segment: 1.0.1
    ChimeraX-Segmentations: 3.5.6
    ChimeraX-SelInspector: 1.0
    ChimeraX-SeqView: 2.14
    ChimeraX-Shape: 1.0.1
    ChimeraX-Shell: 1.0.1
    ChimeraX-Shortcuts: 1.2.0
    ChimeraX-ShowSequences: 1.0.3
    ChimeraX-SideView: 1.0.1
    ChimeraX-SimilarStructures: 1.0.1
    ChimeraX-Smiles: 2.1.2
    ChimeraX-SmoothLines: 1.0
    ChimeraX-SpaceNavigator: 1.0
    ChimeraX-StdCommands: 1.18.1
    ChimeraX-STL: 1.0.1
    ChimeraX-Storm: 1.0
    ChimeraX-StructMeasure: 1.2.1
    ChimeraX-Struts: 1.0.1
    ChimeraX-Surface: 1.0.1
    ChimeraX-SwapAA: 2.0.1
    ChimeraX-SwapRes: 2.5
    ChimeraX-TapeMeasure: 1.0
    ChimeraX-TaskManager: 1.0
    ChimeraX-Test: 1.0
    ChimeraX-Toolbar: 1.2.3
    ChimeraX-ToolshedUtils: 1.2.4
    ChimeraX-Topography: 1.0
    ChimeraX-ToQuest: 1.0
    ChimeraX-Tug: 1.0.1
    ChimeraX-UI: 1.41
    ChimeraX-Umap: 1.0
    ChimeraX-uniprot: 2.3.1
    ChimeraX-UnitCell: 1.0.1
    ChimeraX-ViewDockX: 1.4.4
    ChimeraX-VIPERdb: 1.0
    ChimeraX-Vive: 1.1
    ChimeraX-VolumeMenu: 1.0.1
    ChimeraX-vrml: 1.0
    ChimeraX-VTK: 1.0
    ChimeraX-WavefrontOBJ: 1.0
    ChimeraX-WebCam: 1.0.2
    ChimeraX-WebServices: 1.1.4
    ChimeraX-Zone: 1.0.1
    colorama: 0.4.6
    comm: 0.2.2
    contourpy: 1.3.1
    cxservices: 1.2.3
    cycler: 0.12.1
    Cython: 3.0.10
    debugpy: 1.8.9
    decorator: 5.1.1
    distro: 1.9.0
    docutils: 0.21.2
    executing: 2.1.0
    filelock: 3.15.4
    fonttools: 4.55.3
    funcparserlib: 2.0.0a0
    glfw: 2.8.0
    grako: 3.16.5
    h11: 0.14.0
    h5py: 3.12.1
    html2text: 2024.2.26
    httpcore: 1.0.7
    httpx: 0.28.1
    idna: 3.10
    ihm: 1.3
    imagecodecs: 2024.6.1
    imagesize: 1.4.1
    importlib_metadata: 8.0.0
    importlib_resources: 6.4.0
    inflect: 7.3.1
    ipykernel: 6.29.5
    ipython: 8.26.0
    ipywidgets: 8.1.5
    jaraco.context: 5.3.0
    jaraco.functools: 4.0.1
    jaraco.text: 3.12.1
    jedi: 0.19.1
    Jinja2: 3.1.4
    jupyter_client: 8.6.2
    jupyter_core: 5.7.2
    jupyterlab_widgets: 3.0.13
    kiwisolver: 1.4.7
    line_profiler: 4.1.3
    lxml: 5.2.2
    lz4: 4.3.3
    MarkupSafe: 3.0.2
    matplotlib: 3.9.2
    matplotlib-inline: 0.1.7
    more-itertools: 10.3.0
    msgpack: 1.0.8
    ndindex: 1.9.2
    nest-asyncio: 1.6.0
    netCDF4: 1.6.5
    networkx: 3.3
    nibabel: 5.2.0
    nptyping: 2.5.0
    numexpr: 2.10.2
    numpy: 1.26.4
    openvr: 1.26.701
    ordered-set: 4.1.0
    packaging: 23.2
    packaging: 24.1
    ParmEd: 4.2.2
    parso: 0.8.4
    pep517: 0.13.1
    pexpect: 4.9.0
    pillow: 10.4.0
    pip: 24.2
    pkginfo: 1.11.1
    platformdirs: 4.3.6
    platformdirs: 4.2.2
    prompt_toolkit: 3.0.48
    psutil: 6.0.0
    ptyprocess: 0.7.0
    pure_eval: 0.2.3
    py-cpuinfo: 9.0.0
    pycollada: 0.8
    pydicom: 2.4.4
    pyelftools: 0.31
    Pygments: 2.18.0
    pynmrstar: 3.3.4
    pynrrd: 1.0.0
    PyOpenGL: 3.1.7
    PyOpenGL-accelerate: 3.1.7
    pyopenxr: 1.0.3401
    pyparsing: 3.2.0
    pyproject_hooks: 1.2.0
    PyQt6-commercial: 6.7.1
    PyQt6-Qt6: 6.7.3
    PyQt6-WebEngine-commercial: 6.7.0
    PyQt6-WebEngine-Qt6: 6.7.3
    PyQt6-WebEngineSubwheel-Qt6: 6.7.3
    PyQt6_sip: 13.8.0
    python-dateutil: 2.9.0.post0
    pytz: 2024.2
    pyzmq: 26.2.0
    qtconsole: 5.5.2
    QtPy: 2.4.2
    qtshim: 1.0
    RandomWords: 0.4.0
    rdkit: 2024.9.6
    requests: 2.32.3
    scipy: 1.14.0
    selfies: 2.2.0
    Send2Trash: 1.8.3
    SEQCROW: 1.8.17
    setuptools: 72.1.0
    sfftk-rw: 0.8.1
    six: 1.16.0
    sniffio: 1.3.1
    snowballstemmer: 2.2.0
    sortedcontainers: 2.4.0
    soupsieve: 2.6
    Sphinx: 8.0.2
    sphinx-autodoc-typehints: 2.2.3
    sphinxcontrib-applehelp: 2.0.0
    sphinxcontrib-blockdiag: 3.0.0
    sphinxcontrib-devhelp: 2.0.0
    sphinxcontrib-htmlhelp: 2.1.0
    sphinxcontrib-jsmath: 1.0.1
    sphinxcontrib-qthelp: 2.0.0
    sphinxcontrib-serializinghtml: 2.0.0
    stack-data: 0.6.3
    superqt: 0.6.3
    tables: 3.10.1
    tcia_utils: 1.5.1
    tifffile: 2024.7.24
    tinyarray: 1.2.4
    tomli: 2.0.1
    tornado: 6.4.2
    traitlets: 5.14.3
    typeguard: 4.3.0
    typing_extensions: 4.12.2
    typing_extensions: 4.12.2
    tzdata: 2024.2
    urllib3: 2.2.3
    wcwidth: 0.2.13
    webcolors: 24.6.0
    wheel: 0.43.0
    wheel: 0.43.0
    wheel-filename: 1.4.1
    widgetsnbextension: 4.0.13
    zipp: 3.19.2

Change History (1)

comment:1 by pett, 5 months ago

Component: UnassignedSequence
Owner: set to pett
Platform: all
Project: ChimeraX
Status: newaccepted
Summary: ChimeraX bug report submissionUniProt sequences should retain numbering
Note: See TracTickets for help on using tickets.