This is out of date. We have a better tutorial now
Aim: we will fit a domain of the human respiratory complex I matrix arm starting from an empty (part of the) map.
WARNING::This is non-trivial, a familiarity with Coot is assumed. This is not like the basic tutorial. This is approximately how I use coot. It is tricky because fitting cryo-EM maps is tricky.
Put the scripting files in $HOME/.coot-preferences (make the directory if needed)
Using Coot:
read in the map for EMD-6771
read 5xtb-sans-J.pdb
Fetch the model for accession code 3wj7
from PDBe
Use Copy Fragment to extract the protein A
chain from this
Centre on the unmodelled density for the J
chain
Move the A
fragment from 3wj7
to here
Morph → Jiggle fit this chain (this can take ~60 seconds)
Using the ProSMART
module, generate 4.3A self-restraints for this (now fitted) chain
Chain-refine (Shift E)
Calculate → Scripting → Python:
set_alignment_gap_and_space_penalty(-6, -0.4)
Calculate → Align and mutate the A-chain-fragment to the seqence of the J chain
Use the Validation tools to find problematic regions of the model and fix them.
Hint 1: Density Fit analysis is a good tool to start off.
Hint 2: Use the sequence file J-chain.seq
for alignment
End: You may notice that theres a ligand in the middle of the protein - it's NDP. Get an instance of NDP and fit it to the map using Calculate → Other Modelling Tools → Ligand Fitting