Difference between revisions of "Team:UCopenhagen/Engineering"

 
(29 intermediate revisions by 3 users not shown)
Line 1: Line 1:
 
<!--{{IGEM_TopBar}}-->
 
<!--{{IGEM_TopBar}}-->
{{UCopenhagen/Header3}}
 
 
{{UCopenhagen}}
 
{{UCopenhagen}}
{{UCopenhagen/Back_To_Top_Button}}
+
{{Template: UCopenhagen/Header3}}
{{UCopenhagen/Footnotes}}
+
{{Template: UCopenhagen/Back_To_Top_Button}}
 
+
{{Template: UCopenhagen/Footnotes}}
 
+
 
+
 
<html>
 
<html>
 
   <head>
 
   <head>
 
     <style>
 
     <style>
       .title_piece h2 {
+
       figcaption {font-family: Avenir, Arial; font-size: 14px; text-align: center;}
        font-family: "Montserrat";
+
        font-size: 18px;
+
      }
+
  
       h3 {
+
       .txt-btwn h3 {
         font-family: "Montserrat";
+
        color: #F3557E;
 +
         font-family: Bavro;
 +
        font-size: 35px;
 +
        line-height: 1.3;
 +
        text-align: center;
 +
        margin: 5% 0 2.5% 0;
 
       }
 
       }
  
       .big_wrapper a {
+
       .txt-btwn h4 {
         cursor: pointer;
+
         font-family: Montserrat;
         text-decoration: none;
+
         font-size: 28px;
 +
        line-height: 1.3;
 +
        margin: 5% 0 2.5% 0;
 
       }
 
       }
  
       .big_wrapper a:hover {
+
       .txt-btwn h5, .everything h5 {
         animation: levitating_pills 1.5s infinite;
+
         font-family: Montserrat;
 +
        font-size: 24px;
 +
        line-height: 1.3;
 +
        margin: 5% 0 2.5% 0;
 
       }
 
       }
  
      @keyframes levitating_pills {
+
       .txt-btwn h6 {
        0% {
+
         font-family: Montserrat;
          transform: translate(0,0);
+
         font-size: 20px;
        }
+
        50% {
+
          transform: translate(0,7%);
+
        }
+
        100% {
+
          transform: translate(0,0);
+
        }
+
      }
+
 
+
       .big_wrapper {
+
         display: flex;
+
         flex-direction: column;
+
 
         text-align: center;
 
         text-align: center;
         z-index: -2;
+
         line-height: 1.3;
         width: 100%;
+
         margin: 5% 0 2.5% 0;
 
       }
 
       }
  
       .pill_bottle {
+
       .txt-btwn {
         margin: auto;
+
         font-family: Avenir, Arial;
         margin-top: -165px;
+
         font-size: 16px;
         z-index: -1;
+
         line-height: 1.4;
 
       }
 
       }
  
       .overview_wrapper {
+
       .everything {
        display: flex;
+
        flex-direction: row;
+
        width: 75%;
+
        margin: 0 auto;
+
        justify-content: space-around;
+
        padding: 2.5%;
+
      }
+
 
+
      @keyframes round_guy {
+
        0% {
+
          transform: translate(185%,-100%);
+
        }
+
        100% {
+
          transform: translate(0,0);
+
        }
+
      }
+
 
+
      #Round_Pill {
+
        animation-name: round_guy;
+
        animation-duration: 2s;
+
        z-index: -2;
+
        padding: 6.5% 0% 0% 7%;
+
      }
+
 
+
      @keyframes triangle_guy {
+
        0% {
+
          transform: translate(-50%,-100%);
+
        }
+
        100% {
+
          transform: translate(0,0);
+
        }
+
      }
+
 
+
      #Triangular_Pill {
+
        animation-name: triangle_guy;
+
        animation-duration: 2s;
+
        z-index: -2;
+
        padding: 2.5% 0% 0% 7%;
+
      }
+
 
+
      @keyframes long_guy {
+
        0% {
+
          transform: translate(-270%,-100%);
+
        }
+
        100% {
+
          transform: translate(0,0);
+
        }
+
      }
+
 
+
      #Long_Pill {
+
        animation-name: long_guy;
+
        animation-duration: 2s;
+
        z-index: -2;
+
        padding: 5.5% 0% 0% 0%;
+
      }
+
 
+
      @keyframes text_fade_in {
+
        0% {
+
          opacity:0;
+
        }
+
        10% {
+
          opacity:0;
+
        }
+
        50% {
+
          opacity:0;
+
        }
+
        100% {
+
          opacity:1;
+
        }
+
      }
+
 
+
      .title_piece {
+
        animation-name: text_fade_in;
+
        animation-duration: 3s;
+
        padding: 4% 5% 0% 0%;
+
        width: 100%;
+
      }
+
 
+
      .big_wrapper {
+
 
         display: flex;
 
         display: flex;
 +
        flex-direction: column;
 
         justify-content: center;
 
         justify-content: center;
      }
+
        align-items: center;
 
+
        font-family: Avenir, Arial;
      .tabcontent {
+
        font-size: 16px;
 
         margin: 0 15% 0 15%;
 
         margin: 0 15% 0 15%;
        width: 100%;
 
        height: auto;
 
 
       }
 
       }
  
       .hand {
+
       .everything h2 {
         display: flex;
+
         color: #F3557E;
         justify-content: center;
+
         font-family: Bavro;
         width: 100%;
+
         font-size: 35px;
         overflow: hidden;
+
         line-height: 1.3;
 +
        text-align: center;
 +
        margin: 5% 0 2.5% 0;
 
       }
 
       }
  
       .hand img {
+
       .everything img {
         display: block;
+
         width: 65%;
         margin: auto;
+
         margin: 2.5%;
 
       }
 
       }
 
  
 
       .page-wrapper{
 
       .page-wrapper{
         margin-bottom: -20%;
+
         margin-bottom: -25%;
         margin-top: -2.5%;
+
         margin-top: -5%;
 
         max-width: 100%;
 
         max-width: 100%;
 
         background-image: url("https://static.igem.org/mediawiki/2020/a/a3/T--UCopenhagen--watercolor_bg_smaller.png");
 
         background-image: url("https://static.igem.org/mediawiki/2020/a/a3/T--UCopenhagen--watercolor_bg_smaller.png");
Line 166: Line 79:
 
         background-color: #fcecdf;
 
         background-color: #fcecdf;
 
         padding-top: 5%;
 
         padding-top: 5%;
 +
        padding-bottom: 25%;
 
         text-align: justify;
 
         text-align: justify;
 
       }
 
       }
  
       .titles {
+
       #top-image img {
        background-image: url('https://static.igem.org/mediawiki/2020/3/3b/T--UCopenhagen--watercolor_mix2.png');
+
         width: 100%;
         font-family: Bavro;
+
         height: 94vh;
         color: #F3557E;
+
         position: relative;
         font-size: 42px;
+
         -o-object-fit: cover;
        background-position: top;
+
         object-fit: cover;
         background-repeat: no-repeat;
+
         margin-top: -3em;
        background-size: 100% 100%;
+
         padding: 15.5% 22% 14.5% 19.5%;
+
         margin-bottom: -5%;
+
 
       }
 
       }
  
       .first {
+
       .oversigt {
        margin-top: -15%;
+
      }
+
 
+
      #prob, #SHP, #GHP {
+
 
         display: flex;
 
         display: flex;
        flex-direction: column;
 
 
         align-items: center;
 
         align-items: center;
         margin: 0% 15% 0% 15%;
+
         flex-direction: column;
         text-align: justify;
+
        width: 100%;
 +
         text-align: center;
 +
        justify-content: center;
 
       }
 
       }
  
       .todelt {
+
       .ansigter {
 
         display: flex;
 
         display: flex;
 +
        align-items: center;
 
         flex-direction: row;
 
         flex-direction: row;
         margin: 2% 0 2% 0;
+
         justify-content: space-between;
 
       }
 
       }
  
       .todelt .txt {
+
       .ansigter img {
        display: flex;
+
        padding: 0 1% 0 1%;
+
        flex: 3;
+
      }
+
 
+
 
+
      .txt {
+
        display: flex;
+
        flex-direction: column;
+
        justify-content: center;
+
      }
+
 
+
      .txt img {
+
        object-fit: contain;
+
        width: 45%;
+
        height: 100%;
+
        padding-bottom: 3.5%;
+
        padding-top: 3.5%;
+
        margin: 0 auto;
+
      }
+
 
+
      .txt h2 {
+
        font-family: Montserrat;
+
        font-size: 26px;
+
      }
+
 
+
      .image {
+
        align-items: center;
+
        display: flex;
+
 
         flex: 1;
 
         flex: 1;
        margin: 2%;
 
 
       }
 
       }
  
      .image img {
 
        align-self: center;
 
        width: 100%;
 
      }
 
 
      .format {
 
        font-family: Avenir;
 
        font-size: 18px;
 
        line-height: 1.3;
 
        text-align: justify;
 
      }
 
#top-image img {
 
width: 100%;
 
height: 94vh;
 
margin-top: -3em;
 
position: relative;
 
-o-object-fit: cover;
 
object-fit: cover;
 
}
 
 
     </style>
 
     </style>
 
   </head>
 
   </head>
<body>
+
  <body>
 +
    <!-- placeholders galore!-->
 +
    <div id="top-image">
 +
      <img src="https://static.igem.org/mediawiki/2020/b/b5/T--UCopenhagen--ES-cover.png">
 +
    </div>
 +
    <!-- placeholders galore continue from here but a bit better -->
  
<div id="top-image">
+
    <div class="page-wrapper">
  <img src="https://static.igem.org/mediawiki/2020/b/b5/T--UCopenhagen--ES-cover.png">
+
      <div class="everything">
  </div>
+
        <h2>Introduction</h2>
 
+
        <div class="txt-btwn">
            <div class="format">
+
          This page documents our success with iterative engineering of novel G protein alpha subunits in yeast Saccharomyces cerevisiae. Protein modeling in Rosetta took us most of the way, finally accompanied by wet lab experiments. Our efforts followed the traditional engineering design cycle (see figure 1), with multiple
              As is always the case for engineering work, our dry lab followed the engineering design cycle closely, with going back and forth between extensive research, testing of our models, and going back to the drawing boards multiple times to refine our models and results.
+
          phases of research followed by testing of our models and continuously refining them. Individual steps are noted chronologically as they happened, thus proceeding from relatively simple to more advanced and specialized methods. Finally, the results challenged our initial assumptions about protein interactions.<br>
            </div>
+
        </div>
            <div class="image">
+
        <img src="https://static.igem.org/mediawiki/2020/d/d1/T--UCopenhagen--engineering_cycle.png">
              <img src="https://static.igem.org/mediawiki/2020/d/d1/T--UCopenhagen--engineering_cycle.png">
+
        <figcaption><b>Fig. 1: Expanded engineering cycle.</b><br>
<figcaption style="font-size: 16px;">Expanded engineering cycle that we followed in our engineering processes.
+
          We closely followed the expanded engineering cycle in our engineering process.</figcaption>
</figcaption>
+
            </div>
+
 
+
 
+
 
+
<p>As described on our project design page (link), our most ambitious design called for modification of the alpha subunit of endogenous yeast G protein. The challenge here was that we wanted to maintain its natural functionality while expanding it to serving our purposes.</p>
+
<p>This meant we needed to design a G-alpha that would:</p>
+
<ol>
+
<li>A) have a natural-like affinity for the beta/gamma complex of the yeast G protein to stay bound to the beta/gamma complex when no signal was present, thus inhibiting leaky activation of the pathway.</li>
+
<li>B) have low enough affinity for the beta/gamma complex after cleavage by our TEV-protease so that the fragments dissociated from the beta/gamma complex, thus not interfering with their natural function.</li>
+
</ol>
+
<p>&nbsp;</p>
+
<p>Just to refresh, our work on G-alpha hinged on finding a mutant that preserved natural function when there was no signal (i.e. no cytokine), and also dissociated appropriately in the presence of a signal.</p>
+
<p>&nbsp;</p>
+
<p>This page presents all stages of G-alpha engineering in a chronological order. The progress followed the engineering cycle (design-build-test-learn) as we iterated our designs several times. As we utilized primarily dry-lab experiments, the whole cycle often spans over each stage. Every stage was informed by all the previous ones. FIGURE&nbsp; WITH THE ENGINEERING CYLE</p>
+
<p>&nbsp;</p>
+
<p>1) stage one - design</p>
+
<p>In the first round, we based our design choices on knowledge about the interaction interface of the alpha and beta subunits in bovine G protein, as we found no crystal structure of yeast G protein. Immediately, we noticed that the yeast G protein contained an &lsquo;insertion sequence&rsquo; (highlighted in the figure) that other G proteins did not, meaning that it probably was not essential for the general functionality of G-alpha. This attribute made the insertion sequence a great starting point for our search for a mutation site.</p>
+
<p>&nbsp;</p>
+
 
+
<img src="https://static.igem.org/mediawiki/2020/1/13/T--UCopenhagen--insertion-seq.png">
+
<figcaption style="font-size: 16px;">Sequence of naturally occurring G protein alpha subunit from baker's yeast S. cerevisiae. A unique insertion sequence is highlighted.
+
</figcaption>
+
 
+
<p>2) design-build</p>
+
<p>We received a file with empirically observed binding interface of bovine G-alpha and G-beta through our supervisors. We could translate known characteristics of the bovine G protein into what we expected to see in the yeast G protein because of the similarity of their sequences that we observed after sequence alignment. Our initial analysis was based on our own knowledge of the biophysical properties of different amino acids. As illustrated by the figure, we identified what we believed were most essential residues (red color).</p>
+
<p>&nbsp;</p>
+
 
+
<img src="https://static.igem.org/mediawiki/2020/6/63/T--UCopenhagen--bovine-g.png">
+
<figcaption style="font-size: 16px;">Crystal structure of bovine G protein complex. Red-colored residues are the ones we identified as most important for binding of alpha and beta subunits.
+
</figcaption>
+
 
+
<p>3) test-learn</p>
+
<p>Our dry lab collaboration with iGEM team Aalto-Helsinki started taking shape and we accepted their suggestion to use Rosetta for our protein engineering. Firstly, we supplemented our own analysis with an alanine scan (performed by Robetta server), which provided us with quantitative predictions of selected residues&rsquo; role in G-alpha/G-beta binding interface. Alanine scan was a script that calculated how Gibbs free energy of the two proteins changed upon changing a single residue to alanine. The results confirmed our intuitions to a good extent and offered several valuable insights that led to refinement of our G-alpha designs. As can be seen in the figure, residue 227 was predicted as the single most important residue in the protein-protein interaction. See the dry lab page to learn more about technical details of the method (DRY LAB PAGE LINK).</p>
+
<p>&nbsp;</p>
+
 
+
<img src="https://static.igem.org/mediawiki/2020/1/11/T--UCopenhagen--Alanine_Scan.png">
+
<figcaption style="font-size: 16px;">Change in Gibbs free energy (in kcal/mol) generated by switching a single residue into alanine.
+
</figcaption>
+
 
+
<p>4) design-build-test-learn</p>
+
<p>We designed three versions of G-alpha ranging from more conservative to more extreme modifications. We confirmed that none of the edits lead to unexpected conformational changes by SwissProt modeling. We used Cello to predict protein localization &ndash; correct localization to the membrane was crucial for biosensor functionality. As shown in the figure, our novel designs yielded better localization scores than the wildtype. Cytoplasmic localization was desirable, as certain enzymes in the cytoplasm were known to modify G-alpha for it to go to the membrane.</p>
+
<p>&nbsp;</p>
+
 
+
<img src="https://static.igem.org/mediawiki/2020/0/05/T--UCopenhagen--localization-prediction.png">
+
<figcaption style="font-size: 16px;">Localization prediction with a reliability score (y-axis) for different G-alpha mutants (labels on x-axis).
+
</figcaption>
+
 
+
<p>5) design-build-test-learn</p>
+
<p>Here we faced the reality of laboratory experimentation. Three G-alpha mutants with many diverse mutations could lead to lots of ambiguity as the laboratory results could be hard to interpret due to inability to monitor what mutation had profound impact. Time constraint in the wet lab motivated us to limit the number of possible mutation points. We eliminated point mutations and kept only cleavage sites for TEV protease. To improve the quality of future experiments, we expanded the range of G-alpha candidates from three to seven by introducing all combinations of the three selected mutations. We checked folding predictions for all seven novel G-alphas which resulted in a significant overlap (as visible in the figure). None of the proteins seemed likely to take on a faulty conformation according to SwissProt.</p>
+
<p>&nbsp;</p>
+
 
+
<img src="https://static.igem.org/mediawiki/2020/4/46/T--UCopenhagen--mutants-aligned.png">
+
<figcaption style="font-size: 16px;">Alignment of predicted structures of our G-alpha mutants.
+
</figcaption>
+
 
+
<p>6) test-learn</p>
+
<p>We were still wondering how to predict the behavior of fragments generated by cleavage of a mutant G-alpha by TEV protease. We decided to run such a simulation in Rosetta. Thanks to Kurt V. Mikkelsen, we got the opportunity to use his supercomputer to perform such a computationally demanding analysis. We simulated change in Gibbs free energy of the seven mutant G-alphas (see the figure) and the nine fragments formed after cleavage (see figure). The rational was to find a candidate G-alpha with binding affinity like the wild type, but with lowest possible binding affinity of its fragments after cleavage. (For more technical details about how the simulation was set up, see the dry lab page.) The simulation lead to some unexpected results &ndash; majority of fragments seemed to have higher affinities to G-beta, compared to mutant G-alphas.</p>
+
<p>&nbsp;&nbsp;</p>
+
<p>First, we were startled by the findings but after thorough discussions concluded that the overall G-alpha conformation probably restricted some degrees of freedom of the high-affinity regions. These regions facilitated most interaction between G-alpha and G-beta, as we learned in one of the earlier stages. Our cleavage sites were located between these high-affinity regions, therefore, cleavage effectively freed the high-affinity regions from the previous restrictive conformation and allowed them to bind G-beta more tightly. Obviously, these findings turned our initial intuitions about the protein-protein interaction on its head. We also realized that affinities do not necessarily provided the answer to successful signal transduction, as further analysis showed that conformational changes probably lead to Ste5 recruitment by G-beta (sites where Ste5 bound were not obstructed by G-alpha). Information about affinities without knowledge about conformational changes of G-beta and their causes could not answer our questions. Anyways, mutant M124 was the only whose some fragments had lower affinity to G-beta compared to the mutant itself, so we decided it was the most promising candidate if we followed our reasoning behind the rational of this simulation.</p>
+
<p>&nbsp;</p>
+
 
+
<img src="https://static.igem.org/mediawiki/2020/b/ba/T--UCopenhagen--mutants-affinity.png">
+
<figcaption style="font-size: 16px;">Changes of Gibbs free energy of our G-alpha mutants (compared to the wildtype). Numbers in their labels signify what cleavage sites were included in each mutant.
+
</figcaption>
+
 
+
<img src="https://static.igem.org/mediawiki/2020/f/fe/T--UCopenhagen--fragments-affinity.png">
+
<figcaption style="font-size: 16px;">Change in Gibbs free energy of fragments produced by cleavage of G-alpha mutants. X-axis shows length of each fragment.
+
</figcaption>
+
 
+
<p>7) test-learn</p>
+
<p>Following the conclusion of the previous Rosetta simulation, we designed and ran another simulation to obtain more refined details about fragments from M124 cleavage. The results showed that all generated fragments had a significantly higher affinity to the G-beta compared to the mutant G-alpha M124. Although, exploration of PDB files revealed that spatial assignment of the fragments showed clashes (example in the figure) so that affinity of each fragment could be limited by other presence of the other fragments, ultimately leading to diminished presence at G-beta. Nevertheless, results of our simulations could not inform us about conformational changes of G-beta which protected us from formulating clear conclusions.</p>
+
<p>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</p>
+
<p>PDB (m124-fragments) HERE</p>
+
<p>&nbsp;</p>
+
 
+
<img src="https://static.igem.org/mediawiki/2020/9/99/T--UCopenhagen--m124.png">
+
<figcaption style="font-size: 16px;">Change in Gibbs free energy of fragments produced by cleavage of G-alpha mutant M124. Labels on X-axis signify position of each fragment in the original protein.
+
</figcaption>
+
 
+
<img src="https://static.igem.org/mediawiki/2020/5/5d/T--UCopenhagen--m124-clash.png">
+
<figcaption style="font-size: 16px;">An example of clashing positions of the fragments.
+
</figcaption>
+
 
+
 
+
<p>8) wet lab</p>
+
<p>The wet lab managed to create constructs with four different mutant GPA1s named GPA1m1, GPA1m2, GPA1m3, GPA1m4 . mutant 1, 2 and 4 were the same versions as the ones that were investigated with modeling. The wt GPA1 was replaced with different GPA1 mutants in an adenosine biosensor strain of <i>Saccharomyces Cereviae</i> previously created by Sotirios Kampranis' lab. Furthermore an inducible TEV protease was inserted in the same locus as the rest of the biosensor constructs. In this way it was possible to test out both the mapintainance of the original functions of GPA1 in a GPCR receptor system and to verify that the GPA1 mutants were cleaved på the induced TEV protease and that this would result in higher concentrations of the reporter.
+
 
<br>
 
<br>
The experiments showed that the amount of reporter was increased drastically in all biosensors when the TEV protease was induced by growing the cells in galactose. This included the biosensor with wt GPA1 which made it hard to conclude whether some of the increased signal was caused by cleavage of the GPA1 mutant or that the increase in signaling came from other parts of the system.
 
 
<br>
 
<br>
The experiments further showed a loss-of-function in the sense that none of the biosensors with mutated GPA1 proteins were able to respond to different levels of adenosine whereas we were able to reproduce the A2a adenosine biosensor created in Sotirios Kampranis lab when we used the wt GPA1 and grew the cells in glucose media.
+
        <div class="txt-btwn">
 
+
          As described on our <a href="https://2020.igem.org/Team:UCopenhagen/Design">project design page</a>,
</p>
+
          our most ambitious design called for modification of the alpha subunit of
<p>&nbsp;</p>
+
          endogenous yeast G protein. The challenge here was that we wanted to maintain
 
+
          its natural functionality while expanding it to serving our purposes.
<p>9) future</p>
+
          <br>
<p>The next round of simulations could be based on molecular dynamics for us to see whether the results of Rosetta could be trusted. Apart from that, the focus could shift from affinities to conformations, and particularly G-beta could be more explored.  
+
          This meant we needed to design a G-alpha that:<br>
 +
          <ul>
 +
            <li><b>A)</b> has a natural-like affinity for the beta/gamma complex of the
 +
              yeast G protein to stay bound to the beta/gamma complex when no signal
 +
              was present, thus inhibiting leaky activation of the pathway.</li>
 +
            <li><b>B)</b> has a low enough affinity for the beta/gamma complex after
 +
              cleavage by our TEV-protease to enable dissociation from
 +
              the beta/gamma complex.</li>
 +
          </ul>
 +
          <br>
 +
          Our work on G-alpha hinged on finding a mutant that preserved natural
 +
          function when there was no signal (i.e. no cytokine), and also dissociated appropriately in the presence of a signal.
 +
          <br>
 +
          <br>
 +
          This page presents all stages of G-alpha engineering in a chronological
 +
          order. The progress followed the engineering cycle (design-build-test-learn)
 +
          as we iterated our designs several times. As we utilized primarily dry-lab
 +
          experiments, the whole cycle often spans over each stage. Every stage was
 +
          informed by all the previous ones.
 +
          <h2>Stage One - Design</h2>
 +
          We based our design choices on knowledge about the interaction interface
 +
          of the alpha and beta subunits in bovine G protein, as we found no crystal
 +
          structure of yeast G protein. Immediately, we noticed that the yeast G
 +
          protein contained an ‘insertion sequence’ (highlighted in figure 2) that
 +
          other G proteins did not. This indicated that the sequence was not essential
 +
          for the general functionality of the G-alpha subunit. As such, we chose to
 +
          search this sequence for a mutagenesis site.
 +
          <br>
 
<br>
 
<br>
In the wet lab the next step would be to understand why the signaling also increases in the biosensor with wt GPA1 when the cells are grown in galactose and raffinose media. First we would create a strain of the biosensor without the inducible TEV protease to see whether the spike in signaling is caused by the induction of the TEV protease or something else related to the change of the media. If the TEV protease does not seem to be the cause we would investigate the promoters that we use, some of these are glycolytic promoters like pPGK1 and the different media might cause slightly different expression levels even though the promoters are normally classified as constituent. To further investigate the loss-of-function of the mutants, we would like to create strains of the biosensor with the mutated GPA1 proteins but without the inducible TEV promoter. This would enable us to see whether leaky expression of the TEV protease in the glucose fed cells might have cause the loss-of-function</p>
+
Based on literature review, we identified the residue sequence that the TEV-protease cleaves: <b>EXLYΦQφ</b>, where X is any residue, Φ is any large or medium hydrophobic residue, and φ is any small hydrophobic or polar residue. For the sake of interpretability of experiments, however, it was necessary to introduce a single sequence type at all times, as the TEV protease may display different proteolytic rates depending on the recognition site (Dougherty et al., 1989). Following extensive literature search, we opted for a widely used sequence: <b>ENLYFQG</b>.
<p>&nbsp;</p>
+
  
 
+
          </div>
 
+
          <img src="https://static.igem.org/mediawiki/2020/1/13/T--UCopenhagen--insertion-seq.png">
 
+
          <figcaption>
 
+
            <b>Figure 2: Amino acid sequence of the <em>S. cerevisiae</em> G protein alpha subunit.</b>
</body>
+
            <br>The highlighted sequence was initially selected as an optimal region for mutagenesis.
 +
          </figcaption>
 +
          <div class="txt-btwn">
 +
            <h2>Stage Two - Design-build</h2>
 +
            We received a file with the empirically observed binding interface of
 +
            bovine G-alpha and G-beta through our supervisors. We could translate
 +
            known characteristics of the bovine G protein into what we expected to
 +
            see in the yeast G protein because of the similarity of their sequences.
 +
            Based on the biophysical properties of different amino acids, we identified
 +
            what we believed were the most essential residues for the binding between the
 +
            alpha and beta subunits (see figure 3).
 +
          </div>
 +
          <img src="https://static.igem.org/mediawiki/2020/6/63/T--UCopenhagen--bovine-g.png">
 +
          <figcaption>
 +
            <b>Figure 3: Crystal structure of the G protein complex.</b>
 +
            <br>The red-colored residues are the ones we identified as the most important for binding of the alpha and beta subunits to each other.
 +
          </figcaption>
 +
          <div class="txt-btwn">
 +
            <h2>Stage Three - Test & Learn</h2>
 +
            Based on suggestions from the Sinisens (Aalto-Helsinki iGEM2020 team), we chose to
 +
            use Rosetta for our protein engineering. Rosetta is a powerful tool for
 +
            mapping and characterizing protein structures and interactions, and, we
 +
            theorized, good for mutagenetic analyses of our G alpha protein.
 +
            <br>
 +
            Using Rosetta, we supplemented our own analysis with an alanine scan,
 +
            which provided us with quantitative predictions of the role of selected residues
 +
            in the G-alpha in terms of binding to G-beta. Rosetta's alanine scan is
 +
            a script that calculates how the Gibb's free energy of two proteins change
 +
            upon changing a single amino acid residue to an alanine (see <a href="https://2020.igem.org/Team:UCopenhagen/Model">modeling</a> for more details).<br>
 +
            <br>
 +
            The results of the alanine scan confirmed our intuitions to a good extent,
 +
            and offered several valuable insight points that led to the refinement of
 +
            our G-alpha designs. As illustrated on figure 4, residue number 227 was
 +
            predicted as the single most important residue in the protein-protein interaction
 +
            between the alpha and beta subunits. This result, among mutagenetic data
 +
            found in literature about what the most functionally important amino acid
 +
            residues in G-alpha are<a href="#mutant" aria-describedby="footnote-label" id="mutant-ref"> </a> , guided us in discerning between which ones we
 +
            could change, and which ones were critical for the correct functionality of G-alpha.<br>
 +
            </div>
 +
            <img src="https://static.igem.org/mediawiki/2020/1/11/T--UCopenhagen--Alanine_Scan.png">
 +
            <figcaption>
 +
              <b>Figure 4: Change in Gibbs free energy (in kcal/mol) generated by replacing the given residues with alanine.</b> As Gibbs free energy correlates with a change in affinity, it appears that residue 227 is the most important one for the binding interface of G-alpha and G-beta proteins.
 +
            </figcaption>
 +
            <div class="txt-btwn">
 +
              <h2>Stage Four - Redesigning, Building, Testing & Learning</h2>
 +
              By compiling the information we had received so far, we isolated fitting mutation points in G alpha, where we
 +
              <ol>
 +
                <li><b>1.</b> Avoid changing residues that lead to loss of function in G alpha (Gladue, D. P. 2008).</li>
 +
                <li><b>2.</b> Avoid changing residues identified as very important through Rosetta's alanine scan.</li>
 +
                <li><b>3.</b> Change residues located in the insertion sequence and other non-conserved residues, by comparing the yeast G-alpha to the bovine G-alpha.</li>
 +
                <li><b>4.</b> Change residues located in loops based on visual analysis in Pymol.</li>
 +
              </ol>
 +
              <div style="display: flex; flex-direction:row;">
 +
                <div style="flex: 2; margin:3.5%;">
 +
                  <img style="width:100%;" src="https://static.igem.org/mediawiki/2020/0/05/T--UCopenhagen--localization-prediction.png">
 +
                    <figcaption>
 +
                      <b>Figure 5: Localization prediction of modified alpha subunits.</b>
 +
                      <br>Cello was used to predict the localization of our three GPA1-versions.
 +
                      For all three versions, an increased cytoplasmic reliability score was observed compared to the wildtype.
 +
                    </figcaption>
 +
                </div>
 +
                <div style="display: flex; font-family: Avenir, Arial; font-size: 16px; flex: 2.5; align-items: center;">
 +
                  With this information in mind, we designed three version of G-alpha, each with
 +
                  varying degrees of modifications. Using SwissProt protein modeling, we then
 +
                  confirmed that none of our edits lead to unexpected conformational changes.
 +
                  In addition, we used Cello to predict protein localization to the plasma membrane,
 +
                  as this is crucial for the functionality of our biosensor. As shown on figure 5,
 +
                  our designs yielded even better localization scores than the wildtype.<br>
 +
                  As shown on the figure, some cytoplasmic localization was also predicted,
 +
                  but this was desirable as certain enzymes in the cytoplasm are known to modify
 +
                  G-alpha for it to go to the membrane.
 +
                </div>
 +
              </div>
 +
              <h2>Stage Five - Design, Build, Test & Learn</h2>
 +
              Here we faced the reality of laboratory experimentation. Three G-alpha
 +
              mutants with many diverse mutations could lead to lots of ambiguity
 +
              as the laboratory results could be hard to interpret due to inability
 +
              to monitor what mutation had a profound impact. Time constraints in the
 +
              wet lab motivated us to limit the number of possible mutation points.
 +
              We eliminated point mutations and kept only cleavage sites for TEV
 +
              protease. To improve the quality of future experiments, we expanded
 +
              the range of G-alpha candidates from three to seven by introducing
 +
              all combinations of the three selected mutations. We checked folding
 +
              predictions for all seven novel G-alphas which resulted in a significant
 +
              overlap (as visible in the figure). None of the proteins seemed likely
 +
              to take on a faulty conformation according to the SwissProt analysis.
 +
            </div>
 +
            <img src="https://static.igem.org/mediawiki/2020/4/46/T--UCopenhagen--mutants-aligned.png">
 +
            <figcaption><b>Figure 6: Alignment of predicted structures of our G-alpha mutants.</b> Apart from the intrinsically disordered region on the left, the predicted structures overlap fairly well. This indicates that the conformation of none of the G-alpha mutants was severely impacted by their mutations.
 +
            </figcaption>
 +
            <div class="txt-btwn">
 +
              <h2>Stage Six - Test & Learn</h2>
 +
              <div style="display: flex; flex-direction:row;">
 +
                  <div style="display: flex; font-family: Avenir, Arial; font-size: 16px; flex: 2.5; align-items: center;">
 +
                    In order to predict the behavior of the fragments generated by cleavage
 +
                    of our mutant G-alpha by the TEV protease, we decided to run a simulation
 +
                    in Rosetta (see modeling for more details). However, as this was an extremely demanding computational
 +
                    analysis we could not perform ourselves, we reached out to Assoc. Prof.
 +
                    Kurt V. Mikkelsen of Theoretical Chemistry at the University of Copenhagen,
 +
                    who has a supercomputer capable of such demanding calculations. With a help from Andreas Erbs Hillers-Bendtsen, MSc student of chemistry at Copenhagen University, we were able to simulate the change in Gibbs free energy of
 +
                    the seven mutant G-alphas (see figure 7), and the nine fragments formed
 +
                    after cleavage (figure 8). Our goal was to identify a G-alpha mutant candidate
 +
                    with binding affinity similar to the wildtype, but with the lowest possible
 +
                    binding affinity of its fragments after cleavage. Our simulation lead to
 +
                    some unexpected findings:<br>
 +
                    <br>
 +
                    The majority of the fragments seemed to have a higher affinity towards
 +
                    G-beta than the non-cleaved G-alphas had. Only mutant m124 displayed
 +
                    a higher affinity than its corresponding fragments.
 +
                </div>
 +
                <div style="flex: 2; margin:3.5%;">
 +
                  <img style="width:100%;" src="https://static.igem.org/mediawiki/2020/b/ba/T--UCopenhagen--mutants-affinity.png">
 +
                    <figcaption>
 +
                      <b>Figure 7: Changes in Gibbs free energy of our G-alpha mutants relative to the wildtype.</b> Numbers in their labels signify what cleavage sites were included in each mutant. As all changes in energy are negative, it indicates that all mutants have a higher affinity to G-beta protein.
 +
                    </figcaption>
 +
                  </div>
 +
              </div>
 +
              <div style="display: flex; flex-direction:row;">
 +
                <div style="flex: 2; margin:3.5%;">
 +
                  <img style="width:100%;" src="https://static.igem.org/mediawiki/2020/f/fe/T--UCopenhagen--fragments-affinity.png">
 +
                    <figcaption>
 +
                      <b>Figure 8: Change in Gibbs free energy of fragments produced by cleavage of G-alpha mutants.</b> The x-axis denotes the length of each fragment. There seems to be a slight trend towards a higher affinity in smaller fragments, most likely caused by more degrees of freedom for a binding region constituting a majority of such fragments.
 +
                    </figcaption>
 +
                  </div>
 +
                  <div style="display: flex; font-family: Avenir, Arial; font-size: 16px; flex: 2.5; align-items: center;">
 +
                    Initially, we were startled by these findings, however after thorough
 +
                    discussions we concluded that the overall G-alpha conformation probably
 +
                    restricted some degrees of freedom of the high-affinity regions.
 +
                    These regions facilitated most interaction between G-alpha and G-beta,
 +
                    as we learned in one of the earlier stages. Our cleavage sites were
 +
                    located between these high-affinity regions, therefore, cleavage
 +
                    effectively freed the high-affinity regions from the previous restrictive
 +
                    conformation and allowed them to bind G-beta more tightly. Obviously,
 +
                    these findings turned our initial conclusions about the protein-protein
 +
                    interaction on its head. We also realized that affinities do not
 +
                    necessarily provide the full answer to successful signal transduction.
 +
                    In accordance with this, further analysis showed that conformational changes
 +
                    of G-beta and their causes could not answer our questions.<br> <br>
 +
                    Since mutant m124 was the only mutant whose corresponding fragments had
 +
                    a lower affinity to G-beta compared to the mutant itself, we decided
 +
                    it was the most promising candidate for our biosensor based on the data
 +
                    from the simulation.
 +
                </div>
 +
              </div>
 +
              <h2>Stage Seven - Test & Learn</h2>
 +
              Based on the conclusion from the previous Rosetta simulation, we designed
 +
              and ran another simulation to obtain more refined details about the
 +
              m124 fragments (see <a href="https://2020.igem.org/Team:UCopenhagen/Model">modeling</a> for more details). The results showed that all generated fragments had
 +
              a significantly higher affinity to the G-beta compared to m124. Although
 +
              exploration of the PDB files revealed that spatial assignment of the
 +
              fragments showed clashes (figure 9) so that the affinity of each fragment
 +
              could be limited by the presence of the other fragments, ultimately
 +
              leading to diminished presence at G-beta. Nevertheless, the results of
 +
              our simulations could not inform us about conformational changes of
 +
              G-beta which protected us from formulating clear conclusions.
 +
              <div style="display: flex; flex-direction: row; max-width: 95%; margin: 0 auto; justify-content: space-between; align-items: center;">
 +
                <div style="display: flex; flex-direction: column; justify-content: space-between; max-width: 45%;">
 +
                <img style="width: 100%; align-self: center;" src="https://static.igem.org/mediawiki/2020/9/99/T--UCopenhagen--m124.png">
 +
                  <figcaption>
 +
                    <b>Figure 9: Change in Gibbs free energy of fragments produced by the cleavage of G-alpha mutant m124.</b><br>The labels on the x-axis signify the position of each fragment in the original protein.
 +
                  </figcaption></div>
 +
<div style="display: flex; flex-direction: column; justify-content: space-between; max-width: 45%;">
 +
                <img style="width: 100%; align-self: center;" src="https://static.igem.org/mediawiki/2020/5/5d/T--UCopenhagen--m124-clash.png">
 +
                <figcaption>
 +
                  <b>Figure 10: An example of clashing positions of the fragments.</b> This indicates a situation that cannot happen in the real-world conditions, thus reducing our trust in the prediction.
 +
                </figcaption></div>
 +
              </div>
 +
              <h2>Wet Lab</h2>
 +
              As part of our wet lab, we created four different mutant constructs of GPA1,
 +
              named GPA1m1, GPA1m2, GPA1m3, and GPA1m4. Mutants 1, 2 and 4 correspond to
 +
              the m1, m2, and m4 denotations used in our dry lab modeling. The Wt
 +
              GPA1 was replaced with the different GPA1 mutants in an adenosine biosensor
 +
              strain of <em>Saccharomyces Cerevisiae</em>, previously created by our supervisor Sotirios
 +
              Kampranis' lab. In addition, an inducible TEV protease was inserted into
 +
              the same locus of the <em>S. Cerevisiae</em> chromosome as the rest of
 +
              the biosensor constructs.<br>
 +
              <br>
 +
              As a result, we were able to test out both the maintenance of the original
 +
              functions of GPA1 in a GPCR receptor system, and to verify that the GPA1
 +
              mutants were cleaved by the induced TEV protease, as this would result
 +
              in higher concentrations of the reporter.<br>
 +
              INSERT FIGURE OF A2A.
 +
              <br>
 +
              The experiments showed that the amount of reporter was increased
 +
              drastically in all biosensors when the TEV protease was induced by
 +
              growing the cells in galactose. This included the biosensor with wt
 +
              GPA1 which made it hard to conclude whether some of the increased
 +
              signal was caused by cleavage of the GPA1 mutant or that the increase
 +
              in signaling came from other parts of the system.<br>
 +
              <br>
 +
              The experiments further showed a loss-of-function in the sense that
 +
              none of the biosensors with mutated GPA1 proteins were able to respond
 +
              to different levels of adenosine whereas we were able to reproduce
 +
              the A2A adenosine biosensor created in the Sotirios Kampranis lab when
 +
              we used the wt GPA1 and grew the cells in glucose media.
 +
              <h2>Future Vision</h2>
 +
              <h4 style="margin-top: 0%;">Dry lab</h4>
 +
              The next round of simulations could be based on molecular dynamics
 +
              for us to see whether the results of Rosetta could be trusted. Apart
 +
              from that, the focus could shift from affinities to conformations,
 +
              and particularly G-beta could be more explored.<br>
 +
              <h4>Wet lab</h4>
 +
              In the wet lab, the next step would be to understand why the signaling
 +
              also increases in the biosensor with wt GPA1 when the cells are grown
 +
              in galactose and raffinose media. First, we would create a strain of
 +
              the biosensor without the inducible TEV protease to see whether the
 +
              spike in signaling is caused by the induction of the TEV protease or
 +
              something else related to the change of the media. If the TEV protease
 +
              does not seem to be the cause we would investigate the promoters that
 +
              we use, some of these are glycolytic promoters like pPGK1 and the
 +
              different media might cause slightly different expression levels
 +
              even though the promoters are normally classified as constituent.
 +
              To further investigate the loss-of-function of the mutants, we would
 +
              like to create strains of the biosensor with the mutated GPA1 proteins
 +
              but without the inducible TEV promoter. This would enable us to see
 +
              whether leaky expression of the TEV protease in the glucose fed cells
 +
              might have caused the loss-of-function.
 +
            </div>
 +
      </div>
 +
    </div>
 +
    <!-- inserting our footnotes -->
 +
    <div class="ref_cont">
 +
      <div class="refs">
 +
        <ol>
 +
          <li id="mutant">Douglas P. Gladue, James B. Konopka, Scanning mutagenesis of regions in the Gα protein Gpa1 that are predicted to interact with yeast mating pheromone receptors, FEMS Yeast Research, Volume 8, Issue 1, February 2008, Pages 71–80, https://doi.org/10.1111/j.1567-1364.2007.00311.x<a href="#mutant-ref" aria-label="Back to content">↩</a></li>
 +
        </ol>
 +
      </div>
 +
    </div>
 +
  </body>
 
   <script>
 
   <script>
 
     $(document).ready(function(){
 
     $(document).ready(function(){
Line 395: Line 426:
  
  
 
+
{{Template: UCopenhagen/Footer_Bubbly_sponsors}}
{{UCopenhagen/Footer_Bubbly}}
+

Latest revision as of 02:21, 28 October 2020

Introduction

This page documents our success with iterative engineering of novel G protein alpha subunits in yeast Saccharomyces cerevisiae. Protein modeling in Rosetta took us most of the way, finally accompanied by wet lab experiments. Our efforts followed the traditional engineering design cycle (see figure 1), with multiple phases of research followed by testing of our models and continuously refining them. Individual steps are noted chronologically as they happened, thus proceeding from relatively simple to more advanced and specialized methods. Finally, the results challenged our initial assumptions about protein interactions.
Fig. 1: Expanded engineering cycle.
We closely followed the expanded engineering cycle in our engineering process.


As described on our project design page, our most ambitious design called for modification of the alpha subunit of endogenous yeast G protein. The challenge here was that we wanted to maintain its natural functionality while expanding it to serving our purposes.
This meant we needed to design a G-alpha that:
  • A) has a natural-like affinity for the beta/gamma complex of the yeast G protein to stay bound to the beta/gamma complex when no signal was present, thus inhibiting leaky activation of the pathway.
  • B) has a low enough affinity for the beta/gamma complex after cleavage by our TEV-protease to enable dissociation from the beta/gamma complex.

Our work on G-alpha hinged on finding a mutant that preserved natural function when there was no signal (i.e. no cytokine), and also dissociated appropriately in the presence of a signal.

This page presents all stages of G-alpha engineering in a chronological order. The progress followed the engineering cycle (design-build-test-learn) as we iterated our designs several times. As we utilized primarily dry-lab experiments, the whole cycle often spans over each stage. Every stage was informed by all the previous ones.

Stage One - Design

We based our design choices on knowledge about the interaction interface of the alpha and beta subunits in bovine G protein, as we found no crystal structure of yeast G protein. Immediately, we noticed that the yeast G protein contained an ‘insertion sequence’ (highlighted in figure 2) that other G proteins did not. This indicated that the sequence was not essential for the general functionality of the G-alpha subunit. As such, we chose to search this sequence for a mutagenesis site.

Based on literature review, we identified the residue sequence that the TEV-protease cleaves: EXLYΦQφ, where X is any residue, Φ is any large or medium hydrophobic residue, and φ is any small hydrophobic or polar residue. For the sake of interpretability of experiments, however, it was necessary to introduce a single sequence type at all times, as the TEV protease may display different proteolytic rates depending on the recognition site (Dougherty et al., 1989). Following extensive literature search, we opted for a widely used sequence: ENLYFQG.
Figure 2: Amino acid sequence of the S. cerevisiae G protein alpha subunit.
The highlighted sequence was initially selected as an optimal region for mutagenesis.

Stage Two - Design-build

We received a file with the empirically observed binding interface of bovine G-alpha and G-beta through our supervisors. We could translate known characteristics of the bovine G protein into what we expected to see in the yeast G protein because of the similarity of their sequences. Based on the biophysical properties of different amino acids, we identified what we believed were the most essential residues for the binding between the alpha and beta subunits (see figure 3).
Figure 3: Crystal structure of the G protein complex.
The red-colored residues are the ones we identified as the most important for binding of the alpha and beta subunits to each other.

Stage Three - Test & Learn

Based on suggestions from the Sinisens (Aalto-Helsinki iGEM2020 team), we chose to use Rosetta for our protein engineering. Rosetta is a powerful tool for mapping and characterizing protein structures and interactions, and, we theorized, good for mutagenetic analyses of our G alpha protein.
Using Rosetta, we supplemented our own analysis with an alanine scan, which provided us with quantitative predictions of the role of selected residues in the G-alpha in terms of binding to G-beta. Rosetta's alanine scan is a script that calculates how the Gibb's free energy of two proteins change upon changing a single amino acid residue to an alanine (see modeling for more details).

The results of the alanine scan confirmed our intuitions to a good extent, and offered several valuable insight points that led to the refinement of our G-alpha designs. As illustrated on figure 4, residue number 227 was predicted as the single most important residue in the protein-protein interaction between the alpha and beta subunits. This result, among mutagenetic data found in literature about what the most functionally important amino acid residues in G-alpha are , guided us in discerning between which ones we could change, and which ones were critical for the correct functionality of G-alpha.
Figure 4: Change in Gibbs free energy (in kcal/mol) generated by replacing the given residues with alanine. As Gibbs free energy correlates with a change in affinity, it appears that residue 227 is the most important one for the binding interface of G-alpha and G-beta proteins.

Stage Four - Redesigning, Building, Testing & Learning

By compiling the information we had received so far, we isolated fitting mutation points in G alpha, where we
  1. 1. Avoid changing residues that lead to loss of function in G alpha (Gladue, D. P. 2008).
  2. 2. Avoid changing residues identified as very important through Rosetta's alanine scan.
  3. 3. Change residues located in the insertion sequence and other non-conserved residues, by comparing the yeast G-alpha to the bovine G-alpha.
  4. 4. Change residues located in loops based on visual analysis in Pymol.
Figure 5: Localization prediction of modified alpha subunits.
Cello was used to predict the localization of our three GPA1-versions. For all three versions, an increased cytoplasmic reliability score was observed compared to the wildtype.
With this information in mind, we designed three version of G-alpha, each with varying degrees of modifications. Using SwissProt protein modeling, we then confirmed that none of our edits lead to unexpected conformational changes. In addition, we used Cello to predict protein localization to the plasma membrane, as this is crucial for the functionality of our biosensor. As shown on figure 5, our designs yielded even better localization scores than the wildtype.
As shown on the figure, some cytoplasmic localization was also predicted, but this was desirable as certain enzymes in the cytoplasm are known to modify G-alpha for it to go to the membrane.

Stage Five - Design, Build, Test & Learn

Here we faced the reality of laboratory experimentation. Three G-alpha mutants with many diverse mutations could lead to lots of ambiguity as the laboratory results could be hard to interpret due to inability to monitor what mutation had a profound impact. Time constraints in the wet lab motivated us to limit the number of possible mutation points. We eliminated point mutations and kept only cleavage sites for TEV protease. To improve the quality of future experiments, we expanded the range of G-alpha candidates from three to seven by introducing all combinations of the three selected mutations. We checked folding predictions for all seven novel G-alphas which resulted in a significant overlap (as visible in the figure). None of the proteins seemed likely to take on a faulty conformation according to the SwissProt analysis.
Figure 6: Alignment of predicted structures of our G-alpha mutants. Apart from the intrinsically disordered region on the left, the predicted structures overlap fairly well. This indicates that the conformation of none of the G-alpha mutants was severely impacted by their mutations.

Stage Six - Test & Learn

In order to predict the behavior of the fragments generated by cleavage of our mutant G-alpha by the TEV protease, we decided to run a simulation in Rosetta (see modeling for more details). However, as this was an extremely demanding computational analysis we could not perform ourselves, we reached out to Assoc. Prof. Kurt V. Mikkelsen of Theoretical Chemistry at the University of Copenhagen, who has a supercomputer capable of such demanding calculations. With a help from Andreas Erbs Hillers-Bendtsen, MSc student of chemistry at Copenhagen University, we were able to simulate the change in Gibbs free energy of the seven mutant G-alphas (see figure 7), and the nine fragments formed after cleavage (figure 8). Our goal was to identify a G-alpha mutant candidate with binding affinity similar to the wildtype, but with the lowest possible binding affinity of its fragments after cleavage. Our simulation lead to some unexpected findings:

The majority of the fragments seemed to have a higher affinity towards G-beta than the non-cleaved G-alphas had. Only mutant m124 displayed a higher affinity than its corresponding fragments.
Figure 7: Changes in Gibbs free energy of our G-alpha mutants relative to the wildtype. Numbers in their labels signify what cleavage sites were included in each mutant. As all changes in energy are negative, it indicates that all mutants have a higher affinity to G-beta protein.
Figure 8: Change in Gibbs free energy of fragments produced by cleavage of G-alpha mutants. The x-axis denotes the length of each fragment. There seems to be a slight trend towards a higher affinity in smaller fragments, most likely caused by more degrees of freedom for a binding region constituting a majority of such fragments.
Initially, we were startled by these findings, however after thorough discussions we concluded that the overall G-alpha conformation probably restricted some degrees of freedom of the high-affinity regions. These regions facilitated most interaction between G-alpha and G-beta, as we learned in one of the earlier stages. Our cleavage sites were located between these high-affinity regions, therefore, cleavage effectively freed the high-affinity regions from the previous restrictive conformation and allowed them to bind G-beta more tightly. Obviously, these findings turned our initial conclusions about the protein-protein interaction on its head. We also realized that affinities do not necessarily provide the full answer to successful signal transduction. In accordance with this, further analysis showed that conformational changes of G-beta and their causes could not answer our questions.

Since mutant m124 was the only mutant whose corresponding fragments had a lower affinity to G-beta compared to the mutant itself, we decided it was the most promising candidate for our biosensor based on the data from the simulation.

Stage Seven - Test & Learn

Based on the conclusion from the previous Rosetta simulation, we designed and ran another simulation to obtain more refined details about the m124 fragments (see modeling for more details). The results showed that all generated fragments had a significantly higher affinity to the G-beta compared to m124. Although exploration of the PDB files revealed that spatial assignment of the fragments showed clashes (figure 9) so that the affinity of each fragment could be limited by the presence of the other fragments, ultimately leading to diminished presence at G-beta. Nevertheless, the results of our simulations could not inform us about conformational changes of G-beta which protected us from formulating clear conclusions.
Figure 9: Change in Gibbs free energy of fragments produced by the cleavage of G-alpha mutant m124.
The labels on the x-axis signify the position of each fragment in the original protein.
Figure 10: An example of clashing positions of the fragments. This indicates a situation that cannot happen in the real-world conditions, thus reducing our trust in the prediction.

Wet Lab

As part of our wet lab, we created four different mutant constructs of GPA1, named GPA1m1, GPA1m2, GPA1m3, and GPA1m4. Mutants 1, 2 and 4 correspond to the m1, m2, and m4 denotations used in our dry lab modeling. The Wt GPA1 was replaced with the different GPA1 mutants in an adenosine biosensor strain of Saccharomyces Cerevisiae, previously created by our supervisor Sotirios Kampranis' lab. In addition, an inducible TEV protease was inserted into the same locus of the S. Cerevisiae chromosome as the rest of the biosensor constructs.

As a result, we were able to test out both the maintenance of the original functions of GPA1 in a GPCR receptor system, and to verify that the GPA1 mutants were cleaved by the induced TEV protease, as this would result in higher concentrations of the reporter.
INSERT FIGURE OF A2A.
The experiments showed that the amount of reporter was increased drastically in all biosensors when the TEV protease was induced by growing the cells in galactose. This included the biosensor with wt GPA1 which made it hard to conclude whether some of the increased signal was caused by cleavage of the GPA1 mutant or that the increase in signaling came from other parts of the system.

The experiments further showed a loss-of-function in the sense that none of the biosensors with mutated GPA1 proteins were able to respond to different levels of adenosine whereas we were able to reproduce the A2A adenosine biosensor created in the Sotirios Kampranis lab when we used the wt GPA1 and grew the cells in glucose media.

Future Vision

Dry lab

The next round of simulations could be based on molecular dynamics for us to see whether the results of Rosetta could be trusted. Apart from that, the focus could shift from affinities to conformations, and particularly G-beta could be more explored.

Wet lab

In the wet lab, the next step would be to understand why the signaling also increases in the biosensor with wt GPA1 when the cells are grown in galactose and raffinose media. First, we would create a strain of the biosensor without the inducible TEV protease to see whether the spike in signaling is caused by the induction of the TEV protease or something else related to the change of the media. If the TEV protease does not seem to be the cause we would investigate the promoters that we use, some of these are glycolytic promoters like pPGK1 and the different media might cause slightly different expression levels even though the promoters are normally classified as constituent. To further investigate the loss-of-function of the mutants, we would like to create strains of the biosensor with the mutated GPA1 proteins but without the inducible TEV promoter. This would enable us to see whether leaky expression of the TEV protease in the glucose fed cells might have caused the loss-of-function.
  1. Douglas P. Gladue, James B. Konopka, Scanning mutagenesis of regions in the Gα protein Gpa1 that are predicted to interact with yeast mating pheromone receptors, FEMS Yeast Research, Volume 8, Issue 1, February 2008, Pages 71–80, https://doi.org/10.1111/j.1567-1364.2007.00311.x


Icons made by Freepik from www.flaticon.com