Clipper: density sphere not aligned with structure

[jingjin.fan@…]

The following bug report has been submitted:
Platform:        Windows-10-10.0.19045
ChimeraX Version: 1.7.1 (2024-01-23 01:58:08 UTC)
Description
Dear Sir/Madam, there is an issue when I run ISOLDE. I fit the PDB to the map very well, and then I use clipper assoc #5 to #3(also tried to add the maps to working model), but I think the sphere of density doesn't align with the PDB very well. I tried to delet other distracted molecules like HOH, but doesn't help too much. I also tried to use different PDB(1u77, 1u7c, 1u7g,1xqe), the issue still exists. Could you help me with it? Thanks a lot in advanced!

Log:
Startup Messages  
---  
warning | Replacing fetcher for 'pdb_nmr' and format NMRSTAR from NMRSTAR
bundle with that from NMRSTAR bundle  
  
UCSF ChimeraX version: 1.7.1 (2024-01-23)  
© 2016-2023 Regents of the University of California. All rights reserved.  
How to cite UCSF ChimeraX  

> open 1u77

1u77 title:  
Crystal Structure of Ammonia Channel AmtB from E. Coli [more info...]  
  
Chain information for 1u77 #1  
---  
Chain | Description | UniProt  
A | Probable ammonium transporter | AMTB_ECOLI 1-385  
  
1u77 mmCIF Assemblies  
---  
1| author_and_software_defined_assembly  
  

> sym #1 assembly 1

Made 3 copies for 1u77 assembly 1  

> combine #2

Remapping chain ID 'A' in 1u77 #2.2 to 'B'  
Remapping chain ID 'A' in 1u77 #2.3 to 'C'  

> open
> Q:/__BACKUP__/PROJECTS/TEM/img/raw/0578/cryosparc_P109_J40_004_volume_map.mrc

Opened cryosparc_P109_J40_004_volume_map.mrc as #4, grid size 256,256,256,
pixel 0.85, shown at level 0.0789, step 1, values float32  

> hide #!2 models

> view #4 clip false

> select add #3

9126 atoms, 8589 bonds, 1878 residues, 1 model selected  

> ui mousemode right "translate selected models"

> view matrix models #3,1,0,0,11.599,0,1,0,11.232,0,0,1,47.213

> ui mousemode right "move picked models"

> ui mousemode right "translate selected models"

> view matrix models #3,1,0,0,17.807,0,1,0,3.8859,0,0,1,40.63

> ui mousemode right "rotate selected models"

> view matrix models
> #3,-0.23836,-0.26733,-0.93366,119.36,-0.95448,-0.11301,0.27603,68.057,-0.1793,0.95696,-0.22822,70.817

> volume flip #4

Opened cryosparc_P109_J40_004_volume_map.mrc z flip as #5, grid size
256,256,256, pixel 0.85, shown at step 1, values float32  

> select subtract #3

Nothing selected  

> select add #5

2 models selected  

> view matrix models
> #5,-0.5163,0.70336,0.48858,37.829,0.76469,0.12176,0.63279,-58.271,0.38559,0.70032,-0.60072,56.199

> select subtract #5

Nothing selected  

> select add #3

9126 atoms, 8589 bonds, 1878 residues, 1 model selected  

> view matrix models
> #3,-0.46807,-0.50519,-0.72504,128.45,-0.70119,-0.287,0.65265,46.275,-0.53781,0.81389,-0.2199,91.155

> ui mousemode right "move picked models"

> view matrix models
> #3,-0.46807,-0.50519,-0.72504,161.86,-0.70119,-0.287,0.65265,112.64,-0.53781,0.81389,-0.2199,142.87

> view matrix models
> #3,-0.46807,-0.50519,-0.72504,169.21,-0.70119,-0.287,0.65265,130.42,-0.53781,0.81389,-0.2199,125.81

> ui mousemode right "rotate selected models"

> view matrix models
> #3,-0.35641,-0.0017557,-0.93433,158.45,-0.9315,0.07851,0.35518,142.83,0.072731,0.99691,-0.029617,84.861

> ui tool show "Fit in Map"

> fitmap #3 inMap #5

Fit molecule combination (#3) to map cryosparc_P109_J40_004_volume_map.mrc z
flip (#5) using 9126 atoms  
average map value = 0.1157, steps = 124  
shifted from previous position = 10.2  
rotated from previous position = 15.8 degrees  
atoms outside contour = 2873, contour level = 0.078936  
  
Position of combination (#3) relative to cryosparc_P109_J40_004_volume_map.mrc
z flip (#5) coordinates:  
Matrix rotation and translation  
-0.38553584 0.24719127 0.88896489 87.02072801  
-0.46960714 0.77676146 -0.41965553 125.88615831  
-0.79424885 -0.57925651 -0.18338660 168.82570207  
Axis -0.08690821 0.91656763 -0.39032133  
Axis point 115.92265881 0.00000000 78.68783504  
Rotation angle (degrees) 113.33337857  
Shift along axis 41.92409068  
  

> select subtract #3

Nothing selected  

> delete :HOH

> isolde start

> set selectionWidth 4

Done loading forcefield  

> isolde set simFidelityMode Medium/Medium

ISOLDE: setting sim fidelity mode to Medium/Medium  
nonbonded_cutoff_distance = 0.900000  
use_gbsa = True  
gbsa_cutoff = 1.100000  
Removed all altlocs in #3 and reset associated occupancies to 1.  
ISOLDE: Checking and correcting nomenclature for (pseudo)symmetric side
chains...  
ISOLDE: Corrected atom nomenclature of 48 residues in model #3 to IUPAC-IUB
standards.  

> save Q:/__BACKUP__/PROJECTS/TEM/img/raw/0578/fitting_1u77.cxs

> view #1#3.1,3#2-5#3.2 clip false

> clipper associate #5 toModel #3

Opened cryosparc_P109_J40_004_volume_map.mrc z flip as #3.1.1.1, grid size
256,256,256, pixel 0.85, shown at step 1, values float32  

> ui tool show "Side View"

> isolde sim start #3.2/A-C

Sim termination reason: None  
ISOLDE: stopped sim  

> addh #3.2

Summary of feedback from adding hydrogens to combination #3.2  
---  
warning | Unknown hybridization for atom (SE) of residue type MSE; not adding
hydrogens to it  
notes | Termini for combination (#3.2) chain A determined from SEQRES records  
Termini for combination (#3.2) chain B determined from SEQRES records  
Termini for combination (#3.2) chain C determined from SEQRES records  
Chain-initial residues that are actual N termini:  
Chain-initial residues that are not actual N termini: combination #3.2/A ALA
3, combination #3.2/B ALA 3, combination #3.2/C ALA 3  
Chain-final residues that are actual C termini: combination #3.2/A VAL 385,
combination #3.2/B VAL 385, combination #3.2/C VAL 385  
Chain-final residues that are not actual C termini:  
Missing OXT added to C-terminal residue combination #3.2/A VAL 385  
Missing OXT added to C-terminal residue combination #3.2/B VAL 385  
Missing OXT added to C-terminal residue combination #3.2/C VAL 385  
1161 hydrogen bonds  
Adding 'H' to combination #3.2/A ALA 3  
Adding 'H' to combination #3.2/B ALA 3  
Adding 'H' to combination #3.2/C ALA 3  
8670 hydrogens added  
  

> isolde sim start #3.2/A-C

ISOLDE: started sim  




OpenGL version: 3.3.0 - Build 30.0.100.9864
OpenGL renderer: Intel(R) UHD Graphics 630
OpenGL vendor: Intel

Python: 3.11.2
Locale: en_GB.cp1252
Qt version: PyQt6 6.3.1, Qt 6.3.1
Qt runtime version: 6.3.2
Qt platform: windows

Manufacturer: Dell Inc.
Model: OptiPlex 7070
OS: Microsoft Windows 10 Enterprise (Build 19045)
Memory: 42,721,021,952
MaxProcessMemory: 137,438,953,344
CPU: 8 Intel(R) Core(TM) i7-9700 CPU @ 3.00GHz
OSLanguage: en-GB

Installed Packages:
    alabaster: 0.7.16
    appdirs: 1.4.4
    asttokens: 2.4.1
    Babel: 2.14.0
    backcall: 0.2.0
    beautifulsoup4: 4.11.2
    blockdiag: 3.0.0
    blosc2: 2.0.0
    build: 0.10.0
    certifi: 2023.11.17
    cftime: 1.6.3
    charset-normalizer: 3.3.2
    ChimeraX-AddCharge: 1.5.13
    ChimeraX-AddH: 2.2.5
    ChimeraX-AlignmentAlgorithms: 2.0.1
    ChimeraX-AlignmentHdrs: 3.4.1
    ChimeraX-AlignmentMatrices: 2.1
    ChimeraX-Alignments: 2.12.2
    ChimeraX-AlphaFold: 1.0
    ChimeraX-AltlocExplorer: 1.1.1
    ChimeraX-AmberInfo: 1.0
    ChimeraX-Arrays: 1.1
    ChimeraX-Atomic: 1.49.1
    ChimeraX-AtomicLibrary: 12.1.5
    ChimeraX-AtomSearch: 2.0.1
    ChimeraX-AxesPlanes: 2.3.2
    ChimeraX-BasicActions: 1.1.2
    ChimeraX-BILD: 1.0
    ChimeraX-BlastProtein: 2.1.2
    ChimeraX-BondRot: 2.0.4
    ChimeraX-BugReporter: 1.0.1
    ChimeraX-BuildStructure: 2.10.5
    ChimeraX-Bumps: 1.0
    ChimeraX-BundleBuilder: 1.2.2
    ChimeraX-ButtonPanel: 1.0.1
    ChimeraX-CageBuilder: 1.0.1
    ChimeraX-CellPack: 1.0
    ChimeraX-Centroids: 1.3.2
    ChimeraX-ChangeChains: 1.1
    ChimeraX-CheckWaters: 1.3.2
    ChimeraX-ChemGroup: 2.0.1
    ChimeraX-Clashes: 2.2.4
    ChimeraX-Clipper: 0.22.4
    ChimeraX-ColorActions: 1.0.3
    ChimeraX-ColorGlobe: 1.0
    ChimeraX-ColorKey: 1.5.5
    ChimeraX-CommandLine: 1.2.5
    ChimeraX-ConnectStructure: 2.0.1
    ChimeraX-Contacts: 1.0.1
    ChimeraX-Core: 1.7.1
    ChimeraX-CoreFormats: 1.2
    ChimeraX-coulombic: 1.4.2
    ChimeraX-Crosslinks: 1.0
    ChimeraX-Crystal: 1.0
    ChimeraX-CrystalContacts: 1.0.1
    ChimeraX-DataFormats: 1.2.3
    ChimeraX-Dicom: 1.2
    ChimeraX-DistMonitor: 1.4
    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.1.1
    ChimeraX-Hbonds: 2.4
    ChimeraX-Help: 1.2.2
    ChimeraX-HKCage: 1.3
    ChimeraX-IHM: 1.1
    ChimeraX-ImageFormats: 1.2
    ChimeraX-IMOD: 1.0
    ChimeraX-IO: 1.0.1
    ChimeraX-ISOLDE: 1.7.1
    ChimeraX-ItemsInspection: 1.0.1
    ChimeraX-IUPAC: 1.0
    ChimeraX-Label: 1.1.8
    ChimeraX-ListInfo: 1.2.2
    ChimeraX-Log: 1.1.6
    ChimeraX-LookingGlass: 1.1
    ChimeraX-Maestro: 1.9.1
    ChimeraX-Map: 1.1.4
    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.2
    ChimeraX-MCopy: 1.0
    ChimeraX-MDcrds: 2.6.1
    ChimeraX-MedicalToolbar: 1.0.2
    ChimeraX-Meeting: 1.0.1
    ChimeraX-MLP: 1.1.1
    ChimeraX-mmCIF: 2.12.1
    ChimeraX-MMTF: 2.2
    ChimeraX-Modeller: 1.5.14
    ChimeraX-ModelPanel: 1.4
    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-Neuron: 1.0
    ChimeraX-Nifti: 1.1
    ChimeraX-NRRD: 1.1
    ChimeraX-Nucleotides: 2.0.3
    ChimeraX-OpenCommand: 1.13.1
    ChimeraX-PDB: 2.7.3
    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
    ChimeraX-PubChem: 2.1
    ChimeraX-ReadPbonds: 1.0.1
    ChimeraX-Registration: 1.1.2
    ChimeraX-RemoteControl: 1.0
    ChimeraX-RenderByAttr: 1.1
    ChimeraX-RenumberResidues: 1.1
    ChimeraX-ResidueFit: 1.0.1
    ChimeraX-RestServer: 1.2
    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-SelInspector: 1.0
    ChimeraX-SeqView: 2.11
    ChimeraX-Shape: 1.0.1
    ChimeraX-Shell: 1.0.1
    ChimeraX-Shortcuts: 1.1.1
    ChimeraX-ShowSequences: 1.0.2
    ChimeraX-SideView: 1.0.1
    ChimeraX-Smiles: 2.1.2
    ChimeraX-SmoothLines: 1.0
    ChimeraX-SpaceNavigator: 1.0
    ChimeraX-StdCommands: 1.12.4
    ChimeraX-STL: 1.0.1
    ChimeraX-Storm: 1.0
    ChimeraX-StructMeasure: 1.1.2
    ChimeraX-Struts: 1.0.1
    ChimeraX-Surface: 1.0.1
    ChimeraX-SwapAA: 2.0.1
    ChimeraX-SwapRes: 2.2.2
    ChimeraX-TapeMeasure: 1.0
    ChimeraX-TaskManager: 1.0
    ChimeraX-Test: 1.0
    ChimeraX-Toolbar: 1.1.2
    ChimeraX-ToolshedUtils: 1.2.4
    ChimeraX-Topography: 1.0
    ChimeraX-ToQuest: 1.0
    ChimeraX-Tug: 1.0.1
    ChimeraX-UI: 1.33.3
    ChimeraX-uniprot: 2.3
    ChimeraX-UnitCell: 1.0.1
    ChimeraX-ViewDockX: 1.3.2
    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.3
    ChimeraX-Zone: 1.0.1
    colorama: 0.4.6
    comm: 0.2.1
    comtypes: 1.1.14
    contourpy: 1.2.0
    cxservices: 1.2.2
    cycler: 0.12.1
    Cython: 0.29.33
    debugpy: 1.8.0
    decorator: 5.1.1
    docutils: 0.19
    executing: 2.0.1
    filelock: 3.9.0
    fonttools: 4.47.2
    funcparserlib: 2.0.0a0
    glfw: 2.6.4
    grako: 3.16.5
    h5py: 3.10.0
    html2text: 2020.1.16
    idna: 3.6
    ihm: 0.38
    imagecodecs: 2023.9.18
    imagesize: 1.4.1
    ipykernel: 6.23.2
    ipython: 8.14.0
    ipython-genutils: 0.2.0
    ipywidgets: 8.1.1
    jedi: 0.18.2
    Jinja2: 3.1.2
    jupyter-client: 8.2.0
    jupyter-core: 5.7.1
    jupyterlab-widgets: 3.0.9
    kiwisolver: 1.4.5
    line-profiler: 4.0.2
    lxml: 4.9.2
    lz4: 4.3.2
    MarkupSafe: 2.1.4
    matplotlib: 3.7.2
    matplotlib-inline: 0.1.6
    msgpack: 1.0.4
    nest-asyncio: 1.6.0
    netCDF4: 1.6.2
    networkx: 3.1
    nibabel: 5.0.1
    nptyping: 2.5.0
    numexpr: 2.8.8
    numpy: 1.25.1
    openvr: 1.23.701
    packaging: 23.2
    ParmEd: 3.4.3
    parso: 0.8.3
    pep517: 0.13.0
    pickleshare: 0.7.5
    pillow: 10.2.0
    pip: 23.0
    pkginfo: 1.9.6
    platformdirs: 4.1.0
    prompt-toolkit: 3.0.43
    psutil: 5.9.5
    pure-eval: 0.2.2
    py-cpuinfo: 9.0.0
    pycollada: 0.7.2
    pydicom: 2.3.0
    Pygments: 2.16.1
    pynrrd: 1.0.0
    PyOpenGL: 3.1.7
    PyOpenGL-accelerate: 3.1.7
    pyopenxr: 1.0.2801
    pyparsing: 3.0.9
    pyproject-hooks: 1.0.0
    PyQt6-commercial: 6.3.1
    PyQt6-Qt6: 6.3.2
    PyQt6-sip: 13.4.0
    PyQt6-WebEngine-commercial: 6.3.1
    PyQt6-WebEngine-Qt6: 6.3.2
    python-dateutil: 2.8.2
    pytz: 2023.3.post1
    pywin32: 305
    pyzmq: 25.1.2
    qtconsole: 5.4.3
    QtPy: 2.4.1
    RandomWords: 0.4.0
    requests: 2.31.0
    scipy: 1.11.1
    setuptools: 67.4.0
    sfftk-rw: 0.7.3
    six: 1.16.0
    snowballstemmer: 2.2.0
    sortedcontainers: 2.4.0
    soupsieve: 2.5
    sphinx: 6.1.3
    sphinx-autodoc-typehints: 1.22
    sphinxcontrib-applehelp: 1.0.8
    sphinxcontrib-blockdiag: 3.0.0
    sphinxcontrib-devhelp: 1.0.6
    sphinxcontrib-htmlhelp: 2.0.5
    sphinxcontrib-jsmath: 1.0.1
    sphinxcontrib-qthelp: 1.0.7
    sphinxcontrib-serializinghtml: 1.1.10
    stack-data: 0.6.3
    superqt: 0.5.0
    tables: 3.8.0
    tcia-utils: 1.5.1
    tifffile: 2023.7.18
    tinyarray: 1.2.4
    tomli: 2.0.1
    tornado: 6.4
    traitlets: 5.9.0
    typing-extensions: 4.9.0
    tzdata: 2023.4
    urllib3: 2.1.0
    wcwidth: 0.2.13
    webcolors: 1.12
    wheel: 0.38.4
    wheel-filename: 1.4.1
    widgetsnbextension: 4.0.9
    WMI: 1.5.1

[pett]

Reported by Jingjin Fan

[Tristan Croll]

This is a current limitation of Clipper: the map must be in the absolute
reference frame (i.e. not moved from its starting position after opening).
In other words, you should always move the model to fit the map, and not
the other way around.

On Tue, Sep 3, 2024 at 5:16 PM ChimeraX <ChimeraX-bugs-admin@cgl.ucsf.edu>
wrote:

>
>
>
>