Skip to content

From_sequence_to_model

Jump to bottom
Marta Martinez edited this page May 18, 2023 · 53 revisions

1.- Get the initial model

Basics in homology modeling: We need look for a template to provide the structure to borrow and build a MSA to stablish the best sequence correlation between template and target.

Modeling workflow to get the first model of the alpha subunit of the human haemoglobin using a sequence homologous template (upper part of the alpha subunit modeling workflow):

From sequence to model section

Auxiliary information

Specimen PDB UniProt EMDB
HEMOBLOGIN 1PBX P69905 3488
  1. Step 1b: Import human haemoglobin alpha subunit sequence (details in modeling tutorial and protocol appendix)

Import human haemoglobin alpha subunit sequence (step 1.b)

  1. Step 5a: Import atomic structure to be used as model template (details in modeling tutorial and protocol appendix)

Import atomic structure (step 5.a)

  1. Steps 5b, 5c, 5d and 5e: Import additional sequences to prepare your sequence alignment (details in modeling tutorial and protocol appendix)

Import additional sequences (steps 5.b, 5.c, 5.d and 5.e)

  1. Step 6a: ChimeraX model_from_template protocol to get the first approximation to the atomic structure of the human Hgb alpha subunit (details in modeling tutorial and protocol appendix)

Model_from_template protocol (step 6.a)

Once ChimeraX GUI is opened:

Go to Tools -> Sequence -> Modeller Comparative

  1. Modeller menu to retrieve human Hgb alpha subunit models based on sequence homology (details in modeling tutorial and protocol appendix)

Modeller menu (step 6.a)

At this point wait for a while until the retrieved models appear.

Once you have the retrieved models, select one of them, for instance #3.2. Complete the model with the HEME group from the atomic structure used as template, and save it in Scipion with the next command lines:

  • rename #3.2 id #4
  • save /tmp/chainA.cif format mmcif models #4
  • open /tmp/chainA.cif
  • select #2/A:144
  • save /tmp/HEME.cif format mmcif models #2 selectedOnly true
  • open /tmp/HEME.cif
  • scipioncombine #4,5
  • scipionwrite #6 prefix Hgb_alpha_

ALTERNATIVE

  • rename #3.2 id #4
  • select #4
  • select add #2/A:144
  • delete ~ sel
  • scipioncombine #2,4 modelid 10
  • scipionwrite #10 prefix Hgb_alpha_
  1. Step 6c: Recovering of previous ChimeraX session (details in modeling tutorial and protocol appendix)

Recover session (step 6.c)

Clone this wiki locally