October 2022 - ChimeraX-users - RBVI Mailing Lists

Fit in map, save new coords of map
by Nebojsa Bogdanovic 10 Oct '22

10 Oct '22

Hello, I open two maps and I fit #2 into #1. I need to save map #2 at it's new position, but there is no option "relative to model or map". By doing "volume resample #1 onGrid #1, I match the box sizes but the new position is not maintained after saving that new map. Is there any command/function I could do that with? Thank You and regards, Nebojša Bogdanović, Ph.D Department of Physiology and Biophysics, School of Medicine, Case Western Reserve University

3 3

Re: [chimerax-users] alphafold contacts between 3 and more subunits
by Tom Goddard 10 Oct '22

10 Oct '22

Hi André, I added two options "replace" (default true) and "name" (default "PAE Contacts") to the alphafold contacts command. The replace option determines if existing contact pseudobonds (in the group with the specified or default name) are deleted. If replace is false the new contacts are just added to any existing ones. The name option gives the pseudobond model the specified name. So you could use alphafold contacts /A to /B alphafold contacts /B to /C replace false alphafold contacts /A to /C replace false to put all the contacts in one pseudobond submodel. If instead you wanted them in 3 different models you could use alphafold contacts /A to /B name "AB contacts" alphafold contacts /B to /C name "BC contacts" alphafold contacts /A to /C name "AC contacts" This will be in tomorrow's (Oct 11) ChimeraX daily builds. Tom > On Oct 10, 2022, at 12:38 PM, Michaelis, Andre Clemens <michaelis(a)biochem.mpg.de> wrote: > > Awesome, > Thank you so much! > > Best > André > > > —————————————————— > > >> Am 10.10.2022 um 20:50 schrieb Tom Goddard <goddard(a)sonic.net>: >> >> Hi André, >> >> The "alphafold contacts" command creates a submodel of the atomic structure containing pseudobonds with the submodel named "PAE Contacts". If you rename that submodel then it won't replace it when you do the second alphafold contacts command on the same structure. For instance, if the structure is model #1. >> >> alphafold contacts /A to /B >> rename #1.1 "AB contacts" >> alphafold contacts /B to /C >> rename #1.2 "BC contacts" >> alphafold contacts /A to /C >> >> This should be easier to do. I should add an option to the alphafold contacts command "replace false" that automatically uses new names for each group of pseudobonds. I've made a request for this. >> >> https://www.rbvi.ucsf.edu/trac/ChimeraX/ticket/7770 <https://www.rbvi.ucsf.edu/trac/ChimeraX/ticket/7770> >> >> Tom >> >>> On Oct 10, 2022, at 5:47 AM, Michaelis, Andre Clemens via ChimeraX-users <chimerax-users(a)cgl.ucsf.edu <mailto:chimerax-users@cgl.ucsf.edu>> wrote: >>> >>> Hi, >>> >>> Is it possible to draw contacts (PAE, alphafold) between more than just 2 subunits? >>> e.g.: >>> alphafold contacts /A to /B >>> AND >>> alphafold contacts /B to /C >>> AND >>> alphafold contacts /A to /C >>> >>> Righ now, a new command deletes the old contacts. Not sure about a syntax that would combine/add those >>> >>> Thank you very much! >>> >>> Best André >>> >>> >>> _______________________________________________ >>> ChimeraX-users mailing list >>> ChimeraX-users(a)cgl.ucsf.edu <mailto:ChimeraX-users@cgl.ucsf.edu> >>> Manage subscription: >>> https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users <https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users>

1 0

alphafold contacts between 3 and more subunits
by Michaelis, Andre Clemens 10 Oct '22

10 Oct '22

Hi, Is it possible to draw contacts (PAE, alphafold) between more than just 2 subunits? e.g.: alphafold contacts /A to /B AND alphafold contacts /B to /C AND alphafold contacts /A to /C Righ now, a new command deletes the old contacts. Not sure about a syntax that would combine/add those Thank you very much! Best André

2 1

mrc to pdb
by Tianming Qu 07 Oct '22

07 Oct '22

Hi, is there any way I can convert a .mrc file to .pdb file in chimera? Best, Tianming

5 5

A question about the level option in volume viewer
by Qu Tianming 06 Oct '22

06 Oct '22

Hi ChimeraX-developers, I got a question about the level option in volume viewer. I think the level represent the isovalue for the density map. So when I set the level as 0.03, I got a density map as following : [?? ???????] When I set the level as 0.01, the density map as following: [?? ???????] So my question is that, does level 0.03 or 0.01 means, all grids with a density value lower than 0.03 or 0.01 will be present? Means when I set the level as 0.01, everything left in the density map is the density value lower than 0.01? Best, Tianming

2 1

Moving molecules in VR
by Dr. Eddie 06 Oct '22

06 Oct '22

Hello, I'm trying to move some waters away from a ligand (and eventually some lipids as well) as they overlap too much for minimization to fix it. I can't seem to figure out how to select an entire (water) molecule. The select atom lets me select one atom, but I cannot chain or add more atoms. Some atoms seem to refuse to move even when I try to move one at a time using the "Move atoms" right mouse option. Any help would be appreciated as I can't find much about translating molecules (the bond rotation works perfectly!). Thanks! Eddie

3 2

Problems with Qt6
by Hoover, David (NIH/CIT) [E] 06 Oct '22

06 Oct '22

Hi, I'm trying to create a Ubuntu 22.04 container for ChimeraX v1.4 and no matter what I do I get this error message: NOTE: Traceback (most recent call last): File "/usr/lib/ucsf-chimerax/lib/python3.9/runpy.py", line 197, in _run_module_as_main return _run_code(code, main_globals, None, File "/usr/lib/ucsf-chimerax/lib/python3.9/runpy.py", line 87, in _run_code exec(code, run_globals) File "/usr/lib/ucsf-chimerax/lib/python3.9/site-packages/ChimeraX_main.py", line 1021, in <module> exit_code = init(sys.argv) File "/usr/lib/ucsf-chimerax/lib/python3.9/site-packages/ChimeraX_main.py", line 573, in init from chimerax.ui import initialize_qt File "/usr/lib/ucsf-chimerax/lib/python3.9/site- packages/chimerax/ui/__init__.py", line 21, in <module> from .gui import MainToolWindow, initialize_qt, menu_capitalize File "/usr/lib/ucsf-chimerax/lib/python3.9/site- packages/chimerax/ui/gui.py", line 130, in <module> from Qt.QtWidgets import QApplication File "/usr/lib/ucsf-chimerax/lib/python3.9/site-packages/Qt/__init__.py", line 64, in <module> from PyQt6.QtCore import PYQT_VERSION_STR as PYQT6_VERSION ImportError: libQt6Core.so.6: cannot open shared object file: No such file or directory BUG: ImportError: libQt6Core.so.6: cannot open shared object file: No such file or directory File "/usr/lib/ucsf-chimerax/lib/python3.9/site-packages/Qt/__init__.py", line 64, in from PyQt6.QtCore import PYQT_VERSION_STR as PYQT6_VERSION _See log for complete Python traceback._ I am attempting to install into Ubuntu 22.04 with ever Qt6 package available, but nothing seems to help. The library in question definitely exists at /usr/lib/x86_64-linux-gnu/libQt6Core.so.6, along with a multitude of other libQt6* libraries. Is there something additional that needs to be done or installed to enable ChimeraX v1.4 to run? I have tried export LD_LIBRARY_PATH=/usr/lib/x86_64-linux-gnu:$LD_LIBRARY_PATH but this does not help. Any insight would be appreciated. David -- David Hoover, Ph.D. Computational Biologist High Performance Computing Services, Center for Information Technology, National Institutes of Health 12 South Dr., Rm 2N207 Bethesda, MD 20892, USA TEL: (+1) 301-435-2986 Email: hooverdm(a)hpc.nih.gov

2 1

Re: [chimerax-users] usage for multimer
by Tom Goddard 05 Oct '22

05 Oct '22

Hi Dennis, Are you asking about fetching AlphaFold database models? Those are only of monomer proteins. With AlphaFold predictions energy minimization is done using OpenMM and the Amber force field -- it has nothing to do with binding affinity. I don't know what you mean by a contact map. Maybe you are interested in the ChimeraX interfaces tool that shows contact residues at an interface. Please give enough details in your questions so they are clear, otherwise we will not be able to spend the time to answer. Tom > On Oct 5, 2022, at 7:41 PM, Dennis Poh <pohdennis90(a)gmail.com> wrote: > > Hi Tom > > Thanks, also in the chimerax GUI software, I can fetch my separated sequences delimited with a comma, but does the overall structure that appears the most energy-minimized model or just a superimposition simply based on sequence similarity? > > How is the energy-minimization done (is it based on deltaG = G(bound complex) - G(protein1) -G(protein2)) ? > Could a contact map be produced (if 2 sequences are compared pairwise) too? > > -Dennis > > On Thu, 6 Oct 2022 at 01:43, Tom Goddard <goddard(a)sonic.net <mailto:goddard@sonic.net>> wrote: > Hi Dennis, > > AlphaFold on Google Colab will only work for total sequence length of about 1000 because the old Google Colab GPUs only have 16 Gbytes of memory. Modern GPUs for machine learning like an Nvidia A40 (48 GB) or A100 (40 or 80GB) or A6000 (48 GB) can handle more than 3000 residues but those GPUs cost $5000-$15000. A consumer Nvidia RTX 3090 (24 GB, $1500) can manage up to 2000 residues. Setting up AlphaFold on your own machine is a bit of work. I see from your output file name "af1848" you tried to run total sequence length 1848 which certainly will fail on Google Colab. For examples of AlphaFold runs at different sequence lengths see > > https://www.rbvi.ucsf.edu/chimerax/data/alphafold-jan2022/afspeed.html <https://www.rbvi.ucsf.edu/chimerax/data/alphafold-jan2022/afspeed.html> > > Also AlphaFold has a minimum sequence length for each individual sequence is 16 so your example ACCCC and ALLPAAAA will not work. This minimum length is because AlphaFold uses a sequence alignment made by searching a billion database sequences and searching shorter sequence produces millions of hits that are not evolutionarily related. > > The sequence length limits are described in the ChimeraX documentation > > https://www.rbvi.ucsf.edu/chimerax/docs/user/commands/alphafold.html#caveats <https://www.rbvi.ucsf.edu/chimerax/docs/user/commands/alphafold.html#caveats> > > Tom > > >> On Oct 5, 2022, at 9:35 AM, Dennis Poh via ChimeraX-users <chimerax-users(a)cgl.ucsf.edu <mailto:chimerax-users@cgl.ucsf.edu>> wrote: >> >> Hi Eric >> >> 1. It seems that I don't have sufficient GPU memory in colab. >> >> I have the following error message: >> INFO:colabfold.batch:Running model_1 >> /usr/local/lib/python3.7/dist-packages/haiku/_src/data_structures.py:195: FutureWarning: jax.tree_flatten is deprecated, and will be removed in a future release. Use jax.tree_util.tree_flatten instead. >> leaves, structure = jax.tree_flatten(mapping) >> /usr/local/lib/python3.7/dist-packages/haiku/_src/data_structures.py:203: FutureWarning: jax.tree_unflatten is deprecated, and will be removed in a future release. Use jax.tree_util.tree_unflatten instead. >> self._mapping = jax.tree_unflatten(self._structure, self._leaves) >> /usr/local/lib/python3.7/dist-packages/haiku/_src/stateful.py:457: FutureWarning: jax.tree_leaves is deprecated, and will be removed in a future release. Use jax.tree_util.tree_leaves instead. >> length = jax.tree_leaves(xs)[0].shape[0] >> /usr/local/lib/python3.7/dist-packages/alphafold/model/geometry/struct_of_array.py:136: FutureWarning: jax.tree_flatten is deprecated, and will be removed in a future release. Use jax.tree_util.tree_flatten instead. >> flat_array_like, inner_treedef = jax.tree_flatten(array_like) >> /usr/local/lib/python3.7/dist-packages/alphafold/model/geometry/struct_of_array.py:210: FutureWarning: jax.tree_unflatten is deprecated, and will be removed in a future release. Use jax.tree_util.tree_unflatten instead. >> inner_treedef, data[array_start:array_start + num_array]) >> /usr/local/lib/python3.7/dist-packages/alphafold/model/mapping.py:50: FutureWarning: jax.tree_flatten is deprecated, and will be removed in a future release. Use jax.tree_util.tree_flatten instead. >> values_tree_def = jax.tree_flatten(values)[1] >> /usr/local/lib/python3.7/dist-packages/alphafold/model/mapping.py:54: FutureWarning: jax.tree_unflatten is deprecated, and will be removed in a future release. Use jax.tree_util.tree_unflatten instead. >> return jax.tree_unflatten(values_tree_def, flat_axes) >> /usr/local/lib/python3.7/dist-packages/alphafold/model/mapping.py:129: FutureWarning: jax.tree_flatten is deprecated, and will be removed in a future release. Use jax.tree_util.tree_flatten instead. >> flat_sizes = jax.tree_flatten(in_sizes)[0] >> ERROR:colabfold.batch:Could not predict af1819. Not Enough GPU memory? INTERNAL: cublas error >> INFO:colabfold.batch:Done >> Downloading structure predictions to directory Downloads/ChimeraX/AlphaFold >> cp: cannot stat '*_relaxed_rank_1_model_*.pdb': No such file or directory >> cp: cannot stat '*_unrelaxed_rank_1_model_*_scores.json': No such file or directory >> >> 2. Would it be possible to run this in jupyter? Or are there alternatives? >> >> 3. " Prediction may fail with total sequence length over 1000 residues due to limited GPU memory." - this total sequence length meaning all the sequences in the list to be concatenated? >> >> 4. I seem to also have some pdbxx.m8 and afxxxx.csv files - may I know what these files are for? >> >> Thanks! >> -Dennis >> >> On Thu, 6 Oct 2022 at 00:25, Eric Pettersen <pett(a)cgl.ucsf.edu <mailto:pett@cgl.ucsf.edu>> wrote: >> Also, if you are copying the sequence from a text editor, make sure the editor is displaying the sequence as plain text or that you are copying as plain text, otherwise invisible formatting characters may be embedded in what you paste. >> >> --Eric >> >> Eric Pettersen >> UCSF Computer Graphics Lab >> >> >> > On Oct 5, 2022, at 9:14 AM, Elaine Meng via ChimeraX-users <chimerax-users(a)cgl.ucsf.edu <mailto:chimerax-users@cgl.ucsf.edu>> wrote: >> > >> > Well, the message says there is an invalid character, so all I can say is to make sure that you are pasting plain text, and check to see that you have only standard amino acid codes and commas. >> > >> > WIthout seeing exactly what you pasted, we can't tell which part caused the problem. >> > >> > Elaine >> > >> > >> >> On Oct 5, 2022, at 9:09 AM, Dennis Poh <pohdennis90(a)gmail.com <mailto:pohdennis90@gmail.com>> wrote: >> >> >> >> Hi Elaine, >> >> >> >> Thanks, I was trying to predict a multimer or the overall structure of many subunit chains using individual sequences, each separated with a comma in colab. >> >> But it seems that there was some error and no pdb file was generated; the error message is as follows: >> >> >> >> ERROR:colabfold.batch:Could not generate input features af1848: Invalid character in the sequence: >> >> Traceback (most recent call last): >> >> File "/usr/local/lib/python3.7/dist-packages/colabfold/batch.py", line 1357, in run >> >> model_type, >> >> File "<ipython-input-1-d6881d38b934>", line 122, in generate_input_feature_wrapper >> >> (input_features, domain_names) = batch.generate_input_feature_orig(*args, **kw) >> >> File "/usr/local/lib/python3.7/dist-packages/colabfold/batch.py", line 1018, in generate_input_feature >> >> sequence, input_msa, template_features[sequence_index] >> >> File "/usr/local/lib/python3.7/dist-packages/colabfold/batch.py", line 869, in build_monomer_feature >> >> sequence=sequence, description="none", num_res=len(sequence) >> >> File "/usr/local/lib/python3.7/dist-packages/alphafold/data/pipeline.py", line 43, in make_sequence_features >> >> map_unknown_to_x=True) >> >> File "/usr/local/lib/python3.7/dist-packages/alphafold/common/residue_constants.py", line 580, in sequence_to_onehot >> >> raise ValueError(f'Invalid character in the sequence: {aa_type}') >> >> ValueError: Invalid character in the sequence: >> >> INFO:colabfold.batch:Done >> >> Downloading structure predictions to directory Downloads/ChimeraX/AlphaFold >> >> cp: cannot stat '*_relaxed_rank_1_model_*.pdb': No such file or directory >> >> cp: cannot stat '*_unrelaxed_rank_1_model_*_scores.json': No such file or directory >> >> >> >> -Dennis >> >> >> >> On Wed, 5 Oct 2022 at 23:36, Elaine Meng <meng(a)cgl.ucsf.edu <mailto:meng@cgl.ucsf.edu>> wrote: >> >> Hi Dennis, >> >> Your sequence input is wrong - it should contain only the sequences pasted as plain text, with only a comma between them (NOT the ">description" line because it is not supposed to be in fasta format). How to input sequence(s) is explained in the AlphaFold help page. >> >> >> >> <https://rbvi.ucsf.edu/chimerax/docs/user/tools/alphafold.html#predict <https://rbvi.ucsf.edu/chimerax/docs/user/tools/alphafold.html#predict>> >> >> >> >> "For predicting a complex (multimer), the sequences of all chains in the complex must be given. The same sequence must be given multiple times if it occurs in multiple copies in the complex. The sequences can be specified either collectively as a model number chosen from the menu of currently open models (e.g. when that model contains multiple chains), or individually within a comma-separated list of UniProt identifiers or pasted-in amino acid sequences." >> >> >> >> E.g. something like >> >> >> >> ACCCC,ALLPAAAA >> >> >> >> I hope this helps, >> >> Elaine >> >> ----- >> >> Elaine C. Meng, Ph.D. >> >> UCSF Chimera(X) team >> >> Department of Pharmaceutical Chemistry >> >> University of California, San Francisco >> >> >> >>> On Oct 5, 2022, at 4:12 AM, Dennis Poh via ChimeraX-users <chimerax-users(a)cgl.ucsf.edu <mailto:chimerax-users@cgl.ucsf.edu>> wrote: >> >>> >> >>> Sorry, I couldn't generate the pdb from chimera colab - is there anything that i may have missed? >> >>> >> >>> -Dennis >> >>> >> >>> On Wed, 5 Oct 2022 at 17:11, Dennis Poh <pohdennis90(a)gmail.com <mailto:pohdennis90@gmail.com>> wrote: >> >>> Hi >> >>> >> >>> I encounter a problem in multimer prediction with the sequences I use as input. >> >>> It's always indicating: >> >>> "Missing or invalid "sequences" argument: Sequences argument" >> >>> and " is not a chain specifier, alignment id, UniProt id, or sequence characters" >> >>> >> >>> The input format is always something like a fasta format: >> >>>> seq_id >> >>> ACCCC >> >>> >> >>>> seq_id2 >> >>> ALLPAAAA >> >>> >> >>> May I know how I can rectify this? >> >>> >> >>> Thanks! >> >>> - Dennis >> >> >> > >> > >> > _______________________________________________ >> > ChimeraX-users mailing list >> > ChimeraX-users(a)cgl.ucsf.edu <mailto:ChimeraX-users@cgl.ucsf.edu> >> > Manage subscription: >> > https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users <https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users> >> >> _______________________________________________ >> ChimeraX-users mailing list >> ChimeraX-users(a)cgl.ucsf.edu <mailto:ChimeraX-users@cgl.ucsf.edu> >> Manage subscription: >> https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users <https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users> >

1 0

combining executable html commands
by Krishnan Raman 05 Oct '22

05 Oct '22

I have large presentation with several html executable commands . is there a way I can run them with one click single command for all of them? thanks Krishnan Raman Executive Director, Structural Biology BioCryst Pharmaceuticals CONFIDENTIALITY NOTICE This email, including any attachments, may contain confidential or legally privileged information that is intended only for the individual or entity to whom it is addressed. If you are not the intended recipient, please be advised that any dissemination, distribution or copying of this email and any attachment is strictly prohibited. If you have received this email in error, please reply to the sender so that BioCryst Pharmaceuticals, Inc. can take corrective measures and then permanently delete this email and any attachment, including any printed copies. Thank you.

2 1

usage for multimer
by Dennis Poh 05 Oct '22

05 Oct '22

Hi I encounter a problem in multimer prediction with the sequences I use as input. It's always indicating: "Missing or invalid "sequences" argument: Sequences argument" and " is not a chain specifier, alignment id, UniProt id, or sequence characters" The input format is always something like a fasta format: >seq_id ACCCC >seq_id2 ALLPAAAA May I know how I can rectify this? Thanks! - Dennis

4 7