Hi Ben,

  It would be nice to show templates and sequence alignments used for predicted models from AlphaFold and Modeller.  We could output an html table in the log that lists the templates with a link to show the sequence alignment and a link to load and align the template if it is from the PDB.  

  The AlphaFold models in the EBI AlphaFold database don't appear to say what template structures were used, for instance, I looked at AF-P12004-F1-model_v1.cif

https://alphafold.ebi.ac.uk/entry/P12004

I believe AlphaFold2 finds the 20 best matching structures in the PDB and uses 4 (not sure how they are selected).  I've run AlphaFold many times and the log output says what the 20 matches are but does not appear to say which 4 structures it actually used -- pretty unfortunate.  The AlphaFold per-residue confidence scores are in an mmCIF table _ma_qa_metric_local:

#
loop_
_ma_qa_metric_local.label_asym_id
_ma_qa_metric_local.label_comp_id
_ma_qa_metric_local.label_seq_id
_ma_qa_metric_local.metric_id
_ma_qa_metric_local.metric_value
_ma_qa_metric_local.model_id
_ma_qa_metric_local.ordinal_id
A MET 1   2 91.95 1 1 
A PHE 2   2 96.89 1 1 
A GLU 3   2 98.01 1 1 
A ALA 4   2 98.08 1 1 
A ARG 5   2 97.76 1 1 
A LEU 6   2 96.16 1 1 
..

Currently ChimeraX colors AlphaFold models by confidence using the same scores taken from the bfactor column of the atom site table.

  The Model Archive example you gave as an example (https://www.modelarchive.org/doi/10.5452/ma-bak-cepc-0250) has no templates sequences or alignments in the mmCIF file, and no per-residue scores, but does have some global scores.

  Your ModBase example (https://modbase.compbio.ucsf.edu/modbase-cgi/model_search.cgi?databaseID=Q12321) does have a template sequence and alignment and global scores but no per-residue scores

#
loop_
_ma_template_ref_db_details.template_id
_ma_template_ref_db_details.db_name
_ma_template_ref_db_details.db_accession_code
1 PDB 3nc1

#
loop_
_ma_template_poly.template_id
_ma_template_poly.seq_one_letter_code
_ma_template_poly.seq_one_letter_code_can
1 DMACDTFIKIAQKCRRHFVQVQVGEVMPFIDEILNNINTIICDLQPQQVHTFYEAVGYMIGAQTDQTVQEHLIEKYMLLPNQVWDSIIQQATKNVDILKDPETVKQLGSILKTNVRACKAVGHPFVIQLGRIYLDMLNVYKCLSENISAAIQANGEMVTKQPLIRSMRTVKRETLKLISGWVSRSNDPQMVAENFVPPLLDAVLIDYQRNVPAAREPEVLSTMAIIVNKLGGHITAEIPQIFDAVFECTLNMINKDFEEYPEHRTNFFLLLQAVNSHCFPAFLAIPPAQFKLVLDSIIWAFKHTMRNVADTGLQILFTLLQNVAQEEAAAQSFYQTYFCDILQHIFSVVTDTSHTAGLTMHASILAYMFNLVEEGKISTPLNPN DMACDTFIKIAQKCRRHFVQVQVGEVMPFIDEILNNINTIICDLQPQQVHTFYEAVGYMIGAQTDQTVQEHLIEKYMLLPNQVWDSIIQQATKNVDILKDPETVKQLGSILKTNVRACKAVGHPFVIQLGRIYLDMLNVYKCLSENISAAIQANGEMVTKQPLIRSMRTVKRETLKLISGWVSRSNDPQMVAENFVPPLLDAVLIDYQRNVPAAREPEVLSTMAIIVNKLGGHITAEIPQIFDAVFECTLNMINKDFEEYPEHRTNFFLLLQAVNSHCFPAFLAIPPAQFKLVLDSIIWAFKHTMRNVADTGLQILFTLLQNVAQEEAAAQSFYQTYFCDILQHIFSVVTDTSHTAGLTMHASILAYMFNLVEEGKISTPLNPN

#
loop_
_ma_alignment.ordinal_id
_ma_alignment.alignment_id
_ma_alignment.target_template_flag
_ma_alignment.sequence
1 1 2 DMACDTFIKIAQKCRRHFVQVQVGEVMPFIDEILNNINTIICDLQPQQVHTFYEAVGYMIGAQTDQTVQEHLIEKYMLLPNQVWDSIIQQATKNVDILKDPETVKQLGSILKTNVRACKAVGHPFVIQLGRIYLDMLNVYKCLSENISAAIQANGEMVTKQPLIRSMRTVKRETLKLISGWVSRSNDPQMVAENFVPPLLDAVLI---------DYQRNVPAAREPEVLSTMAIIVNKLGGHITAEIPQIFDAVFECTLNMINKDFEE---------YPEHRTNFFLLLQAVNSHCFPAFLAIPPAQ---FKLVLDSIIWAFKHTMRNVADTGLQILFTLLQNVAQEEAAAQSFYQTYFCDILQHIFSVVTDTSHTAGLTMHASILAYMFNLVEEGKISTPLNPN
2 1 1 DSYVETLDSMIELFKDYKPGSITLENITRLCQTL-GLESFTEELSNELSR--LSTASKIIVIDVDYNKKQDRIQDVKLVLASNFDNFDYFNQRDGEHEKSNILLNSLTKYPDLKAFHNNLKFLYLLDAYSHIESDSTSHNNGSSDKSLDSSNASFNNQGKLDLFKYFTELSHYIRQCFQDNCCDFKVRTNLNDKFGIYILTQGINGKEVPLAKIYLEENKSDSQYRFYEYIYSQETKSWINESAENFSNGISLVMEIVANAKESNYTDLIWFPEDFISPELIIDKVTCSSNSSSSPPIIDLFSNNNYNSRIQLMNDFTTKLINIKKFDISNDNLDLISEILKWV------------QWSRIVLQNVFKLVSTPSSNSNSSELEPDYQAPFSTSTKDKNSSTSNTE

#
loop_
_ma_qa_metric.id
_ma_qa_metric.name
_ma_qa_metric.description
_ma_qa_metric.type
_ma_qa_metric.mode
_ma_qa_metric.other_details
_ma_qa_metric.software_group_id
1 MPQS 'ModPipe Quality Score' other global
'composite score, values >1.1 are considered reliable' 1
2 zDOPE 'Normalized DOPE' zscore global . 2
3 'TSVMod RMSD' 'TSVMod predicted RMSD (MSALL)' distance global . .
4 'TSVMod NO35' 'TSVMod predicted native overlap (MSALL)' other global . .

#
loop_
_ma_qa_metric_global.ordinal_id
_ma_qa_metric_global.model_id
_ma_qa_metric_global.metric_id
_ma_qa_metric_global.metric_value
1 1 1 0.665346
2 1 2 -0.11
3 1 3 14.527
4 1 4 0.036

  So it looks like only ModBase would currently benefit from ChimeraX reading template sequences and alignments.  I do not think it would be too hard to implement it.  I've made a ChimeraX feature request for that

https://www.rbvi.ucsf.edu/trac/ChimeraX/ticket/5601

  Tom

On Nov 12, 2021, at 10:20 AM, Ben Webb via ChimeraX-users <chimerax-users@cgl.ucsf.edu> wrote:

Do you have any plans to extend ChimeraX's mmCIF reader to parse and display metadata on theoretical models, such as quality scores or the alignments to template structures?

The folks at PDB have recently done a lot of work to standardize this metadata in the MA mmCIF dictionary:
https://mmcif.wwpdb.org/dictionaries/mmcif_ma.dic/Index/

The dictionary has already been adopted by ModelArchive (e.g. AlphaFold2 models) and by ModBase (Modeller models) and I believe that other repositories such as SwissModel are also moving in that direction. See e.g. mmCIF downloads at
https://www.modelarchive.org/doi/10.5452/ma-bak-cepc-0250
https://modbase.compbio.ucsf.edu/modbase-cgi/model_search.cgi?databaseID=Q12321

(My ulterior motive: we've previously built Chimera web data files to download a ModBase model and the accompanying alignment, and display them in Chimera; now that this data is embedded in the mmCIF file, in principle ChimeraX could do this itself in a less clunky and not ModBase-specific fashion.)

Ben
--
ben@salilab.org                      https://salilab.org/~ben/
"It is a capital mistake to theorize before one has data."
- Sir Arthur Conan Doyle
_______________________________________________
ChimeraX-users mailing list
ChimeraX-users@cgl.ucsf.edu
Manage subscription:
https://www.rbvi.ucsf.edu/mailman/listinfo/chimerax-users