Opened 5 months ago
Last modified 5 months ago
#17899 accepted defect
UniProt sequences should retain numbering
| Reported by: | Tristan Croll | Owned by: | Eric Pettersen |
|---|---|---|---|
| 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 , 5 months ago
| Component: | Unassigned → Sequence |
|---|---|
| Owner: | set to |
| Platform: | → all |
| Project: | → ChimeraX |
| Status: | new → accepted |
| Summary: | ChimeraX bug report submission → UniProt sequences should retain numbering |
Note:
See TracTickets
for help on using tickets.