Difference between revisions of "Part:BBa K4905006"

 
(33 intermediate revisions by 4 users not shown)
Line 11: Line 11:
 
<html>
 
<html>
 
<body>
 
<body>
 +
<style>
 +
.image-container {
 +
    width: 100%;
 +
    margin: auto;
 +
    text-align: center;
 +
}
 +
.img123 {
 +
    width: 40%;
 +
    height: auto;
 +
}
 +
.row{
 +
    text-align: center;
 +
}
 +
table td, table th {
 +
    border: 1px solid black;
 +
}
 +
table {
 +
    margin: auto;
 +
}
 +
figcaption{
 +
  text-align:center;
 +
  font-size: 12px;
 +
}
  
<html>
+
</style>
<body>
+
 
<h1>Information</h1>
 
<h1>Information</h1>
 
<p>
 
<p>
 
This part is made up of the basic parts: Leucine zipper Z1 (<a href="https://parts.igem.org/Part:BBa_K4905004">BBa_K4905004</a>), Leucine zipper Z2 (<a href="https://parts.igem.org/Part:BBa_K4905005">BBa_K4905005</a>), and two times Elastin-Like Polypeptide (ELP) sequence A<sub>[60]</sub>I<sub>[60]</sub> (<a href="https://parts.igem.org/Part:BBa_K4905001">BBa_K4905001</a>]). This results in the sequence Z1-I<sub>[60]</sub>-A<sub>[120]</sub>-I<sub>[60]</sub>-Z2. With A<sub>[5]</sub> the sequence (VPGAG<sub>[3]</sub>VPGGG<sub>[2]</sub>). </p>
 
This part is made up of the basic parts: Leucine zipper Z1 (<a href="https://parts.igem.org/Part:BBa_K4905004">BBa_K4905004</a>), Leucine zipper Z2 (<a href="https://parts.igem.org/Part:BBa_K4905005">BBa_K4905005</a>), and two times Elastin-Like Polypeptide (ELP) sequence A<sub>[60]</sub>I<sub>[60]</sub> (<a href="https://parts.igem.org/Part:BBa_K4905001">BBa_K4905001</a>]). This results in the sequence Z1-I<sub>[60]</sub>-A<sub>[120]</sub>-I<sub>[60]</sub>-Z2. With A<sub>[5]</sub> the sequence (VPGAG<sub>[3]</sub>VPGGG<sub>[2]</sub>). </p>
<p>
+
<div class="image-container">
This construct was used by the TU Eindhoven 2023 team to form a hydrogel outside as well as inside <i>E.coli</i> BL21 cells. The goal of this was to stop bacteria from dividing, however many more applications can be thought of for this part. Inspiration for this part was taken from Gradišar et al.<sup>[1]</sup> where they used Leucine zippers form heterodimers. We wanted to use these zippers in a similar way that Fernández‐Colino et al.<sup>[2]</sup> used leucine zippers together with Elastin-Like Polypeptides (ELPs) to form a reversible, injectable hydrogel. Furthermore, we noticed that Contreras-Llano et al.<sup>[3]</sup> were able to form a hydrogel inside of cells, however not with a protein based approach. We wanted to combine these ideas and form a protein based-hydrogel inside of <i>E.coli</i> and use these bacteria for therapeutic applications. However, with the hydrogel itfels, many more applications can be thought of so we wanted to characterize it as well as possible so that other teams can make use of it in the future too. Below you can read more about the experiments we did to characterize this part and about the molecular dynamics simulations that we did.
+
<figure>
 +
<img src="https://static.igem.wiki/teams/4905/wiki/simulatie/elp-gif-25.gif" style="width:100%">
 +
 
 +
<figcaption><p align="justify">
 +
<b>Figure 1: </b> <i>Video with one of the MD simulations at 25°C. The red residues represent the Leucine zipper, the blue residues are the hydrophobic part, and the light-blue residues are the hydrophilic part. </i>
 
</p>
 
</p>
<p>
+
</figcaption></figure>
A schematic overview of what the part looks like is shown in figure 1. It is also easy to ligate other pieces of sequences together to make different constructs, with either different lengths of ELPs, different amino acids or other gBlocks at their end.  
+
<p align="justify">
 +
This construct was used by the TU Eindhoven 2023 team to form a hydrogel extracellular as well as intracellular in <i>E.coli</i> BL21 cells. The goal of this was to stop bacteria from dividing, however many more applications can be thought of for this part. Inspiration for this part was taken from Gradišar et al.<sup>[1]</sup> where they used Leucine zippers to form heterodimers. We wanted to use these zipper domains in a similar way that Fernández‐Colino et al.<sup>[2]</sup> used Leucine zippers together with ELPs to form a reversible, injectable hydrogel. Furthermore, we noticed that Contreras-Llano et al.<sup>[3]</sup> were able to form a hydrogel inside of cells, however not with a protein based approach. We wanted to combine these ideas and form a protein based-hydrogel inside of <i>E.coli</i> and use these bacteria for therapeutic applications. However, with the hydrogel itself, many more applications can be thought of that we wanted to explore as much as possible so that other teams can make use of it in the future too. Below you can read more about the experiments we did to characterize this part and about the molecular dynamics simulations that we did.
 +
</p>
 +
<p align="justify">
 +
A schematic overview of what the part looks like is shown in figure 2. It is simple to ligate other pieces of sequences together to make different constructs, with either different lengths of ELPs, different amino acids or other gBlocks at their end.  
 
</p>
 
</p>
 +
<div class="image-container">
 
<figure><img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/part6-1.png" width="640px">
 
<figure><img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/part6-1.png" width="640px">
  
<figcaption,width="640px">
+
<figcaption>
<p><b>Figure 1 |</b> Schematic overview of the sequence of this construct. (VPGAG<sub>[3]</sub>VPGGG<sub>[2]</sub>) is from now on referred to as A<sub>[5]</sub> and VPGIG will be referred to as I. This created the final name Z1-I<sub>[60]</sub>A<sub>[120]</sub>I<sub>[60]</sub-Z2. This name will besimplified to Z1-A120-Z2. </p>
+
<p><b>Figure 2 |</b> <i> Schematic overview of the sequence of this construct. (VPGAG<sub>[3]</sub>VPGGG<sub>[2]</sub>) is from now on referred to as A<sub>[5]</sub> and VPGIG will be referred to as I. This created the final name Z1-I<sub>[60]</sub>A<sub>[120]</sub>I<sub>[60]</sub>-Z2. This name will besimplified to Z1-A120-Z2. </i></p>
  
 
</figcaption>
 
</figcaption>
</figure><be>
+
</figure>
 +
</div>
  
  
 
<h2>General applications</h2>
 
<h2>General applications</h2>
<p>
+
<p align="justify">
 
ELPs are protein polymers derived from human tropoelastin. One of their key features is that they exhibit a phase separation that is often reversible whereby samples remain soluble below their transition temperature (T<sub>t</sub>) but form coacervates above T<sub>t</sub><sup>[4]</sup>. They have many possible applications in purification, sensing, activation, and nano assembly. Furthermore, they are non-immunogenic, substrates for proteolytic biodegradation, and can be decorated with pharmacologically active peptides, proteins, and small molecules. Recombinant synthesis additionally allows precise control over ELP architecture and molecular weight, resulting in protein polymers with uniform physicochemical properties suited to the design of multifunctional biologics. As such, ELPs have been employed for various uses including as anti-cancer agents, ocular drug delivery vehicles, and protein trafficking modulators<sup>[4]</sup>.
 
ELPs are protein polymers derived from human tropoelastin. One of their key features is that they exhibit a phase separation that is often reversible whereby samples remain soluble below their transition temperature (T<sub>t</sub>) but form coacervates above T<sub>t</sub><sup>[4]</sup>. They have many possible applications in purification, sensing, activation, and nano assembly. Furthermore, they are non-immunogenic, substrates for proteolytic biodegradation, and can be decorated with pharmacologically active peptides, proteins, and small molecules. Recombinant synthesis additionally allows precise control over ELP architecture and molecular weight, resulting in protein polymers with uniform physicochemical properties suited to the design of multifunctional biologics. As such, ELPs have been employed for various uses including as anti-cancer agents, ocular drug delivery vehicles, and protein trafficking modulators<sup>[4]</sup>.
 
</p>
 
</p>
<p>
+
<p align="justify">
Currently ELPs are also being explored in the formation of hydrogels and tissue engineering applications. We are the first iGEM team to use ELPs for the purpose of formation of a hydrogel and have made many different parts while investigating the possibilities. This part turned out to be the best. The reversible behaviour of ELPs is even suitable for injectable hydrogels for example<sup>[2]</sup>. We think this is a part that can have many applications inside as well as outside of cells and we hope that me other iGEM teams will make use of it in the future.
+
Currently ELPs are also being explored in the formation of hydrogels and tissue engineering applications. We are the first iGEM team to use ELPs for the purpose of formation of a hydrogel and have made many different parts while investigating the possibilities. This part turned out to be the most sucessfull. The reversible behaviour of ELPs is even suitable for injectable hydrogels for example<sup>[2]</sup>. We think this is a part that can have many applications inside as well as outside of cells and we hope that many other iGEM teams will make use of it in the future.
  
 
<h2>Construct design</h2>
 
<h2>Construct design</h2>
<p>
+
<p align="justify">
The construct consists of ELPs and two different leucine zippers that have affinity for each other. In general, ELPs have hydrophilic and hydrophobic domains that exhibit reversible phase separation behavior that is temperature-dependent. They are made from a repeating VPGXG sequence, with X replaced by specific amino acids. This results in specific properties of the ELPs, especially related to the T<sub>t</sub> at which the ELPs will interact with each other on the hydrophobic sites<sup>[2]</sup>. When the temperature is below T<sub>t</sub>, the water molecules surrounding the hydrophobic parts will go into the bulk water phase which gains the solvent entropy. This makes it possible to form interactions with other ELP molecules<sup>[3]</sup>.  
+
The construct consists of ELPs and two different Leucine zippers that have affinity for each other. In general, ELPs have hydrophilic and hydrophobic domains that exhibit reversible phase separation behavior that is temperature-dependent. They are made from a repeating VPGXG sequence, with X replaced by specific amino acids. This results in specific properties of the ELPs, especially related to the T<sub>t</sub> at which the ELPs will interact with each other on the hydrophobic sites<sup>[2]</sup>. When the temperature is below T<sub>t</sub>, the water molecules surrounding the hydrophobic parts will go into the bulk water phase which increases the solvent entropy. This makes it possible to form interactions with other ELP molecules<sup>[3]</sup>.  
 
</p>
 
</p>
<p>
+
<p align="justify">
As shown in figure 2, this construct has a hydrophilic region in the middle (A<sub>[120]</sub>) and a hydrophobic region on each side of it (I<sub>[60]</sub>). On the ends the leucine zippers Z1 and Z2 are located for stronger interactions between the ELPs. Leucine zippers consist of a repeating unit that forms an alpha helix. Two leucine zippers together form ion pairs between the helices, which causes association<sup>[4]</sup>. These stronger and reversible interactions make them useful in the formation of a hydrogel at a specific T<sub>t</sub>. In the end, the hydrogel is formed with electrostatic and hydrophobic interactions between the ELPs.  
+
As shown in figure 3, this construct has a hydrophilic region in the middle (A<sub>[120]</sub>) and a hydrophobic region on each side of it (I<sub>[60]</sub>). At the ends the Leucine zippers Z1 and Z2 are located for stronger interactions between the ELPs. Leucine zippers consist of a repeating unit that forms an alpha helix. Two Leucine zippers together form ion pairs between the helices, which causes association<sup>[4]</sup>. These stronger and reversible interactions make them useful in the formation of a hydrogel at a specific T<sub>t</sub>. In the end, the hydrogel is formed with electrostatic and hydrophobic interactions between the ELPs.  
 
</p>
 
</p>
 
+
<div class="image-container">
 
<figure><img src="https://static.igem.wiki/teams/4905/wiki/notebook/parts-goede.png" width="300px">
 
<figure><img src="https://static.igem.wiki/teams/4905/wiki/notebook/parts-goede.png" width="300px">
  
 
<figcaption>
 
<figcaption>
<p><b>Figure 2 |</b> Schematic representation of the composite part, an ELP with leucine zippers on the ends. It has a triblock ELP sequence that can form a thermoresponsive hydrogel. </p>
+
<p><b>Figure 3 |</b> <i> Schematic representation of the composite part, an ELP with Leucine zippers on the ends. It has a triblock ELP sequence that can form a thermoresponsive hydrogel. </i></p>
 
</figcaption>
 
</figcaption>
</figure><be>
+
</figure>
 +
</div>
  
<p>
+
<p align="justify">
As soon as the hydrogel is formed inside <i>E.coli</i> BL21 cells, the cells are prevented from dividing. However, the cells remain functional. So they can still be used to express therapeutic agents, like Interleukin 10 in the TU Eindhoven 2023 teams project. As mentioned before, this is just one of the many applications of this part.  
+
The goal of our project specifically was to use this part to stop the cell from dividing and to co-express a therapeutic protein (Interleukin 10).  
 
</p>
 
</p>
<p>
+
<p align="justify">
Another part that is often mentioned on this page is <a style="color:#F6B227" href="https://parts.igem.org/Part:BBa_K4905008"> BBa_K4905008</a>. This part was often used as a control. The ELP domains of the protein is the same as for the part on this page, however, part <a style="color:#F6B227" href="https://parts.igem.org/Part:BBa_K4905008"> BBa_K4905008</a> has two Leucine zippers that are not complmentary, which means that the interactions between the proteins will be weaker and this way, it could be studies what the effect of only the ELP triblocks was.
+
Another part that is often mentioned on this page is <a style="color:#F6B227" href="https://parts.igem.org/Part:BBa_K4905008"> BBa_K4905008</a>. This part was often used as a control. The ELP domains of the protein is the same as for the part on this page, however, part <a style="color:#F6B227" href="https://parts.igem.org/Part:BBa_K4905008"> BBa_K4905008</a> has two Leucine zippers that are not complementary, which means that the interactions between the proteins will be weaker and this way, it could be studied what the effect of only the ELP triblocks was. This control part will be reffered to as Z2-A120-Z2 on the rest of this page. Lastly, another protein that can be seen a lot on this page is abbreviated as A120, which is the same sequence as the part on this page, but without the Leucine zippers. So it is part <a style="color:#F6B227" href="https://parts.igem.org/Part:BBa_K4905001">BBa_K4905001</a> ligated to itself. This protein also serves as a control group in many experiments.
  
 
<h1>Results</h1>
 
<h1>Results</h1>
 
<h3>Protein expression and purification</h3>
 
<h3>Protein expression and purification</h3>
<p> This part was expressed in BL2 <i>E.coli</i>. The full plasmid that the DNa was cloned into can be seen below in figure 3. To check its expression we ran an SDS-PAGE gel of which the results can be seen in figure 4. The bands on the gel showed up where we expected, so this means that protein expression was succesfull.  
+
<p> This part was expressed in BL21 <i>E.coli</i>. The full plasmid that the DNA was cloned into can be seen below in figure 4. To check its expression we ran an SDS-PAGE gel of which the results can be seen in figure 5. The bands on the gel showed up where we expected, so this means that protein expression was succesfull.  
<figure><img src="https://static.igem.wiki/teams/4905/wiki/notebook/compositepartprizeplasmid.png" width="300px">
+
</p>
 +
 
 +
<div class="image-container">
 +
<figure><img src="https://static.igem.wiki/teams/4905/wiki/notebook/compositepartprizeplasmid.png" width="400px">
  
 
<figcaption>
 
<figcaption>
<p><b>Figure 3 |</b> Full plasmid containing the sequence for this part in a pET24(+) vector that was transformed into <i>E.coli</i> BL21 cells for protein expression.</p>
+
<p><b>Figure 4 |</b> <i>Full plasmid containing the sequence for this part in a pET24(+) vector that was transformed into <i>E.coli</i> BL21 cells for protein expression.</i></p>
 
</figcaption>
 
</figcaption>
</figure><be>
+
</figure>
 +
</div>
  
 +
<div class="image-container">
 
<figure><img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/z1-a120-z2-sds-page.png" width="300px">
 
<figure><img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/z1-a120-z2-sds-page.png" width="300px">
  
 
<figcaption>
 
<figcaption>
<p><b>Figure 4 |</b> SDS page gel run with this part. The band can be seen at the expected weight of 105.1 kDa.</p>
+
<p><b>Figure 5 |</b> <i>SDS page gel run with this part. Samples were taken after 12, 20 and 36 hours of inducing protein expression. The band can be seen at the expected weight of 105.1 kDa.</i></p>
 
</figcaption>
 
</figcaption>
</figure><be>
+
</figure>
 
+
</div>
  
 
<h2>Characterization</h2>
 
<h2>Characterization</h2>
 
<h3>Hydrogel formation</h3>
 
<h3>Hydrogel formation</h3>
 
<p align="justify">
 
<p align="justify">
All gels were formed by resuspending the ELPs in MilliQ (MQ) water at different w/v%.  The leucine zipper-containing ELPs were expected to instantaneously form a gel upon increasing the temperature above T<sub>t</sub>, as has also been reported in the literature<sup>[5]</sup>.
+
First, we wanted to find out if we could use our part to form a hydrogel. All gels were formed by resuspending the ELPs in MilliQ (MQ) at different weight/volume%(w/v%).  The Leucine zipper containing ELPs were expected to instantaneously form a gel upon increasing the temperature above T<sub>t</sub>.
 
</p>
 
</p>
 
<p align="justify">
 
<p align="justify">
We saw that for the constructs containing the leucine zippers, there formed a viscous fluid upon dissolving the ELPs. Subsequently, after warming the solution to room temperature, it became very turbid and a gel started to form. This is supported by the images below in figure 5.
+
We saw that for both Z1-A120-Z2 and Z2-A120-Z2, a viscous fluid formed upon dissolving the ELPs in MQ water at 4°C. Subsequently, after warming the solution to room temperature, the liquid became very turbid and a gel started to form. Images of the gels can be seen in figure 6.
 
</p><br>
 
</p><br>
  
 
+
<div class="image-container">
 
<div class="row">
 
<div class="row">
  <div class="col">
+
 
     <img src="https://static.igem.wiki/teams/4905/wiki/5-gel-black.png" style="width:50vw;" class="results-img enlarge-image" alt="Image 1">
+
     <img src="https://static.igem.wiki/teams/4905/wiki/5-gel-black.png" style="width:49%;" class="results-img enlarge-image" alt="Image 1">
  </div>
+
 
  <div class="col">
+
 
     <img src="https://static.igem.wiki/teams/4905/wiki/10-gels-black.png" style="width:50vw;" class="results-img enlarge-image" alt="Image 2">
+
     <img src="https://static.igem.wiki/teams/4905/wiki/10-gels-black.png" style="width:49%;" class="results-img enlarge-image" alt="Image 2">
  </div>
+
 
 
</div><br>
 
</div><br>
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 5 |</b> Hydrogel solutions at 5 and 10 weight/volume % (w/v%). It can be seen that their solubility changes when the transition temperature is surpassed. Above this temperature a hydrogel is formed and the substance in the tube becomes cloudy and gel-like.</figcaption><br>
+
<figcaption><b>Figure 6 |</b> <i> Hydrogel solutions at 5 and 10 (w/v%). It can be seen that their solubility changes when the transition temperature is surpassed. Above this temperature a hydrogel is formed and the substance in the tube becomes cloudy and gel-like (w/v%).</i></figcaption><br>
 +
</div>
  
 
<p align="justify">
 
<p align="justify">
The construct containing Z1-A120-Z2 formed a gel at 5% and 10% (w/v%) as described above. Upon cooling the gel down to 4°C again, the gel remained intact.   As a control, Z2-A120-Z2 that had non-complementary leucine zippers was also dissolved to form solutions of 5% and 10% (w/v%). This construct too, was able to form a gel instantly upon heating to RT. However, in contrast to the leucine zipper-ELPs, the gel disassembled again when the solution was cooled to 4°C. This indicates that the latter of the two constructs shows reversibility, which is expected when the ELPs are brought below their T<sub>t</sub>, while the former construct does not show this reversibility.  
+
The construct containing Z1-A120-Z2 formed a gel at 5% and 10% (w/v%) as described above. Upon cooling the gel down to 4°C again, the gel remained intact. As a control, Z2-A120-Z2 (<a style="color:#F6B227" href="https://parts.igem.org/Part:BBa_K4905008"> BBa_K4905008</a>) was also dissolved to form solutions of 5% and 10% (w/v%). This construct too, was able to form a gel instantly upon heating to RT. However, in contrast to the Leucine zipper-ELPs, the gel disassembled again when the solution was cooled to 4°C. This indicates that the latter of the two constructs shows reversibility, which is expected when the ELPs are brought below their T<sub>t</sub>, while the former construct does not show this reversibility.
 
</p>
 
</p>
  
 
<h3>Transition temperature determination with UV-Vis</h3>
 
<h3>Transition temperature determination with UV-Vis</h3>
As mentioned before, the transition temperature of the hydrogel is influenced by the amino acids that are used in the ELP, furthermore we expected that the leucine zippers would strengthen the interaction between the peptides and cause a lower transition temperature. The transition temperatures of our constructs Z1-A120-Z2 and, as a control to see the transition temperature without the effect of leucine zippers, Z2-A120-Z2 were determined with the use of UV-visible spectroscopy. Solutions of 5 and 20 μM were made in MQ water. Absorbances were measured at 350 nm. T<sub>t</sub>s were determined to be around 23.5°C and 21°C for both constructs, respectively. This is remarkable, since the expectation was that the construct containing Z1 and Z2, two leucine zippers that are be able to dimerize, would have a lower T<sub>t</sub>, with the reasoning that the interaction of these domains would cause the ELPs to become insoluble sooner. However, it is possible that the Leucine zippers already bring the hydrophobic domains in a different conformation which causes the transition temperature to actually become higher.
+
<p align="justify">
 +
The transition temperatures (T<sub>t</sub>s) Z1-A<sub>[120]</sub>-Z2 and control Z2-A<sub>[120]</sub>-Z2 were determined with the use of UV-visible spectroscopy. Solutions of 5 and 20 μM were made in MQ water. Absorbances were measured at 350 nm. T<sub>t</sub>s were determined to be 23.5°C and 21°C for both constructs, respectively (Figure 7).  
 +
</p>
 +
<p>
 +
When the control protein Z2-A<sub>[120]</sub>-Z2 crosses the T<sub>t</sub> at 21°C, aggregates form leading to a sharp increase in absorbance. The absorbance then drops, which might be explained by the non-complementary two zipper domains repelling each other, leading to disassembly of the aggregates. Upon further increasing the temperature, the hydrophobic interactions become stronger, which further increases the absorbance, until finally even the hydrophilic block becomes insoluble above 50°C.
 +
</p>
 +
<p>
 +
Z1-A<sub>[120]</sub>-Z2 shows a different temperature dependent behavior compared to Z2-A<sub>[120]</sub>-Z2, in that it shows a smaller increase in absorbance when it crosses T<sub>t</sub> of 23.5°C. We think the interaction between Z1 and Z2 stabilizes the aggregates, forming into a branched network of ELPs, even at low concentrations of 5 μM.
 
</p>
 
</p>
  
  
<head>
 
    <meta charset="UTF-8">
 
    <meta name="viewport" content="width=device-width, initial-scale=1.0">
 
    <title>Two Images Side by Side</title>
 
    <style>
 
        .image-container {
 
            display: flex;
 
        }
 
  
        .image-container img {
 
            width: 50%; /* Adjust this value to control the width of the images */
 
            height: auto; /* Maintain aspect ratio */
 
        }
 
    </style>
 
</head>
 
 
     <div class="image-container">
 
     <div class="image-container">
         <img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/uv-vis-5-um-350-nm.png" alt="Image 1">
+
         <img style="width:49%;" src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/uv-vis-5-um-350-nm.png" alt="Image 1">
         <img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/uv-vis-20-um-350-nm.png" alt="Image 2">
+
         <img style="width:49%;" src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/uv-vis-20-um-350-nm.png" alt="Image 2">
 
     </div>
 
     </div>
  
 
<br>
 
<br>
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 6 |</b> UV-Vis spectra for constructs Z1-A120-Z2 and Z2-A120-Z2 measured at 350 nm. The solution was heated at a rate of 0.5 degrees celcius per minute. It can be seen that the transition temperature lays around 23.5 and 21 degrees Celcius for the two proteins respectively. This is the case for 5 uM as well as when a higher concentration of 20 uM is used, which inducated that this behaviour is independent of the concentration of protein used. 
+
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 7 |</b> <i>UV-Vis spectra for constructs Z1-A120-Z2 and Z2-A120-Z2 measured at 350 nm. The solution was heated at a rate of 0.5 degrees Celcius per minute. It can be seen that the transition temperature lays at 23.5 and 21 degrees Celcius for the two proteins respectively. </i>
 
</figcaption><br>
 
</figcaption><br>
  
Line 134: Line 166:
  
 
<p align="justify">
 
<p align="justify">
Absorbances were also measured upon cooling the solution back down to find out whether the interaction was reversible. Once the solution goes below T<sub>t</sub>, the absorbances of returned back to base levels, indicating that the ELPs can be reversibly precipitated and resuspended again. Furthermore, the spectra below look identical, which further confirms that the interactions are all reversible.
+
Absorbances were also measured upon cooling the solution back down. Once the solution goes below T<sub>t</sub>, the absorbances of returned back to base levels, indicating that the ELPs can be reversibly precipitated and resuspended again (Figure 8) at these low concentrations.  
 
</p><br>
 
</p><br>
  
<head>
 
    <meta charset="UTF-8">
 
    <meta name="viewport" content="width=device-width, initial-scale=1.0">
 
    <title>Two Images Side by Side</title>
 
    <style>
 
        .image-container {
 
            display: flex;
 
        }
 
 
        .image-container img {
 
            width: 50%; /* Adjust this value to control the width of the images */
 
            height: auto; /* Maintain aspect ratio */
 
        }
 
    </style>
 
</head>
 
 
     <div class="image-container">
 
     <div class="image-container">
         <img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/5-um-cooling-down-350nm.png" alt="Image 3">
+
         <img style="width:49%;" src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/5-um-cooling-down-350nm.png" alt="Image 3">
         <img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/20-um-cooling-down-350nm.png" alt="Image 4">
+
         <img style="width:49%;" src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/20-um-cooling-down-350nm.png" alt="Image 4">
 
     </div>
 
     </div>
  
 
<br>
 
<br>
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 7 |</b> UV-Vis spectra for constructs Z1-A120-Z2 and Z2-A120-Z2 measured at 350 nm. The solution was cooled at a rate of 2 degrees celcius per minute. It can be seen that the transition temperature lays around 23.5 and 21 degrees Celcius for the two proteins respectively. This is the case for 5 uM as well as when a higher concentration of 20 uM is used, which inducated that this behaviour is independent of the concentration of protein used. 
+
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 8 |</b> <i>UV-Vis spectra for constructs Z1-A120-Z2 and Z2-A120-Z2 measured at 350 nm. The solution was cooled at a rate of 2 degrees celcius per minute. It can be seen that the transition temperature lays around 23.5 and 21 degrees Celcius for the two proteins respectively. </i>
 
</figcaption><br>
 
</figcaption><br>
  
 
<h3>Dye release from hydrogel</h3>
 
<h3>Dye release from hydrogel</h3>
 
<p align="justify">
 
<p align="justify">
To determine the relative ability for diffusion of small molecules through our hydrogels, we made gels at different percentages containing the fluorescent dye rhodamine B. In the future, apart from dyes, it could also be interesting to look at drug diffusion out of the hydrogel.
+
To determine the relative ability for the diffusion of small molecules through our hydrogels, we made gels at two different weight-to-volume percentages containing the fluorescent dye rhodamine B. Finding out whether small molecules can gradually diffuse out of the hydrogel might be interesting for future applications, a drug could by incorporated for instance.
First, a calibration curve was made to determine at which concentration rhodamine B should be added to the gel. Figure 8 shows a linear relationship up to 26 μM of Rhodamine B, with equation: y = 2*106x + 963.16 (R² = 0.993).  We decided to make the gels with a Rhodamine B concentration of 104 μM, since the amount of dye diffusing out of the gels would be significantly lower than this amount.  
+
</p>
The gels were made at a volume of 100 μL. They were then brought to room temperature (RT) and washed with 1.5 mL of warm MQ. Then, 1 mL of warm MQ was pipetted on top of the gels, and they were incubated up to 8 hours at RT, with samples of 100 μL being taken every 2 hours. After measuring the fluorescence intensity, the values were corrected for the reduction in volume upon taking the samples.
+
<p align="justify">
 +
First, a calibration curve was made to determine at which concentration rhodamine B should be added to the gel. Figure 9 shows a linear relationship up to 26 μM of Rhodamine B. We decided to make the gels with a Rhodamine B concentration of 104 μM (0.05 mg/mL), since the amount of dye diffusing out of the gels was expected to be far lower than this amount. A 5% gel as well as a 10% gel were tested. Our hypothesis was that diffusion through the 5% gel would be faster than when a higher weight percentage of the ELPs was used,
 +
 
</p><br>
 
</p><br>
  
<head>
 
    <meta charset="UTF-8">
 
    <meta name="viewport" content="width=device-width, initial-scale=1.0">
 
    <title>Two Images Side by Side</title>
 
    <style>
 
        .image-container {
 
            display: flex;
 
        }
 
  
        .image-container img {
+
</p><br>
            width: 50%; /* Adjust this value to control the width of the images */
+
 
            height: auto; /* Maintain aspect ratio */
+
        }
+
    </style>
+
</head>
+
 
     <div class="image-container">
 
     <div class="image-container">
         <img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/rhodamine-b-calibration-curve.png" alt="Image 5">
+
         <img style="width:49%;" src="https://static.igem.wiki/teams/4905/wiki/rhodb-calibration.png" alt="Image 3">
         <img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/rhodamine-gels.png" alt="Image 6">
+
         <img style="width:49%;" src="https://static.igem.wiki/teams/4905/wiki/linear-chart.png" alt="Image 4">
 
     </div>
 
     </div>
 
+
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 9:</b><i>  Left: Calibration curve for Rhodamine B, a linear relationship can be seen up to 26 μM of Rhodamine B, with equation: y = 2*106x + 963.16 (R² = 0.993) Right: The hydrogels made with Rhodamine B at a concentration of 104 μM. </i>
<br>
+
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 8 |</b> Left: Calibration curve for Rhodamine B, a linear relationship can be seen up to 26 μM of Rhodamine B, with equation: y = 2*106x + 963.16 (R² = 0.993) Right: The hydrogels made with Rhodamine B at a concentration of 104 μM.  
+
 
</figcaption><br>
 
</figcaption><br>
  
<p align="justify">
 
The hydrogel consisting with a w/v% of 10% shows a linear relationship between fluorescence intensity over time between t=0 and t=8, with a linear equation of: y = 3310.8x + 5222.4.  Therefore the diffusion rate, as calculated from the calibration curve, is 3.7*10-3 mg/mL/h. The 5% gel shows a faster diffusion rate , which appears to be non-linear between t=0 and t=8. This can be attributed to the fact that the supernatant on top of the gel came close to being saturated, which is likely also the case for the 10% gel, but at a later timepoint. These results can be seen in figure 9 below.
 
</p>
 
  
 
<p align="justify">
 
<p align="justify">
The faster diffusion rate of the 5% gel is in line with our hypothesis that a lower weight percentage gel is able to release small molecules faster than a 10% gel. Since Rhodamine B is a much smaller molecule than a cytokine like IL-10, it would be interesting to see how the diffusion of such a larger molecule would be affected by our gels. The weight percentage of the gel might therefore not only affect the division rate of the bacteria, but it can also play a role in the release rate of the therapeutic.  
+
The gels were made at a volume of 100 μL. They were then brought to RT and washed with 1.5 mL of warm MQ. Then, 1 mL of warm MQ was pipetted on top of the gels, and they were incubated for up to 8 hours at RT, with samples of 100 μL being taken every 2 hours. After measuring the fluorescence intensity, the values were corrected for the reduction in volume upon taking the samples. From figure 10, it becomes clear that the dye diffuses faster from the 5% gel than from the 10% gel. However, more experiments are needed to truly find the rate of diffusion, it is hard to draw a conclusion with limited data points.
 +
The faster diffusion rate of the 5% gel is in line with our hypothesis that a lower weight percentage gel is able to release small molecules faster than a 10% gel.  
 
</p><br>
 
</p><br>
  
<img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/rhodamine-b-5-vs-10.png" style="width:40vw;" class="results-img enlarge-image" alt="Image 11"><br>
+
    <div class="image-container">
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 9 |</b> Fluorescence intensity measured from the supernatant on top of the gel over a period of time.The hydrogel consisting with a w/v% of 10% shows a linear relationship between fluorescence intensity over time between t=0 and t=8, with a linear equation of: y = 3310.8x + 5222.4.  Therefore the diffusion rate, as calculated from the calibration curve, is 3.7*10-3 mg/mL/h. The 5% gel shows a faster diffusion rate , which appears to be non-linear between t=0 and t=8.</figcaption><br>
+
        <img style="width:49%;" src="https://static.igem.wiki/teams/4905/wiki/dye-diff-10-vs-5-g.png" alt="Image 3">
 +
        <img style="width:49%;" src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/rhodamine-gels.png" alt="Image 4">
 +
    </div>
 +
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 10 |</b> <i>Fluorescence intensity measured from the supernatant on top of the gel over a period of time.</i></figcaption><br>
  
  
 
<h3>CCK8</h3>
 
<h3>CCK8</h3>
 
<p align="justify">  
 
<p align="justify">  
Cell counting kit 8 (CCK8) is an assay which is commonly used to determine the number of alive cells in a sample. The kit uses a highly water-soluble tetrazolium salt, WST-8, which produces a formazan dye upon reduction in the presence of an electron mediator.   The amount of the formazan generated by dehydrogenases is directly in proportion to the number of living cells. Yang et al. (2021) have demonstrated that this assay is also suitable for the detection of living E.coli<sup>[6]</sup>.
+
We wanted to find a way that we could measure the amount of cells that we can stop from dividing with our hydrogel and whether they would stay alive. We came across Cell counting kit 8 (CCK8). This is an assay that is commonly used to determine the number of alive cells in a sample. The kit uses a highly water-soluble tetrazolium salt, WST-8, which produces a formazan dye upon reduction in the presence of an electron mediator. The amount of formazan generated by dehydrogenases is directly in proportion to the number of living cells. Yang et al. (2021) have demonstrated that this assay is also suitable for the detection of living <i>E. coli</i><sup>[5]</sup>.
 
</p>
 
</p>
 
<p align="justify">  
 
<p align="justify">  
We used this assay to test whether our hydrogelated cells would be more resistant to several conditions, including incubation with ampicillin, H<sub>2</sub>O<sub>2</sub>, and freeze-drying. Ampicillin prevents the synthesis of the bacterium’s cell wall. It does this by binding to enzymes necessary for the formation of the cell wall<sup>[7]</sup>. Our hypothesis is that if our hydrogel prevents the bacteria from dividing, they should be affected less by the presence of this antibiotic, and therefore show a higher survival rate than for bacteria which do not contain the hydrogel.  
+
We used this assay to test whether our hydrogelated cells would be more resistant to several conditions, including incubation with ampicillin and freeze-drying. Our biggest goal was to find out if cells would stop dividing, but stay functional with our hydrogel. To test this, we came up with the idea of combining the CCK8 assay described above with the antibiotic ampicillin. Ampicillin prevents the synthesis of the bacterium’s cell wall. It does this by binding to enzymes necessary for the formation of the cell wall <sup>[6]</sup>. This means it only kills bacteria when they divide. Our hypothesis is that if our hydrogel prevents the bacteria from dividing, they should be affected less by the presence of this antibiotic, and therefore show a higher survival rate than bacteria that do not contain the hydrogel since these will still divide and be killed by the antibiotic. A schematic representation of the experiment can be seen below in figure 11.
 
</p>
 
</p>
<p align="justify">
+
 
We also wanted to test whether our bacteria would have a higher survival rate when exposed to oxidizing conditions, since inflammatory bowel disease (IBD) is known to cause oxidative stress in the gut<sup>[8,9]</sup>.
+
<div class="image-container">
</p>
+
<img src="https://static.igem.wiki/teams/4905/wiki/cck8zwart.png" class="pd-img enlarge-image-width" style="width:90%"><br>
<p align="justify">
+
</div>
Lastly, since the concept of our treatment is to put them in a capsule which can be orally administered, the bacteria need to be able to survive freeze-drying, so that they can produce the therapeutic locally when they are reconstituted in the GI tract.
+
 
 +
<figcaption style="text-align:center; font-size: 12px;"> <b>Figure 11: </b><i>A schematic representation of the CCK8 assay protocol. Samples were taken and incubated with or without ampicillin. After 1 hour, the cell pellets were collected and resuspended in PBS. Then the samples were incubated again with the substrate WST-8. The amount of produced fluorescent product that is measured via the <sub>OD450</sub> with a plate reader is directly proportional to the amount of cells in the sample.</i></figcaption><br>
 +
 
 +
<p>
 +
For this assay, it was nessecary to measure the OD<sub>450</sub> with a plate reader. This can directly be related to the OD<sub>600</sub><sup>[5]</sup>. For this relationship, the calibration curve in figure 12 was made.
 
</p>
 
</p>
 +
<div class="image-container">
 +
<img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/cck8-calibration-curve.png" style="width:49%;" class="results-img enlarge-image" alt="Image 7"><br>
 +
<figcaption style="text-align:center; font-size: 12px;"> <b>Figure 12: </b><i> Calibration curve relating the OD<sub>600</sub> to the OD<sub>450</sub>.</i></figcaption><br>
 +
</div>
  
 
<h3>Ampicillin resistance</h3>
 
<h3>Ampicillin resistance</h3>
  
 
<p align="justify">
 
<p align="justify">
Three different conditions were tested over a period of time. In one sample, the ELPs without any crosslinker (A120) were expressed and samples were taken at 12h, 20h and 36h of protein expression. This same procedure was followed for the other sample, in which the ELPs including the leucine zippers were expressed (Z1-A120-Z2). As a control, a culture of bacteria was also made in which no proteins were expressed (Figure 10). Before incubation, all samples started at a baseline of OD<sub>600</sub> = 0.5. They were then either incubated with ampicillin (1 mg/mL) or without for 1h at 37°C.
+
Three different conditions were tested over a period of time. Namely; no (protein) expression, ELPs with no cross linkers (A120), and ELPs with Leucine zippers (Z1-A120-Z2). Samples were taken at 12h, 20h, and 36h of protein expression. Before incubation ampicillin, samples were diluted to at a baseline of OD<sub>600</sub> = 0.5. They were then either incubated with ampicillin (1 mg/mL) or without for 1h at 37°C. The results are summarized in Table 1 and Figure 13.
 
</p>
 
</p>
 +
<center><table>
 +
  <tr>
 +
    <th>WST-8 conversion Amp/No amp</th><th>No expression</th><th>A120</th> <th>Z1-A120-Z2</th>
 +
</tr>
 +
<tr> 
 +
<td> 12h</td><td>0.71 ± 0.05 </td><td>1.02 ± 1.00</td> <td>3.05 ± 0.33</td>
 +
</tr>
 +
<tr>
 +
<td>20h</td><td>0.60 ± 0.21</td><td>5.81 ± 0.22</td> <td>20.85 ± 4.14</td>
 +
</tr>
 +
<tr>
 +
<td>36h</td><td>0.28 ± 0.08</td><td>6.79 ± 0.17</td> <td>17.05 ± 2.13</td> 
 +
</tr>
 +
 +
</tr>
 +
</table></center><br>
 +
<figcaption style="text-align:center; font-size: 12px;"><b>Table 1:</b><i> Measured WST-8 conversion for samples with no expression, A120 and Z1-A120-Z2 at different timepoints. </i></figcaption><br>
 +
<p align="justify">
 +
 
<p align="justify">
 
<p align="justify">
<i>E.coli</i> where no protein expression was induced, and which were treated with ampicillin, showed at most a conversion of 0.71% ± 0.05% compared to cells which were not treated. E.coli where expression of A120 was induced, resulted in a maximum conversion of 6.80% ± 0.17% after 36h of protein expression. This is likely due to the fact that this construct is also able to form a gel, although less strong compared to the Z1-A120-Z2 construct, leading to an attenuation in cell division. Finally, Z1-A120-Z2 was expressed for the three timepoints. A significant increase was observed in the conversion of substrate, most notably after 20h of expression (20.86% ± 4.14%). We had expected to see a higher percentage after 36h of expression, since at this timepoint, more bacteria should be gelated and thus unable to divide. That being said, The standard deviation of the samples is quite high, with this small number of measurements, it is hard to say whether the values would be similar given a larger sample size. Other explanations could be that some of the cells were over-gelated, which might have led to a premature death unrelated to the incubation with ampicillin. Additionally, cell death could also have occurred due to nutrients in the culturing medium running out.  
+
Ampicillin-treated <i>E. coli</i> that were not induced to express proteins, and which were treated with ampicillin, showed at most a WST-8 conversion of 0.71% ± 0.05% compared to cells which were not treated, as measured by absorbance at 450 nm. In contrast, <i>E. coli</i>  where expression of A120 was induced, showed a maximum conversion of 6.80% ± 0.17% after 36h of protein expression. This is likely due to the fact that this construct is also able to form a gel, although less strong compared to the Z1-A120-Z2 construct, leading to an attenuation in cell division. Finally, Z1-A120-Z2 was expressed for the three-time points. An increase was observed in the conversion of substrate, most notably after 20h of expression (20.86% ± 4.14%). We expected to see a higher percentage after 36h of expression, since at this time point, more bacteria should be gelated and thus unable to divide. The standard deviation of the samples is quite high, with this small number of measurements and large deviation, it is hard to tell whether there is really a difference between these two time points.  Overall, a larger sample size could help to reduce the variability of this assay in future tests. Other explanations could be that some of the cells were over-gelated, which might have led to a premature death unrelated to the incubation with ampicillin. Additionally, cell death could also have occurred due to nutrients in the culturing medium running out.  
 
</p>
 
</p>
 
<p align="justify">
 
<p align="justify">
For comparison, OD<sub>600</sub> of all samples was also measured after incubation with or without ampicillin. When the samples were incubated with ampicillin, they all decreased in OD<sub>600</sub>.  The samples incubated without ampicillin in which ELPs were expressed show a relatively stable OD<sub>600</sub> compared to the sample in which no expression was induced. This was the only sample in which a significant increase in OD<sub>600</sub> could be observed. </p><br>
+
For comparison, OD<sub>600</sub> of all samples was also measured after incubation with and without ampicillin. When the samples were incubated with ampicillin, they all decreased in OD<sub>600</sub>.  The samples incubated without ampicillin in which ELPs were expressed show a relatively stable OD<sub>600</sub> compared to the sample in which no expression was induced. This was the only sample in which a substantial increase in OD<sub>600</sub> could be observed. </p><br>
<head>
+
 
    <meta charset="UTF-8">
+
    <meta name="viewport" content="width=device-width, initial-scale=1.0">
+
    <title>Two Images Side by Side</title>
+
    <style>
+
        .image-container {
+
            display: flex;
+
        }
+
  
        .image-container img {
 
            width: 50%; /* Adjust this value to control the width of the images */
 
            height: auto; /* Maintain aspect ratio */
 
        }
 
    </style>
 
</head>
 
 
     <div class="image-container">
 
     <div class="image-container">
         <img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/cck8-measurement-4.png" alt="Image 7">
+
         <img style="width:49%;" src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/cck8-measurement-4.png" alt="Image 7">
         <img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/od600.png" alt="Image 8">
+
         <img style="width:49%;" src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/od600.png" alt="Image 8">
 
     </div>
 
     </div>
  
  
 
<br>
 
<br>
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 10 |</b> Left: CCK8 conversion of bacteria treated with ampicillin versus untreated bacteria. Right: corresponding OD<sub>600</sub> values. Measurements were done at differrent timepoints (12h, 24h and 36h). It can be seen that When the samples were incubated with ampicillin, they all decreased in OD<sub>600</sub>. The samples incubated without ampicillin in which ELPs were expressed show a relatively stable OD<sub>600</sub> compared to the sample in which no expression was induced.
+
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 13 |</b> <i>Left) WST-8 conversion of bacteria treated with ampicillin versus untreated bacteria. right) corresponding OD<sub>600</sub> values. Measurements were done at different timepoints (12h, 20h, and 36h). It can be seen that when the samples were incubated with ampicillin, they all decreased in OD<sub>600</sub>. The samples incubated without ampicillin in which ELPs were expressed show a relatively stable OD<sub>600</sub> compared to the sample in which no expression was induced.</i>
 
</figcaption><br>
 
</figcaption><br>
  
 
<p align="justify">
 
<p align="justify">
The OD600 was also measured over time for each of the three samples. After 20h of expression, these were diluted to OD = 0.1 and then measured every 30 minutes. From the graph, it becomes clear that when no expression of ELPs was induced, the bacteria divided more rapidly then when expression was induced. This can be seen in figure 11. These findings are in line with the results from the ampicillin treatment. It seems that the presence of the ELPs does not altogether inhibit division, but it does slow it down, which seems to be the reason why they are killed slower upon addition of ampicillin. That being said, it is hard to draw a conclusion as to whether the attenuation of bacterial division is driven by the formation of a gel inside the cytosol, or simply because of the crowded environment inside the cell, caused by the presence of a high amount of protein. Therefore, additional controls could be added, like samples in which a mono or diblock ELP is expressed, in order to account for molecular crowding.  
+
The OD<sub>600</sub> was also measured without addition of ampicillin over time for each of the three samples (No expression, A120 and Z1-A120-Z2). After 20h of expression, these were diluted to OD = 0.1 and then measured every 30 minutes. From figure 14, it becomes clear that when no expression of ELPs was induced, the bacteria divided more rapidly than when expression was induced. These findings are in line with the results of the ampicillin treatment. It seems that the presence of the ELPs does not altogether inhibit division, but it does slow it down, which seems to be the reason why they are killed slower upon the addition of ampicillin. However, it is hard to draw a conclusion as to whether the attenuation of bacterial division is driven by the formation of a gel inside the cytosol, or simply because of the crowded environment inside the cell, caused by the presence of a high amount of protein. Therefore, additional controls need to be added in future experiments, like samples in which a mono-or diblock ELP is expressed, in order to account for molecular crowding. Furthermore, a sample where a completely different protein is also interesting, to see what the effect of protein expression is on cell division.  
 
</p><br>
 
</p><br>
  
<img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/od600-over-time.png" style="width:25vw;" class="results-img enlarge-image" alt="Image 6"><br>
+
<div class="image-container">
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 11 |</b> OD600 measurements of the three samples over a period of time. Protein expression in two of the three samples was induced for 20h, after which the samples were diluted to OD = 0.1 and measured every 30 min. It can be seen that no expression of ELPs was induced, the bacteria divided more rapidly then when expression was induced.</figcaption><br>
+
<img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/od600-over-time.png" style="width:49%;" class="results-img enlarge-image" alt="Image 6"><br>
 +
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 14 |</b> <i>measurements of three samples over a period of time. Protein expression in two of the three samples was induced for 20h, after which the samples were diluted to OD = 0.1 and measured every 30 min. Expression of gel-forming ELPs greatly reduces the rate of replication compared to cells without expression.</i></figcaption><br>
 +
</div>
  
 
<p align="justify">
 
<p align="justify">
These results were compared to the classic colony counting method. From each sample, a diltion of 105 was made and then plated out on an LB-agar plate. The plates were incubated at 37°C overnight. Images of the plates were taken (Figure 12) and the number of colonies was counted using the colony counter tool from ImageJ. The results of which are summarized in Table 1 As expected, the wild-type <i>E.coli</i> treated with ampicillin did not form any colonies, whereas untreated bacteria formed several hundreds of colonies. Samples containing A120 and Z1-A120-Z2 both showed a significantly smaller amount of colonies.</p><br>
+
These results were compared to the classic colony counting method. From each sample, a dilution of 10<sup>5</sup> was made and then plated out on an LB-agar plate. The plates were incubated at 37°C overnight. Images of the plates were taken (Figure 15) and the number of colonies was counted using the colony counter tool from ImageJ. The results of which are summarized in Figure 15. As expected, the wild-type <i>E. coli</i> treated with ampicillin did not form any colonies, whereas untreated bacteria formed several hundreds of colonies. Samples containing A120 and Z1-A120-Z2 both showed a significantly smaller amount of colonies.</p><br>
 
+
<img src="https://static.igem.wiki/teams/4905/wiki/all-plates-black1.png" style="display: block; margin-left: auto; margin-right: auto; width: 50vw;"><br>
+
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 12 |</b> Photos of all LB-agar plates after incubation at 37°C overnight at different conditions.</figcaption><br>
+
  
 +
<div class="image-container">
 +
<img src="https://static.igem.wiki/teams/4905/wiki/all-plates-black1.png" style="display: block; margin-left: auto; margin-right: auto; width: 90%;"><br>
 +
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 15 |</b> <i>Photos of all LB-agar plates after incubation at 37°C overnight, including a a graph and table in which the amount of colonies is summarized. It can be seen that the bacteria which do not contain the ELPs form a greater amount of colonies than when expression of these proteins was induced.</i></figcaption><br>
 +
</div>
 +
<div class="image-container">
 
<table>
 
<table>
 
   <tr>
 
   <tr>
Line 285: Line 309:
 
</tr>
 
</tr>
 
</table><br>
 
</table><br>
 +
</div>
  
<figcaption style="text-align:center; font-size: 12px;"><b>Table 1 |</b> Table with the amount of colonies on each plate</figcaption><br>
+
<figcaption style="text-align:center; font-size: 12px;"><b>Table 2 |</b> <i>Table with the amount of colonies on each plate</i></figcaption><br>
  
 
<h3>Microscopy (FRAP)</h3>
 
<h3>Microscopy (FRAP)</h3>
 
<p align="justify">
 
<p align="justify">
Characterization of our gelated cells with the use of microscopy required us to design a model protein which could interact with the hydrogel, and which we could visualize under the microscope. We decided to use VPGIG<sub>[60]</sub>-mNeonGreen (part <a style="color:#F6B227" href="https://parts.igem.org/Part:BBa_K4905016"> BBa_K4905016</a>), since the hydrophobic ELP part of the protein is able to make hydrophobic interactions with the gel. Because of these interactions, we expected that the proteins would al localize in the same places in the cell. We cloned the gene fragment into a pBad vector under the control of an arabinose promotor, to allow for orthogonality between this vector and the pET24a(+) vector encoding the zipper-ELP construct. Success of cloning was verified by a double digestion of the plasmids as well as sequencing.  
+
Characterization of the ability of our part to form a gell inside the cell with the use of microscopy required us to design a model protein which could interact with the hydrogel, and which we could visualize under the microscope. We decided to use VPGIG<sub>[60]</sub>-mNeonGreen (part <a style="color:#F6B227" href="https://parts.igem.org/Part:BBa_K4905016"> BBa_K4905016</a>), since the hydrophobic ELP part of the protein is able to make hydrophobic interactions with the gel. Because of these interactions, we expected that the proteins would all localize in the same places in the cell. We cloned the gene fragment into a pBad vector under the control of an arabinose promotor, to allow for orthogonality between this vector and the pET24a(+) vector encoding the zipper-ELP construct. Success of cloning was verified by a double digestion of the plasmids as well as sequencing.  
 
</p>
 
</p>
  
Line 297: Line 322:
 
</p><br>
 
</p><br>
  
<img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/mneongreen.png" style="width:25vw;" class="results-img enlarge-image" alt="Image 12"><br>
+
<p align="justify">  
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 13 |</b> Double digestion of pBad vector containing [A3G2]<sub>[60]</sub>-mNeonGreen. The band can be seen at the expected weight.<</figcaption><br>
+
Figure 15 shows two samples before and after photobleaching. Following bleach of <i>E.coli</i> which co-expressed Z1-A120-Z2 and and [A3G2]<sub>12</sub>-mNeonGreen, a gradual recovery of fluorescence could be observed up to 37.5 seconds post-bleach (Fig. 16A, Fig. 17B). The control, containing only [A3G2]<sub>12</sub>-mNeonGreen, showed a virtually instantaneous recovery, as evidenced by Fig 16B, Fig 17C.
 
+
<p align="justify">
+
Z1-A120-Z2 and [A3G2]<sub>[12]</sub>-mNeonGreen were co-expressed in E.coli BL21(DE3). SDS-PAGE was then used to confirm the presence of both proteins. After some trial and error using different conditions for the co-expression, we managed to express both the ELPs, as well as the fluorescent protein.
+
 
</p><br>
 
</p><br>
<head>
+
<div class="image-container">
    <meta charset="UTF-8">
+
<img src="https://static.igem.wiki/teams/4905/wiki/frap-4-timepoints-black-g.png" style="display: block; margin-left: auto; margin-right: auto; width: 90%;"><br>
    <meta name="viewport" content="width=device-width, initial-scale=1.0">
+
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 16:</b><i> A) sample containing Z1-A120-Z2 and [A3G2]<sub>12</sub>-mNeonGreen. B) Control containing only [A3G2]<sub>12</sub>-mNeonGreen. Both samples were first measured at three timepoints (1s each)  pre-bleach, then bleached for 5 seconds and finally measured for another 30 timepoints (1s each). The ROI that was bleached is indicated with the red circles.</i></figcaption><br>
    <title>Two Images Side by Side</title>
+
</div>
    <style>
+
        .image-container {
+
            display: flex;
+
        }
+
 
+
        .image-container img {
+
            width: 50%; /* Adjust this value to control the width of the images */
+
            height: auto; /* Maintain aspect ratio */
+
        }
+
    </style>
+
</head>
+
    <div class="image-container">
+
        <img src="https://static.igem.wiki/teams/4905/wiki/co-expression-mng-z1z2-g.png" alt="Image 9">
+
        <img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/img-2903.jpeg" alt="Image 10">
+
    </div>
+
 
+
<br>
+
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 14 | </b> A) SDS-PAGE gel showing the co-expressed proteins (left three lanes) versus a negative control in which no arabinose was added and therefore no expression of [A3G2]60-mNeonGreen was induced (right three lanes). B) Bacterial cell pellet excited by a blue light transilluminator showing the fluorescence emitted by mNeonGreen.
+
</figcaption><br>
+
  
 
<p align="justify">
 
<p align="justify">
The bacteria containing Z1-A120-Z2 and [A3G2]<sub>[12]</sub>-mNeonGreen were then used to determine if diffusion throughout the cell was affected by the formation of the gel. This was done by using Fluorescence recovery after photobleaching (FRAP). We hypothesized that once the gel is formed intracellularly, the diffusion rate is negatively affected due to an increase of viscosity of the cytosol.
 
Figure 15 shows two samples before and after photobleaching. Following bleach of <i>E.coli</i> which co-expressed Z1-A120-Z2 and and [A3G2]<sub>[12]</sub>-mNeonGreen, a gradual recovery of fluorescence could be observed up to 30 seconds post-bleach (Fig. 15A, Fig. 16A). The control, containing only [A3G2]<sub>[12]</sub>-mNeonGreen, showed a virtually instantaneous recovery, as evidenced by Fig 15B, Fig 16B.
 
 
</p><br>
 
  
<img src="https://static.igem.wiki/teams/4905/wiki/frap-4-timepoints-black.png" style="display: block; margin-left: auto; margin-right: auto; width:50vw;"><br>
+
The raw data was corrected for photobleaching during post-bleach image acquisition and normalized to the intensity pre-bleach. The data was then fitted with equation 1, Where I(t) stands for the normalized fluorescence intensity at time t, A is the mobile fraction, t is time and τ is the time constant<sup>[7]</sup>. This was done for each of the measurements, resulting in four graphs for the co-expressed sample and two graphs for the control. Using the equations, the recovery half-life (τ<sub>1/2</sub>) was calculated, using t<sub>1/2</sub> = ln(2)*τ for each data set which is summarized in Figure 15. This experiment shows that proteins are distributed throughout the whole cell and that the diffusion rate becomes lower when our ELP constructs are expressed. It is impossible to exactly determine the intracellular concentration of ELPs, and thus to make a comparison with the gels made extracellularly. However, it would be worthwhile doing these FRAP measurements on more samples in the future where, just as was done for the CCK8 assays, protein expression is induced for a series of different timepoints to see the effect on diffusion.
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 15 | </b> A) sample containing Z1-A120-Z2 and [A3G2]<sub>[12]</sub>-mNeonGreen. B) Control containing only [A3G2]<sub>[12]</sub>-mNeonGreen. Both samples were first measured at three timepoints (1s each) pre-bleach, then bleached for 5 seconds and finally measured for another 30 timepoints (1s each). The ROI that was bleached is indicated with the red circles.<</figcaption><br>
+
</p>
 +
<div style="text-align: center;">
 +
<p>
 +
I(t) = A(1-e<sup>(-t/τ)</sup>)+c  (1)
 +
</p>
 +
</div>
  
  
<p align="justify">
 
  
The data was fitted with the following equation:
 
y = a∗(1-exp(-t/τ))+c (4)
 
Where A is the mobile fraction, t is time and τ is the time constant.  This was done for each of the measurements, resulting in four graphs for the co-expressed sample and two graphs for the control. Using the equations, the recovery half-life (τ1/2) was calculated for each data set which is summarized in Table 2 below.
 
  
<div class="row">
+
    <div class="image-container">
   <div class="col">
+
<div class="row" style="display: flex;">
     <table>
+
   <div class="col" style="width:30%;">
 +
     <center><table>
 
     <tr>
 
     <tr>
  
         <th>Measurement</th><th>τ1/2 (s)</th>
+
         <th>Measurement</th><th>τ<sub>1/2</sub> (s)</th>
 
     </tr>
 
     </tr>
 
     <tr>   
 
     <tr>   
     <td>1</td><td>11.89</td>  
+
     <td>1</td><td>23.9</td>  
 
     </tr>
 
     </tr>
 
     <tr>
 
     <tr>
     <td>3</td><td>7.08</td>
+
     <td>2</td><td>7.0</td>
 
     </tr>
 
     </tr>
 
     <tr>
 
     <tr>
     <td>4</td><td>10.87</td>  
+
     <td>3</td><td>16.5</td>  
 
     </tr>
 
     </tr>
 
     <tr>
 
     <tr>
     <td>5</td><td>8.43</td>  
+
     <td>4</td><td>11.7</td>  
 
     </tr>
 
     </tr>
 
     <tr>
 
     <tr>
     <td>Average</td><td>9.57+/-1.91</td>  
+
     <td>Average + stdev</td><td>14.8 ± 6.2</td>  
 
     </tr>
 
     </tr>
 
   
 
   
  
 
     </tr>
 
     </tr>
     </table>
+
     </table></center>
 
   </div>
 
   </div>
 
+
        <img style="width:40%;" src="https://static.igem.wiki/teams/4905/wiki/z1z2-final-g.png" alt="Image 7">
</p><br>
+
         <img style="width:40%;" src="https://static.igem.wiki/teams/4905/wiki/mng-final.png" alt="Image 8">
 
+
<head>
+
    <meta charset="UTF-8">
+
    <meta name="viewport" content="width=device-width, initial-scale=1.0">
+
    <title>Two Images Side by Side</title>
+
    <style>
+
        .image-container {
+
            display: flex;
+
        }
+
 
+
        .image-container img {
+
            width: 50%; /* Adjust this value to control the width of the images */
+
            height: auto; /* Maintain aspect ratio */
+
        }
+
    </style>
+
</head>
+
    <div class="image-container">
+
        <img src="https://static.igem.wiki/teams/4905/wiki/frap-z1z2.png" alt="Image 11">
+
         <img src="https://static.igem.wiki/teams/4905/wiki/frap-mng.png" alt="Image 12">
+
 
     </div>
 
     </div>
 +
 
 +
</div><br>
 +
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 17:<i></b> A) Table with the average recovery half-lives. B) Fluorescence recovery of samples induced to express Z1-A120-Z2 and [A3G2]<sub>12</sub>-mNeonGreen, normalized to the pre-bleach intensity. All samples show a gradual recovery post-bleach with a mean recovery half-life of 14.8 ± 6.2 s. C) Fluorescence recovery of bacteria induced to express only [A3G2]<sub>12</sub>-mNeonGreen, normalized to the pre-bleach intensity. The fluorescense recovery was almost instantaneous.</i>
 +
</figcaption>
  
 
<br>
 
<br>
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 16 | </b> Measured intensities from the FRAP measurement. Top: Measurements from the samples containing both Z1-A120-Z2 and [A3G2]<sub>[12]</sub>-mNeonGreen. Bottom: Measurements from samples containing only [A3G2]<sub>[12]</sub>-mNeonGreen. It can be seen that the fluorescence slowly recovers over time for the samples that contain the zipper proteins, and recovery is almost instant in the control group. This means that the visocity in the samples with our ELP-zipper proteins is higher.
+
<p align="justify">
</figcaption>
+
To further characterize the gelated versus non-gelated bacteria, we used a DNA staining, Hoechst (Figure 18). We saw that a large number of gelated bacteria showed a clustering of DNA in the center of the cell, whereas for the control sample, the DNA was more homogenously spread throughout the cell. This might indicate that the DNA in gelated cells is not as accessible to replication machinery and that therefore division of the cells is stopped or slowed down.
 
+
 
</p>
 
</p>
  
 
+
<div class="image-container">
<h3>Freeze-drying</h3>
+
<img src="https://static.igem.wiki/teams/4905/wiki/live-or-dye-black.png" style="display: block; margin-left: auto; margin-right: auto; width: 80%;">
 
+
<p align="justify">
+
For our application, it was important that we made the storage conditions of our bacteria as easy as possible. We came up with he idea of freeze-drying our bacteria to see if they would survive that. Since our ML-I Biolab did not contain a freeze-drying setup, we had to make one ourselves. We did this in one of the fumehoods in the biolab and we attempted to freeze-dry six different types of sample (Table 3). Three for each of the different conditions also tested previously (no expression, A120 and Z1-A120-Z2) and these .same three, but with 100 mM threhalose added. This is a sugar which protects the bacterial membrane during freeze-drying<sup>[10]</sup>, and was therefore meant to be used in the samples for a positive control. However, even after running the setup for over 8h, some samples did not sublimate unfortunately. And since this setup required us to refill the cold trap with liquid nitrogen every couple of hours, we could not leave it overnight. Therefore, only two of the samples, Wt and Z1-A120-Z2, were dried relatively well. These were then reconstituted and subseqeuntly assayed using CCK8. Using the calibration curve made for this assay, the OD<sub>600</sub> values were calculated from the measured OD<sub>450</sub> values. These were then divided by 0.5, the original OD600 at which the samples were made, to obtain the final survival rate of the bacteria.</p><br>
+
 
+
<div class="row">
+
  <div class="col">
+
    <table>
+
    <tr>
+
 
+
        <th>Sample</th><th>Success?</th>
+
    </tr>
+
    <tr> 
+
    <td>No expression</td><td>Yes</td>
+
    </tr>
+
    <tr>
+
    <td>No expression + 100 mM trehalose</td><td>No</td>
+
    </tr>
+
    <tr>
+
    <td>A120</td><td>No</td>
+
    </tr>
+
    <tr>
+
    <td>A120 + 100 mM trehalose</td><td>No</td>
+
    </tr>
+
    <tr>
+
    <td>Z1-A120-Z2</td><td>Yes</td>
+
    </tr>
+
    <tr>
+
    <td>Z1-A120-Z2 + 100 mM trehalose</td><td>No</td>
+
    </tr>
+
 
+
    </tr>
+
    </table>
+
  </div>
+
  <div class="col">
+
    <img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/cck8-calibration-curve.png" style="width:25vw;" class="results-img enlarge-image" alt="Image 7">
+
 
   </div>
 
   </div>
 
</div><br>
 
</div><br>
 +
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 18:</b><i> A) bacteria gelated with Z1-A120-Z2 stained with Hoechst. The DNA is densely clustered in the center of the cells. B) non-gelated cells stained with Hoechst. The DNA is homogenously spread throughout the cell.</i>
 +
</figcaption><br>
  
<figcaption style="text-align:center; font-size: 12px;"><b>igure 17 |</b> Left: Freeze-drying samples that were made and whether they succeeded or not. Right: CCK8 calibration curve, used to calculate the OD<sub>600</sub> from the OD<sub>450</sub>.</figcaption><br>
+
</p>
 
+
<img src="https://static.igem.wiki/teams/4905/wiki/bba-k4905006/cck8-freeze-dry.png" style="width:25vw;" class="results-img enlarge-image" alt="Image 8"><br>
+
<figcaption style="text-align:center; font-size: 12px;"><b>Figure 18 |</b> Survival rate of E.coli after freeze-drying and reconstitution. It can be seen that the survival rate of the samples that contained the Z1-A120-Z2 proteins was higher than the control sample.</figcaption><br>
+
 
+
  
  
Line 450: Line 398:
 
<h2>Molecular Dynamics simulation</h2>  
 
<h2>Molecular Dynamics simulation</h2>  
 
<h3>Introduction</h3>
 
<h3>Introduction</h3>
<p>
+
<p align="justify">
 +
 
 
Molecular dynamics (MD) simulations are commonly used to see how proteins or other molecules move over time. This can contribute to the understanding of molecular systems, which cannot be seen with the eye or even with the best microscopes.
 
Molecular dynamics (MD) simulations are commonly used to see how proteins or other molecules move over time. This can contribute to the understanding of molecular systems, which cannot be seen with the eye or even with the best microscopes.
In our application of forming a hydrogel with Elastin-Like Polypeptides (ELPs), it can be useful to see how the ELPs will behave at different temperatures. Above a certain transition temperature (Tt), there is a disruption in the thermal energy of the water molecules around the ELPs, which causes the water molecules to move away from the ELPs. This makes room for interactions between the ELPs. It has been found that there will be conformational changes in another ELP sequence when the temperature rises above Tt, which makes them more ordered and increases the hydrophobicity<sup>[11]</sup>.  
+
In our application of forming a hydrogel with ELPs, it can be useful to see how the ELPs will behave at different temperatures. Above a certain transition temperature (T<sub>t</sub>), there is a disruption in the thermal energy of the water molecules around the ELPs, which causes the water molecules to move away from the ELPs. This makes room for interactions between the ELPs. It has been found that there will be conformational changes in another ELP sequence when the temperature rises above T<sub>t</sub>, which makes them more ordered and increases the hydrophobicity<sup>[8]</sup>.  
 
</p>
 
</p>
 
<p>
 
<p>
It is interesting to see how this applies to our own ELP. The MD simulation will give a better understanding of the differences in the structure below and above the T<sub>t</sub> and it will help with understanding how the ELP actually looks like. We will analyze different properties of the proteins, like the secondary structure formation, compactness, and the surface area that is accessible for the solvent, to see how the conformation differs at specific temperatures below and above the Tt.  
+
It is interesting to see how this applies to our own ELP. The MD simulation will give a better understanding of the differences in the structure below and above the T<sub>t</sub> and it will help with understanding what the ELP structurally looks like. We have analyzed different properties of the proteins, like the secondary structure formation, compactness, and the surface area that is accessible for the solvent, to see how the conformation differs at specific temperatures below and above the T<sub>t</sub>.  
 
</p>
 
</p>
 
<h3>Methods</h3>
 
<h3>Methods</h3>
<p>
+
<p align="justify">
MD simulations are typically computationally expensive. This made it necessary to shorten our ELP sequence to make the computational time manageable. We chose to shorten the hydrophobic parts to I<sub>[10]</sub> instead of I<sub>[60]</sub> and to shorten the hydrophilic part to A<sub>[50]</sub> instead of A<sub>[120]</sub> or A<sub>[100]</sub>, as shown in Figure 19. The leucine zippers stay the same. In this length reduction, the proportions are not the same as in the original ELP. However, we wanted to make sure that the hydrophobic parts were not shorter than the leucine zippers. This maintains a greater distance between the Leucine zippers and the hydrophilic part, and it prevents unrealistic folding. The hydrophilic part is kept long since this part gives the ELP its length and it might be interesting to see how this part behaves during the folding. The choices made for this might affect the outcome, but the combination of regions is still the same and the difference between the temperatures can still be compared.
+
MD simulations are typically computationally expensive. This made it necessary to shorten our ELP sequence to make the computational time manageable. We chose to shorten the hydrophobic parts to I<sub>[9]</sub> instead of I<sub>[60]</sub> and to shorten the hydrophilic part to A<sub>[50]</sub> instead of A<sub>[120]</sub> or A<sub>[100]</sub>, as shown in Figure 19. The Leucine zippers stay the same. In this length reduction, the proportions are not the same as in the original ELP. However, we wanted to make sure that the hydrophobic parts were not shorter than the leucine zippers. This maintains a greater distance between the Leucine zippers and the hydrophilic part, and it prevents unrealistic folding. The hydrophilic part is kept long, since this part gives the ELP its length and it will provide relevant information on how the protein behaves when it folds. The choice to change the sequence might affect the outcome a bit, however we still expect to find relevant information on how the ELP behaves at different temperatures.  
 
</p>
 
</p>
  
 +
<div class="image-container">
 +
<figure>
 
<img src="https://static.igem.wiki/teams/4905/wiki/partsconstructs/partsconstructs/simulatieseqnieuw.png" style="width:100%">  
 
<img src="https://static.igem.wiki/teams/4905/wiki/partsconstructs/partsconstructs/simulatieseqnieuw.png" style="width:100%">  
 +
 +
<figcaption>
 
<p align="justify">
 
<p align="justify">
<b>Figure 19: </b> <i>ELP construct used for the MD simulation. In this construct, the hydrophobic VPGIG parts are repeated 10 times instead of 60 times and the hydrophilic part in the middle is repeated 10 times instead of 24 or 20 times. At the ends, the two different leucine zippers are located, which remain the same.</i>
+
<b>Figure 19: </b> <i>ELP construct used for the MD simulation. In this construct, the hydrophobic VPGIG parts are repeated 10 times instead of 60 times and the hydrophilic part in the middle is repeated 10 times instead of 24 or 20 times. At the ends, the two different Leucine zippers are located, which remain the same.</i>
</p>  
+
</p>
 +
</figcaption>
 +
</figure>
 +
</div>
  
<p>  
+
<p align="justify">
A generalized born implicit solvent is used to fold the protein. In an implicit solvent, the solvent (in this case water) molecules are not individually modeled, but with an approximation of the mean force from the solvent on the protein. Different theories can be used for this, and we chose to use the generalized born model. Implicit solvents are faster than an explicit solvent with simulated water molecules and more detailed than a simulation in vacuum[12]. The MD simulation is performed at temperatures 3 °C, 15 °C, 25 °C, and 37 °C, since the transition temperature was found to be 20 °C in the conducted UV-VIS experiments. The simulations are run in triplet to reduce the chance of coincidence in the results. This is necessary, since the initial velocities are chosen by the NAMD software and thus differ every time. Different velocities can cause differences in the final stable conformations, which can affect the results of the analytical methods.
+
A generalized born implicit solvent is used to fold the protein. In an implicit solvent, the solvent (in this case water) molecules are not individually modeled, but with an approximation of the mean force from the solvent on the protein. Different theories can be used for this, and we chose to use the generalized born model. Implicit solvents are faster than an explicit solvent with simulated water molecules and more detailed than a simulation in vacuum<sup>[10]</sup>. The MD simulation is performed at temperatures 3 °C, 15 °C, 25 °C, and 37 °C, since the transition temperature was found to be 20 °C in the conducted UV-VIS experiments. The simulations are run in triplet to reduce the chance of coincidence in the results. This is necessary, since the initial velocities are chosen by the NAMD software and thus differ every time. Different velocities can cause differences in the final stable conformations, which can affect the results of the analytical methods.
 
</p>
 
</p>
<p>
+
<p align="justify">
Just like in López Barreiro et al., first energy minimalization is done for 20 000 time steps. Then, Langevin dynamics is used with Generalized Born implicit solvent. The simulation step was 2 fs with a simulation time of 50 ns. The cutoff distance is 18 Å, the switch distance 16 Å<sup>[13]</sup>.
+
Just like in López Barreiro et al., first energy minimalization is done for 20 000 time steps. Then, Langevin dynamics is used with Generalized Born implicit solvent. The simulation step was 2 fs with a simulation time of 50 ns. The cutoff distance is 18 Å, the switch distance 16 Å<sup>[11]</sup>.
 
</p>
 
</p>
  
 
<b>Analysis methods</b>
 
<b>Analysis methods</b>
<p>
+
<p align="justify">
The results are analyzed with Root Mean Square Distance (RMSD), Radius of Gyration and relative shape anisotropy, Ratio of Secondary Structures, and Solvent Accessible Surface Area (SASA) to show the different behaviors of the ELP at different temperatures. They are calculated with the MDTraj Python library. For the analysis, the analysis methods are calculated for the 1000 frames of the last 5 ns, which has a stable conformation at all temperatures. Even when the conformation is stable, there are still some small differences. These differences are taken into account with this distribution of the results of the last 5 ns.
+
The results are analyzed with Root Mean Square Distance (RMSD), Radius of Gyration and relative shape anisotropy, Ratio of Secondary Structures, and Solvent Accessible Surface Area (SASA) to show the different behaviors of the ELP at different temperatures. They are calculated with the MDTraj Python library. For the analysis, methods are calculated for the last 5 ns, which contain 1000 frames. This is chosen based on the RMSD results, where a stable conformation is shown at all temperatures. Even when the conformation is stable, there are still some small conformation differences. These differences are taken into account with this distribution of the results of the last 5 ns.
 
</p>
 
</p>
  
  
 
<h3>Results</h3>
 
<h3>Results</h3>
<p>
+
<p align="justify">
In Figure 20, a video is shown from one of the simulations at 25°C. Not all videos are included, since the process is comparable for the other simulations, especially when comparing them by eye. Some frames of those simulations can be found in the file below. When comparing the simulations at the four temperatures, they all start with folding at the ends of the ELP. This was expected since those parts are the most hydrophobic. Some of the simulations also show the folding of the ELP in the middle at the more hydrophilic part. However, there is no clear relationship between temperature and this folding behavior, so it is probably caused by luck in the simulation.
+
In Figure 20, a video is shown from one of the simulations at 25°C. Not all videos are included, since the process is comparable for the other simulations, especially when comparing them visually. Some frames of those simulations can be found in the file below. When comparing the simulations at the four temperatures, they all start with folding at the ends of the ELP. This was expected since those parts are the most hydrophobic. Some of the simulations also show the folding of the ELP in the middle at the more hydrophilic part. However, there is no clear relationship between temperature and this folding behavior, so it is probably caused by luck in the simulation.
 
</p>
 
</p>
  
<img src="https://static.igem.wiki/teams/4905/wiki/simulatie/video-elp-folding-1.gif" style="width:100%">
+
<div class="image-container">
 +
<figure>
 +
<img src="https://static.igem.wiki/teams/4905/wiki/simulatie/elp-gif-25.gif" style="width:100%">
 +
<figcaption>
 
<p align="justify">
 
<p align="justify">
 
<b>Figure 20: </b> <i>Video with one of the MD simulations at 25°C. The red residues represent the Leucine zipper, the blue residues are the hydrophobic part, and the light-blue residues are the hydrophilic part. </i>
 
<b>Figure 20: </b> <i>Video with one of the MD simulations at 25°C. The red residues represent the Leucine zipper, the blue residues are the hydrophobic part, and the light-blue residues are the hydrophilic part. </i>
 
</p>  
 
</p>  
 +
</figcaption>
 +
</figure>
 +
</div>
  
 
<b>RMSD</b>
 
<b>RMSD</b>
<p>
+
<p align="justify">
RMSD plots are made for each of the three simulations at the four temperatures. They are made for the complete 50 ns, but also for the last 5 ns to show how stable the conformation is related to frame 9000 (45 ns). This is important since there can still be a lot of changes in the structure while the moved distance of the atoms stays the same. Since there are a lot of figures made, they are not visualized here. However, they can be found on the wiki of the TU Eindhoven 2023 team.   
+
RMSD plots are made for each of the three simulations at the four temperatures. They are made for the complete 50 ns, but also for the last 5 ns to show how stable the conformation is related to frame 9000 (45 ns). This is important since there can still be a lot of changes in the structure while the moved distance of the atoms stays the same. Since there are a lot of figures made, they are not visualized here. However, they can be found on the wiki of the TU-Eindhoven 2023 team.   
 
</p>
 
</p>
  
  
<p>
+
<p align="justify">
For all of the measurements, the equilibrium seems to be reached. When looking at the plots of the last 5 ns, most RMSD plots have a variability of a maximum of 1 nm. Only the third simulation at 37 °C has a slightly larger difference. An RMSD of 0 describes a stable conformation, but there are always some small deviations. So it can be said that an RMSD of 0.25 nm has a high stability, but also 0.2-0.5 nm will give a high stability<sup>[14]</sup>. With this information and when looking at the course of the curves, it can still be concluded that the conformations are relatively stable since the RMSD is still very small.  
+
For all of the measurements, the equilibrium seems to be reached. When looking at the plots of the last 5 ns, most RMSD plots have a variability of a maximum of 1 nm. Only the third simulation at 37 °C has a slightly larger difference. An RMSD of 0 describes a stable conformation, but there are always some small deviations. So it can be said that an RMSD of 0.25 nm has a high stability, but also 0.2-0.5 nm will give a high stability<sup>[12]</sup>. With this information and when looking at the course of the plots, it can still be concluded that the conformations are relatively stable since the RMSD is still very small.  
 
</p>
 
</p>
  
 
<b>Secondary structure</b>
 
<b>Secondary structure</b>
<p>
+
<p align="justify">
 
In Figures 21a, b, and c, the boxplots for the secondary structure analysis of the individual simulations are shown. This shows some variability between the simulations of a single temperature, but also overlapping deviations. This happens especially for the α-helix and β-sheet plots. The simulation results are taken together for each temperature to make them easier to analyze, which is shown in Figure 21d, e, and f.  
 
In Figures 21a, b, and c, the boxplots for the secondary structure analysis of the individual simulations are shown. This shows some variability between the simulations of a single temperature, but also overlapping deviations. This happens especially for the α-helix and β-sheet plots. The simulation results are taken together for each temperature to make them easier to analyze, which is shown in Figure 21d, e, and f.  
 
</p>
 
</p>
  
<img src="https://static.igem.wiki/teams/4905/wiki/simulatie/secondary-structure3.png" style="width:60%">  
+
<div class="image-container">
 +
<figure>
 +
<img src="https://static.igem.wiki/teams/4905/wiki/simulatie/secondary-structure3.png" style="width:90%">
 +
<figcaption>
 
<p align="justify">
 
<p align="justify">
 
<b>Figure 21: </b> <i>Secondary structure analysis with <b>a, b, and c</b> for respectively the coils, α-helices, and β-sheets three individual simulations for each of the four temperatures and <b>d, e, and f</b> the overall results with the different simulations per temperature taken together.</i>
 
<b>Figure 21: </b> <i>Secondary structure analysis with <b>a, b, and c</b> for respectively the coils, α-helices, and β-sheets three individual simulations for each of the four temperatures and <b>d, e, and f</b> the overall results with the different simulations per temperature taken together.</i>
 
</p>  
 
</p>  
 +
</figcaption>
 +
</figure>
 +
</div>
  
<p>
+
<p align="justify">
Overall, almost no secondary structures seem to be formed. The rate of the coils is much higher than the rates for alpha-helices and beta-sheets. It is possible that more secondary structures make the ELP more stable, but that it is stuck in a local minimum. This can be proved with more simulations and some other simulation techniques.  
+
Overall, almost no secondary structures seem to be formed. The ratio of the coils is much higher than the ratios for alpha-helices and beta-sheets. It is possible that more secondary structures make the ELP more stable, but that it is stuck in a local minimum. This can be proved with more simulations with for example other force fields used.  
 
</p>
 
</p>
<p>
+
<p align="justify">
 
When comparing the boxplots of the different temperatures, there does not seem to be a convincing difference in the secondary structure formation. There are some small differences in the medians and distributions, but there is no clear relation between those differences. We can thus conclude that there is no difference in the secondary structure between the temperatures.  
 
When comparing the boxplots of the different temperatures, there does not seem to be a convincing difference in the secondary structure formation. There are some small differences in the medians and distributions, but there is no clear relation between those differences. We can thus conclude that there is no difference in the secondary structure between the temperatures.  
 
</p>
 
</p>
  
 
<b>Radius of gyration</b>
 
<b>Radius of gyration</b>
<p>
+
<p align="justify">
 
The radius of gyration results of the individual simulations are shown in Figure 22a. The individual simulation results of 3 and 37°C show large differences in the results. This can be caused by the different velocities that are chosen for the single simulations. Since the individual simulations of 3 and 37°C differ a lot, the total boxplot in Figure 22b shows a large distribution of the results.  
 
The radius of gyration results of the individual simulations are shown in Figure 22a. The individual simulation results of 3 and 37°C show large differences in the results. This can be caused by the different velocities that are chosen for the single simulations. Since the individual simulations of 3 and 37°C differ a lot, the total boxplot in Figure 22b shows a large distribution of the results.  
 
</p>
 
</p>
  
<img src="https://static.igem.wiki/teams/4905/wiki/simulatie/radiusofgyration3.png" style="width:60%">  
+
<div class="image-container">
 +
<figure>
 +
<img src="https://static.igem.wiki/teams/4905/wiki/simulatie/radiusofgyration3.png" style="width:90%">
 +
<figcaption>
 
<p align="justify">
 
<p align="justify">
 
<b>Figure 22: </b> <i>Radius of gyration results with <b>a</b> three individual simulations for each of the four temperatures and <b>b</b> the overall results with the different simulations per temperature taken together.</i>
 
<b>Figure 22: </b> <i>Radius of gyration results with <b>a</b> three individual simulations for each of the four temperatures and <b>b</b> the overall results with the different simulations per temperature taken together.</i>
 
</p>  
 
</p>  
 +
</figcaption></figure>
 +
</div>
 +
 +
<p align="justify">
  
<p>
 
 
There is a clear difference between the radius of gyration at 15 °C and 25 °C, which are just below and just above the transition temperature. When looking at the median at 37 °C, also these results show a lower radius of gyration than 15 °C. At 3 °C, the radius of gyration is for two of the simulations very low compared to the other temperatures. When looking at the final conformations of the simulations, this can be caused by the shape of the stable configurations. When the shape is more cylindrical instead of spherical, which seems especially the case for the simulations at 15°C, the radius of gyration will increase. As mentioned in the methods section, we also computed the relative shape anisotropy.
 
There is a clear difference between the radius of gyration at 15 °C and 25 °C, which are just below and just above the transition temperature. When looking at the median at 37 °C, also these results show a lower radius of gyration than 15 °C. At 3 °C, the radius of gyration is for two of the simulations very low compared to the other temperatures. When looking at the final conformations of the simulations, this can be caused by the shape of the stable configurations. When the shape is more cylindrical instead of spherical, which seems especially the case for the simulations at 15°C, the radius of gyration will increase. As mentioned in the methods section, we also computed the relative shape anisotropy.
 
</p>
 
</p>
  
<img src="https://static.igem.wiki/teams/4905/wiki/simulatie/relativeshape3.png" style="width:60%">  
+
<div class="image-container">
 +
<figure>
 +
<img src="https://static.igem.wiki/teams/4905/wiki/simulatie/relativeshape3.png" style="width:90%">  
 +
<figcaption>
 +
<p align="justify">
 +
<b>Figure 23: </b> <i>Relative shape anisotropy results with <b>a</b> three individual simulations for each of the four temperatures and <b>b</b> the overall results with the different simulations per temperature taken together.</i></p></figcaption></figure></div>
 +
 
 
<p align="justify">
 
<p align="justify">
<b>Figure 23: </b> <i>Relative shape anisotropy results with <b>a</b> three individual simulations for each of the four temperatures and <b>b</b> the overall results with the different simulations per temperature taken together.</i>
 
</p>
 
  
<p>
 
 
In Figure 23a, the relative shape anisotropy is shown for the individual simulations. There is clearly a difference in the shape between the different simulations, also in those with the same temperature. For example 25°C, which has a clear variability between the simulations. Some of the differences were already expected after the analysis of the radius of gyration. The second simulation at 3°C is clearly less spherical than the other two simulations at that temperature. Overall, when looking at Figure 23b, the shape of the ELP at 15°C is less spherical than at 25°C. Because of the large variabilities between the measurements at 3 and 37°C, it cannot be concluded how this is affected at lower and higher temperatures.  
 
In Figure 23a, the relative shape anisotropy is shown for the individual simulations. There is clearly a difference in the shape between the different simulations, also in those with the same temperature. For example 25°C, which has a clear variability between the simulations. Some of the differences were already expected after the analysis of the radius of gyration. The second simulation at 3°C is clearly less spherical than the other two simulations at that temperature. Overall, when looking at Figure 23b, the shape of the ELP at 15°C is less spherical than at 25°C. Because of the large variabilities between the measurements at 3 and 37°C, it cannot be concluded how this is affected at lower and higher temperatures.  
 
</p>
 
</p>
<p>
+
<p align="justify">
This analysis suggests that there is a difference in the compactness of the ELP below and above the transition temperature. Below the transition temperature, the ELP is less compact. At a temperature above the transition temperature, the ELP will become more hydrophobic with fewer solvent molecules around it. This makes the folded ELP more compact, but also more likely to induce interactions with other ELPs. However, for harder conclusions with the role of solvent molecules included, explicit simulations are needed.  
+
This analysis suggests that there is a difference in the compactness of the ELP below and above the T<sub>t</sub>. Below the T<sub>t</sub>, the ELP is less compact. At a temperature above the T<sub>t</sub>, the ELP will become more hydrophobic with fewer solvent molecules around it. This makes the folded ELP more compact, but also more likely to induce interactions with other ELPs. However, for harder conclusions with the role of solvent molecules included, explicit simulations are needed.  
 
</p>
 
</p>
  
 
<b>SASA</b>
 
<b>SASA</b>
<p>
+
<p align="justify">
 
In Figure 24a, the individual results show some deviations between the simulations at specific temperatures, especially at 3 and 37 °C. When taking the individual simulations together for each temperature in Figure 24b, some small deviations are shown. SASA at 25 °C seems to be lower than SASA at 15 °C, while the median SASA at 3 °C is also lower than the rest of the temperatures, but with a large deviation. The results for 37 °C also have a large deviation, with a median higher than the median of 15 °C. When comparing the results with the relative shape anisotropy, there is a relationship between more spherical simulation results and a lower SASA. There seems to be a smaller difference than the difference found in the radius of gyration and the anisotropy, which might be because ends that stick out probably have more effect on SASA than the relative shape anisotropy.  
 
In Figure 24a, the individual results show some deviations between the simulations at specific temperatures, especially at 3 and 37 °C. When taking the individual simulations together for each temperature in Figure 24b, some small deviations are shown. SASA at 25 °C seems to be lower than SASA at 15 °C, while the median SASA at 3 °C is also lower than the rest of the temperatures, but with a large deviation. The results for 37 °C also have a large deviation, with a median higher than the median of 15 °C. When comparing the results with the relative shape anisotropy, there is a relationship between more spherical simulation results and a lower SASA. There seems to be a smaller difference than the difference found in the radius of gyration and the anisotropy, which might be because ends that stick out probably have more effect on SASA than the relative shape anisotropy.  
 
</p>
 
</p>
  
<img src="https://static.igem.wiki/teams/4905/wiki/simulatie/sasa3.png" style="width:60%>  
+
<div class="image-container">
 +
<figure>
 +
<img src="https://static.igem.wiki/teams/4905/wiki/simulatie/sasa3.png" style="width:90%">
 +
<figcaption>
 
<p align="justify">
 
<p align="justify">
<b>Figure 24: </b> <i>SASA results with <b>a</b> three individual simulations for each of the four temperatures and <b>b</b> the overall results with the different simulations per temperature taken together. </sup>.</i>
+
<b>Figure 24: </b> <i>SASA results with <b>a</b> three individual simulations for each of the four temperatures and <b>b</b> the overall results with the different simulations per temperature taken together.</i>
 
</p>  
 
</p>  
 
+
</figcaption></figure>
<p>
+
</div>
The conclusion will be that there is a small but significant difference in SASA at 15 and 25 °C, but that SASA at 37°C becomes higher again. So there is a small effect of the transition temperature on SASA, but not as clear as the radius of gyration. It was expected that SASA would become much lower when passing the transition temperature<sup>[11]</sup>, but in our ELPs it seems more subtle. This might thus be caused by the ends of the ELP that stick out of the spherical shape of the folded ELP.
+
</p>
+
  
 
<h3>Conclusion</h3>
 
<h3>Conclusion</h3>
<p>
+
<p align="justify">
In the analysis of the MD simulation of a shortened version of our ELP, different analysis methods were applied. These bring us to the conclusion that the shape of the folded ELP will be altered when increasing the temperature. At a higher temperature, the results suggest a more compact ELP with a more spherical shape. When the temperature increases even more to 37°C, the difference compared to 15°C will decrease. A temperature of 3°C generally results in more compact conformations but with a large variability between the simulations. The explanation for the more compact ELP at higher temperatures is that the residues become more hydrophobic and thus have more interactions with each other. When multiple ELPs are brought together, it is expected that this property results in increased interactions between the ELPs, especially between the hydrophobic ends and between the Leucine zippers.
+
The MD simulation of the folding of a shortened version of our ELP taught us that an increase in temperature above the T<sub>t</sub> results in a more compact and more spherical ELP due to increased hydrophobicity. This causes more and stronger interactions between the ELPs, which is an important driver of the hydrogel formation<sup>[8]</sup>.  
 
</p>
 
</p>
<p>
+
<p align="justify">
The radius of gyration and relative shape anisotropy shows the most convincing difference between temperatures below and above the transition temperature, while the secondary structure analysis did not find a clear difference. For SASA, there seems to be a difference, but it is not as clear as the radius of gyration and relative shape anisotropy. For more convincing results, it would probably be better to use explicit solvent simulations and to perform more simulations at more temperatures. However, this and other more extensive experiments on the ELP interactions will be done in further research.
+
The greatest difference in results is found for the radius of gyration and relative shape anisotropy between 15°C and 25°C. SASA also shows a difference, but this is less convincing. There is no difference found in the secondary structure between the different temperatures. When comparing the triplicate simulations, there are large differences for especially 3°C but also 37°C. This makes these results less reliable, but there still seems to be a difference in properties between 15 and 37°C. In further research, explicit solvent simulations can be used with more simulated temperatures. This would give more details for better and possibly more convincing results. Also, more extensive research can be done on the interactions between the ELPs at different temperatures and conditions.  
 
</p>
 
</p>
 +
  
 
<h1>References</h1>
 
<h1>References</h1>
Line 579: Line 558:
 
</p>
 
</p>
 
<p>
 
<p>
[5] Christensen, T., Hassouneh, W., Trabbic-Carlson, K., & Chilkoti, A. (2023). Predicting Transition Temperatures of Elastin-Like Polypeptide Fusion Proteins. https://doi.org/10.1021/bm400167h
+
[5] X. Yang, Y. Zhong, D. Wang, and Z. Lu, “A simple colorimetric method for viable bacteria detection based on cell counting Kit-8,” Anal. Methods, vol. 13, no. 43, pp. 5211–5215, Nov. 2021, doi: 10.1039/D1AY01624E
 
</p>
 
</p>
 
<p>
 
<p>
[6] Yang, X., Zhong, Y., Dong, W., & Lu, Z. (2021). A simple colorimetric method for viable bacteria detection based on cell counting Kit-8. Analytical Methods, 13(43), 5211–5215. https://doi.org/10.1039/d1ay01624e
+
[6] E. L. Foltz, J. W. West, I. H. Breslow, and H. Wallick, “Clinical pharmacology of pivampicillin.,” Antimicrob. Agents Chemother., vol. 10, no. 3, pp. 442–454, 1970, doi: 10.47363/jprsr/2022(3)129.
 
</p>
 
</p>
<p>
+
[7] Poudyal, R. R., Guth-Metzler, R., Veenis, A. J., Frankel, E. A., Keating, C. D., & Bevilacqua, P. C. (2019). Template-directed RNA polymerization and enhanced ribozyme catalysis inside membraneless compartments formed by coacervates. Nature Communications, 10(1). https://doi.org/10.1038/s41467-019-08353-4
[7] Bereda, G. (2022) Clinical Pharmacology of Ampicillin. Journal of Pharmaceutical Research & Reports Volume 3(3): 3-3
+
</p>
+
<p>
+
[8] Israeli, E., Shaffer, B. T., & Lighthart, B. (1993). Protection of Freeze-Dried Escherichia coli by trehalose upon exposure to environmental conditions. Cryobiology, 30(5), 519–523. https://doi.org/10.1006/cryo.1993.1052
+
</p>
+
<p>
+
[9] Eils, R., & Kappel, C. (2004). Fluorescence Recovery after Photobleaching with the Leica TCS SP2. Biology. https://www.leica-microsystems.com/science-lab/fluorescence-recovery-after-photobleaching-with-the-leica-tcs-sp2/
+
 
</p>
 
</p>
 
<p>
 
<p>
[10] Pille, J., Van Lith, S. A. M., Van Hest, J. C. M., & Leenders, W. P. J. (2017). Self-Assembling VHH-Elastin-Like peptides for photodynamic nanomedicine. Biomacromolecules, 18(4), 1302–1310. https://doi.org/10.1021/acs.biomac.7b00064
+
[8] N. K. Li, F. G. Quiroz, C. K. Hall, A. Chilkoti, and Y. G. Yingling, “Molecular description of the lcst behavior of an elastin-like polypeptide,” Biomacromolecules, vol. 15, no. 10, pp. 3522–3530, Oct. 2014, doi: 10.1021/BM500658W/SUPPL_FILE/BM500658W_SI_001.PDF.
 
</p>
 
</p>
 
<p>
 
<p>
[11] N. K. Li, F. G. Quiroz, C. K. Hall, A. Chilkoti, and Y. G. Yingling, “Molecular description of the lcst behavior of an elastin-like polypeptide,” Biomacromolecules, vol. 15, no. 10, pp. 3522–3530, Oct. 2014, doi: 10.1021/BM500658W/SUPPL_FILE/BM500658W_SI_001.PDF.
+
[9] J. Kleinjung and F. Fraternali, “Design and application of implicit solvent models in biomolecular simulations,” Curr Opin Struct Biol, vol. 25, no. 100, p. 126, 2014, doi: 10.1016/J.SBI.2014.04.003.
 
</p>
 
</p>
 
<p>
 
<p>
[12] J. Kleinjung and F. Fraternali, “Design and application of implicit solvent models in biomolecular simulations,” Curr Opin Struct Biol, vol. 25, no. 100, p. 126, 2014, doi: 10.1016/J.SBI.2014.04.003.
+
[10] J. Iwata and T. Ando, “Molecular Dynamics Study on Behavior of Resist Molecules in UV-Nanoimprint Lithography Filling Process,” Nanomaterials, vol. 12, no. 15, Aug. 2022, doi: 10.3390/NANO12152554/S1.
 
</p>
 
</p>
 
<p>
 
<p>
[13] D. López Barreiro, A. Folch-Fortuny, I. Muntz, J. C. Thies, C. M. J. Sagt, and G. H. Koenderink, “Sequence Control of the Self-Assembly of Elastin-Like Polypeptides into Hydrogels with Bespoke Viscoelastic and Structural Properties,” Biomacromolecules, vol. 24, no. 1, pp. 489–501, Jan. 2023, doi: 10.1021/ACS.BIOMAC.2C01405/ASSET/IMAGES/LARGE/BM2C01405_0003.JPEG.
+
[11] D. López Barreiro, A. Folch-Fortuny, I. Muntz, J. C. Thies, C. M. J. Sagt, and G. H. Koenderink, “Sequence Control of the Self-Assembly of Elastin-Like Polypeptides into Hydrogels with Bespoke Viscoelastic and Structural Properties,” Biomacromolecules, vol. 24, no. 1, pp. 489–501, Jan. 2023, doi: 10.1021/ACS.BIOMAC.2C01405/ASSET/IMAGES/LARGE/BM2C01405_0003.JPEG.
 
</p>
 
</p>
 
<p>
 
<p>
[14] M. Arnittali, A. N. Rissanou, and V. Harmandaris, “ScienceDirect Structure Of Biomolecules Through Molecular Dynamics Simulations,” 2019, doi: 10.1016/j.procs.2019.08.181.
+
[12] M. Arnittali, A. N. Rissanou, and V. Harmandaris, “ScienceDirect Structure Of Biomolecules Through Molecular Dynamics Simulations,” 2019, doi: 10.1016/j.procs.2019.08.181.
 
</p>
 
</p>
  

Latest revision as of 13:09, 12 October 2023


Elastin-Like Polypeptide Triblock with Leucine Zippers


Sequence and Features

Molecular weight: 105.1 kDa

Codon optimized for E.coli BL21 cells.

Assembly Compatibility:
  • 10
    INCOMPATIBLE WITH RFC[10]
    Illegal EcoRI site found at 2023
    Illegal EcoRI site found at 3949
    Illegal XbaI site found at 140
    Illegal XbaI site found at 2066
  • 12
    INCOMPATIBLE WITH RFC[12]
    Illegal EcoRI site found at 2023
    Illegal EcoRI site found at 3949
    Illegal NheI site found at 4077
  • 21
    INCOMPATIBLE WITH RFC[21]
    Illegal EcoRI site found at 2023
    Illegal EcoRI site found at 3949
    Illegal XhoI site found at 2040
    Illegal XhoI site found at 3966
  • 23
    INCOMPATIBLE WITH RFC[23]
    Illegal EcoRI site found at 2023
    Illegal EcoRI site found at 3949
    Illegal XbaI site found at 140
    Illegal XbaI site found at 2066
  • 25
    INCOMPATIBLE WITH RFC[25]
    Illegal EcoRI site found at 2023
    Illegal EcoRI site found at 3949
    Illegal XbaI site found at 140
    Illegal XbaI site found at 2066
    Illegal NgoMIV site found at 197
    Illegal NgoMIV site found at 377
    Illegal NgoMIV site found at 467
    Illegal NgoMIV site found at 647
    Illegal NgoMIV site found at 2123
    Illegal NgoMIV site found at 2303
    Illegal NgoMIV site found at 2393
  • 1000
    COMPATIBLE WITH RFC[1000]

Information

This part is made up of the basic parts: Leucine zipper Z1 (BBa_K4905004), Leucine zipper Z2 (BBa_K4905005), and two times Elastin-Like Polypeptide (ELP) sequence A[60]I[60] (BBa_K4905001]). This results in the sequence Z1-I[60]-A[120]-I[60]-Z2. With A[5] the sequence (VPGAG[3]VPGGG[2]).

Figure 1: Video with one of the MD simulations at 25°C. The red residues represent the Leucine zipper, the blue residues are the hydrophobic part, and the light-blue residues are the hydrophilic part.

This construct was used by the TU Eindhoven 2023 team to form a hydrogel extracellular as well as intracellular in E.coli BL21 cells. The goal of this was to stop bacteria from dividing, however many more applications can be thought of for this part. Inspiration for this part was taken from Gradišar et al.[1] where they used Leucine zippers to form heterodimers. We wanted to use these zipper domains in a similar way that Fernández‐Colino et al.[2] used Leucine zippers together with ELPs to form a reversible, injectable hydrogel. Furthermore, we noticed that Contreras-Llano et al.[3] were able to form a hydrogel inside of cells, however not with a protein based approach. We wanted to combine these ideas and form a protein based-hydrogel inside of E.coli and use these bacteria for therapeutic applications. However, with the hydrogel itself, many more applications can be thought of that we wanted to explore as much as possible so that other teams can make use of it in the future too. Below you can read more about the experiments we did to characterize this part and about the molecular dynamics simulations that we did.

A schematic overview of what the part looks like is shown in figure 2. It is simple to ligate other pieces of sequences together to make different constructs, with either different lengths of ELPs, different amino acids or other gBlocks at their end.

Figure 2 | Schematic overview of the sequence of this construct. (VPGAG[3]VPGGG[2]) is from now on referred to as A[5] and VPGIG will be referred to as I. This created the final name Z1-I[60]A[120]I[60]-Z2. This name will besimplified to Z1-A120-Z2.

General applications

ELPs are protein polymers derived from human tropoelastin. One of their key features is that they exhibit a phase separation that is often reversible whereby samples remain soluble below their transition temperature (Tt) but form coacervates above Tt[4]. They have many possible applications in purification, sensing, activation, and nano assembly. Furthermore, they are non-immunogenic, substrates for proteolytic biodegradation, and can be decorated with pharmacologically active peptides, proteins, and small molecules. Recombinant synthesis additionally allows precise control over ELP architecture and molecular weight, resulting in protein polymers with uniform physicochemical properties suited to the design of multifunctional biologics. As such, ELPs have been employed for various uses including as anti-cancer agents, ocular drug delivery vehicles, and protein trafficking modulators[4].

Currently ELPs are also being explored in the formation of hydrogels and tissue engineering applications. We are the first iGEM team to use ELPs for the purpose of formation of a hydrogel and have made many different parts while investigating the possibilities. This part turned out to be the most sucessfull. The reversible behaviour of ELPs is even suitable for injectable hydrogels for example[2]. We think this is a part that can have many applications inside as well as outside of cells and we hope that many other iGEM teams will make use of it in the future.

Construct design

The construct consists of ELPs and two different Leucine zippers that have affinity for each other. In general, ELPs have hydrophilic and hydrophobic domains that exhibit reversible phase separation behavior that is temperature-dependent. They are made from a repeating VPGXG sequence, with X replaced by specific amino acids. This results in specific properties of the ELPs, especially related to the Tt at which the ELPs will interact with each other on the hydrophobic sites[2]. When the temperature is below Tt, the water molecules surrounding the hydrophobic parts will go into the bulk water phase which increases the solvent entropy. This makes it possible to form interactions with other ELP molecules[3].

As shown in figure 3, this construct has a hydrophilic region in the middle (A[120]) and a hydrophobic region on each side of it (I[60]). At the ends the Leucine zippers Z1 and Z2 are located for stronger interactions between the ELPs. Leucine zippers consist of a repeating unit that forms an alpha helix. Two Leucine zippers together form ion pairs between the helices, which causes association[4]. These stronger and reversible interactions make them useful in the formation of a hydrogel at a specific Tt. In the end, the hydrogel is formed with electrostatic and hydrophobic interactions between the ELPs.

Figure 3 | Schematic representation of the composite part, an ELP with Leucine zippers on the ends. It has a triblock ELP sequence that can form a thermoresponsive hydrogel.

The goal of our project specifically was to use this part to stop the cell from dividing and to co-express a therapeutic protein (Interleukin 10).

Another part that is often mentioned on this page is BBa_K4905008. This part was often used as a control. The ELP domains of the protein is the same as for the part on this page, however, part BBa_K4905008 has two Leucine zippers that are not complementary, which means that the interactions between the proteins will be weaker and this way, it could be studied what the effect of only the ELP triblocks was. This control part will be reffered to as Z2-A120-Z2 on the rest of this page. Lastly, another protein that can be seen a lot on this page is abbreviated as A120, which is the same sequence as the part on this page, but without the Leucine zippers. So it is part BBa_K4905001 ligated to itself. This protein also serves as a control group in many experiments.

Results

Protein expression and purification

This part was expressed in BL21 E.coli. The full plasmid that the DNA was cloned into can be seen below in figure 4. To check its expression we ran an SDS-PAGE gel of which the results can be seen in figure 5. The bands on the gel showed up where we expected, so this means that protein expression was succesfull.

Figure 4 | Full plasmid containing the sequence for this part in a pET24(+) vector that was transformed into E.coli BL21 cells for protein expression.

Figure 5 | SDS page gel run with this part. Samples were taken after 12, 20 and 36 hours of inducing protein expression. The band can be seen at the expected weight of 105.1 kDa.

Characterization

Hydrogel formation

First, we wanted to find out if we could use our part to form a hydrogel. All gels were formed by resuspending the ELPs in MilliQ (MQ) at different weight/volume%(w/v%). The Leucine zipper containing ELPs were expected to instantaneously form a gel upon increasing the temperature above Tt.

We saw that for both Z1-A120-Z2 and Z2-A120-Z2, a viscous fluid formed upon dissolving the ELPs in MQ water at 4°C. Subsequently, after warming the solution to room temperature, the liquid became very turbid and a gel started to form. Images of the gels can be seen in figure 6.


Image 1 Image 2

Figure 6 | Hydrogel solutions at 5 and 10 (w/v%). It can be seen that their solubility changes when the transition temperature is surpassed. Above this temperature a hydrogel is formed and the substance in the tube becomes cloudy and gel-like (w/v%).

The construct containing Z1-A120-Z2 formed a gel at 5% and 10% (w/v%) as described above. Upon cooling the gel down to 4°C again, the gel remained intact. As a control, Z2-A120-Z2 ( BBa_K4905008) was also dissolved to form solutions of 5% and 10% (w/v%). This construct too, was able to form a gel instantly upon heating to RT. However, in contrast to the Leucine zipper-ELPs, the gel disassembled again when the solution was cooled to 4°C. This indicates that the latter of the two constructs shows reversibility, which is expected when the ELPs are brought below their Tt, while the former construct does not show this reversibility.

Transition temperature determination with UV-Vis

The transition temperatures (Tts) Z1-A[120]-Z2 and control Z2-A[120]-Z2 were determined with the use of UV-visible spectroscopy. Solutions of 5 and 20 μM were made in MQ water. Absorbances were measured at 350 nm. Tts were determined to be 23.5°C and 21°C for both constructs, respectively (Figure 7).

When the control protein Z2-A[120]-Z2 crosses the Tt at 21°C, aggregates form leading to a sharp increase in absorbance. The absorbance then drops, which might be explained by the non-complementary two zipper domains repelling each other, leading to disassembly of the aggregates. Upon further increasing the temperature, the hydrophobic interactions become stronger, which further increases the absorbance, until finally even the hydrophilic block becomes insoluble above 50°C.

Z1-A[120]-Z2 shows a different temperature dependent behavior compared to Z2-A[120]-Z2, in that it shows a smaller increase in absorbance when it crosses Tt of 23.5°C. We think the interaction between Z1 and Z2 stabilizes the aggregates, forming into a branched network of ELPs, even at low concentrations of 5 μM.

Image 1 Image 2

Figure 7 | UV-Vis spectra for constructs Z1-A120-Z2 and Z2-A120-Z2 measured at 350 nm. The solution was heated at a rate of 0.5 degrees Celcius per minute. It can be seen that the transition temperature lays at 23.5 and 21 degrees Celcius for the two proteins respectively.

Absorbances were also measured upon cooling the solution back down. Once the solution goes below Tt, the absorbances of returned back to base levels, indicating that the ELPs can be reversibly precipitated and resuspended again (Figure 8) at these low concentrations.


Image 3 Image 4

Figure 8 | UV-Vis spectra for constructs Z1-A120-Z2 and Z2-A120-Z2 measured at 350 nm. The solution was cooled at a rate of 2 degrees celcius per minute. It can be seen that the transition temperature lays around 23.5 and 21 degrees Celcius for the two proteins respectively.

Dye release from hydrogel

To determine the relative ability for the diffusion of small molecules through our hydrogels, we made gels at two different weight-to-volume percentages containing the fluorescent dye rhodamine B. Finding out whether small molecules can gradually diffuse out of the hydrogel might be interesting for future applications, a drug could by incorporated for instance.

First, a calibration curve was made to determine at which concentration rhodamine B should be added to the gel. Figure 9 shows a linear relationship up to 26 μM of Rhodamine B. We decided to make the gels with a Rhodamine B concentration of 104 μM (0.05 mg/mL), since the amount of dye diffusing out of the gels was expected to be far lower than this amount. A 5% gel as well as a 10% gel were tested. Our hypothesis was that diffusion through the 5% gel would be faster than when a higher weight percentage of the ELPs was used,



Image 3 Image 4
Figure 9: Left: Calibration curve for Rhodamine B, a linear relationship can be seen up to 26 μM of Rhodamine B, with equation: y = 2*106x + 963.16 (R² = 0.993) Right: The hydrogels made with Rhodamine B at a concentration of 104 μM.

The gels were made at a volume of 100 μL. They were then brought to RT and washed with 1.5 mL of warm MQ. Then, 1 mL of warm MQ was pipetted on top of the gels, and they were incubated for up to 8 hours at RT, with samples of 100 μL being taken every 2 hours. After measuring the fluorescence intensity, the values were corrected for the reduction in volume upon taking the samples. From figure 10, it becomes clear that the dye diffuses faster from the 5% gel than from the 10% gel. However, more experiments are needed to truly find the rate of diffusion, it is hard to draw a conclusion with limited data points. The faster diffusion rate of the 5% gel is in line with our hypothesis that a lower weight percentage gel is able to release small molecules faster than a 10% gel.


Image 3 Image 4
Figure 10 | Fluorescence intensity measured from the supernatant on top of the gel over a period of time.

CCK8

We wanted to find a way that we could measure the amount of cells that we can stop from dividing with our hydrogel and whether they would stay alive. We came across Cell counting kit 8 (CCK8). This is an assay that is commonly used to determine the number of alive cells in a sample. The kit uses a highly water-soluble tetrazolium salt, WST-8, which produces a formazan dye upon reduction in the presence of an electron mediator. The amount of formazan generated by dehydrogenases is directly in proportion to the number of living cells. Yang et al. (2021) have demonstrated that this assay is also suitable for the detection of living E. coli[5].

We used this assay to test whether our hydrogelated cells would be more resistant to several conditions, including incubation with ampicillin and freeze-drying. Our biggest goal was to find out if cells would stop dividing, but stay functional with our hydrogel. To test this, we came up with the idea of combining the CCK8 assay described above with the antibiotic ampicillin. Ampicillin prevents the synthesis of the bacterium’s cell wall. It does this by binding to enzymes necessary for the formation of the cell wall [6]. This means it only kills bacteria when they divide. Our hypothesis is that if our hydrogel prevents the bacteria from dividing, they should be affected less by the presence of this antibiotic, and therefore show a higher survival rate than bacteria that do not contain the hydrogel since these will still divide and be killed by the antibiotic. A schematic representation of the experiment can be seen below in figure 11.


Figure 11: A schematic representation of the CCK8 assay protocol. Samples were taken and incubated with or without ampicillin. After 1 hour, the cell pellets were collected and resuspended in PBS. Then the samples were incubated again with the substrate WST-8. The amount of produced fluorescent product that is measured via the OD450 with a plate reader is directly proportional to the amount of cells in the sample.

For this assay, it was nessecary to measure the OD450 with a plate reader. This can directly be related to the OD600[5]. For this relationship, the calibration curve in figure 12 was made.

Image 7
Figure 12: Calibration curve relating the OD600 to the OD450.

Ampicillin resistance

Three different conditions were tested over a period of time. Namely; no (protein) expression, ELPs with no cross linkers (A120), and ELPs with Leucine zippers (Z1-A120-Z2). Samples were taken at 12h, 20h, and 36h of protein expression. Before incubation ampicillin, samples were diluted to at a baseline of OD600 = 0.5. They were then either incubated with ampicillin (1 mg/mL) or without for 1h at 37°C. The results are summarized in Table 1 and Figure 13.

WST-8 conversion Amp/No ampNo expressionA120 Z1-A120-Z2
12h0.71 ± 0.05 1.02 ± 1.00 3.05 ± 0.33
20h0.60 ± 0.215.81 ± 0.22 20.85 ± 4.14
36h0.28 ± 0.086.79 ± 0.17 17.05 ± 2.13

Table 1: Measured WST-8 conversion for samples with no expression, A120 and Z1-A120-Z2 at different timepoints.

Ampicillin-treated E. coli that were not induced to express proteins, and which were treated with ampicillin, showed at most a WST-8 conversion of 0.71% ± 0.05% compared to cells which were not treated, as measured by absorbance at 450 nm. In contrast, E. coli where expression of A120 was induced, showed a maximum conversion of 6.80% ± 0.17% after 36h of protein expression. This is likely due to the fact that this construct is also able to form a gel, although less strong compared to the Z1-A120-Z2 construct, leading to an attenuation in cell division. Finally, Z1-A120-Z2 was expressed for the three-time points. An increase was observed in the conversion of substrate, most notably after 20h of expression (20.86% ± 4.14%). We expected to see a higher percentage after 36h of expression, since at this time point, more bacteria should be gelated and thus unable to divide. The standard deviation of the samples is quite high, with this small number of measurements and large deviation, it is hard to tell whether there is really a difference between these two time points. Overall, a larger sample size could help to reduce the variability of this assay in future tests. Other explanations could be that some of the cells were over-gelated, which might have led to a premature death unrelated to the incubation with ampicillin. Additionally, cell death could also have occurred due to nutrients in the culturing medium running out.

For comparison, OD600 of all samples was also measured after incubation with and without ampicillin. When the samples were incubated with ampicillin, they all decreased in OD600. The samples incubated without ampicillin in which ELPs were expressed show a relatively stable OD600 compared to the sample in which no expression was induced. This was the only sample in which a substantial increase in OD600 could be observed.


Image 7 Image 8

Figure 13 | Left) WST-8 conversion of bacteria treated with ampicillin versus untreated bacteria. right) corresponding OD600 values. Measurements were done at different timepoints (12h, 20h, and 36h). It can be seen that when the samples were incubated with ampicillin, they all decreased in OD600. The samples incubated without ampicillin in which ELPs were expressed show a relatively stable OD600 compared to the sample in which no expression was induced.

The OD600 was also measured without addition of ampicillin over time for each of the three samples (No expression, A120 and Z1-A120-Z2). After 20h of expression, these were diluted to OD = 0.1 and then measured every 30 minutes. From figure 14, it becomes clear that when no expression of ELPs was induced, the bacteria divided more rapidly than when expression was induced. These findings are in line with the results of the ampicillin treatment. It seems that the presence of the ELPs does not altogether inhibit division, but it does slow it down, which seems to be the reason why they are killed slower upon the addition of ampicillin. However, it is hard to draw a conclusion as to whether the attenuation of bacterial division is driven by the formation of a gel inside the cytosol, or simply because of the crowded environment inside the cell, caused by the presence of a high amount of protein. Therefore, additional controls need to be added in future experiments, like samples in which a mono-or diblock ELP is expressed, in order to account for molecular crowding. Furthermore, a sample where a completely different protein is also interesting, to see what the effect of protein expression is on cell division.


Image 6
Figure 14 | measurements of three samples over a period of time. Protein expression in two of the three samples was induced for 20h, after which the samples were diluted to OD = 0.1 and measured every 30 min. Expression of gel-forming ELPs greatly reduces the rate of replication compared to cells without expression.

These results were compared to the classic colony counting method. From each sample, a dilution of 105 was made and then plated out on an LB-agar plate. The plates were incubated at 37°C overnight. Images of the plates were taken (Figure 15) and the number of colonies was counted using the colony counter tool from ImageJ. The results of which are summarized in Figure 15. As expected, the wild-type E. coli treated with ampicillin did not form any colonies, whereas untreated bacteria formed several hundreds of colonies. Samples containing A120 and Z1-A120-Z2 both showed a significantly smaller amount of colonies.



Figure 15 | Photos of all LB-agar plates after incubation at 37°C overnight, including a a graph and table in which the amount of colonies is summarized. It can be seen that the bacteria which do not contain the ELPs form a greater amount of colonies than when expression of these proteins was induced.

No expression A120 Z1-A120-Z2
Amp No Amp Amp No Amp Amp No Amp
12h0 614 0 4 2 16
20h0 597 0 166 21 242
36h0 1610 35 169 142 297

Table 2 | Table with the amount of colonies on each plate

Microscopy (FRAP)

Characterization of the ability of our part to form a gell inside the cell with the use of microscopy required us to design a model protein which could interact with the hydrogel, and which we could visualize under the microscope. We decided to use VPGIG[60]-mNeonGreen (part BBa_K4905016), since the hydrophobic ELP part of the protein is able to make hydrophobic interactions with the gel. Because of these interactions, we expected that the proteins would all localize in the same places in the cell. We cloned the gene fragment into a pBad vector under the control of an arabinose promotor, to allow for orthogonality between this vector and the pET24a(+) vector encoding the zipper-ELP construct. Success of cloning was verified by a double digestion of the plasmids as well as sequencing.

However, to our surprise, the sequencing results indicated that instead of a VGPIG[60] fragment, the mNeonGreen was attached to a VPGAG[3]VPGGG[2][12] fragment, which is hydrophilic, relative to the VPGIG[60] fragment. This is likely due to a labeling mistake by a previous user of the plasmid. Due to a lack of time, we decided to continue using this construct in our microscopy experiments as a model protein to monitor the diffusion of proteins within the gelated bacteria.


Figure 15 shows two samples before and after photobleaching. Following bleach of E.coli which co-expressed Z1-A120-Z2 and and [A3G2]12-mNeonGreen, a gradual recovery of fluorescence could be observed up to 37.5 seconds post-bleach (Fig. 16A, Fig. 17B). The control, containing only [A3G2]12-mNeonGreen, showed a virtually instantaneous recovery, as evidenced by Fig 16B, Fig 17C.



Figure 16: A) sample containing Z1-A120-Z2 and [A3G2]12-mNeonGreen. B) Control containing only [A3G2]12-mNeonGreen. Both samples were first measured at three timepoints (1s each) pre-bleach, then bleached for 5 seconds and finally measured for another 30 timepoints (1s each). The ROI that was bleached is indicated with the red circles.

The raw data was corrected for photobleaching during post-bleach image acquisition and normalized to the intensity pre-bleach. The data was then fitted with equation 1, Where I(t) stands for the normalized fluorescence intensity at time t, A is the mobile fraction, t is time and τ is the time constant[7]. This was done for each of the measurements, resulting in four graphs for the co-expressed sample and two graphs for the control. Using the equations, the recovery half-life (τ1/2) was calculated, using t1/2 = ln(2)*τ for each data set which is summarized in Figure 15. This experiment shows that proteins are distributed throughout the whole cell and that the diffusion rate becomes lower when our ELP constructs are expressed. It is impossible to exactly determine the intracellular concentration of ELPs, and thus to make a comparison with the gels made extracellularly. However, it would be worthwhile doing these FRAP measurements on more samples in the future where, just as was done for the CCK8 assays, protein expression is induced for a series of different timepoints to see the effect on diffusion.

I(t) = A(1-e(-t/τ))+c (1)

Measurementτ1/2 (s)
123.9
27.0
316.5
411.7
Average + stdev14.8 ± 6.2
Image 7 Image 8

Figure 17: A) Table with the average recovery half-lives. B) Fluorescence recovery of samples induced to express Z1-A120-Z2 and [A3G2]12-mNeonGreen, normalized to the pre-bleach intensity. All samples show a gradual recovery post-bleach with a mean recovery half-life of 14.8 ± 6.2 s. C) Fluorescence recovery of bacteria induced to express only [A3G2]12-mNeonGreen, normalized to the pre-bleach intensity. The fluorescense recovery was almost instantaneous.

To further characterize the gelated versus non-gelated bacteria, we used a DNA staining, Hoechst (Figure 18). We saw that a large number of gelated bacteria showed a clustering of DNA in the center of the cell, whereas for the control sample, the DNA was more homogenously spread throughout the cell. This might indicate that the DNA in gelated cells is not as accessible to replication machinery and that therefore division of the cells is stopped or slowed down.


Figure 18: A) bacteria gelated with Z1-A120-Z2 stained with Hoechst. The DNA is densely clustered in the center of the cells. B) non-gelated cells stained with Hoechst. The DNA is homogenously spread throughout the cell.

Molecular Dynamics simulation

Introduction

Molecular dynamics (MD) simulations are commonly used to see how proteins or other molecules move over time. This can contribute to the understanding of molecular systems, which cannot be seen with the eye or even with the best microscopes. In our application of forming a hydrogel with ELPs, it can be useful to see how the ELPs will behave at different temperatures. Above a certain transition temperature (Tt), there is a disruption in the thermal energy of the water molecules around the ELPs, which causes the water molecules to move away from the ELPs. This makes room for interactions between the ELPs. It has been found that there will be conformational changes in another ELP sequence when the temperature rises above Tt, which makes them more ordered and increases the hydrophobicity[8].

It is interesting to see how this applies to our own ELP. The MD simulation will give a better understanding of the differences in the structure below and above the Tt and it will help with understanding what the ELP structurally looks like. We have analyzed different properties of the proteins, like the secondary structure formation, compactness, and the surface area that is accessible for the solvent, to see how the conformation differs at specific temperatures below and above the Tt.

Methods

MD simulations are typically computationally expensive. This made it necessary to shorten our ELP sequence to make the computational time manageable. We chose to shorten the hydrophobic parts to I[9] instead of I[60] and to shorten the hydrophilic part to A[50] instead of A[120] or A[100], as shown in Figure 19. The Leucine zippers stay the same. In this length reduction, the proportions are not the same as in the original ELP. However, we wanted to make sure that the hydrophobic parts were not shorter than the leucine zippers. This maintains a greater distance between the Leucine zippers and the hydrophilic part, and it prevents unrealistic folding. The hydrophilic part is kept long, since this part gives the ELP its length and it will provide relevant information on how the protein behaves when it folds. The choice to change the sequence might affect the outcome a bit, however we still expect to find relevant information on how the ELP behaves at different temperatures.

Figure 19: ELP construct used for the MD simulation. In this construct, the hydrophobic VPGIG parts are repeated 10 times instead of 60 times and the hydrophilic part in the middle is repeated 10 times instead of 24 or 20 times. At the ends, the two different Leucine zippers are located, which remain the same.

A generalized born implicit solvent is used to fold the protein. In an implicit solvent, the solvent (in this case water) molecules are not individually modeled, but with an approximation of the mean force from the solvent on the protein. Different theories can be used for this, and we chose to use the generalized born model. Implicit solvents are faster than an explicit solvent with simulated water molecules and more detailed than a simulation in vacuum[10]. The MD simulation is performed at temperatures 3 °C, 15 °C, 25 °C, and 37 °C, since the transition temperature was found to be 20 °C in the conducted UV-VIS experiments. The simulations are run in triplet to reduce the chance of coincidence in the results. This is necessary, since the initial velocities are chosen by the NAMD software and thus differ every time. Different velocities can cause differences in the final stable conformations, which can affect the results of the analytical methods.

Just like in López Barreiro et al., first energy minimalization is done for 20 000 time steps. Then, Langevin dynamics is used with Generalized Born implicit solvent. The simulation step was 2 fs with a simulation time of 50 ns. The cutoff distance is 18 Å, the switch distance 16 Å[11].

Analysis methods

The results are analyzed with Root Mean Square Distance (RMSD), Radius of Gyration and relative shape anisotropy, Ratio of Secondary Structures, and Solvent Accessible Surface Area (SASA) to show the different behaviors of the ELP at different temperatures. They are calculated with the MDTraj Python library. For the analysis, methods are calculated for the last 5 ns, which contain 1000 frames. This is chosen based on the RMSD results, where a stable conformation is shown at all temperatures. Even when the conformation is stable, there are still some small conformation differences. These differences are taken into account with this distribution of the results of the last 5 ns.

Results

In Figure 20, a video is shown from one of the simulations at 25°C. Not all videos are included, since the process is comparable for the other simulations, especially when comparing them visually. Some frames of those simulations can be found in the file below. When comparing the simulations at the four temperatures, they all start with folding at the ends of the ELP. This was expected since those parts are the most hydrophobic. Some of the simulations also show the folding of the ELP in the middle at the more hydrophilic part. However, there is no clear relationship between temperature and this folding behavior, so it is probably caused by luck in the simulation.

Figure 20: Video with one of the MD simulations at 25°C. The red residues represent the Leucine zipper, the blue residues are the hydrophobic part, and the light-blue residues are the hydrophilic part.

RMSD

RMSD plots are made for each of the three simulations at the four temperatures. They are made for the complete 50 ns, but also for the last 5 ns to show how stable the conformation is related to frame 9000 (45 ns). This is important since there can still be a lot of changes in the structure while the moved distance of the atoms stays the same. Since there are a lot of figures made, they are not visualized here. However, they can be found on the wiki of the TU-Eindhoven 2023 team.

For all of the measurements, the equilibrium seems to be reached. When looking at the plots of the last 5 ns, most RMSD plots have a variability of a maximum of 1 nm. Only the third simulation at 37 °C has a slightly larger difference. An RMSD of 0 describes a stable conformation, but there are always some small deviations. So it can be said that an RMSD of 0.25 nm has a high stability, but also 0.2-0.5 nm will give a high stability[12]. With this information and when looking at the course of the plots, it can still be concluded that the conformations are relatively stable since the RMSD is still very small.

Secondary structure

In Figures 21a, b, and c, the boxplots for the secondary structure analysis of the individual simulations are shown. This shows some variability between the simulations of a single temperature, but also overlapping deviations. This happens especially for the α-helix and β-sheet plots. The simulation results are taken together for each temperature to make them easier to analyze, which is shown in Figure 21d, e, and f.

Figure 21: Secondary structure analysis with a, b, and c for respectively the coils, α-helices, and β-sheets three individual simulations for each of the four temperatures and d, e, and f the overall results with the different simulations per temperature taken together.

Overall, almost no secondary structures seem to be formed. The ratio of the coils is much higher than the ratios for alpha-helices and beta-sheets. It is possible that more secondary structures make the ELP more stable, but that it is stuck in a local minimum. This can be proved with more simulations with for example other force fields used.

When comparing the boxplots of the different temperatures, there does not seem to be a convincing difference in the secondary structure formation. There are some small differences in the medians and distributions, but there is no clear relation between those differences. We can thus conclude that there is no difference in the secondary structure between the temperatures.

Radius of gyration

The radius of gyration results of the individual simulations are shown in Figure 22a. The individual simulation results of 3 and 37°C show large differences in the results. This can be caused by the different velocities that are chosen for the single simulations. Since the individual simulations of 3 and 37°C differ a lot, the total boxplot in Figure 22b shows a large distribution of the results.

Figure 22: Radius of gyration results with a three individual simulations for each of the four temperatures and b the overall results with the different simulations per temperature taken together.

There is a clear difference between the radius of gyration at 15 °C and 25 °C, which are just below and just above the transition temperature. When looking at the median at 37 °C, also these results show a lower radius of gyration than 15 °C. At 3 °C, the radius of gyration is for two of the simulations very low compared to the other temperatures. When looking at the final conformations of the simulations, this can be caused by the shape of the stable configurations. When the shape is more cylindrical instead of spherical, which seems especially the case for the simulations at 15°C, the radius of gyration will increase. As mentioned in the methods section, we also computed the relative shape anisotropy.

Figure 23: Relative shape anisotropy results with a three individual simulations for each of the four temperatures and b the overall results with the different simulations per temperature taken together.

In Figure 23a, the relative shape anisotropy is shown for the individual simulations. There is clearly a difference in the shape between the different simulations, also in those with the same temperature. For example 25°C, which has a clear variability between the simulations. Some of the differences were already expected after the analysis of the radius of gyration. The second simulation at 3°C is clearly less spherical than the other two simulations at that temperature. Overall, when looking at Figure 23b, the shape of the ELP at 15°C is less spherical than at 25°C. Because of the large variabilities between the measurements at 3 and 37°C, it cannot be concluded how this is affected at lower and higher temperatures.

This analysis suggests that there is a difference in the compactness of the ELP below and above the Tt. Below the Tt, the ELP is less compact. At a temperature above the Tt, the ELP will become more hydrophobic with fewer solvent molecules around it. This makes the folded ELP more compact, but also more likely to induce interactions with other ELPs. However, for harder conclusions with the role of solvent molecules included, explicit simulations are needed.

SASA

In Figure 24a, the individual results show some deviations between the simulations at specific temperatures, especially at 3 and 37 °C. When taking the individual simulations together for each temperature in Figure 24b, some small deviations are shown. SASA at 25 °C seems to be lower than SASA at 15 °C, while the median SASA at 3 °C is also lower than the rest of the temperatures, but with a large deviation. The results for 37 °C also have a large deviation, with a median higher than the median of 15 °C. When comparing the results with the relative shape anisotropy, there is a relationship between more spherical simulation results and a lower SASA. There seems to be a smaller difference than the difference found in the radius of gyration and the anisotropy, which might be because ends that stick out probably have more effect on SASA than the relative shape anisotropy.

Figure 24: SASA results with a three individual simulations for each of the four temperatures and b the overall results with the different simulations per temperature taken together.

Conclusion

The MD simulation of the folding of a shortened version of our ELP taught us that an increase in temperature above the Tt results in a more compact and more spherical ELP due to increased hydrophobicity. This causes more and stronger interactions between the ELPs, which is an important driver of the hydrogel formation[8].

The greatest difference in results is found for the radius of gyration and relative shape anisotropy between 15°C and 25°C. SASA also shows a difference, but this is less convincing. There is no difference found in the secondary structure between the different temperatures. When comparing the triplicate simulations, there are large differences for especially 3°C but also 37°C. This makes these results less reliable, but there still seems to be a difference in properties between 15 and 37°C. In further research, explicit solvent simulations can be used with more simulated temperatures. This would give more details for better and possibly more convincing results. Also, more extensive research can be done on the interactions between the ELPs at different temperatures and conditions.

References

[1] Gradišar, H., & Jerala, R. (2010). De novodesign of orthogonal peptide pairs forming parallel coiled-coil heterodimers. Journal of Peptide Science, 17(2), 100–106. https://doi.org/10.1002/psc.1331

[2] Fernández‐Colino, A., Arias, F. J., Alonso, M., & Rodríguez-Cabello, J. C. (2015). Amphiphilic Elastin-Like Block Co-Recombinamers Containing Leucine Zippers: Cooperative Interplay between Both Domains Results in Injectable and Stable Hydrogels. Biomacromolecules, 16(10), 3389–3398. https://doi.org/10.1021/acs.biomac.5b01103

[3] L. E. Contreras-Llano et al., “Engineering Cyborg Bacteria Through Intracellular Hydrogelation,” Adv. Sci., Mar. 2023, doi: 10.1002/ADVS.202204175.

[4] Alber, T. (1992). Structure of the leucine zipper. Current Opinion in Genetics and Development, 2, 205–210

[5] X. Yang, Y. Zhong, D. Wang, and Z. Lu, “A simple colorimetric method for viable bacteria detection based on cell counting Kit-8,” Anal. Methods, vol. 13, no. 43, pp. 5211–5215, Nov. 2021, doi: 10.1039/D1AY01624E

[6] E. L. Foltz, J. W. West, I. H. Breslow, and H. Wallick, “Clinical pharmacology of pivampicillin.,” Antimicrob. Agents Chemother., vol. 10, no. 3, pp. 442–454, 1970, doi: 10.47363/jprsr/2022(3)129.

[7] Poudyal, R. R., Guth-Metzler, R., Veenis, A. J., Frankel, E. A., Keating, C. D., & Bevilacqua, P. C. (2019). Template-directed RNA polymerization and enhanced ribozyme catalysis inside membraneless compartments formed by coacervates. Nature Communications, 10(1). https://doi.org/10.1038/s41467-019-08353-4

[8] N. K. Li, F. G. Quiroz, C. K. Hall, A. Chilkoti, and Y. G. Yingling, “Molecular description of the lcst behavior of an elastin-like polypeptide,” Biomacromolecules, vol. 15, no. 10, pp. 3522–3530, Oct. 2014, doi: 10.1021/BM500658W/SUPPL_FILE/BM500658W_SI_001.PDF.

[9] J. Kleinjung and F. Fraternali, “Design and application of implicit solvent models in biomolecular simulations,” Curr Opin Struct Biol, vol. 25, no. 100, p. 126, 2014, doi: 10.1016/J.SBI.2014.04.003.

[10] J. Iwata and T. Ando, “Molecular Dynamics Study on Behavior of Resist Molecules in UV-Nanoimprint Lithography Filling Process,” Nanomaterials, vol. 12, no. 15, Aug. 2022, doi: 10.3390/NANO12152554/S1.

[11] D. López Barreiro, A. Folch-Fortuny, I. Muntz, J. C. Thies, C. M. J. Sagt, and G. H. Koenderink, “Sequence Control of the Self-Assembly of Elastin-Like Polypeptides into Hydrogels with Bespoke Viscoelastic and Structural Properties,” Biomacromolecules, vol. 24, no. 1, pp. 489–501, Jan. 2023, doi: 10.1021/ACS.BIOMAC.2C01405/ASSET/IMAGES/LARGE/BM2C01405_0003.JPEG.

[12] M. Arnittali, A. N. Rissanou, and V. Harmandaris, “ScienceDirect Structure Of Biomolecules Through Molecular Dynamics Simulations,” 2019, doi: 10.1016/j.procs.2019.08.181.