You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: education/HADDOCK24/HADDOCK24-protein-protein-basic/index.md
+14-10Lines changed: 14 additions & 10 deletions
Display the source diff
Display the rich diff
Original file line number
Diff line number
Diff line change
@@ -16,14 +16,15 @@ This tutorial was last updated on 12-03-2026
16
16
<hr>
17
17
## Introduction
18
18
19
-
This tutorial will demonstrate the use of HADDOCK for predicting the structure of a protein-protein complex from NMR chemical shift perturbation (CSP) data. Namely, we will dock two E. coli proteins involved in glucose transport: the glucose-specific enzyme IIA (E2A) and the histidine-containing phosphocarrier protein (HPR).
19
+
This tutorial will demonstrate the use of HADDOCK for predicting the structure of a protein-protein complex from NMR chemical shift perturbation (CSP) data. Namely, we will dock two E. coli proteins involved in glucose transport: the [glucose-specific enzyme IIA](https://www.uniprot.org/uniprotkb/P69783/){:target="_blank"} (E2A) and the [histidine-containing phosphocarrier protein](https://www.uniprot.org/uniprotkb/P0AA04/){:target="_blank"} (HPR).
20
20
21
-
Bacteria use a specific mechanism to import glucose from outside the cell.
21
+
Bacteria use a specific mechanism to import glucose from outside the cell.
22
22
As glucose enters the cell, a phosphate group is attached to it, i.e. glucose becomes phosphorylated.
23
23
This phosphorylation prevents glucose from diffusing back out of the cell and at the same time marks it for further metabolism.
24
24
The phosphate group used for the glucose transport process originates from phosphoenolpyruvate (PEP) and is transferred through a cascade of proteins.
25
25
It first moves from PEP to enzyme I, then to HPR, next to E2A, and finally to enzyme IIB. Enzyme IIB is located on the cytoplasmic side of the membrane, where the phosphate group is ultimately transferred to glucose as it crosses the membrane.
26
-
This animation provides a simple visualisation of the entire process:
26
+
More information can be found in [Jeckelmann *et al*, Eur J Physiol (2020)](https://doi.org/10.1007/s00424-020-02379-0){:target="_blank"}.
27
+
In the mean time, this animation provides a simple visualisation of the entire process:
27
28
28
29
<divstyle="width:60%; margin:auto;">
29
30
<iframe
@@ -99,8 +100,11 @@ The second stage of the docking protocol introduces flexibility to the interacti
99
100
</details>
100
101
<br>
101
102
102
-
**3. Refinement in Cartesian space with explicit solvent (water)**
103
-
**_Note_** that this stage was part of the standard HADDOCK protocol up to (and including) v2.2. As of v2.4 it is no longer performed by default but the user still has the option of enabling it. In its place, a short energy minimisation is performed instead. The final stage of the docking protocol immerses the complex in a solvent shell so as to improve the energetics of the interaction. HADDOCK currently supports water (TIP3P model) and DMSO environments. The latter can be used as a membrane mimic. In this short explicit solvent refinement the models are subjected to a short molecular dynamics simulation at 300K, with position restraints on the non-interface heavy atoms. These restraints are later relaxed to allow all side chains to be optimized.
103
+
**3. Refinement in Cartesian space with explicit solvent (itw)**
104
+
The final stage of the docking protocol immerses the complex in a solvent shell so as to improve the energetics of the interaction. HADDOCK currently supports water (TIP3P model) and DMSO environments. The latter can be used as a membrane mimic. In this short explicit solvent refinement the models are subjected to a short molecular dynamics simulation at 300K, with position restraints on the non-interface heavy atoms. These restraints are later relaxed to allow all side chains to be optimized. <br>
105
+
**_Note_** that as of v2.4, it is no longer performed by default, which used to be the case up to (and including) v2.2.
106
+
Instead, a short energy minimisation in Cartesian space is performed.
107
+
Users can still opt for refinement with explicit solvent as an alternative to this energy minimisation.
104
108
105
109
<details >
106
110
<summarystyle="bold">
@@ -126,7 +130,7 @@ hide lines<br>
126
130
</a>
127
131
128
132
You should see a cartoon representation of the protein.
129
-
It is known from the literature that a phosphate group interacts with the side chain of a histidine residue.
133
+
It is known from the literature that a phosphate group interacts and can form a covalent bond with the side chain of a histidine residue.
130
134
Let us first check whether histidine residues are present in this structure:
131
135
<aclass="prompt prompt-pymol">
132
136
show sticks, resn HIS<br>
@@ -169,7 +173,7 @@ Inspect the surface.
169
173
<aclass="prompt prompt-question">Do the identified residues form a well defined patch on the surface?</a>
170
174
<aclass="prompt prompt-question">Do they form a contiguous surface?</a>
171
175
172
-
The answer to the last question should be no: we can observe residue in the center of the patch that do not seem significantly affected while still being in the middle of the defined interface. This is the reason why in HADDOCK we also define "*passive*" residues that correspond to surface neighbors of active residues. These can be selected manually, or more conveniently you can let the HADDOCK server do it for you (see [Setting up the docking run](#setting-up-the-docking-run) below).
176
+
The answer to the last question should be **no**: we can observe residue in the center of the patch that do not seem significantly affected while still being in the middle of the defined interface. This is the reason why in HADDOCK we also define "*passive*" residues that correspond to surface neighbors of active residues. These can be selected manually, or more conveniently you can let the HADDOCK server do it for you (see [Setting up the docking run](#setting-up-the-docking-run) below).
173
177
174
178
As final step save the molecule as a new PDB file which we will call: *e2a_1F3G.pdb*<br>
175
179
For this in the PyMOL menu on top select:
@@ -200,7 +204,7 @@ hide lines<br>
200
204
201
205
<aclass="prompt prompt-question">Are there any histidines present in this structure?</a>
202
206
<aclass="prompt prompt-question">Is there a phosphate group present in this structure?</a>
203
-
<aclass="prompt prompt-question">Are there any irrelevant molecules present in this structure?</a>
207
+
<aclass="prompt prompt-question">Are there any irrelevant (for the docking) molecules present in this structure?</a>
204
208
205
209
Let's vizualize the residues affected by binding as identified by NMR. From [Wang *et al*, EMBO J (2000)](https://doi.org/10.1093/emboj/19.21.5635){:target="_blank"} the following residues were identified has having significant chemical shift perturbations:
206
210
@@ -229,11 +233,11 @@ To stop the playback:
229
233
230
234
HADDOCK is able to handle such ensembles by using each conformer in turn as a starting point for docking.
231
235
We generally recommend limiting the number of conformers used.
232
-
Otherwise, the number of possible combinations between the input molecules can quickly explode (i.e. become very large).
236
+
Otherwise, the number of possible combinations between the input molecules can quickly escalate (i.e. become very large).
233
237
For example, if both partners contain 10 conformers, this results in 100 possible starting combinations. If 1000 rigid-body models are generated (see [HADDOCK general concepts](#haddock-general-concepts) above), each combination would then be sampled only 10 times!
234
238
235
239
In case if limiting number of input conformers is an unreasonable choice, it is possible to increase the number of models generated in the rigid-body docking stage (it0).
236
-
However, this requires elevated access privileges on the HADDOCK 2.4 server.
240
+
However, this requires elevated permissions level on the HADDOCK 2.4 server, which you can request via "[User Dashboard](https://wenmr.science.uu.nl/dashboard){:target="_blank"}".
237
241
238
242
Now let's display all models of this NMR ensemble simultaneously in ribbon representation.
239
243
This representation is handy for visualizing backbone conformation:
0 commit comments