Difference between revisions of "Part:BBa K1330000"

Line 47: Line 47:
  
 
===2022 Team IISER-Tirupati_India's Contribution===
 
===2022 Team IISER-Tirupati_India's Contribution===
<h1>Contribution to BBa_K1330000:</h1>
 
<h2>&nbsp;Potential <a href="https://parts.igem.org/Part:BBa_K1330000">Spinach 2.1</a> mutants using Molecular Dynamic Simulations and MFE calculations with higher fluorescence intensity as compared to BBa_K1330000</h2>
 
<p>Towards the modification of parts from the Registry of Standard Biological Parts (RSBP), we aimed to design investigations that would enhance the fluorescence intensity of Spinach2.1, a light-up aptamer [17]. In order to achieve this, we altered the sequence by introducing mutations (single, double), inversions, and deletions at different regions of the Spinach2.1 aptamer sequence to increase its relative fluorescence intensity than the wild-type. We employed two approaches in which, firstly, the sequence in the tetraloop (UUCG or TTCG) was mutated and inversion near this tetraloop sequence [10, 12]. Secondly, point mutations (A12G and U86C) were made at different places in the Spinach2.1 sequence. These mutations were chosen by a thorough inspection of the Spinach2.1 3D structure with the help of Dr. Hussain Bhukya, IISER Tirupati, and Dr. Harikrishna S, Senior Scientist at Syngene.</p>
 
<div class = 'w3-container w3-center'>
 
  <img src = 'https://static.igem.wiki/teams/4438/wiki/contribution/modifications/parts-mod-fig1.jpg' width="100%" height="auto" class = 'w3-center'>
 
</div>
 
<p><strong>Fig.1 </strong>Aptamer sequence and its RNAfold structure of Spinach2 and Spinach2.1 from the parts registry made by DTU, Denmark. The sequences highlighted in bold are the aptamer sequence, and those flanking are the tRNA scaffold sequence, tRNALys3. The sequences displayed are the DNA sequences of the RNA aptamers (Spinach2 and Spinach2.1) as given in the RSBP. The tetraloop region is indicated by a box colored blue.</p>
 
<p>All possible single-point mutations (SPMs) were made in the tetraloop using combinatorics and inversions. The free energy of the thermodynamic ensembles, frequency of the minimum free energy (MFE) structure, and the ensemble diversity were predicted using RNAfold. Based on these predictions, we chose a few sequences that displayed the best parameters for the above predictions. The mutated sequences and their corresponding RNAfold structures are given below.&nbsp;</p>
 
<h3><strong>Mutations and Inversion: </strong>Inversion near tetraloop and SPM in the loop:</h3>
 
<p><em>Mutant 1</em> - <strong>Spinach2.2</strong></p>
 
<p>The free energy of the thermodynamic ensemble for the sequence below is predicted to be <strong>-34.44</strong> kcal/mol.</p>
 
<div class = 'w3-container w3-center'>
 
  <img src = 'https://static.igem.wiki/teams/4438/wiki/contribution/modifications/parts-mod-fig2.jpg' width="100%" height="auto" class = 'w3-center'>
 
</div>
 
<p><strong>Fig.2 </strong>Mutated aptamer sequence, Spinach2.2 and its RNAfold structure. The tetraloop is highlighted in red text and the inversions are underlined. The sequences highlighted in bold are the aptamer sequence and those flanking are the tRNA scaffold sequence, tRNALys3. The mutation and inversions are indicated in the RNAfold structures using blue circles.</p>
 
<p><em>Mutant 2</em> - <strong>Spinach2.3</strong></p>
 
<p>The free energy of the thermodynamic ensemble for the sequence below is predicted to be <strong>-34.44</strong> kcal/mol.</p>
 
<div class = 'w3-container w3-center'>
 
  <img src = 'https://static.igem.wiki/teams/4438/wiki/contribution/modifications/parts-mod-fig3.jpg' width="100%" height="auto" class = 'w3-center'>
 
</div>
 
<p><strong>Fig.3 </strong>Mutated aptamer sequence, Spinach2.3, and its RNAfold structure. The tetraloop is highlighted in red text, and the inversions are underlined. The sequences highlighted in bold are the aptamer sequence, and those flanking are the tRNA scaffold sequence, tRNALys3. The mutation and inversions are indicated in the RNAfold structures using blue circles.</p>
 
<p><em>Mutant 3</em> - <strong>Spinach2.4</strong></p>
 
<p>The free energy of the thermodynamic ensemble for the sequence below is predicted to be <strong>-34.80</strong> kcal/mol.</p>
 
<div class = 'w3-container w3-center'>
 
  <img src = 'https://static.igem.wiki/teams/4438/wiki/contribution/modifications/parts-mod-fig4.jpg' width="100%" height="auto" class = 'w3-center'>
 
</div>
 
<p><strong>Fig. 4 </strong>Mutated aptamer sequence, Spinach2.4 and its RNAfold structure. The tetraloop is highlighted in red text and the inversions are underlined. The sequences highlighted in bold are the aptamer sequence and those flanking are the tRNA scaffold sequence, tRNALys3. The mutation and inversions are indicated in the RNAfold structures using blue circles.</p>
 
<p><em>Mutant 4</em> - <strong>Spinach2.5</strong></p>
 
<p>The free energy of the thermodynamic ensemble for the sequence below is predicted to be <strong>-35.50</strong> kcal/mol.</p>
 
<div class = 'w3-container w3-center'>
 
  <img src = 'https://static.igem.wiki/teams/4438/wiki/contribution/modifications/parts-mod-fig5.jpg' width="100%" height="auto" class = 'w3-center'>
 
</div>
 
<p><br /><strong>Fig. 5 </strong>Mutated aptamer sequence, Spinach2.5 and its RNAfold structure. The tetraloop is highlighted in red text and the inversions are underlined. The sequences highlighted in bold are the aptamer sequence and those flanking are the tRNA scaffold sequence, tRNALys3. The mutation and inversions are indicated in the RNAfold structures using blue circles.</p>
 
  
 +
<h2>&nbsp;Potential Spinach 2.1 mutants using Molecular Dynamic Simulations and MFE calculations with higher fluorescence intensity as compared to BBa_K1330000</h2>
 +
<p>Towards the modification of parts from the Registry of Standard Biological Parts (RSBP), we aimed to design investigations that would enhance the fluorescence intensity of Spinach 2.1(BBa_K130000), a light-up aptamer. In order to achieve this, we altered the sequence by introducing mutations (single, double), inversions, and deletions at different regions of the Spinach2.1 aptamer sequence to increase its relative fluorescence intensity than the wild-type.</p>
 +
<p>All possible single-point mutations (SPMs) were made in the tetraloop using combinatorics and inversions. The free energy of the thermodynamic ensembles, frequency of the minimum free energy (MFE) structure, and the ensemble diversity were predicted using RNAfold. Based on these predictions, we chose a few sequences that displayed the best parameters for the above predictions. The mutated sequences and their corresponding RNAfold structures are given Figure 1.</p>
  
 +
<h3><strong>Mutations and Inversions: </strong></h3>
 +
<p>The free energy of the thermodynamic ensemble are predicted to be -34.44 kcal/mol, -34.44 kcal/mol, -34.80 kcal/mol and -35.50 kcal/mol for the sequences (modified Spinach2.1) of Spinach2.2(<partinfo>BBa_K4438200</partinfo>), Spinach2.3(<partinfo>BBa_K4438201</partinfo>), Spinach2.4(<partinfo>BBa_K4438202</partinfo>) and Spinach2.5(<partinfo>BBa_K4438203</partinfo>) respectively.</p>
  
 +
[[File:Parts-mod-Figure1.jpg|center]]
  
 +
<p>Using these sequences we ran GROMACS simulation to find the energy minimization results. These results were plotted using the steepest descent algorithm and it was observed that the native structure converged in almost 200 steps (Figure 4a) while for Spinach2.6 it converged in nearly 180 steps (Figure 4d). An NVT equilibration was performed and trajectory files were generated to visualise the simulation. It was found that the system reached equilibrium at less than 20 ps (Figure 4b, e). This equilibrated system was subjected to the Steer MD simulations to dislocate the ligand from the G-Quadruplex binding site with spring force 10000 KJ/mol/nm2, see the Figure 4c, d. The force applied was able to pull the ligand apart from the binding site of the RNA aptamer as seen in the animation below. Overall, the energy minimization plot suggests that the Spinach2.6 complex structure is relatively stable when compared with the native structure. Moreover, the Steer MD plots for both these complexes are significantly different suggesting the binding affinities of the ligand to RNA aptamer structure are different.</p>
  
<h2>B. In silico Analysis: Molecular Dynamics Simulation of RNA-ligand interactions to obtain energy</h2>
+
<div class = 'w3-content w3-center w3-padding-16'> <img class = 'w3-center' src = "https://static.igem.wiki/teams/4438/wiki/model/parts/final/parts-mod-figure1.jpg" width="100%" style = 'overflow-x:auto'>
<p>We modeled one PDB called 4TS2 [15] to show our idea is working. Because 4TS2 has maximum similarity with our desired part, Spinach 2.1 [17]. (90% match, verified through NCBI BLAST). We used Gromacs <strong>Umbrella Sampling MD </strong>[13] to obtain the free energy of the wild-type DFHBI-4TS2 complex. Later we replaced the 12th position A with G and 86th position U with C to generate a modified 4TS2 PDB file. We did the same Umbrella sampling for this mutated 4TS2 and compared the results between 4TS2 WT and 4TS2 Mutated.&nbsp;</p>
+
</div><br>
<p>The umbrella sampling method will tell the user about the binding energy between two species. The binding energy (&Delta;Gbind) is obtained from the potential of mean force (PMF), extracted from a series of umbrella sampling simulations.&nbsp;</p>
+
<p>Sequence and structures of the RNA Aptamers. The DNA sequence and the RNAfold structures of Spinach2 and Spinach2.1 are adopted from the parts registry made by DTU, Denmark. The Spinach2 sequences highlighted in bold are the aptamer sequence, and those flanking are the tRNA scaffold sequence, tRNALys3. The yellow colored box highlights the tetraloop sequence (colored red) and the bases adjacent (colored blue) to the tetraloop. The bases colored green are mutated from A/T and T/A to give a modified Spinach2 version, Spinach2.1. The other Spinach aptamer versions are derived from the Spinach2 sequence by base mutations in the tetraloop and the region adjacent to it (highlighted in the yellow colored boxes) and keeping the rest of the sequence unchanged. These mutants correspond to Spinach2.2, Spinach2.3, Spinach2.4 and Spinach2.5. The mutations in the DNA sequence are displayed on the RNAfold structures of the RNA aptamers (Spinach). The tetraloop region is indicated by a box colored blue for Spinach2 and the other mutations are circled in blue for other modified Spinach aptamers.</p>
<p><strong>4TS2 Wild Type : </strong>[16]</p>
+
<p><strong>ACGCGACCGAATGAAATGGTGAAGGACGGGTCCAGCCGGCTGCGCAGCCGGCTTGTTGAGTAGAGTGTGAGCTCCGTAACTGGTCGCGTC</strong></p>
+
<p><strong>GACGCGACCGAAUGAAAUGGUGAAGGACGGGUCCAGCCGGCUGCGCAGCCGGCUUGUUGAGUAGAGUGUGAGCUCCGUAACUGGUCGCGUC</strong></p>
+
<div class = 'w3-container w3-center w3-padding-16'>
+
  <img src = 'https://static.igem.wiki/teams/4438/wiki/contribution/modifications/wildtype.png' width="80%" height="auto" class = 'w3-center'>
+
</div>
+
<p><strong>4TS2 Mutated: </strong>[16]</p>
+
<p><strong>GACGCGACCGA</strong><strong>G</strong><strong>TGAAATGGTGAAGGACGGGTCCAGCCGGCTGCGCAGCCGGCTTGTTGAGTAGAGTGTGAGCTCCGTAAC</strong><strong>C</strong><strong>GGTCGCGTC</strong></p>
+
<p><strong>GACGCGACCGA</strong><strong>G</strong><strong>UGAAAUGGUGAAGGACGGGUCCAGCCGGCUGCGCAGCCGGCUUGUUGAGUAGAGUGUGAGCUCCGUAAC</strong><strong>C</strong><strong>GGUCGCGUC</strong></p>
+
<div class = 'w3-container w3-center w3-padding-16'>
+
  <img src = 'https://static.igem.wiki/teams/4438/wiki/contribution/modifications/mutated.png' width="80%" height="auto" class = 'w3-center'>
+
</div>
+
  
<p><strong>TECHNIQUE FOR SIMULATION OF 4TS2 WILD TYPE OR MUTATED PDB IN GROMACS</strong></p>
+
<p>Using these sequences we ran GROMACS simulation to find the energy minimization results. These results were plotted using the steepest descent algorithm and it was observed that the native structure converged in almost 200 steps (Figure 4a) while for Spinach2.6 it converged in nearly 180 steps (Figure 4d). An NVT equilibration was performed and trajectory files were generated to visualise the simulation. It was found that the system reached equilibrium at less than 20 ps (Figure 4b, e). This equilibrated system was subjected to the Steer MD simulations to dislocate the ligand from the G-Quadruplex binding site with spring force 10000 KJ/mol/nm2, see the Figure 4c, d. The force applied was able to pull the ligand apart from the binding site of the RNA aptamer as seen in the animation below. Overall, the energy minimization plot suggests that the Spinach2.6 complex structure is relatively stable when compared with the native structure. Moreover, the Steer MD plots for both these complexes are significantly different suggesting the binding affinities of the ligand to RNA aptamer structure are different.</p>
<p>We simulated the Spinach X-ray crystal structure (PDB entry: 4TS2) to validate the forcefield and other parameters for the MD as it was the closest (90% identity using NCBI BLAST) to our desired <a href="https://parts.igem.org/Part:BBa_K1330000">part</a> sequence. We employed Gromacs <strong>Umbrella Sampling MD</strong> to obtain the free energy of the wild-type DFHBI-4TS2 complex. Later, we mutated 12AG and U86C to generate the modified Spinach structure. The two structures (4TS2 and modified 4TS2) were then subjected to Umbrella Sampling using identical simulation parameters to compare the results.&nbsp;</p>
+
<p>Umbrella sampling method will tell the user about binding energy between two species. The binding energy (&Delta;Gbind) is obtained from the potential of mean force (PMF), extracted from a series of umbrella sampling simulations.</p>
+
  
 +
[[File:Parts-mod-Figure4.jpg|center]]
  
<h2><strong>RESULT: 4TS2 WILD TYPE PDB</strong></h2>
+
<div class = 'w3-content w3-center w3-padding-16'> <img class = 'w3-center' src = "https://static.igem.wiki/teams/4438/wiki/model/parts/final/parts-mod-figure1.jpg" width="100%" style = 'overflow-x:auto'>
<p>At first, we did energy minimization and plotted it as the steepest descent algorithm, which converged in almost 200 steps. Then we did NVT equilibration and generated a Trajectory (XTC) and .gro file; from it, we can visualize how the simulation is occurring. We also plotted temperature Vs. Time and found fluctuations are significantly less, proving our equilibration is good enough. We did Steer MD Simulation (SMD) to remove the ligand from the pocket of G-Quadruplex with force 10000 KJ/mol/nm, and while visualizing its trajectory file in PyMol, we found that it is unbinding. We have also plotted SMD results as Spring force vs. time and get an umbrella-like graph. Then we generated configurations and put a window to start umbrella sampling using a bash script get_distances.sh.Further simulation will generatepotential of Mean Force Vs. Reaction Coordinate graph, from there we could conclude the binding energy.</p>
+
</div><br>
<div class = 'w3-container w3-center w3-padding-16'>
+
<p>Results of MD simulations of native and modified RNA aptamers:a, b) Energy minimization and equilibration plots of native (4TS2) complexes respectively. c) Steer MD simulation plot for the native structure. d, e) Energy minimization and equilibration plots of modified Spinach RNA aptamer, Spinach2.6 complex structure respectively and f) Steer MD simulation plot for the Spinach2.6 complex structure</p>
  <img src = 'https://static.igem.wiki/teams/4438/wiki/contribution/modifications/energy-minimization-of-wild-4ts2.png' width="80%" height="auto" class = 'w3-center'>
+
</div>
+
<p><strong>Energy minimization for 4TS2-Wild type</strong></p>
+
<p><strong>( generated from potential.xvg file)</strong></p>
+
<div class = 'w3-container w3-center w3-padding-16'>
+
  <img src = 'https://static.igem.wiki/teams/4438/wiki/contribution/modifications/4ts2-wild-type-equillibration-graph.png' width="80%" height="auto" class = 'w3-center'>
+
</div>
+
<p><strong>Equilibration of 4TS2-Wild type</strong></p>
+
<p><strong>(Generated from temperature.xvg file)</strong><br /><br /></p>
+
<div class = 'w3-container w3-center w3-padding-16'>
+
  <img src = 'https://static.igem.wiki/teams/4438/wiki/contribution/modifications/4ts2-wild-type-smd.png' width="80%" height="auto" class = 'w3-center'>
+
</div>
+
<p><strong>Graph of SMD (Steered MD) of 4TS2-Wild type</strong></p>
+
<p>(<strong>Generated from pullf.xvg)</strong><br /><br /></p>
+
  
<p><strong>RESULT: 4TS2 Mutated PDB</strong></p>
 
 
<p>At first, we did Energy Minimization and plotted it as the steepest descent algorithm, and it converged in almost 177 steps. Then we did NVT equilibration and generated Trajectory(XTC) and .gro file; from it we can visualize how the simulation is occurring. We also plotted temperature Vs. Time and found fluctuations are significantly less, proving our equilibration is good enough. We did Steer MD Simulation(SMD) to remove the ligand from the pocket of G-Quadruplex with force 10000 KJ/mol/nm, and while visualizing its trajectory file in PyMol we found that it is unbinding. We have also plotted SMD results as Spring force vs. time and get an umbrella-like graph. Then we generated configurations and put a window to start umbrella sampling using a bash script get_distances.sh. Further simulation will generatepotential of Mean Force Vs. Reaction Coordinate graph, from there we could conclude the binding energy.</p>
 
<p><strong>Below graphs are showing results for our simulation.&nbsp;</strong></p>
 
<div class = 'w3-container w3-center w3-padding-16'>
 
  <img src = 'https://static.igem.wiki/teams/4438/wiki/contribution/modifications/energy-minimization-of-mutated-4ts2.png' width="80%" height="auto" class = 'w3-center'>
 
</div>
 
<p><strong>Energy minimization for 4TS2-Mutated (generated from potential.xvg file)</strong><br /><br /></p>
 
 
<div class = 'w3-container w3-center w3-padding-16'>
 
  <img src = 'https://static.igem.wiki/teams/4438/wiki/contribution/modifications/4ts2-mutated-equillibration.png' width="80%" height="auto" class = 'w3-center'>
 
</div>
 
<p><strong>Equilibration of 4TS2-Mutated</strong></p>
 
<p><strong>(Generated from temperature.xvg file)&nbsp;</strong></p>
 
 
<div class = 'w3-container w3-center w3-padding-16'>
 
  <img src = 'https://static.igem.wiki/teams/4438/wiki/contribution/modifications/4ts2-mutated-smd.png' width="80%" height="auto" class = 'w3-center'>
 
</div>
 
<p><strong>Graph of SMD (Steered MD) of 4TS2-Mutated&nbsp;</strong></p>
 
<p>(<strong>Generated from pullf.xvg)&nbsp;</strong></p>
 
<p><strong>FINAL RESULT</strong>: However we couldn&rsquo;t make it to PMF vs Reaction coordinate plot, energy minimization plot and other plots are indicative that 4TS2 mutated can have higher binding energy than the wild type&rsquo;</p>
 
<p>The 4TS2 Wild type converges at 200 steps after energy minimization, while 4TS2 Mutated converges at 177 steps, hence it is more stable and predicted to give higher fluorescence.</p>
 
  
  

Revision as of 00:27, 14 October 2022

Spinach2.1 flanked by tRNALys3

This BioBrick contains the gene coding for the Spinach2.1 RNA, flanked by a tRNA scaffold sequence.

Usage and Biology

Spinach 2.1 RNA can bind and activate the fluorophore [http://pubs.acs.org/doi/abs/10.1021/ja410819x DFHBI-1T], both in vivo and in vitro. RNA molecules can thus be fluorescently tagged, by gene fusion with this BioBrick. The flanking tRNA comes from human tRNALys3 and increases stability and folding efficiency of the RNA, leading to higher fluorescent signal. It is possible to detect signal from Spinach2.1 when expressed from moderate strength promoters, e.g. the Anderson promoter library.

This part does not contain or need a Ribosome Binding Site.

Characterization

Spinach2.1 is derived from the Spinach2 sequence described in [http://www.nature.com/nmeth/journal/v10/n12/full/nmeth.2701.html Strack et al. (2013)]. Spinach2.1 has had 2 nucleotides swapped to remove a SpeI site, present in Spinach2.

A folding-prediction is shown below, demonstrating that Spinach2.1 is expected adopt the same fold as Spinach2.

Spinach2.1 folding

To confirm that Spinach2.1 functions as well as Spinach2, fluorescence measurements were conducted on Spinach2 and Spinach2.1 RNA in vitro. Fluorescence was measured on samples with different concentrations of DFHBI-1T and excess RNA, to compare the brightness of the Spinach-DFHBI-1T complexes. The resulting fluorescence standard curves are shown below:

Spinach2.1 fluorescence standard curve.

The slopes of the curves represent the brightness of the RNA-DFHBI-1T complex. The slopes are similar, although Spinach2.1 might be moderately less bright. More experiments need to be conducted to establish whether the difference is significant.

To assess folding efficiency of Spinach2.1 fluorescence was measured on samples with excess DFHBI-1T. Fluorescence normalized to RNA concentration is shown below:

Spinach bar charts

As seen from the chart, there was no significant difference between fluorescence from the two RNA's. This means that equal fractions of RNA was bound to DFHBI-1T, indicating equal folding efficiency.

See [http://2014.igem.org/Team:DTU-Denmark/Achievements/Experimental_Results#lab-comparison-div our wiki] for more details.

2018 Team Hong_Kong-CUHK's Improvement

Spinach2.1 was originally constructed to prevent the illegal site in the superfolding Spinach2. However, its activity might be moderately less bright. On the other hand, Spinach2 has the melting temperature of ~38 degrees Celsius according to previous literature, which might not favor the measurement of heat-induced promoter activity. Comparing the 3 RNA reporters we constructed and the pre-existed Spinach2.1, iSpinach-D5 is moderately brighter.

Fluorescence data obtained from plate reader. DH5a transformed with pSB1C3-lpp-Spinach2.1 or other Spinach variants were incubated overnight at 37 degrees Celsius. Overnight culture was normalized to OD=0.2 and incubated in 200uM DFHBI for 45 minutes. N=1. More replicates are needed.

Collaborating with NUS-A team, we also observed that iSpinach-D5 is moderately more resistant to 37->45 degrees Celsius change than Spinach2.1. However, it was also more vulnerable to 37->30 degrees Celsius change.

(A) Fluorescence data obtained from plate reader. DH5a transformed with pSB1C3-lpp-Spinach2.1 or iSpinach-D5 were incubated overnight at 30, 37 or 45 degrees Celsius. Overnight culture was normalized to OD=0.2 and incubated in 200uM DFHBI for 45 minutes at their incubation temperatures. (B) Signal-to-noise ratio calculated by dividing fluorescence level of Spinach-expressing E. coli with that of untransformed E. coli. (C) Multiple t-tests across temperatures. 37 vs 45 P-value of iSpinach is larger than that of Spinach2.1. N=3.

2022 Team IISER-Tirupati_India's Contribution

 Potential Spinach 2.1 mutants using Molecular Dynamic Simulations and MFE calculations with higher fluorescence intensity as compared to BBa_K1330000

Towards the modification of parts from the Registry of Standard Biological Parts (RSBP), we aimed to design investigations that would enhance the fluorescence intensity of Spinach 2.1(BBa_K130000), a light-up aptamer. In order to achieve this, we altered the sequence by introducing mutations (single, double), inversions, and deletions at different regions of the Spinach2.1 aptamer sequence to increase its relative fluorescence intensity than the wild-type.

All possible single-point mutations (SPMs) were made in the tetraloop using combinatorics and inversions. The free energy of the thermodynamic ensembles, frequency of the minimum free energy (MFE) structure, and the ensemble diversity were predicted using RNAfold. Based on these predictions, we chose a few sequences that displayed the best parameters for the above predictions. The mutated sequences and their corresponding RNAfold structures are given Figure 1.

Mutations and Inversions:

The free energy of the thermodynamic ensemble are predicted to be -34.44 kcal/mol, -34.44 kcal/mol, -34.80 kcal/mol and -35.50 kcal/mol for the sequences (modified Spinach2.1) of Spinach2.2(BBa_K4438200), Spinach2.3(BBa_K4438201), Spinach2.4(BBa_K4438202) and Spinach2.5(BBa_K4438203) respectively.

Parts-mod-Figure1.jpg

Using these sequences we ran GROMACS simulation to find the energy minimization results. These results were plotted using the steepest descent algorithm and it was observed that the native structure converged in almost 200 steps (Figure 4a) while for Spinach2.6 it converged in nearly 180 steps (Figure 4d). An NVT equilibration was performed and trajectory files were generated to visualise the simulation. It was found that the system reached equilibrium at less than 20 ps (Figure 4b, e). This equilibrated system was subjected to the Steer MD simulations to dislocate the ligand from the G-Quadruplex binding site with spring force 10000 KJ/mol/nm2, see the Figure 4c, d. The force applied was able to pull the ligand apart from the binding site of the RNA aptamer as seen in the animation below. Overall, the energy minimization plot suggests that the Spinach2.6 complex structure is relatively stable when compared with the native structure. Moreover, the Steer MD plots for both these complexes are significantly different suggesting the binding affinities of the ligand to RNA aptamer structure are different.

<img class = 'w3-center' src = "parts-mod-figure1.jpg" width="100%" style = 'overflow-x:auto'>

Sequence and structures of the RNA Aptamers. The DNA sequence and the RNAfold structures of Spinach2 and Spinach2.1 are adopted from the parts registry made by DTU, Denmark. The Spinach2 sequences highlighted in bold are the aptamer sequence, and those flanking are the tRNA scaffold sequence, tRNALys3. The yellow colored box highlights the tetraloop sequence (colored red) and the bases adjacent (colored blue) to the tetraloop. The bases colored green are mutated from A/T and T/A to give a modified Spinach2 version, Spinach2.1. The other Spinach aptamer versions are derived from the Spinach2 sequence by base mutations in the tetraloop and the region adjacent to it (highlighted in the yellow colored boxes) and keeping the rest of the sequence unchanged. These mutants correspond to Spinach2.2, Spinach2.3, Spinach2.4 and Spinach2.5. The mutations in the DNA sequence are displayed on the RNAfold structures of the RNA aptamers (Spinach). The tetraloop region is indicated by a box colored blue for Spinach2 and the other mutations are circled in blue for other modified Spinach aptamers.

Using these sequences we ran GROMACS simulation to find the energy minimization results. These results were plotted using the steepest descent algorithm and it was observed that the native structure converged in almost 200 steps (Figure 4a) while for Spinach2.6 it converged in nearly 180 steps (Figure 4d). An NVT equilibration was performed and trajectory files were generated to visualise the simulation. It was found that the system reached equilibrium at less than 20 ps (Figure 4b, e). This equilibrated system was subjected to the Steer MD simulations to dislocate the ligand from the G-Quadruplex binding site with spring force 10000 KJ/mol/nm2, see the Figure 4c, d. The force applied was able to pull the ligand apart from the binding site of the RNA aptamer as seen in the animation below. Overall, the energy minimization plot suggests that the Spinach2.6 complex structure is relatively stable when compared with the native structure. Moreover, the Steer MD plots for both these complexes are significantly different suggesting the binding affinities of the ligand to RNA aptamer structure are different.

Parts-mod-Figure4.jpg
<img class = 'w3-center' src = "parts-mod-figure1.jpg" width="100%" style = 'overflow-x:auto'>

Results of MD simulations of native and modified RNA aptamers:a, b) Energy minimization and equilibration plots of native (4TS2) complexes respectively. c) Steer MD simulation plot for the native structure. d, e) Energy minimization and equilibration plots of modified Spinach RNA aptamer, Spinach2.6 complex structure respectively and f) Steer MD simulation plot for the Spinach2.6 complex structure


Sequence and Features


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