Opened 4 years ago

Closed 4 years ago

#5165 closed defect (fixed)

AlphaFold second prediction overwrites first prediction

Reported by: goddard@… Owned by: Tom Goddard
Priority: normal Milestone:
Component: Structure Analysis Version:
Keywords: Cc:
Blocked By: Blocking:
Notify when closed: Platform: all
Project: ChimeraX

Description

The following bug report has been submitted:
Platform:        macOS-10.16-x86_64-i386-64bit
ChimeraX Version: 1.3.dev202108251821 (2021-08-25 18:21:58 UTC)
Description
Running a second AlphaFold prediction without closing the Colab window put the results in the same Downloads directory prediction_12 and loaded the wrong file for the second prediction since it was named best_model (2).pdb.

Log:
UCSF ChimeraX version: 1.3.dev202108251821 (2021-08-25)  
© 2016-2021 Regents of the University of California. All rights reserved.  
How to cite UCSF ChimeraX  

> open 7jtl

Summary of feedback from opening 7jtl fetched from pdb  
---  
note | Fetching compressed mmCIF 7jtl from
http://files.rcsb.org/download/7jtl.cif  
  
7jtl title:  
Structure of SARS-CoV-2 ORF8 accessory protein [more info...]  
  
Chain information for 7jtl #1  
---  
Chain | Description | UniProt  
A B | ORF8 protein | NS8_SARS2  
  
Non-standard residues in 7jtl #1  
---  
NA — sodium ion  
  

> ui tool show "Blast Protein"

> blastprotein /A database pdb cutoff 1e-3 matrix BLOSUM62 maxSeqs 100 name
> bp1

ChimeraX REST job id: job_3wikysds  
BlastProtein finished.  
Parsing BLAST results.  

> blastprotein /A database pdb cutoff 1e-1 matrix BLOSUM62 maxSeqs 100 name
> bp1

ChimeraX REST job id: job_1tomv5gm  
BlastProtein finished.  
Parsing BLAST results.  

> open pdb:7JX6

Summary of feedback from opening 7JX6 fetched from pdb  
---  
note | Fetching compressed mmCIF 7jx6 from
http://files.rcsb.org/download/7jx6.cif  
  
7jx6 title:  
STRUCTURE OF THE SARS-CoV-2 ORF8 ENCODED ACCESSORY PROTEIN [more info...]  
  
Chain information for 7jx6 #2  
---  
Chain | Description | UniProt  
A B | ORF8 protein | NS8_SARS2  
  
Non-standard residues in 7jx6 #2  
---  
NA — sodium ion  
  
7jx6 mmCIF Assemblies  
---  
1| author_and_software_defined_assembly  
2| software_defined_assembly  
  

> matchmaker #2/B to #1/A

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 7jtl, chain A (#1) with 7jx6, chain B (#2), sequence alignment
score = 509.5  
RMSD between 82 pruned atom pairs is 0.471 angstroms; (across all 88 pairs:
1.033)  
  

> close

> open /Users/goddard/Desktop/rotavirus_nsp5.fasta

Summary of feedback from opening /Users/goddard/Desktop/rotavirus_nsp5.fasta  
---  
note | Alignment identifier is rotavirus_nsp5.fasta  
  
Opened 1 sequences from rotavirus_nsp5.fasta  

> ui dockable true "Sequence Viewer"

> blastprotein rotavirus_nsp5.fasta:1

ChimeraX REST job id: job_nt3nk849  
BlastProtein finished.  
Parsing BLAST results.  

> usage blastprotein

blastprotein atoms [database database] [cutoff a number] [matrix matrix]
[maxSeqs an integer] [log true or false] [name a text string]  
— Search PDB/NR using BLAST  
atoms: an atom specifier or [alignment-id:]sequence-name-or-number  
database: one of alphafold, nr, or pdb  
matrix: one of BLOSUM45, BLOSUM50, BLOSUM62, BLOSUM80, BLOSUM90, IDENTITY,
PAM250, PAM30, or PAM70

blastprotein mav [name] [selected true or false]  
— Display BLAST results in MultAlign Viewer  
name: a text string  

> help blastprotein

> blastprotein rotavirus_nsp5.fasta:1 cutoff 1e-1

ChimeraX REST job id: job_ai_aq1xu  
BlastProtein finished.  
Parsing BLAST results.  

> ui tool show AlphaFold

> alphafold match
> MSLSIDVTSLPSIPSTIYKNESSSTTSTLSGKSIGRSEQYISPDAEAFNKYMLSKSPEDIGPSDSASNDPLTSFSIRSNAVKTNADAGVSMDSSAQSRPSSNVGCDQVDFSLNKGLKVKANLDSSISISTDTKKEKSKQNHKSRKHYPRIEAESDSDDYVLDDSDSDDGKCKNCKYKKKYFALRMRMKQVAMQLIEDL

No AlphaFold model with similar sequence for 1 sequences  

Opened 0 AlphaFold model  

> alphafold search
> MSLSIDVTSLPSIPSTIYKNESSSTTSTLSGKSIGRSEQYISPDAEAFNKYMLSKSPEDIGPSDSASNDPLTSFSIRSNAVKTNADAGVSMDSSAQSRPSSNVGCDQVDFSLNKGLKVKANLDSSISISTDTKKEKSKQNHKSRKHYPRIEAESDSDDYVLDDSDSDDGKCKNCKYKKKYFALRMRMKQVAMQLIEDL

ChimeraX REST job id: job_og1yw3d6  
BlastProtein finished.  
Parsing BLAST results.  

> alphafold predict
> MSLSIDVTSLPSIPSTIYKNESSSTTSTLSGKSIGRSEQYISPDAEAFNKYMLSKSPEDIGPSDSASNDPLTSFSIRSNAVKTNADAGVSMDSSAQSRPSSNVGCDQVDFSLNKGLKVKANLDSSISISTDTKKEKSKQNHKSRKHYPRIEAESDSDDYVLDDSDSDDGKCKNCKYKKKYFALRMRMKQVAMQLIEDL

Running AlphaFold prediction  
Fetching A2T3S0 UniProt info from https://www.uniprot.org/uniprot/A2T3S0.xml  

> alphafold search RDRP_ROTSH

ChimeraX REST job id: job_g9n13bed  
BlastProtein finished.  
Parsing BLAST results.  
Fetching A2T3M4 UniProt info from https://www.uniprot.org/uniprot/A2T3M4.xml  

> alphafold search NSP1_ROTSH

ChimeraX REST job id: job_qv33bcre  
BlastProtein finished.  
Parsing BLAST results.  
AlphaFold prediction finished  
Results in /Users/goddard/Downloads/ChimeraX/AlphaFold/prediction_12  
Chain information for best_model.pdb #1  
---  
Chain | Description  
A | No description available  
  

> show atoms

> hide atoms

Fetching P06179 UniProt info from https://www.uniprot.org/uniprot/P06179.xml  

> alphafold match FLIC_SALTY

1 UniProt id does not have an AlphaFold database model: P06179 (UniProt
FLIC_SALTY)  

Fetching compressed AlphaFold P04949 from
https://alphafold.ebi.ac.uk/files/AF-P04949-F1-model_v1.cif  
1 AlphaFold model found using sequence similarity searches: P04949 (UniProt
FLIC_SALTY)  
Opened 1 AlphaFold model  

> hide #1 models

> open FLIC_SALTY fromDatabase uniprot

Summary of feedback from opening FLIC_SALTY fetched from uniprot  
---  
notes | UniProt identifier FLIC_SALTY maps to entry P06179  
Fetching FLIC_SALTY UniProt info from
https://www.uniprot.org/uniprot/FLIC_SALTY.xml  
Alignment identifier is FLIC_SALTY  
  
Opened UniProt FLIC_SALTY  

> sequence associate #2/A FLIC_SALTY:1

Associated UniProt P04949 chain A to FLIC_SALTY with 238 mismatches and/or
gaps  

> alphafold search FLIC_SALTY

ChimeraX REST job id: job_woy004sw  
BlastProtein finished.  
Parsing BLAST results.  

> alphafold fetch P04949

Chain information for AlphaFold P04949 #3  
---  
Chain | Description | UniProt  
A | Flagellin | FLIC_ECOLI  
  

> select clear

Fetching A2T3R0 UniProt info from https://www.uniprot.org/uniprot/A2T3R0.xml  

> alphafold predict NSP6_ROTSH

Running AlphaFold prediction  

> hide #3 models

> hide #2 models

> show #1 models

AlphaFold prediction finished  
Results in /Users/goddard/Downloads/ChimeraX/AlphaFold/prediction_12  
Chain information for best_model.pdb #4  
---  
Chain | Description  
A | No description available  
  

> hide #1 models

> open "/Users/goddard/Downloads/ChimeraX/AlphaFold/prediction_12/best_model
> (2).pdb"

Chain information for best_model (2).pdb #5  
---  
Chain | Description  
A | No description available  
  




OpenGL version: 4.1 ATI-4.6.20
OpenGL renderer: AMD Radeon Pro Vega 20 OpenGL Engine
OpenGL vendor: ATI Technologies Inc.Hardware:

    Hardware Overview:

      Model Name: MacBook Pro
      Model Identifier: MacBookPro15,3
      Processor Name: 8-Core Intel Core i9
      Processor Speed: 2.4 GHz
      Number of Processors: 1
      Total Number of Cores: 8
      L2 Cache (per Core): 256 KB
      L3 Cache: 16 MB
      Hyper-Threading Technology: Enabled
      Memory: 32 GB
      System Firmware Version: 1554.140.20.0.0 (iBridge: 18.16.14759.0.1,0)

Software:

    System Software Overview:

      System Version: macOS 11.5.2 (20G95)
      Kernel Version: Darwin 20.6.0
      Time since boot: 2 days 9:04

Graphics/Displays:

    Intel UHD Graphics 630:

      Chipset Model: Intel UHD Graphics 630
      Type: GPU
      Bus: Built-In
      VRAM (Dynamic, Max): 1536 MB
      Vendor: Intel
      Device ID: 0x3e9b
      Revision ID: 0x0002
      Automatic Graphics Switching: Supported
      gMux Version: 5.0.0
      Metal Family: Supported, Metal GPUFamily macOS 2

    Radeon Pro Vega 20:

      Chipset Model: Radeon Pro Vega 20
      Type: GPU
      Bus: PCIe
      PCIe Lane Width: x8
      VRAM (Total): 4 GB
      Vendor: AMD (0x1002)
      Device ID: 0x69af
      Revision ID: 0x00c0
      ROM Revision: 113-D2060I-087
      VBIOS Version: 113-D20601MA0T-016
      Option ROM Version: 113-D20601MA0T-016
      EFI Driver Version: 01.01.087
      Automatic Graphics Switching: Supported
      gMux Version: 5.0.0
      Metal Family: Supported, Metal GPUFamily macOS 2
      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: Yes
          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
    chardet: 4.0.0
    charset-normalizer: 2.0.4
    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.1
    ChimeraX-AlphaFold: 1.0
    ChimeraX-AltlocExplorer: 1.0
    ChimeraX-AmberInfo: 1.0
    ChimeraX-Arrays: 1.0
    ChimeraX-Atomic: 1.28.1
    ChimeraX-AtomicLibrary: 4.1
    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.5.3
    ChimeraX-Bumps: 1.0
    ChimeraX-BundleBuilder: 1.1
    ChimeraX-ButtonPanel: 1.0
    ChimeraX-CageBuilder: 1.0
    ChimeraX-CellPack: 1.0
    ChimeraX-Centroids: 1.1
    ChimeraX-ChemGroup: 2.0
    ChimeraX-Clashes: 2.1
    ChimeraX-ColorActions: 1.0
    ChimeraX-ColorGlobe: 1.0
    ChimeraX-ColorKey: 1.3.3
    ChimeraX-CommandLine: 1.1.4
    ChimeraX-ConnectStructure: 2.0
    ChimeraX-Contacts: 1.0
    ChimeraX-Core: 1.3.dev202108251821
    ChimeraX-CoreFormats: 1.1
    ChimeraX-coulombic: 1.3
    ChimeraX-Crosslinks: 1.0
    ChimeraX-Crystal: 1.0
    ChimeraX-CrystalContacts: 1.0
    ChimeraX-DataFormats: 1.2
    ChimeraX-Dicom: 1.0
    ChimeraX-DistMonitor: 1.1.4
    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
    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
    ChimeraX-MDcrds: 2.4
    ChimeraX-MedicalToolbar: 1.0.1
    ChimeraX-Meeting: 1.0
    ChimeraX-MLP: 1.1
    ChimeraX-mmCIF: 2.3
    ChimeraX-MMTF: 2.1
    ChimeraX-Modeller: 1.1
    ChimeraX-ModelPanel: 1.1
    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.5.1
    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.1
    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.12
    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
    comtypes: 1.1.10
    cxservices: 1.0.1
    cycler: 0.10.0
    Cython: 0.29.23
    decorator: 5.0.9
    distlib: 0.3.1
    docutils: 0.17.1
    filelock: 3.0.12
    funcparserlib: 0.3.6
    grako: 3.16.5
    h5py: 3.3.0
    html2text: 2020.1.16
    idna: 2.10
    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: 2.11.3
    jupyter-client: 6.1.12
    jupyter-core: 4.7.1
    kiwisolver: 1.3.1
    lxml: 4.6.3
    lz4: 3.1.3
    MarkupSafe: 1.1.1
    matplotlib: 3.4.2
    matplotlib-inline: 0.1.2
    msgpack: 1.0.2
    netCDF4: 1.5.6
    networkx: 2.6.2
    numexpr: 2.7.3
    numpy: 1.21.0
    numpydoc: 1.1.0
    openvr: 1.16.801
    packaging: 21.0
    pandas: 1.3.2
    ParmEd: 3.2.0
    parso: 0.8.2
    pexpect: 4.8.0
    pickleshare: 0.7.5
    Pillow: 8.2.0
    pip: 21.1.1
    pkginfo: 1.7.0
    prompt-toolkit: 3.0.19
    psutil: 5.8.0
    ptyprocess: 0.7.0
    pycollada: 0.7.1
    pydicom: 2.1.2
    Pygments: 2.9.0
    PyOpenGL: 3.1.5
    PyOpenGL-accelerate: 3.1.5
    pyparsing: 2.4.7
    PyQt5: 5.15.2
    PyQt5-sip: 12.8.1
    PyQtWebEngine: 5.15.2
    python-dateutil: 2.8.1
    pytz: 2021.1
    pyzmq: 22.2.1
    qtconsole: 5.1.0
    QtPy: 1.9.0
    RandomWords: 0.3.0
    requests: 2.25.1
    scipy: 1.6.3
    setuptools: 57.0.0
    sfftk-rw: 0.7.0.post1
    six: 1.16.0
    snowballstemmer: 2.1.0
    sortedcontainers: 2.4.0
    Sphinx: 4.0.1
    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.0.5
    urllib3: 1.26.6
    wcwidth: 0.2.5
    webcolors: 1.11.1
    wheel: 0.36.2
    wheel-filename: 1.3.0

Change History (2)

comment:1 by Tom Goddard, 4 years ago

Component: UnassignedStructure Analysis
Owner: set to Tom Goddard
Platform: all
Project: ChimeraX
Status: newassigned
Summary: ChimeraX bug report submissionAlphaFold second prediction overwrites first prediction

comment:2 by Tom Goddard, 4 years ago

Resolution: fixed
Status: assignedclosed

Fixed.

AlphaFold Run tool was a single instance and remembered the download directory from the first run. Made it forget the directory once the results are downloaded and shown.

Note: See TracTickets for help on using tickets.