|(1999) Chapter 18: Determination of Protein Structures in Introduction to Protein Structure, Second Edition.|
|(1985) Interpretation of Electron Density Maps in Methods in Enzymology, 115: 189-206. PMID: 3841180|
The two aims are:
There are 4 parts, for different fitting tasks at varying resolutions. Download and install KiNG for standalone use. See the instructions on the KiNG software page.
The 1mjf map file is available from the class web site if needed, but you should try downloading it from a very useful web site called the Upssala Electron Density Server (http://eds.bmc.uu.se/eds/). Enter the 1mjf code, go to "Maps" on the left-hand menu, and ask to generate the 2Fo-Fc map as O format; with the right button, ask to save link target, to download the map file.
Launch KiNG with file 1mjfHa.kin.gz (106KB). With "Find point" on the Edit menu, type in "cg 263"; zoom in (drag down with the right mouse button), and clip to remove most of the obscuring things in front of the His ring - or, just choose view 2. Under the Tools menu, choose "Open electron density maps" and browse to find file 1mjf.omap.gz (1.6MB). You should see gray contours appear, and a small control window that should be moved off to the side. From this view above the ring, the fit looks good, meaning that chi1 is good. Choose view2, or rotate to view the ring in-plane from the sidechain end; the density is rounded, but about twice as wide as thick - the model is not correctly aligned to the plane of the density. Under the Tools menu, choose "Sidechain rotator"; select 1mjfH.pdb as the model file; move the Model manager dialog off to one side, and turn on Probe dots. Click an atom in His 263 with the middle mouse button (with a one-button mouse, add the Control key) to activate fitting of His 263 - you should get an orange, idealized sidechain and a rotation-control window, which can be moved to a convenient position that doesn't block the graphics. Try the nearest rotamer in the hypertext list (which would be t60), and/or drag across the bottom of the chi2 circle to refit the angle to match the density. (Note that the left button makes a rapid angle change and the right button a slow one.)
What value of chi2 fits? _________
Are the contacts and H-bonds for the His ring better or worse after your change? _________
Choose "Release" and "Yes", to accept your refitting.
Explore the neighborhood in this map at 1.8Å resolution. (Note: KiNG's "Find point" function is on the Edit menu.) For the nearby Phe 267, look for a dimple in the center of the ring; look end-on, to see that the cross-section of its 6-membered ring is more flattened than for the 5-membered His 263 ring.
For Pro 63, look edge-on to the ring:
does the density clearly support the direction of Cγ pucker? _______
Quit out of KiNG without saving anything.
Launch KiNG with file 1BKRH-contact.kin.gz (139KB) (which you could have made in MolProbity and downloaded). Move it around and locate the biggest clump of red clash spikes. Center on it with a right-button click (or shift-click, for a one-button mouse), zoom in to view the Arg 32 guanidinium, and narrow the zclip. Under the Tools menu, choose "Open electron density maps" and browse to find file 1bkr_2Fo-Fc.omap.gz (1.3MB); move the map dialog to one side. There are bad clashes with waters, and the sidechain occupies only part of the density. It looks as though one of the waters was fit where the sidechain should be. To set up for refitting, choose "Sidechain rotator" under the Tools menu, choose the 1BKRH.pdb.gz (70KB) model file, and enable Probe dots in the Model manager box. Then right-click on an atom in Arg 32 to enable rotations. Look at Arg Cγ to confirm that chi1 should be plus (opposite to the alpha H). Try out all the p rotamers, to find the one positioned closest to the sidechain density.
Shift the 4 dihedral angles to come as close as you can to matching the density; the values are: chi1 ______, chi2 _______, chi3 _______, chi4 _______.
Release and accept.
A refit Arg similar to yours was put back into crystallographic refinement. Append file 1bkr_later.kin.gz (25KB) to see the position it ended up in.
Does it fit the electron density convincingly better than the original model? ________
Quit out of KiNG
Launch KiNG with file cabenolH.kin.gz (70KB). Find trp 179, zoom way in, center on the ring, and clip down around the ring flat in the screen plane.
Measure (under Tools) the chi1 angle: _________.
That is too close to eclipsed; we initially found this problem on MolProbity's list of bad rotamers. Open the electron density map called cabenol_2FoFc.omap.gz (410KB).
Does the Trp ring fit the density fairly well? _______
Now open the difference density map ( cabenol_FoFc.omap.gz (400KB)); set the lower color to sky and turn it on; set the upper color to orange and move its value to -3.0 sigma; move the dialog out of the way. Orange contours mean there is too much model in that region, and blue ones that there is not enough - in this case they strongly suggest that the ring is elongated in the wrong direction. Choose the Sidechain rotator tool with cabenolH.pdb.gz (61KB) as the model, turn on Probe dots, and move the box out of the way.
Center-click an atom in the Trp, try rotamers (you know chi1 should be m) to find the nearest, and adjust the angles to get a good fit in the density and with good Probe contacts (but the other way around - it should have more atoms in the blue difference density than before). Chi1: _______, chi2: ______.
Find the ne1 atom; is there something that might provide an H-bond partner for it? _____________________
Release, yes, to keep your change.
Discard the difference map, and zoom out. Center on a nearby helix and study the appearance of its density (perhaps adjust the map scope a bit larger with the top slider in its dialog box).
Look end-on down the helix axis; is there a hole down the density? ______
Try turning the model on & off to get a feeling for the density shape of helical backbone; which backbone atom sticks out the most from the general backbone worm shape? ___________
Does it point N-terminal or C-terminal along the helix direction? _________
Find Tyr 180; compare the planarity of the ring, OH, and Cbeta atoms with the pronounced "elbow" at the _________ atom.
Find a Met, to see what the sidechain sulfur looks like; turn on the 3sigma purple density, and slide it to higher values until it almost disappears: at what value? _______ for which Met? ____
Do a find on cg for lysines 23, 53, and 81, and evaluate their density; which one has strong density all the way to nz? _____; which has density that gradually peters out? _____;
which disappears abruptly _____, last enclosing atom ______?
For this last case, adjust the density contour lower, to see why they fitted the sidechain where they did; at what value can you see to the end? _______
Sometimes the density for an exposed Lys is invisible beyond Cbeta, because it is free to adopt so many conformations.
Launch KiNG with file 301dH_mcCon.kin.gz (41KB), noting that there are some bad clashes along the backbone. Zoom in on the tetraloop at the top (or choose View2), turn on the bases ("sidechain"), and open electron density file 301d.omap.gz (115KB). Move around in it to see the shapes of bases and backbone; try changing both contour levels, to see if you can find a gray contour that encloses the bases and a purple one that is just the phosphates; then try it in the crowded 151-152 region. Fill in the table below with your estimates of yes, no, or maybe for whether the listed feature can reliably be placed in density at this resolution.
|overall fold or trace|
|position of base|
|orientation of base ring|
|backbone dihedral angles|