<div dir="ltr">Hi Elaine,<div><br></div><div>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.</div><div>But it seems that there was some error and no pdb file was generated; the error message is as follows:</div><div><br></div><div>ERROR:colabfold.batch:Could not generate input features af1848: Invalid character in the sequence:  <br>Traceback (most recent call last):<br>  File "/usr/local/lib/python3.7/dist-packages/colabfold/batch.py", line 1357, in run<br>    model_type,<br>  File "<ipython-input-1-d6881d38b934>", line 122, in generate_input_feature_wrapper<br>    (input_features, domain_names) = batch.generate_input_feature_orig(*args, **kw)<br>  File "/usr/local/lib/python3.7/dist-packages/colabfold/batch.py", line 1018, in generate_input_feature<br>    sequence, input_msa, template_features[sequence_index]<br>  File "/usr/local/lib/python3.7/dist-packages/colabfold/batch.py", line 869, in build_monomer_feature<br>    sequence=sequence, description="none", num_res=len(sequence)<br>  File "/usr/local/lib/python3.7/dist-packages/alphafold/data/pipeline.py", line 43, in make_sequence_features<br>    map_unknown_to_x=True)<br>  File "/usr/local/lib/python3.7/dist-packages/alphafold/common/residue_constants.py", line 580, in sequence_to_onehot<br>    raise ValueError(f'Invalid character in the sequence: {aa_type}')<br>ValueError: Invalid character in the sequence:  <br>INFO:colabfold.batch:Done<br>Downloading structure predictions to directory Downloads/ChimeraX/AlphaFold<br>cp: cannot stat '*_relaxed_rank_1_model_*.pdb': No such file or directory<br>cp: cannot stat '*_unrelaxed_rank_1_model_*_scores.json': No such file or directory<br></div><div><div><br></div><div>-Dennis</div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, 5 Oct 2022 at 23:36, Elaine Meng <<a href="mailto:meng@cgl.ucsf.edu">meng@cgl.ucsf.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hi Dennis,<br>
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.<br>
<br>
<<a href="https://rbvi.ucsf.edu/chimerax/docs/user/tools/alphafold.html#predict" rel="noreferrer" target="_blank">https://rbvi.ucsf.edu/chimerax/docs/user/tools/alphafold.html#predict</a>><br>
<br>
"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."<br>
<br>
E.g. something like<br>
<br>
ACCCC,ALLPAAAA<br>
<br>
I hope this helps,<br>
Elaine<br>
-----<br>
Elaine C. Meng, Ph.D.                       <br>
UCSF Chimera(X) team<br>
Department of Pharmaceutical Chemistry<br>
University of California, San Francisco<br>
<br>
> On Oct 5, 2022, at 4:12 AM, Dennis Poh via ChimeraX-users <<a href="mailto:chimerax-users@cgl.ucsf.edu" target="_blank">chimerax-users@cgl.ucsf.edu</a>> wrote:<br>
> <br>
> Sorry, I couldn't generate the pdb from chimera colab - is there anything that i may have missed?<br>
> <br>
> -Dennis<br>
> <br>
> On Wed, 5 Oct 2022 at 17:11, Dennis Poh <<a href="mailto:pohdennis90@gmail.com" target="_blank">pohdennis90@gmail.com</a>> wrote:<br>
> Hi <br>
> <br>
> I encounter a problem in multimer prediction with the sequences I use as input.<br>
> It's always indicating:<br>
> "Missing or invalid "sequences" argument: Sequences argument"<br>
> and " is not a chain specifier, alignment id, UniProt id, or sequence characters"<br>
> <br>
> The input format is always something like a fasta format:<br>
> >seq_id<br>
> ACCCC<br>
> <br>
> >seq_id2<br>
> ALLPAAAA<br>
> <br>
> May I know how I can rectify this?<br>
> <br>
> Thanks!<br>
> - Dennis<br>
<br>
</blockquote></div>