-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathsoftware-results.html
160 lines (144 loc) · 7.76 KB
/
software-results.html
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
{% extends "layout.html" %} {% block title %}Software Results{% endblock %}
{% block lead %}{% endblock %} {% block page_content %}
<style>
body {
background-image: linear-gradient(to bottom,
#f1b9e7 0%,
#e9c4ed 20%,
#e2cef3 40%,
#dad9f9 60%,
#d2e4ff 100%);
}
#float1 {
background-image: url("https://static.igem.wiki/teams/5346/doodle/408-signing-up-removebg-preview.png");
}
#float2 {
background-image: url("https://static.igem.wiki/teams/5346/doodle/algorithms.png");
}
</style>
<div class="col-md-3 sidebar">
</div>
<div class="col-md-9 content">
<h1>IMPROViSeD Pipeline</h1>
<figure>
<img src="https://static.igem.wiki/teams/5346/eng-images/code11.png" id="code11" alt="" />
<figcaption>
Pseudocode for IMPROViSeD. Modelling three dimensional structures for protein complexes using crosslinking is
broken down into two steps: Localization and Registration. The subunits are oriented in multiple
iterations. Each set of results obtained from a subset of crosslinks is followed by localization and registration. The structures with low crosslink violations and those that do not have
clashes make it to as final result.
</figcaption>
</figure>
<h1>Flip Correction Refinement</h1>
<p>
Recall from the "Engineering cycles", that although geometrically mirror images have similar inter-atomic
distances. They are valid solutions to the registration problem. However, one of the instances do not occur in
nature. We therefore to verify the same and discard one of the choice of orientations.
</p>
<p>
This was achieved by the flip correction algorithm. The operation was encoded as a reflection across the origin
in the Ramachandran map, as can be seen in the "Engineering cycles" section. As a result we were successful in
obtaining the correct orientation for the subunits in the protein complex. Speaking in a general sense, the step
can be seen as an integral component for modelling protein structures based on the distance constraints obtained
from experiments and the covalent geometry of the protein.
Below is the visual representation of the flip correction refinement for Ribonuclease Inhibitor Complexed With Ribonuclease A(PDB ID - 1DFJ) . IMPROViSeD was run starting with
just 6 out of 12 crosslinks. Note that after the flip correction, the three dimensional structure is computed
even starting with just half the available crosslinks. This also demonstrates the robustness of the IMPROViSeD
method.
</p>
<figure>
<img src="https://static.igem.wiki/teams/5346/exp-images/1dfj-flip.png" id="noflip" alt="" />
<figcaption>
Registered Structure without flip correction.
</figcaption>
</figure>
<figure>
<img src="https://static.igem.wiki/teams/5346/exp-images/1dfj-proof-2.png" id="flip" alt="" />
<figcaption>
Registered Structure with flip correction.
</figcaption>
</figure>
<h1>
Result for Ribonuclease Inhibitor Complexed With Ribonuclease A (PDB 1DFJ)
</h1>
<p>
We first present the result obtained for the complex with PDB ID 1DFJ. The experimental data was obtained from
<a href='https://github.com/isblab/Integrative_docking_benchmark/tree/main'> link </a>.
Following the localizaton, and having performed the flip correction, as required, the final results for the
complex with PDB ID 1DFJ, is shown below. The results of multiple iterations was displayed.
The backbone RMSD after aligning the obtained result with that availabe in PDB is shown underneath each of the
figure. The number beside represents the number of violated crosslinks out of the total number of crosslinks.
</p>
<p>
Note that we were able to achieve multiple orientations, which are distinct from the original structure in PDB,
but satisfy most of the crosslinks. This not only demonstrates the robustness of the IMPROViSeD method, but also
highlights its efficiently to model new interfaces.
</p>
<figure>
<img src="https://static.igem.wiki/teams/5346/exp-images/whatsapp-image-2024-10-02-at-10-17-38.jpeg" id="1dfj"
alt="" />
<figcaption>
Complex with PDB ID 1DFJ structure modelled using a subset of crosslinks, backbone RMSD with respect to the
structure in PDB,
and crosslink violations (denoted as number of violations/total number of crosslinks) are shown.
</figcaption>
</figure>
<figure>
<img src="https://static.igem.wiki/teams/5346/exp-images/1dfj-movie.gif" id="1dfj" alt="" />
<figcaption>
More new orientations that satisfy most of the crosslinks. The top left image is the original structure in
PDB. The others are due to the IMPROViSeD method. All the structures are shown in similar orientation
</figcaption>
</figure>
<h1>
Result for Human Neutrophil Gelatinase-associated Lipocalin (HNGAL)(PDB ID LCN2)- Hydrolase/Hydrolase Inhibitor
Complex(PDB ID MMP9)
</h1>
<p>
Next we present the result obtained for the complex with LCN2-MMP9. The crosslinking data was obtained
artificially, as descibed in the <a href="https://2024.igem.wiki/iisc-software/software-experiments">Experiments</a> section.
Similar to the previous case, we have performed the flip correction, as required, and the final results for the
complex with LCN2-MMP9, is shown below. The results of multiple iterations was displayed.
</p>
<p>
The backbone RMSD after aligning the obtained result with that availabe in PDB is shown underneath each of the
figure. The number beside represents the number of violated crosslinks out of the total number of crosslinks.
As before, multple results from IMPROViSeD demonstrate its efficiency to model orientations which differ from
those available in the public repository (PDB). This is a significant advantage as it allows us to model new
interfaces.
</p>
<figure>
<img src="https://static.igem.wiki/teams/5346/exp-images/lcn2-drawing.png" id="lcn2-mmp9" alt="" />
<figcaption>
LCN2-MMP9 structure modelled using a subset of crosslinks, backbone RMSD with respect to the structure in
PDB, and crosslink violations (denoted as number of violations/total number of crosslinks) are shown.
</figcaption>
</figure>
<figure>
<img src="https://static.igem.wiki/teams/5346/exp-images/lcn-movie.gif" id="lcn2-mmp9" alt="" />
<figcaption>
Comparison of orientations for LCN2-MMP9 complex with PDB. The one downloaded from PDB is shown in top left.
The ones modelled by IMPROViSeD are shown in similar orientation.
</figcaption>
</figure>
<h1>
Analysis of Results
</h1>
<ul>
<li>
IMPROViSeD was able to model the structure of the protein complex, shown here for two instances, both from the
experimental and artificial data.
</li>
<li>
The time taken for execution is less than 1.5 minutes. This is a significant improvement over the existing
methods, which take hours. The user certainly has the choice of running molecular dynamics simulations to refine
the structure further.
</li>
<li>
The method is robust and can model multiple orientations. It is a unique advantage as it allows an end-user to
model new interfaces. This has the potential to offer insights into the physiological functions of the protein
complexes.
</li>
</ul>
</div>
{% endblock %}