Difference between revisions of "Part:BBa K4165030"

 
(22 intermediate revisions by 6 users not shown)
Line 3: Line 3:
 
<partinfo>BBa_K4165030 short</partinfo>
 
<partinfo>BBa_K4165030 short</partinfo>
  
This composite part consists of T7 promoter (BBa_K3633015), lac operator (BBa_K4165062), pGS-21a RBS (BBa_K4165016), 6x His-tag (BBa_K4165020), SPINK8 Inhibitor (BBa_K4165010), TD28rev (BBa_K4165006), WWW (BBa_K4165007), H1A peptide (BBa_K4165000) and T7 terminator (BBa_K731721).
+
This composite part consists of T7 promoter (BBa_K3633015), lac operator (BBa_K4165062), pGS-21a RBS (BBa_K4165016), 6x His-tag (BBa_K4165020), H1A peptide (BBa_K4165000), GGGGSG (BBa_K4165017), TD28rev (BBa_K4165006), GGSGGGG (BBa_K4165018), WWW (BBa_K4165007), SPINK8 inhibitor (BBa_K4165010) and T7 terminator (BBa_K731721).
  
 
===Usage and Biology===
 
===Usage and Biology===
Switch 10 is used to mediate the activity of HTRA1. Activating HTRA1 requires a conformational change in the linker, eliminating the attached inhibitor from the active site. The conformational rearrangement can be mediated through the affinity clamp for tau and beta-amyloid binding. These clamps are used for stabilizing the inhibitor away from the active site. These two domains (inhibitor and affinity clamp connected with linker1). Additionally, (H1A) binding peptide bound to the PDZ domain and connected to the affinity clamp on the other side with linker3.
+
Switch 10 is used to mediate the activity of HTRA1. It is composed of 3 parts connected by different linkers; an HtrA1 PDZ peptide, a clamp of two targeting peptides for tau or amyloid beta, and a catalytic domain inhibitor. Activating HTRA1 requires a conformational change in the linker, eliminating the attached inhibitor from the active site. The conformational rearrangement can be mediated through the binding of affinity clamp to tau or beta-amyloid. This binding will result in a tension that detaches the inhibitor from the active site.
  
===Modeling===
+
The TD28REV and WWW peptides which are considered as tau binding peptides are proved experimentally to bind with tau inhibit the aggregations of tau aggregations respectively. The H1A peptide was also proven to bind with the PDZ of HtrA1 experimentally. The last part which is the inhibitor which is mainly a serine protease inhibitor, and since our protease is a serine protease, so it will act and inhibit the Protein. The whole construction was similarly proved from literature.
TRrosetta models this composite part with a score of 5 out of 6 according to our quality assessment code (you can find the python script file on the programming club page with further explanation of how you can optimize it to your needs).
+
  
 +
<!-- -->
 +
<span class='h3bb'><p style=" font-weight: bold; font-size:17px;"> Sequence and Features</p></span>
 +
<partinfo>BBa_K4165030 SequenceAndFeatures</partinfo>
 +
 +
===Dry Lab===
 +
The top model is TrRoseeta model and the workflow to reach the model will be described below.
 
<html>
 
<html>
<p><img src="https://static.igem.wiki/teams/4165/wiki/parts-registry/11.png" style="margin-left:200px;" alt="" width="500" /></p>
+
<p><img src="https://static.igem.wiki/teams/4165/wiki/parts-registry/pep10.png" style="margin-left:200px;" alt="" width="500" /></p>
 
</html>
 
</html>
  
                            Figure 1. The 3D structure of switch 10 modeled by TRrosetta
+
                      Figure 1. The 3D structure of switch 10 protein Visualized by Pymol. Red: Tau binding peptides,
   
+
                                      blue: H1A peptide, cyan: inhibitor, and green: linkers.
<!-- -->
+
 
<span class='h3bb'>Sequence and Features</span>
+
The pipeline for creating this model is discussed in detail in the section below:
<partinfo>BBa_K4165030 SequenceAndFeatures</partinfo>
+
 
 +
<h1>Switch construction Pipeline</h1>
 +
 
 +
<html>
 +
<img src="https://static.igem.wiki/teams/4165/wiki/registry/dry-lab-modelling-pipeline.png" style="margin-left:200px;" alt="" width="500" /> <br>
 +
 
 +
 
 +
 
 +
<p style="text-align:center;">Figure 2. A figure which describes our Dry-Lab Modelling Pipeline. By team CU_Egypt 2022.</p>
 +
 
 +
<p style=" font-weight: bold; font-size:14px;"> 1) Modelling </p>
 +
<p> Since our Switch parts (HTRA1 binding peptide, TAU, and Beta-amyloid Binding peptide) do not have experimentally acquired structures, we modeled each one of them separately. This approach is done using both denovo modeling (ab initio) and template-based modeling. For modeling small peptides of our system, we used AppTest and Alphafold.</p>
 +
 
 +
<p style=" font-weight: bold; font-size:14px;"> 2) Structure Assessment </p>
 +
<p>In order to assess the quality of generated structures, we used the Swiss-Model tool, which gives an overall quality of any 3D structure (For more information, please check our 
 +
<a href="https://2022.igem.wiki/cu-egypt/ProteinModelling.html">Modeling page</a>.</p>
 +
 
 +
 
 +
<p style=" font-weight: bold; font-size:14px;"> 3) Quality Assessment </p>
 +
<p>Using the code created by us (CU_Egypt 2022), we use the JSON files created from the structure assessment step in Swiss-Model to rank all the models out of score 6. For more information: <a href="https://2022.igem.wiki/cu-egypt/ProgrammingClub.html">Programming club page code under the name of Modric.</a>.</p>
 +
 
 +
<p style=" font-weight: bold; font-size:14px;">4) Filtering</p>
 +
<p>We take the top-ranked models from the previous steps that have either a score of 5 or 6 </p>
 +
 
 +
<p style=" font-weight: bold; font-size:14px;">5) Docking</p>
 +
<p>The top models of inhibitor and HTRA Binding Peptide are docked with HtrA1, and the top models of the clamps are docked with the Target protein, that is, in our case is Beta-amyloid (BBa_K4165004).</p>
 +
 
 +
 
 +
<p style=" font-weight: bold; font-size:14px;"> switch 10 vs HtrA1 trimer: </p>
 +
 
 +
<p style=" font-weight: bold; color: red; font-size:14px;"> ΔG = -20.041 </p>
 +
 
 +
<html>
 +
<p><img src="https://static.igem.wiki/teams/4165/wiki/parts-registry/pep10htra1-dock.png" style="margin-left:200px;" alt="" width="500" /></p>
 +
</html>
 +
 
 +
                          Figure 3. The 3D structure of switch 10 docked to HtrA1 Visualized by Pymol
 +
 
 +
 
 +
 
 +
 
 +
<p style=" font-weight: bold; font-size:14px;">6) Ranking</p>
 +
<p>The docking results are ranked according to their PRODIGY results. For more information: (Link Docking page).</p>
 +
 
 +
<p style=" font-weight: bold; font-size:14px;">7) Top Models</p>
 +
<p>The results that came out from PRODIGY are ranked and the top models are chosen to proceed with to the next step. For more information: (Link Docking page).</p>
 +
 
 +
<p style=" font-weight: bold; font-size:14px;">8) Alignment</p>
 +
<p>Docked structures are aligned. This means that the HtrA1- binding peptide complex is aligned with the second complex which
 +
  is the HtrA1-inhibitor complex to check whether they bound to the same site or not.</p>
 +
<html>
 +
<p><img src="https://static.igem.wiki/teams/4165/wiki/parts-registry/switches/switch31/picture10.png" style="margin-left:200px;" alt="" width="500" /></p>
 +
</html>
 +
 
 +
 
 +
      Figure 5. Aligned structures of HtrA1 binding peptide 1 docked to HtrA1 and inhibitor docked to HtrA1. </p>
 +
 
 +
<p style=" font-weight: bold; font-size:14px;">9) Linker length</p>
 +
<p>The linker lengths are acquired by seeing the distance between the inhibitor and the HtrA1 binding peptide, which is between both C terminals, N terminals, C- and N- terminal, and N- and C-terminals.</p>
 +
 
 +
 
 +
<html>
 +
<style>
 +
table, th, td {
 +
  border:1px solid black; margin-left:auto;margin-right:auto;
 +
}
 +
</style>
 +
<body>
 +
<table style="width:65%">
 +
<table>
 +
  <tr>
 +
    <th>CN</th>
 +
    <th>CC</th>
 +
    <th>NC</th>
 +
    <th>NN</th>
 +
  </tr>
 +
  <tr>
 +
    <td>6.073</td>
 +
    <td>6.476</td>
 +
    <td>6.555</td>
 +
    <td>6.869</td>
 +
  </tr>
 +
</table>
 +
</body>
 +
</html>
 +
 
 +
                            Table 1: The average linker amino acid length of all possibilities.
 +
 
 +
<p style=" font-weight: bold; font-size:14px;">10) Assembly</p>
 +
<p>After settling on the linkers lengths, now we will proceed to the assembly step of the whole system, which is done using TRrosetta, AlphaFold, RosettaFold, and Modeller.</p>
 +
 
 +
<html>
 +
<p><img src="https://static.igem.wiki/teams/4165/wiki/h1a.jpg" style="margin-left:70px;" alt="" width="200" /><img src="https://static.igem.wiki/teams/4165/wiki/inhibitor.png" style="margin-left:70px;" alt="" width="200" /></p>
 +
</html>
 +
<html>
 +
<p style=" font-weight: bold; font-size:14px;">11) Structure Assessment</p>
 +
<p>In order to assess the quality of our structures, we used the Swiss-Model tool, which gives an overall quality of any 3D structure (For more information, please check our <a href="https://2022.igem.wiki/cu-egypt/ProteinModelling.html">Modeling page</a>.</p>
 +
 
 +
 
 +
<html>
 +
<p style=" font-weight: bold; font-size:14px;">12) Quality Assessment </p>
 +
<p>Using the code created by us (CU_Egypt 2022), we use the JSON files created from the structure assessment step in Swiss-Model to rank all the models For more information, please proceed to our <a href="https://2022.igem.wiki/cu-egypt/ProgrammingClub.html">Programming club</a> under the name of Modric.</p>
 +
 
 +
 
 +
<html>
 +
<style>
 +
table, th, td {
 +
  border:1px solid black; margin-left:auto;margin-right:auto;
 +
}
 +
</style>
 +
<body>
 +
<table style="width:65%">
 +
<table>
 +
  <tr>
 +
    <th>cbeta_deviations</th>
 +
    <th>clashscore</th>
 +
    <th>molprobity</th>
 +
    <th>ramachandran_favored</th>
 +
    <th>ramachandran_outliers</th>
 +
    <th>Qmean_4</th>
 +
    <th>Qmean_6</th>
 +
  </tr>
 +
  <tr>
 +
    <td>0</td>
 +
    <td>0</td>
 +
    <td>1.49</td>
 +
    <td>100</td>
 +
    <td>0</td>
 +
    <td>2.023586</td>
 +
    <td>1.78381</td>
 +
  </tr>
 +
</table>
 +
</body>
 +
</html>
 +
               
 +
                              Table 2: Quality assessment parameters of Switch 10 model.
 +
<p style=" font-weight: bold; font-size:14px;">13) Alignment</p>
 +
<p>The docked structures are then aligned and compared to the basic parts, which are docked with the protein of interest (HtrA1). The structures with the least RMSD are chosen following the recommended range provided by CAPRI protocol.</p>
 +
 
 +
 
 +
<html>
 +
<style>
 +
table, th, td {
 +
  border:1px solid black; margin-left:auto;margin-right:auto;
 +
}
 +
</style>
 +
<body>
 +
<table style="width:65%">
 +
<table>
 +
  <tr>
 +
    <th>RMSD Before Docking</th>
 +
    <th>RMSD After Docking</th>
 +
 
 +
  </tr>
 +
  <tr>
 +
    <td>1.739</td>
 +
    <td>1.894</td>
 +
 
 +
  </tr>
 +
</table>
 +
</body>
 +
</html>
 +
 
 +
===Conclusion===
 +
The top model was HtrA1 switch 10 (BBa_K4165030) since it was the best switch fulfilling the criteria of structure assessment, docking, and RMSD.
 +
 
 +
 
 +
 
 +
<p style=" font-weight: bold; font-size:14px;"> Mathematical modeling </p>
 +
<p style=" font-weight: bold; font-size:14px;">Transcription rate and translation rate </p>
 +
the mathematical modeling was based on our code for the calculation of transcription and translation (you can find it in the code section) besides the estimated results from the wet lab.
 +
 
 +
<html>
 +
<p><img src="https://static.igem.wiki/teams/4165/wiki/parts-registry/peptide-1o2.png" style="margin-left:200px;" alt="" width="500" /></p>
 +
</html>
 +
 
 +
              Figure 4 shows the results from the transcription and translation code showing the
 +
                                      variation of mRNA and protein concentrations.
 +
===References===
 +
1. Goedert, M., & Spillantini, M. G. (2017). Propagation of Tau aggregates. Molecular Brain, 10. https://doi.org/10.1186/s13041-017-0298-7
 +
 
 +
2. Etienne, M. A., Edwin, N. J., Aucoin, J. P., Russo, P. S., McCarley, R. L., & Hammer, R. P. (2007). Beta-amyloid protein aggregation. Methods in molecular biology (Clifton, N.J.), 386, 203–225. https://doi.org/10.1007/1-59745-430-3_7
 +
 
 +
4. Seidler, P., Boyer, D., Rodriguez, J., Sawaya, M., Cascio, D., Murray, K., Gonen, T., & Eisenberg, D. (2018). Structure-based inhibitors of tau aggregation. Nature chemistry, 10(2), 170. https://doi.org/10.1038/nchem.2889
 +
 
 +
5. Romero-Molina, S., Ruiz-Blanco, Y. B., Mieres-Perez, J., Harms, M., Münch, J., Ehrmann, M., & Sanchez-Garcia, E. (2022). PPI-Affinity: A Web Tool for the Prediction and Optimization of Protein–Peptide and Protein–Protein Binding Affinity. Journal of Proteome Research.
 +
 
 +
6. Stein, V., & Alexandrov, K. (2014). Protease-based synthetic sensing and signal amplification. Proceedings of the National Academy of Sciences, 111(45), 15934-15939.
 +
 
  
  

Latest revision as of 03:47, 14 October 2022


HtrA1 switch 10

This composite part consists of T7 promoter (BBa_K3633015), lac operator (BBa_K4165062), pGS-21a RBS (BBa_K4165016), 6x His-tag (BBa_K4165020), H1A peptide (BBa_K4165000), GGGGSG (BBa_K4165017), TD28rev (BBa_K4165006), GGSGGGG (BBa_K4165018), WWW (BBa_K4165007), SPINK8 inhibitor (BBa_K4165010) and T7 terminator (BBa_K731721).

Usage and Biology

Switch 10 is used to mediate the activity of HTRA1. It is composed of 3 parts connected by different linkers; an HtrA1 PDZ peptide, a clamp of two targeting peptides for tau or amyloid beta, and a catalytic domain inhibitor. Activating HTRA1 requires a conformational change in the linker, eliminating the attached inhibitor from the active site. The conformational rearrangement can be mediated through the binding of affinity clamp to tau or beta-amyloid. This binding will result in a tension that detaches the inhibitor from the active site.

The TD28REV and WWW peptides which are considered as tau binding peptides are proved experimentally to bind with tau inhibit the aggregations of tau aggregations respectively. The H1A peptide was also proven to bind with the PDZ of HtrA1 experimentally. The last part which is the inhibitor which is mainly a serine protease inhibitor, and since our protease is a serine protease, so it will act and inhibit the Protein. The whole construction was similarly proved from literature.

Sequence and Features


Assembly Compatibility:
  • 10
    COMPATIBLE WITH RFC[10]
  • 12
    COMPATIBLE WITH RFC[12]
  • 21
    COMPATIBLE WITH RFC[21]
  • 23
    COMPATIBLE WITH RFC[23]
  • 25
    COMPATIBLE WITH RFC[25]
  • 1000
    COMPATIBLE WITH RFC[1000]

Dry Lab

The top model is TrRoseeta model and the workflow to reach the model will be described below.

                     Figure 1. The 3D structure of switch 10 protein Visualized by Pymol. Red: Tau binding peptides, 
                                     blue: H1A peptide, cyan: inhibitor, and green: linkers.

The pipeline for creating this model is discussed in detail in the section below:

Switch construction Pipeline


Figure 2. A figure which describes our Dry-Lab Modelling Pipeline. By team CU_Egypt 2022.

1) Modelling

Since our Switch parts (HTRA1 binding peptide, TAU, and Beta-amyloid Binding peptide) do not have experimentally acquired structures, we modeled each one of them separately. This approach is done using both denovo modeling (ab initio) and template-based modeling. For modeling small peptides of our system, we used AppTest and Alphafold.

2) Structure Assessment

In order to assess the quality of generated structures, we used the Swiss-Model tool, which gives an overall quality of any 3D structure (For more information, please check our Modeling page.

3) Quality Assessment

Using the code created by us (CU_Egypt 2022), we use the JSON files created from the structure assessment step in Swiss-Model to rank all the models out of score 6. For more information: Programming club page code under the name of Modric..

4) Filtering

We take the top-ranked models from the previous steps that have either a score of 5 or 6

5) Docking

The top models of inhibitor and HTRA Binding Peptide are docked with HtrA1, and the top models of the clamps are docked with the Target protein, that is, in our case is Beta-amyloid (BBa_K4165004).

switch 10 vs HtrA1 trimer:

ΔG = -20.041

                          Figure 3. The 3D structure of switch 10 docked to HtrA1 Visualized by Pymol



6) Ranking

The docking results are ranked according to their PRODIGY results. For more information: (Link Docking page).

7) Top Models

The results that came out from PRODIGY are ranked and the top models are chosen to proceed with to the next step. For more information: (Link Docking page).

8) Alignment

Docked structures are aligned. This means that the HtrA1- binding peptide complex is aligned with the second complex which is the HtrA1-inhibitor complex to check whether they bound to the same site or not.


     Figure 5. Aligned structures of HtrA1 binding peptide 1 docked to HtrA1 and inhibitor docked to HtrA1. </p>

9) Linker length

The linker lengths are acquired by seeing the distance between the inhibitor and the HtrA1 binding peptide, which is between both C terminals, N terminals, C- and N- terminal, and N- and C-terminals.


CN CC NC NN
6.073 6.476 6.555 6.869

                            Table 1: The average linker amino acid length of all possibilities. 

10) Assembly

After settling on the linkers lengths, now we will proceed to the assembly step of the whole system, which is done using TRrosetta, AlphaFold, RosettaFold, and Modeller.

11) Structure Assessment

In order to assess the quality of our structures, we used the Swiss-Model tool, which gives an overall quality of any 3D structure (For more information, please check our Modeling page.

12) Quality Assessment

Using the code created by us (CU_Egypt 2022), we use the JSON files created from the structure assessment step in Swiss-Model to rank all the models For more information, please proceed to our Programming club under the name of Modric.

cbeta_deviations clashscore molprobity ramachandran_favored ramachandran_outliers Qmean_4 Qmean_6
0 0 1.49 100 0 2.023586 1.78381

                              Table 2: Quality assessment parameters of Switch 10 model. 

13) Alignment

The docked structures are then aligned and compared to the basic parts, which are docked with the protein of interest (HtrA1). The structures with the least RMSD are chosen following the recommended range provided by CAPRI protocol.


RMSD Before Docking RMSD After Docking
1.739 1.894

Conclusion

The top model was HtrA1 switch 10 (BBa_K4165030) since it was the best switch fulfilling the criteria of structure assessment, docking, and RMSD.


Mathematical modeling

Transcription rate and translation rate

the mathematical modeling was based on our code for the calculation of transcription and translation (you can find it in the code section) besides the estimated results from the wet lab.

              Figure 4 shows the results from the transcription and translation code showing the 
                                      variation of mRNA and protein concentrations.

References

1. Goedert, M., & Spillantini, M. G. (2017). Propagation of Tau aggregates. Molecular Brain, 10. https://doi.org/10.1186/s13041-017-0298-7

2. Etienne, M. A., Edwin, N. J., Aucoin, J. P., Russo, P. S., McCarley, R. L., & Hammer, R. P. (2007). Beta-amyloid protein aggregation. Methods in molecular biology (Clifton, N.J.), 386, 203–225. https://doi.org/10.1007/1-59745-430-3_7

4. Seidler, P., Boyer, D., Rodriguez, J., Sawaya, M., Cascio, D., Murray, K., Gonen, T., & Eisenberg, D. (2018). Structure-based inhibitors of tau aggregation. Nature chemistry, 10(2), 170. https://doi.org/10.1038/nchem.2889

5. Romero-Molina, S., Ruiz-Blanco, Y. B., Mieres-Perez, J., Harms, M., Münch, J., Ehrmann, M., & Sanchez-Garcia, E. (2022). PPI-Affinity: A Web Tool for the Prediction and Optimization of Protein–Peptide and Protein–Protein Binding Affinity. Journal of Proteome Research.

6. Stein, V., & Alexandrov, K. (2014). Protease-based synthetic sensing and signal amplification. Proceedings of the National Academy of Sciences, 111(45), 15934-15939.