22layout : page
33title : " PowerFit Tutorial"
44excerpt : " A small tutorial on PowerFit for automatic rigid body fitting"
5- tags : [PowerFit, Cryo-EM, HADDOCK, Ribosome, Chimera , rigid body fitting]
5+ tags : [PowerFit, Cryo-EM, HADDOCK, Ribosome, ChimeraX , rigid body fitting]
66image :
77 feature : pages/banner_education-thin.jpg
88---
@@ -18,11 +18,12 @@ This tutorial consists of the following sections:
1818
1919PowerFit is a software application developed to fit atomic resolution
2020structures of biomolecules to cryo-electron microscopy (cryo-EM) density maps.
21- It is open-source and available for download on [ Github] [ link-powerfit ] .
21+ It is open-source and available for download on [ Github] [ link-powerfit ] . In addition
22+ you can find a manual on usage of PowerFit [ here] [ link-manual ] .
2223
2324This tutorial will show you how to utilize PowerFit by applying it to an E.coli
2425ribosome case. To follow this tutorial, you need, in addition to PowerFit, the
25- [ UCSF Chimera ] [ link-chimera ] visualization software, a popular tool in the
26+ [ UCSF ChimeraX ] [ link-chimerax ] visualization software, a popular tool in the
2627cryo-electron microscopy community for its volume visualization capabilities.
2728We will further discuss the limits of rigid body fitting, and how HADDOCK can
2829alleviate some of the shortcomings. We provide the data necessary to run this
@@ -41,13 +42,13 @@ _AIMS Biophysics_. *2*, 73-87 (2015).
4142_ Structure._ * 23* , 949-960 (2015).
4243
4344Throughout the tutorial, colored text will be used to refer to questions or
44- instructions, Linux and/or Chimera commands.
45+ instructions, Linux and/or ChimeraX commands.
4546
4647<a class =" prompt prompt-question " >This is a question prompt: try answering
4748it!</a >
4849<a class =" prompt prompt-info " >This an instruction prompt: follow it!</a >
49- <a class =" prompt prompt-pymol " >This is a Chimera prompt: write this in the
50- Chimera command line prompt !</a >
50+ <a class =" prompt prompt-pymol " >This is a ChimeraX prompt: write this in the
51+ ChimeraX command line!</a >
5152<a class =" prompt prompt-cmd " >This is a Linux prompt: insert the commands in the
5253terminal!</a >
5354
@@ -74,7 +75,7 @@ ribosome has already been properly fitted in the density.
7475</a >
7576
7677* In case you might run this tutorial on your own* , make sure to have the required
77- software installed ([ UCSF Chimera ] [ link-chimera ] and [ PowerFit] [ link-powerfit ] ),
78+ software installed ([ UCSF ChimeraX ] [ link-chimerax ] and [ PowerFit] [ link-powerfit ] ),
7879and download the data to run this tutorial from our GitHub
7980data repository [ here] [ link-data ] or clone it from the command line
8081
@@ -89,26 +90,24 @@ data repository [here][link-data] or clone it from the command line
8990Let us first inspect the data we have available, namely the cryo-EM density map
9091and the structures we will attempt to fit.
9192
92- Using Chimera , we can easily visualize and inspect the density and models,
93+ Using ChimeraX , we can easily visualize and inspect the density and models,
9394mostly through a few mouse clicks.
9495
9596<a class =" prompt prompt-info " >
9697 Open the density map together with the ribosome and KsgA.
9798</a >
9899<a class =" prompt prompt-cmd " >
99- chimera ribosome-KsgA.map ribosome.pdb KsgA.pdb
100+ chimerax ribosome-KsgA.map ribosome.pdb KsgA.pdb
100101</a >
101102
102103In the ` Volume Viewer ` window, the middle slide bar provides control on the
103104value at which the isosurface of the density is shown. At high values, the
104105envelope will sink while lower values might even display the noise in the map.
105106We will first make the density transparent, to see the fitted structure inside:
106107
107- * Within the ` Volume Viewer ` click on the gray box next to ` Color ` , which opens
108- the ` Color Editor ` window.
109- * In there, check the ` Opacity ` box. An extra slider bar appears in the box
110- called ` A ` , for the alpha channel.
111- * Set the alpha channel value to around 0.6.
108+ <a class =" prompt prompt-pymol " >
109+ transparency #1 60
110+ </a >
112111
113112Notice that the density becomes transparent providing a better view of the fit
114113of the ribosome model. On closer inspection, you can also discern a region of
@@ -133,7 +132,7 @@ on the map and calculate a cross-correlation score for each of them.
133132 map.
134133</a >
135134<a class =" prompt prompt-cmd " >
136- powerfit ribosome-KsgA.map 13 KsgA.pdb -d run-KsgA -a 20 -p 2 -l
135+ powerfit ribosome-KsgA.map 13 KsgA.pdb -d run-KsgA -a 20 -p 2
137136</a >
138137
139138While performing the search, PowerFit will update you on the progress of the
@@ -161,10 +160,8 @@ option specifies where the results will be stored while the `-p` option
161160specifies the number of processors that PowerFit can use during the search, to
162161leverage available CPU resources.
163162
164- Finally, the ` -l ` flag applies a Laplace pre-filter on the density data, which
165- increases the cross-correlation sensitivity by enhancing edges in the density.
166163In this example scenario, all other options are left at their default values
167- but feel free to explore them.
164+ but feel free to explore [ them] [ link-manual ] .
168165
169166
170167## Analyzing the results
@@ -186,14 +183,14 @@ with date and timing information.
186183
187184<a class =" prompt prompt-info " >
188185 Open the density map, the * lcc.mrc* cross-correlation map, and the 10
189- best-ranked solutions in Chimera .
186+ best-ranked solutions in ChimeraX .
190187</a >
191188<a class =" prompt prompt-cmd " >
192- chimera ribosome-KsgA.map run-KsgA/lcc.mrc ribosome.pdb run-KsgA/fit_ * .pdb
189+ chimerax ribosome-KsgA.map run-KsgA/lcc.mrc ribosome.pdb run-KsgA/fit_ * .pdb
193190</a >
194191
195- Make the density map transparent again, by adjusting the alpha channel value to
196- 0.6. The values of the ` lcc.mrc ` slider bar correspond to the cross-correlation
192+ Make the density map transparent again to 60.
193+ The values of the ` lcc.mrc ` slider bar correspond to the cross-correlation
197194score found. In this way, you can selectively visualize regions of high or low
198195cross-correlation values: i.e., pushing the slider to the right (higher cutoff)
199196shows only regions of the grid with high cross-correlation scores.
@@ -202,10 +199,10 @@ As you can see, PowerFit found quite some local optima, one of which stands out
202199(if the rotational search was tight enough). Further, the 10 best-ranked
203200solutions are centered on regions corresponding to local cross-correlation maxima.
204201
205- To view each fitted solution individually, in the main panel, go to ` Favorites `
206- → ` Model Panel ` to open the ` Model Panel ` window. The window shows each model
207- and its associated color that Chimera has processed. To show or hide a specific
208- model you can click the box in the ` S ` column.
202+ To view each fitted solution individually, in the main panel, go to ` Tools `
203+ → ` Models ` to open the ` Model Panel ` window. The window shows each model
204+ and its associated color that ChimeraX has processed. To show or hide a specific
205+ model you can click the box in the column with an eye .
209206
210207<a class =" prompt prompt-info " >
211208 Go through the 10 solutions one by one to appreciate their goodness-of-fit
@@ -214,12 +211,13 @@ model you can click the box in the `S` column.
214211<a class =" prompt prompt-question " >
215212 Do you agree with what PowerFit proposes as the best solution?
216213</a >
214+
217215<a class =" prompt prompt-info " >
218- In a new Chimera session, reopen the density map and the fit that you find
216+ In a new ChimeraX session, reopen the density map and the fit that you find
219217best. Replace * ?* by the appropriate solution number.
220218</a >
221219<a class =" prompt prompt-cmd " >
222- chimera ribosome-KsgA.map ribosome.pdb run-KsgA/fit_ ?.pdb
220+ chimerax ribosome-KsgA.map ribosome.pdb run-KsgA/fit_ ?.pdb
223221</a >
224222
225223You now have combined the ribosome structure with the rigid-body fit of KsgA
@@ -228,29 +226,30 @@ experiments performed on this complex indicate three charged residues of KsgA -
228226` R221 ` , ` R222 ` , and ` K223 ` - that are of special importance for the
229227interaction.
230228
231- In the same session of Chimera where you have your chosen fitted KsgA
232- structure, go to ` Favorites ` → ` Command Line ` . A command line is now present
233- below the main viewing window. In the command line of Chimera, type the
229+ In the same session of ChimeraX where you have your chosen fitted KsgA
230+ structure, typ in the command line of ChimeraX the
234231following instructions to center your view on these residues and highlight
235232their interactions:
236233
237234<a class =" prompt prompt-pymol " >
238- show #2 :221-223 zr<5 & #1 || #2 :221-223
239- center #2 :221-223 zr<5 & #1 || #2 :221-223
235+ sel #3 :221-223 <br >
236+ view sel <br >
237+ contacts sel distanceOnly 5.0 makePseudobonds true reveal true <br >
238+ </a >
240239</a >
241240<a class =" prompt prompt-info " >
242241 Take some time to inspect the model, paying particular attention to these three
243242 residues and their spatial neighbors.
244243</a >
245244<a class =" prompt prompt-question " >
246- Are there any clashes between the ribosome and KsgA chains?
245+ Are there any clashes between the ribosome and KsgA chains? Show the selection as spheres to visualize this better.
247246</a >
248247<a class =" prompt prompt-question " >
249248 Is the mutagenesis data explained by the model, i.e. are the three charged
250249amino acids involved in strong interactions?
251250</a >
252251
253- Chimera also includes a tool to locally optimize the fit of a rigid structure
252+ ChimeraX also includes a tool to locally optimize the fit of a rigid structure
254253against a given density map, which can be an additional help on top of the
255254PowerFit calculations. Make the main display window active by clicking on it,
256255then go to ` Tools ` → ` Volume data ` → ` Fit in Map ` . In the newly opened `Fit in
@@ -259,14 +258,14 @@ Map` window, select the best-fitted structure of PowerFit (`fit_?.pdb`) as
259258Press ` Fit ` to start the optimization.
260259
261260<a class =" prompt prompt-question " >
262- Does the Chimera local fit optimization tool improve the results of PowerFit?
261+ Does the ChimeraX local fit optimization tool improve the results of PowerFit?
263262</a >
264263
265- The scoring function used by Chimera to estimate the quality of the fit makes
264+ The scoring function used by ChimeraX to estimate the quality of the fit makes
266265our model worse, increasing the number of clashes between the ribosomal RNA and
267266KsgA. Click ` Undo ` in the ` Fit in Map ` window to undo the optimization.
268267
269- Next, we will try to optimize the fit using the cross-correlation that Chimera
268+ Next, we will try to optimize the fit using the cross-correlation that ChimeraX
270269provides. Click ` Options ` and check the `Use map simulated from atoms,
271270resolution` box and fill in ` 13` for resolution. Check the ` correlation` radio
272271button and uncheck the ` Use only data above contour level from first map ` .
@@ -304,10 +303,10 @@ HADDOCK model, generated by combining the cryo-EM map, the PowerFit centroid
304303positions, and the mutagenesis data, in the tutorial data folder.
305304
306305<a class =" prompt prompt-info " >
307- Open the density map in Chimera and load the best-ranked HADDOCK model.
306+ Open the density map in ChimeraX and load the best-ranked HADDOCK model.
308307</a >
309308<a class =" prompt prompt-cmd " >
310- chimera ribosome-KsgA.map HADDOCK-ribosome.pdb HADDOCK-KsgA.pdb
309+ chimerax ribosome-KsgA.map HADDOCK-ribosome.pdb HADDOCK-KsgA.pdb
311310</a >
312311<a class =" prompt prompt-question " >
313312Does HADDOCK improve the quality of the model, i.e. are the number of clashes
@@ -367,6 +366,7 @@ suggestions, feel free to contact us via email or by submitting an issue in the
367366appropriate Github repository.
368367
369368[ link-powerfit ] : https://github.com/haddocking/powerfit " PowerFit "
370- [ link-chimera ] : https://www.cgl.ucsf.edu/chimera / " UCSF Chimera "
369+ [ link-chimerax ] : https://www.cgl.ucsf.edu/chimerax / " UCSF ChimeraX "
371370[ link-data ] : https://github.com/haddocking/powerfit-tutorial " PowerFit tutorial data "
372371[ link-density ] : https://www.ebi.ac.uk/pdbe/entry/emdb/EMD-2017 " Ribosome KsgA density "
372+ [ link-manual ] : https://bonvinlab.org/powerfit/manual.html " PowerFit manual "
0 commit comments