Part 2: Rebuilding in KiNG
For this lesson, you will need downloaded files from Part 1: MolProbity Analysis :
1bkrFH-multi.kin.gz and 1bkrFH.pdb
and 1bkr.omap.gz from the EDS.
A. Thr 77
Fire-up KiNG and read in the model
For this tutorial you will need KiNG installed on your system. KiNG is available directly from the Kinemage website or as part of the PHENIX system.
KiNG is a Java-based kinemage graphics viewer, portable to all platforms capable of running a Java virtual machine and amenable to the addition of new plug-in modules. You've already seen the stripped-down web version, KiNGlet, in Part 1.
KiNG has numerous tools available; we'll use the Structural Biology subset to rebuild a poorly-fit threonine in 1BKR. In Part 1, you assessed the quality of this model using the MolProbity web service and downloaded a Reduce-protonated PDB file and a multi-criterion kinemage file. You also downloaded the 2Fo-Fc map for 1BKR from the EDS.
In this rebuilding demo, we are using a very high-resolution structure for didactic purposes -- problems are rare in such cases, but when they do occur their nature and their correction can be seen and understood extremely clearly.
In your working directory, type "king 1bkrFH-multi.kin.gz" or
type "phenix.king 1bkrFH-multi.kin.gz" (if using the PHENIX-embedded version), or open the kinemage from the File menu inside KiNG. Conveniently, the file need not be unzipped. [On a GUI system, you can use drag-n-drop.]
Adjust the kinemage view and read in the map
The multi-crit kinemage contains much information, most of it abstracted in the default view. Use the kinemage display buttons on the right panel of the KiNG window to prepare for the rebuild: turn off the Cαs and turn on mainchain, sidechain, H's, and water.
Go to Thr 77 using Edit > Find point... Use either the Zoom slider or right-mouse vertical-drag to zoom in. (Adjust the slab with the Clipping slider or right-mouse horizontal drag.) Note the big ball that marks the 0.349Å Cβ deviation, and then turn off the "Cbeta dev" button.
Read in the 2Fo-Fc ED map (1bkr.omap.gz) using Tools > Structural Biology > Electron density maps... The map will be displayed in the KiNG graphics window, along with the map's settings widget shown at left (which you'll want to move to a better location). The 15Å default size is OK and the default gray for the first map contour is a decent color. You may want to increase the contour level to 1.4.
Choose a new rotamer
MolProbity identified Thr 77 as an outier both for rotamer and for Cβ deviation. You may be able to tell for yourself that some of the bond angles are distorted and chi1 is eclipsed (look down on the sidechain from its end). Turn the contact dots off (the "dots" button on the panel) and look at the beautiful density versus the model placement, to see that there's definitely a problem with Thr 77. Which atom seems to be fit worst? Since we know that the current model has a non-rotameric conformer, let's start the refitting process by seeing if any of the 3 real rotamers will match the density and give decent contacts.
KiNG uses our Penultimate rotamer library (Lovell, 2000). One way of gaining access to rotamer selection is via Tools > Structural biology > Sidechain Rotator. In the resulting file-selection dialog, navigate to and open the 1bkrFH.pdb file. The kinemage will update with some new buttons: 2 near the top of the list represent the "frozen" and "molten" models used in the rotation. The Model Manager Widget is floated on your screen; in it, turn on "Probe dots", to enable interactive contact display. Select the Thr77 residue for rotation by middle-clicking on one of its atoms.
This will pop up another widget, for control of the sidechain dihedral rotations. On this widget, you'll find a dial for each torsion angle, a listing of rotamers for this residue (with % occurrence), a rotamericity score for the current conformer, and a checkbox for sidechain idealization (toggle to see how this setting affects the model, but leave it turned on).
The dials work by dragging around the edge of the circle, left-mouse for fast and right-mouse for slow change. Move one slightly, to trigger update of the contact dots for this starting position (or click the ">" in the Model manager next to the Probe command line).
Click on each listed rotamer (p, t, or m, in the sidechain-control window) to evaluate its match to the electron density and its interactive probe dots (red clashes and pale green H-bods). [Make sure the "dots" button for the original multi-crit model is turned off, to allow you to evaluate the new contacts.] Right-click on the Cβ and narrow the clipping to isolate the sidechain model and density clearly. Look down the Cβ-Cα bond, but rotate your view as needed to get a feel for the fit of each rotamer. [At this stage, matching the shape and orientation of the density is much more important than being centered on it.]
The m rotamer gives the best agreement to the electron density shape. But it needs some improvement: the whole sidechain is shifted to one side, and the CG methyl clashes with the Lys 78 N amide. Adjustment of chi angles will not correct this problem (try it). We need some way to recenter the Cbeta back into its density, which in turn will alleviate the clashes. For this, we will use the "Backrub".
Backrub the Thr 77 mainchain
We've now got what appears to be the right rotamer for Thr 77, but the entire residue needs to be tilted back into the density. For this, KiNG has a unique and effective backbone adjustment tool, the "Backrub". [For more information about the backrub motion, see Davis et al. (2006).]
Click the "BACKRUB mainchain" button on the A 77 THR rotator widget to invoke the backrub tool for this residue. The BACKRUB widget will pop up; take a look at the information and controls that it displays. Beneath the three rotation dials are reports for the Tau deviation, Ramachandran status, and phi,psi angles for the subject residues. Fundamental to the backrub tool are the three labeled dials that control movement of the central residue (Thr77) and its two neighboring peptides. The primary backrub motion rotates the entire dipeptide as a rigid unit around the axis between the two Cαs at i-1 and i+1; a small shift of the backbone is leveraged into a much larger change of the central sidechain position. The two secondary motions (less often needed, and at smaller amplitudes) each rotate one of the flanking peptides around the axis between its end Cαs, which helps to maintain H-bonding and carbonyl O position.
You may want to get a feel for the effect of the three motions of backrub, before trying to get a good fit. To do so at this time, turn off the map and the frozen model in order to concentrate on the movement of the orange sidechains and brown mainchain of the molten model. Play with the three dials of the Backrub widget (especially the primary, center one) to get a feel for their connection with the model. Return the dials to their original values when done.
Now, let's use the Backrub and Sidechain rotator dials to make the model better fit the density. First, make sure that the dials are in their original positions by matching the black and gray numbers for each dial. Be sure the frozen model and *1bkrFH model are on. (Optional: Center on Cβ of Thr77 (right-mouse click on the Cβ) and then decrease the map coverage to 4Å at top of the map Presets widget.)
Start out by pulling the Cβ of Thr77 back into its density by using the A 77 THR dial (middle dial) on the Backrub widget. Something around 7 to 9° is a reasonable movement.
Turn on the purple 3σ map contours, to see that the Og1 is now in the higher-density arm as it should be. Now use the chi1 dial on the sidechain rotator to move both Og1 and Cg2 more toward the centers of their lobes of density. A chi1 angle near -50 to -55° seems close.
In a real refinement situation, we would quit here - this is close enough to know we've found the right local minimum, and the refinement program will twiddle the rest. If you want, you can continue to optimize the rebuild. Extend the map coverage to about 8Å, and try the two individual-peptide adjustments to center the backbone COs in their density and maintain H-bonding. Rotate the OH to point at the nearby water peak (no atom shown, hence no dots, since it's a symmetry related one).
Apply and save the changes
Backrub tool: Click the "Finished" button. Answer "Yes" to keeping the changes.
Sidechain rotator: Click the "Finished" button. Answer "Yes" to keeping the changes.
In the Model manager window: Choose File > Save PDB file... and in the file-selection dialog give your new PDB file a name and click save.
B. His 42
To see density for the His 42 flip from Part 1:
turn off both frozen model and molten model in button panel on right side of screen,
then enter "42 cg" in Edit/Find point, get a good view, and turn on the "H-bonds" button to remember its contacts in this conformation that was preferred by Reduce. The ring electron density is quite clean, with no evidence of disorder. In the map dialog, move the purple contour level up to 4.7σ, where you should see a difference between peak levels for N vs C atoms that corroborates Reduce's flip choice (in spite of model bias from having refined it the other way around).