Difference between revisions of "Part:BBa K5396001"

 
(29 intermediate revisions by the same user not shown)
Line 1: Line 1:
  
__NOTOC__
+
__TOC__
 +
 
 
<partinfo>BBa_K5396001 short</partinfo>
 
<partinfo>BBa_K5396001 short</partinfo>
  
<p>BARBIE1 is a synthetic protein derived from BaCBM2 (BBa_K5396000) through a process of reverse engineering. It has the increased ability to bind to plastics when compared to BaCBM2.</p><p>The BARBIE1 protein is fused with the red fluorescent protein (RFP)[ ], which exhibits an excitation maximum at 558 nm and an emission maximum at 583 nm. This fusion enhances the visualization of BARBIE1 by fluorescence-based methods.</p>
+
Barbie1 is a synthetic protein derived from BaCBM2 through a process of reverse engineering. It has the increased ability to bind to plastics when compared to BaCBM2.
  
===Protein Design===
+
The Barbie1 protein is fused with the red fluorescent protein (miRFP670). This fusion enhances the visualization of Barbie1 by fluorescence-based methods.
 +
 
 +
This part was used as template to construct <partinfo>BBa_K5396004</partinfo>
 +
 
 +
= Usage and Biology =
 +
 
 +
===Barbie1 design===
 +
 
 +
Starting from the BaCBM2 structure model generated by the AlphaFold2 software, we performed docking assays with six types of plastic: polypropylene (PP),  polyethylene (PE), polyethylene terephthalate (PET), nylon (NY), polyvinyl chloride (PVC) and polystirene (PS). We made the docking using Gnina software with relaxed parameters to screen many proteins and features for plastic affinity.
 +
 
 +
Thereafter, the produced overlaps were removed by the docking assays using the ChimeraX software, as well as used for visualization and sequence manipulation. A reverse folding was then performed with the protein output from the docking using the LigandMPNN tool. The original protein set generated from the docking was filtered to maintain just unique positions, considering the associated score , without overlap between them.
 +
 
 +
By doing that, 6.000 sequences were generated for each ligand, totalizing 6 plastics x 6.000 sequences = 36.000 sequences, as illustrated in Figure 1. The consensus sequence from the 36000 sequence  originated our most optimized protein sequence sensitive for several plastics types was named as Barbie1!
  
<p>Starting from the BaCBM2 structure model generated by the AlphaFold2 software, we performed docking assays with six types of plastic: polypropylene (PP),  polyethylene (PE), polyethylene terephthalate (PET), nylon (NY), polyvinyl chloride (PVC) and polystirene (PS). We made the docking using Gnina software with relaxed parameters to screen many proteins and features for plastic affinity, which was calculated as shown below, where P stands for Protein and L for Ligand </p>
 
<p>PL \xrightarrow{} P + L \quad \xrightarrow{} \quad K_D = \frac{[P][L]}{[PL]} </p>
 
<p>
 
Thereafter, the produced overlaps were removed by the docking assays using the ChimeraX software, as well as used for visualization and sequence manipulation. A reverse folding was then performed with the protein output from the docking using the LigandMPNN tool. The original protein set generated from the docking was filtered to maintain just unique positions, considering the associated score , without overlap between them. </p>
 
<p>
 
By doing that, 6.000 sequences were generated for each ligand, totalizing 6 plastics x 6.000 sequences = 36.000 sequences, as illustrated in Figure 1. The consensus sequence from the 36000 sequence  originated our most optimized protein sequence sensitive for several plastics types named as BARBIE1! </p>
 
 
https://static.igem.wiki/teams/5396/registry/barbie-docking-mps.jpg
 
https://static.igem.wiki/teams/5396/registry/barbie-docking-mps.jpg
<p style="font-size: 11px;"><b>Figure 1.</b> Protein-ligand docking representation of the BARBIE1 protein docked with  PP, PE, PET, NY, PVC.</p>
 
<p>The BARBIE1 protein was redocked with the same plastics as before, once more using Gnina software. The result comparing its affinity with BaCBM2 in silico assays performed are described in the barplot of Figure 5. Comparing the predicted affinity between the original and the modified protein, it is notable a substantial increase for all plastics, in particular for PE, PP, and PS, highlighting the effectiveness of the processing pipeline.
 
</p>
 
<p>Besides the monomers tests, we also wanted to test the affinity using different sizes of plastic in order to guarantee that this could be a valuable parameter to future analysis and experiments. Therefore, the tested ligands were PE and PET plastics with 50 and 25 repeating units, respectively. As a result, the previous behavior at maintaining a higher KD for BARBIE1 when compared to BaCBM2 was preserved.</p>
 
  
=== Computational Modeling===
+
'''Figure 1.''' Protein-ligand docking representation of the Barbie1 protein docked with  PP, PE, PET, NY, PVC.
<p>Proteins with carbohydrate-binding modules (CBM) can not only be found as single units but also in more units or as part of larger multi-domain proteins. In light of the importance of understanding the thermodynamics basis for structural composition, the newest version of Alpha Fold 3 (AF3) was used [doi.org/10.1038/s41586-024-07487-w]. The model was used to optimize protein-protein interaction and test high-order oligomers. To effectively model the water filter system closer to reality, the first step was to predict the state of the proteins.</p>
+
 
 +
https://static.igem.wiki/teams/5396/registry/imagem-2024-09-20-141552930.png
 +
 
 +
'''Figure 2.''' AlphaFold 3 3D simulation of Barbie1 protein with miRFP and three Mad10 tags.
 +
 
 +
===miRFP670===
 +
 
 +
miRFP670 is a monomeric near-infrared fluorescent protein (NIR FP) derived from the bacterial photoreceptor ''Rhodopseudomonas palustris''. It exhibits excitation and emission maxima at 643 nm and 670 nm, respectively, making it suitable for applications in multicolor microscopy. With a molecular weight of approximately 34.5 kDa, miRFP670 is characterized by its high brightness and photostability in various cells, outperforming many other NIR FPs in terms of fluorescence intensity. The protein utilizes '''biliverdin''' as its chromophore, which enhances its effective brightness and allows for efficient labeling in live-cell imaging. The miRFP670 demonstrates low cytotoxicity and stability across a range of pH levels. Its structural properties enable the development of biosensors and facilitate protein-protein interaction studies. [https://www.fpbase.org/protein/mirfp670/]
 +
 
 +
===Mad10 peptide===
 +
 
 +
Mad10 (Magnetosome-associated deep-branching 10) is a protein derived from the magnetotactic bacterium ''Desulfamplus magnetovallimortis'' BW-1, known for its ability to synthesize magnetite crystals, which are essential for its navigation in aquatic environments. This protein is important for the formation and stabilization of magnetosomes—organelles that contain magnetic particles. [https://pubs.acs.org/doi/10.1021/acs.nanolett.9b03560]
 +
 
 +
The Mad10 tag refers to a peptide derived from Mad10, designed to include the most conserved amino acids within the protein sequence. The Mad10-derived peptide exhibits a strong affinity for magnetite, allowing it to effectively bind to magnetic columns during purification processes.
 +
 
 +
= Protein Expression and Purification =
 +
 
 +
The sequence for this part was chemically synthesized by Genscript and arrived pre-cloned in the pET-15b(+) vector. This vector includes a T7 promoter, lacI, a ribosome binding site (RBS), a 6xHis tag, and a T7 terminator. After the arrival, we proceeded with transforming the plasmid into the ''Escherichia coli'' BL21(DE3) strain.
 +
 
 +
==Optimization of protein expression==
 +
 
 +
We inoculated ''E. coli'' cells expressing Barbie1-RFP-3xMad10 into LB medium with ampicillin, incubating them at 37°C until an OD600 of 0.5 was reached. At this point, 1 mM of IPTG was added to each culture to induce protein expression, and the cultures were left shaking at 37°C for different durations: 3, 4, and 5 hours.
 +
 
 +
After harvesting the pellets, we extracted the proteins using a protein extraction buffer and sonication, followed by protein quantification using the BSA assay. The samples were then prepared for SDS-PAGE by denaturing them at 95°C and loaded into the gel for electrophoresis.
 +
 
 +
The electrophoresis results showed a clear and significant difference in protein expression when comparing the IPTG-induced samples (3, 4, and 5 hours) to the control at 0 hours, which had no visible expression bands. However, when comparing the 3, 4, and 5-hour induction times to each other, there was no substantial visual difference in the intensity of the bands, indicating that protein expression was already strong at the 3-hour mark and remained consistent through the later induction times.
 +
 
 +
https://static.igem.wiki/teams/5396/registry/sds-barbie1-mad10.png
 +
 
 +
'''Figure 3.''' Analysis of protein expression of Barbie1-RFP-3xMad10 over different induction times (0h, 3h, 4h, 5h) via SDS-PAGE. Lanes: 1- Ladder (molecular weight marker), 2- Barbie1-RFP-3xMad10 (0h), 3- Barbie1-RFP-3xMad10 (3h), 4- Barbie1-RFP-3xMad10 (4h), 5- Barbie1-RFP-3xMad10 (5h), 6- BLS1 (untransformed control).
 +
 
 +
==Final expression and purification==
 +
 
 +
During the sample preparation, we observed a high level of viscosity, which complicated the process. To address this, we centrifuged the sample multiple times and filtered it twice before loading it onto the equipment. Unfortunately, a significant portion of the sample was lost during this additional processing. Despite these challenges, the chromatogram still showed a clear peak corresponding to Barbie1-RFP-Mad10.
 +
 
 +
After purification, we sought to understand the cause of the excessive viscosity in our sample. To investigate, we performed a computational analysis using Aggrescan4D (A4D) to predict the aggregation propensities of the protein fold state. The results indicate that Barbie1-RFP-3xMad10 exhibits a higher aggregation propensity compared to BaCBM2-RFP-3xMad10. Figure 4 shows the scores for each residue in both proteins, where residues with negative values are considered "soluble residues", and those with positive values are classified as "aggregation-prone residues". The difference between the two proteins is particularly noticeable in the initial residues, corresponding to BaCBM2 and Barbie1.
 +
 
 +
https://static.igem.wiki/teams/5396/registry/imagem-2024-10-01-091300553.png
 +
 
 +
'''Figure 4.''' Scores for each residue in both proteins (BaCBM2-RFP-3xMad10 and Barbie1-RFP-3xMad10). Residues with negative values are identified as "soluble residues", while those with positive values are considered "aggregation-prone residues”.
 +
 
 +
Based on this finding, we decided to perform new rounds of expression and purification, this time incorporating detergents to improve protein solubilization. In subsequent cycles, we adapted our strategy to account for the aggregation tendency of Barbie1-RFP-Mad10, aiming to minimize this characteristic and improve the purification process.
 +
 
 +
https://static.igem.wiki/teams/5396/registry/chromatogram-barbie1-mad10.png
 +
 
 +
'''Figure 5.''' Chromatogram of Barbie1-RFP-Mad10 purification using IMAC (Immobilized Metal Affinity Chromatography) on a Ni-column. The peak corresponding to the BARBIE1-RFP-Mad10 protein is highlighted, indicating successful, although partial, purification.
 +
 
 +
=Characterization=
 +
==SEC-MALS==
 +
 
 +
After purifying Barbie1-RFP-3xMad10, SDS-PAGE analysis revealed that the sample was eluted along with other E. coli proteins. To better characterize the sample and improve its purity, a 30 kDa concentrator was used to separate the larger proteins from the smaller ones. With the sample more concentrated, we proceeded with SEC-MALS analysis. Unlike the results observed for BaCBM2-RFP-3xMad10, the SEC-MALS outcome for Barbie1-RFP-3xMad10 displayed three distinct peaks in the curve. However, none of the three peaks had an estimated molecular mass consistent with the expected mass of Barbie1-RFP-3xMad10 (46.47 kDa).
 +
 
 +
As noted earlier, during sample preparation for purification, we detected the formation of large aggregates. Additionally, the aggregation score for Barbie1-RFP-3xMad10, calculated by Aggrescan4D, is significantly higher than that of BaCBM2-RFP-3xMad10. Considering the SEC-MALS results along with the SDS-PAGE findings, which confirmed the presence of Barbie1-RFP-3xMad10 in the collected sample, we believe that the peak with a molecular mass of 157.7 kDa corresponds to a protein aggregate.
 +
 
 +
https://static.igem.wiki/teams/5396/registry/sec-mals-barbie1.png
 +
 
 +
'''Figure 6.''' SEC-MALS result of 6His-BARBIE1-RFP-3xMAD10.
 +
 
 +
==Circular Dichroism (CD)==
 +
 
 +
CD can provide crucial insights due to the presence of beta sheets in their structures. Typically, β-sheets exhibit a characteristic negative band around 218 nanometers (nm) and a positive band near 195 nm, while α-helices show a negative band around 222 nm and a positive band around 190 nm. However, a limitation of the equipment used restricts the analyzed wavelength range to between 200 and 250 nm.
 +
 
 +
The pipeline acquisition of the spectrum was followed by a 2 seconds of integration for each point in each 6 spectrum calculated at room temperature (20 ºC). From the collected data, it is calculated its average, subtracted from the water absorbance, normalized, and smoothed using a Savitzky-Golay filter.
 +
 
 +
https://static.igem.wiki/teams/5396/registry/imagem-2024-09-23-150429757.png
 +
 
 +
'''Figure 7.''' CD acquisition of Barbie1-RFP-3xMad10 protein at 20 ºC.
 +
 
 +
===Barbie1-RFP-3xMad10 Heating Ramp===
 +
 
 +
An essential evaluation for the protein is the heating ramp, in which we can evaluate how the protein acts in different temperatures and their melting point. For doing so, the same procedure was followed, in which each protein was heated from 20 to 95ºC, as shown in Figure 8. The data acquisition was made in steps of 5ºC.
 +
 
 +
The CD absorbance is shown as the figure’s color map. It is evident that the protein show absorbance in the secondary structure regions, indicating a thermal resistivity.
 +
 
 +
https://static.igem.wiki/teams/5396/registry/imagem-2024-09-23-151406532.png
 +
 
 +
'''Figura 8.''' Heating ramp for Barbie1-RFP-3xMad10.
 +
 
 +
===Barbie1-RFP-3xMad10 Pontual Spectra===
 +
 
 +
With this result, it is interesting to observe that, even after the protein engineering, the protein could still maintain great results for tests like the described. To further explore it, we also ran a protein pontual spectra, in which we measured CD absorbance for only 215 nanometers, but each measure made at 1ºC increase.
 +
 
 +
Figure 9 shows the result of the experiment, which emphasizes the decrease in the CD absorbance for Barbie1-RFP-3xMad10 as already expected from Figure 8. Although Barbie1-RFP-3xMad10 loses part of its absorbance with the increase of temperature, it is not possible to confirm its denaturation. A viable hypothesis is some reduction in the beta-sheets structures.
 +
 
 +
https://static.igem.wiki/teams/5396/registry/imagem-2024-09-23-151807308.png
 +
 
 +
'''Figure 9.''' Pontual spectra for Barbie1-RFP-3xMad10 at 215 nm.
 +
 
 +
In the first instant, this behavior might show a bad characteristic of Barbie1-RFP-3xMad10. However, the loss of function at high temperatures can be used to properly design the water filter. Once saturated with plastics, the system can be heated to high temperatures making Barbie1-RFP-3xMad10 lose part of its structure. Thus, knowing the function is correlated to the form, the protein might lose its plastic affinity, allowing the removal of microplastics and further processing for reutilizing it.
 +
 
 +
===Heating and Cooling Barbie1-RFP-3xMad10 Ramp===
 +
 
 +
With this possibility in our sights, the following experiment is heating and cooling the protein, in a way we can observe how its structures will behave. Illustrated in Figure 10, the plot was rotated for better visualizing the CD absorbance in function of temperature variation. From left to right, it is possible to firstly see how Barbie1-RFP-3xMad10behaves when heated and then cooled.
 +
 
 +
As a result, protein’s CD absorbance in the cooling ramp still presented a signal, indicating it maintained its structure and a consequently renovelation of the possible denoveled parts.
 +
 
 +
https://static.igem.wiki/teams/5396/registry/heating-cooling-barbie.png
 +
 
 +
'''Figure 10.''' Heating and cooling ramp CD absorbance of Barbie1-RFP-3xMad10.
 +
 
 +
===Specific Wavelength Analysis===
 +
 
 +
If we deepen our analysis into 210 and 222 nm, which are in the region of secondary structure absorbance, we can better understand the variancy. In Figure 11, the graphs are related to the 210 nm wavelength, whereas the graphs from below to the 222 nm.
 +
 
 +
In a dashed gray line, an average initial absorbance for Barbie1-RFP-3xMad10 was added to highlight the loss of the original absorbance. Additionally, both wavelengths show an interesting behavior, that is an important trend change between 55 and 60 ºC, which may be a critical point for the protein. From the cooling graphs, it is notable how the protein loses part of the original signal, indicating a conformational alteration from the temperature change.
 +
 
 +
For future adaptation of Barbie1-RFP-3xMad10, it is fundamental to assess and design the protein in a way it can be further reused for multiple filtration.
 +
 
 +
https://static.igem.wiki/teams/5396/registry/imagem-2024-09-23-160301522.png
 +
 
 +
'''Figure 11.''' Two specific wavelengths, 210 nm above and the 222 nm below, plots for the heating and cooling ramp for Barbie1-RFP-3xMad10.
 +
 
 +
==Barbie1-RFP-3xMad10 Melting Temperature==
 +
 
 +
From the previous result of analyzing a specific wavelength, it is possible to note the curve’s behavior, which is usual for a protein denaturation. Thus, a viable way to uncover the protein’s melting temperature is by fitting the Boltzmann curve in the collected data, which is given by:
 +
 
 +
'''$$A2 + (A1 - A2) / (1 + np.exp((x - x0) / dx)),$$'''
 +
 
 +
Where:
 +
*A1 is the protein’s native state;
 +
*A2 is the completely unfolded protein’s state;
 +
*Tm is the protein’s melting temperature;
 +
*k is a constant that represents the curve inclination.
 +
 
 +
After fitting the curve in Figure 12, the melting point found for Barbie1-RFP-3xMad10 is approximately 58.81ºC. This information is fundamental for building the filter, since it represents stability in a higher temperature. In particular, it is known that in water purifiers the temperatures do not exceed 30 to 35ºC, making the protein appropriate to the system.
 +
 
 +
https://static.igem.wiki/teams/5396/registry/fitting-curve-barbie.png
 +
 
 +
'''Figure 12.''' Boltzmann fit in the absorbance values of Barbie1-RFP-3xMad10 in a heating ramp.
 +
 
 +
==Plastic Interaction Analysis==
 +
 
 +
Through secondary structures absorbance, it was possible to observe that Barbie1-RFP-3xMad10 structure do not change in the presence of plastic nanoparticles
 +
 
 +
Finally, the last experimental evaluation at CD was with the polystyrene nanoplastics. We used a 100 nanometer particle which was inserted into a solution with Barbie1-RFP-3xMad10.
 +
 
 +
Nevertheless, as reported from the computational simulations, we do not expect any conformational change. Since CD is only evaluating the absorbance in secondary structures, there may not happen major differences in the experiment.
  
<p>The proposed proteins to be evaluated in Alpha Pulldown are the B1-CBM (pink), BaCBM2 (blue), and 1A3N (green), which is a reference protein that forms the deoxy human hemoglobin.</p>
+
Consequently, Figure 13 shows as projected. When the Barbie1-RFP-3xMad10a spectra is plotted with the nanoparticle solution spectra, minor changes can be observed. Specifically, this contrast might be a result of the polystyrene nanoparticle scattering, creating noise in the measurement.
https://static.igem.wiki/teams/5396/registry/barbie-oligomeric-state.png
+
<p style="font-size: 11px;"><b>Figure X.</b> Protein comparative of the tested proteins. In (a), it is shown the BARBIE1 structure, in (b) the BaCBM2, and in (c) the 1A3N protein.
+
</p>
+
<p>With 4 subunits as shown in Figure X+1, the 1A3N multicomplex protein forms a pocket in the middle region to store oxygen molecules. Since this reference protein was resolved with more subunits, it can serve as a baseline for the prediction of the other proteins.
+
</p>
+
https://static.igem.wiki/teams/5396/registry/1a3n.jpg
+
<p style="font-size: 11px;"><b>Figure X+1.</b> 1A3N protein representation with four repeating units.</p>
+
  
<p>About the predicted template modeling (pTM) values:</p>
+
https://static.igem.wiki/teams/5396/registry/barbie-np-cd.png
https://static.igem.wiki/teams/5396/registry/imagem-2024-09-09-161834616.png
+
<p>Concerning the interface predicted template modeling (iPTM) values:</p>
+
https://static.igem.wiki/teams/5396/registry/imagem-2024-09-09-162505500.png
+
https://static.igem.wiki/teams/5396/registry/imagem-2024-09-09-162714537.png
+
<p style="font-size: 11px;"><b>Figure X+2.</b> pTM and iPTM values calculated by AlphaFold3 for BARBIE1, BaCBM2, and 1A3N with different repeating units.
+
</p>
+
<p>As it can be seen in Figure X+2, both BARBIE1, BaCBM2, and 1A3N achieved high scores as monomers. On the other hand, only the 1A3N as a multimer achieved higher values, such as its two and three subunits iPTM, as expected. Therefore, it is possible to affirm that it is very unlikely for both B1-CBM and BaCBM2 to form a multimer, which can be assured by the deoxy human hemoglobin results. This result allows the modeling of B1-CBM as a monomer, which simplifies the system.
+
</p>
+
=== Interaction Properties===
+
<p>In order to further our knowledge in the protein properties, we compared our resulting sequence BARBIE1 with the original BaCBM2 in different tests. Firstly, we compared each protein's tertiary structure, as shown in Figure X+3. As a result, it is notable the similarity between them, specially in the secondary structures, in which is notable the presence of the same amount of beta sheets. In the upper left part of the BARBIE1 structure, however, there is a notable difference between them.
+
</p>
+
https://static.igem.wiki/teams/5396/registry/imagem-2024-09-09-163621557.png
+
<p style="font-size: 11px;"><b>Figure X+3.</b> Tertiary structure of BARBIE1 on the left (pink) and BaCBM2 on the right (blue).
+
</p>
+
<p>This subtable contrast between the original and modified protein is shown in Figure X+4, in which we aligned the structures using VMD and calculated the root mean square deviation (RMSD). As it is possible to analyze in the aligned structures, it is very similar visibly, with a resulting RMSD of 0.56 Å (RMSDs < 2 Å indicates high similarity).
+
</p>
+
https://static.igem.wiki/teams/5396/registry/barbie-alligned-w-cbm.png
+
<p style="font-size: 11px;"><b>Figure X+4.</b> Alignment of BARBIE1 (magenta) and BaCBM2 (blue) tertiary structures resulting in a RMSD of 0.56.</p>
+
<p>After that, we assessed the electrostatic surface of each protein using the ChimeraX tool, which can be valuable for identifying binding sites for ligands, as well as stability and its behavior when solvated.</p>
+
<p>In Figure X+5, both proteins are represented with the electrostatic surface. While on one hand the red parts stand for the negative electrostatic, the blue parts stand for the positive electrostatic. Therefore, it is notable a major presence of negative values for BARBIE1 when compared to BaCBM2, which may indicate a higher affinity with positive ions or positive charged ligands.</p>
+
https://static.igem.wiki/teams/5396/registry/electrostatic-barbie-cbm.png
+
<p style="font-size: 11px;"><b>Figure X+5.</b> Electrostatic surface represented in the left for BARBIE1 and in the right for BaCBM2. The red regions indicate positive electrostatic and the negative are represented as blue.</p>
+
<p>Following on, we calculated the proteins hydrophobicities also using ChimeraX. In this way, we can at the same time understand the water interacting regions, corresponding to the more hydrophilic parts, and possible binding sites, generally indicated by hydrophobic parts. In Figure X+6, the hydrophilic regions are represented as blue and the hydrophobic regions as yellow.</p>
+
https://static.igem.wiki/teams/5396/registry/hydrophobicity-barbie-cbm.png
+
<p style="font-size: 11px;"><b>Figure X+6.</b> Hydrophobicity surface represented in the left for BARBIE1 and in the right for BaCBM2. The blue regions indicate hydrophilic and the hydrophobicity is represented as yellow.</p>
+
<p>In general, while on one hand hydrophilic regions are exposed to aqueous solvents in the exterior of the structure, hydrophobic regions are buried in its interior. When confronting both structures, it is notable an inverse behavior on the BARBIE1 structure compared to BaCBM2: the hydrophobic parts are not buried, but exposed to the surface.
+
Since more hydrophobic sites are exposed to the surface, it may point to a better understanding of our results. The better scores were achieved by choosing more hydrophobic amino acids in the protein primary structure, which enabled the creation of more pockets and the subsequent increase of KD. Therefore, the choice is notable for a higher plastic affinity, but also a lesser water solubility.
+
</p>
+
===BARBIE1 molecular dynamics===
+
  
<p>Following on, with a deeper understanding of BARBIE1 protein, it is now possible to study its interaction with plastic molecules. As a proof of concept of the system, only the polystyrene (PS) is analyzed in the experiments, as in other parts of the project. To better comprehend the binding between BARBIE1 with PS, we used Charmm-GUI to create the solution system for running on Gromacs.</p>
+
'''Figure 13.''' Plastic binding proteins Barbie1-RFP-3xMad10and BaCBM2-RFP-3xMad10 spectra acquisition with 100 nm polystyrene nanoparticles.
<p>Considering that building a microplastic particle is a very expensive task in computation terms, we used a plastic molecule instead. In most of the tests, the PS molecule was built with 10 repeating units, as shown in Figure X+7. For this kind of plastic in specific, there are three principal taticities formation, that is the extent to where the functional group is.
+
</p>
+
<p>For this particular case, polystyrene is composed of a hydrocarbon chain with phenyl groups, making it possible the creation of the particular taticites: atactic, syndiotactic, and isotactic. Since atactic is the only important for commercial purposes, as consequence the majority of PS microplastics must be atactic, which justified the creation of the ligand as atactic.
+
</p>
+
https://static.igem.wiki/teams/5396/registry/ps.png
+
<p style="font-size: 11px;"><b>Figure X+7.</b> On the left, it is represented a single molecule of styrene. On the right, it is represented a molecule with 10 repeating units.
+
</p>
+
<p>The created system simulated the one used in the wet lab experiment, which is in ambient temperature (300 K), neutral pH (pH=7), and 0.15 mM as salt concentration of NaCl. As represented in Figure X+8, the system is colored with its solvent as blue, BARBIE1 as magenta, the styrene molecule in gray and white, and the ions in green. The system configuration is setted with periodic boundary conditions (PBCs) with an edge distance of 15 Å.
+
</p>
+
https://static.igem.wiki/teams/5396/registry/barbie-ps.png
+
<p style="font-size: 11px;"><b>Figure X+8.</b> Representation of the solvated system of BARBIE1 with polystyrene.
+
</p>
+
<p>Aiming to keep on the comparison between BARBIE1 and BaCBM2, another similar system was generated for BaCBM2. It is fundamental to note that since calculating MDs for large molecules, such as proteins, it is not viable using Quantum Mechanics. For this reason, we are not able to check the actual bind between protein and ligand, but it is possible to check its energies and behavior through classical dynamics.</p>
+
<p>With the calculated systems, it was notable how rapidly both proteins get close to the plastic as a result of their affinity for them. For quantifying this proximity, we calculated the minimum distance between plastic and protein for each simulation, as shown in Figure X+9. It is relevant to realize that the peak generated in BaCBM2 simulation is given by a periodic boundary condition - i.e. the plastic is getting close to the borders and traveling to the other side of the box.
+
</p>
+
https://static.igem.wiki/teams/5396/registry/minimum-distances.png
+
<p style="font-size: 11px;"><b>Figure X+9.</b> Minimum distances between each protein and the polystyrene molecule.
+
</p>
+
<p>As it is possible to see in the lower part of the graph, the proximity between each component of the simulations are very stable, with a minimum distance average of 2.17 Å for BARBIE1. This is a very interesting result, since a hydrogen bond has an average distance between 2.6 and 3.3 Å, which is very relevant for understanding the interaction type we expect.</p>
+
  
<p>A possible complementary analysis to the minimum distance is the root mean square deviation, which calculates the deviation trajectory through the simulation. Essentially, it allows us to understand a molecule or atom's stability through time, which can be useful for protein-ligand interaction. In Figure X+10, we calculated the RMSD for polystyrene molecule.</p>
+
Considering we were not able to assess the plastic-protein binding, further studies might be done. One important form to evaluate this is the infrared spectroscopy, which allows studying the vibrational frequencies of molecules.
https://static.igem.wiki/teams/5396/registry/ps-trajectory.png
+
<p style="font-size: 11px;"><b>Figure X+10.</b> Root mean square deviation for polystyrene trajectory.
+
</p>
+
<p></p>
+
<p></p>
+
  
 
<span class='h3bb'>Sequence and Features</span>
 
<span class='h3bb'>Sequence and Features</span>

Latest revision as of 20:28, 1 October 2024

Barbie1_RFP_3xMad10

Barbie1 is a synthetic protein derived from BaCBM2 through a process of reverse engineering. It has the increased ability to bind to plastics when compared to BaCBM2.

The Barbie1 protein is fused with the red fluorescent protein (miRFP670). This fusion enhances the visualization of Barbie1 by fluorescence-based methods.

This part was used as template to construct BBa_K5396004

Usage and Biology

Barbie1 design

Starting from the BaCBM2 structure model generated by the AlphaFold2 software, we performed docking assays with six types of plastic: polypropylene (PP), polyethylene (PE), polyethylene terephthalate (PET), nylon (NY), polyvinyl chloride (PVC) and polystirene (PS). We made the docking using Gnina software with relaxed parameters to screen many proteins and features for plastic affinity.

Thereafter, the produced overlaps were removed by the docking assays using the ChimeraX software, as well as used for visualization and sequence manipulation. A reverse folding was then performed with the protein output from the docking using the LigandMPNN tool. The original protein set generated from the docking was filtered to maintain just unique positions, considering the associated score , without overlap between them.

By doing that, 6.000 sequences were generated for each ligand, totalizing 6 plastics x 6.000 sequences = 36.000 sequences, as illustrated in Figure 1. The consensus sequence from the 36000 sequence originated our most optimized protein sequence sensitive for several plastics types was named as Barbie1!

barbie-docking-mps.jpg

Figure 1. Protein-ligand docking representation of the Barbie1 protein docked with PP, PE, PET, NY, PVC.

imagem-2024-09-20-141552930.png

Figure 2. AlphaFold 3 3D simulation of Barbie1 protein with miRFP and three Mad10 tags.

miRFP670

miRFP670 is a monomeric near-infrared fluorescent protein (NIR FP) derived from the bacterial photoreceptor Rhodopseudomonas palustris. It exhibits excitation and emission maxima at 643 nm and 670 nm, respectively, making it suitable for applications in multicolor microscopy. With a molecular weight of approximately 34.5 kDa, miRFP670 is characterized by its high brightness and photostability in various cells, outperforming many other NIR FPs in terms of fluorescence intensity. The protein utilizes biliverdin as its chromophore, which enhances its effective brightness and allows for efficient labeling in live-cell imaging. The miRFP670 demonstrates low cytotoxicity and stability across a range of pH levels. Its structural properties enable the development of biosensors and facilitate protein-protein interaction studies. [1]

Mad10 peptide

Mad10 (Magnetosome-associated deep-branching 10) is a protein derived from the magnetotactic bacterium Desulfamplus magnetovallimortis BW-1, known for its ability to synthesize magnetite crystals, which are essential for its navigation in aquatic environments. This protein is important for the formation and stabilization of magnetosomes—organelles that contain magnetic particles. [2]

The Mad10 tag refers to a peptide derived from Mad10, designed to include the most conserved amino acids within the protein sequence. The Mad10-derived peptide exhibits a strong affinity for magnetite, allowing it to effectively bind to magnetic columns during purification processes.

Protein Expression and Purification

The sequence for this part was chemically synthesized by Genscript and arrived pre-cloned in the pET-15b(+) vector. This vector includes a T7 promoter, lacI, a ribosome binding site (RBS), a 6xHis tag, and a T7 terminator. After the arrival, we proceeded with transforming the plasmid into the Escherichia coli BL21(DE3) strain.

Optimization of protein expression

We inoculated E. coli cells expressing Barbie1-RFP-3xMad10 into LB medium with ampicillin, incubating them at 37°C until an OD600 of 0.5 was reached. At this point, 1 mM of IPTG was added to each culture to induce protein expression, and the cultures were left shaking at 37°C for different durations: 3, 4, and 5 hours.

After harvesting the pellets, we extracted the proteins using a protein extraction buffer and sonication, followed by protein quantification using the BSA assay. The samples were then prepared for SDS-PAGE by denaturing them at 95°C and loaded into the gel for electrophoresis.

The electrophoresis results showed a clear and significant difference in protein expression when comparing the IPTG-induced samples (3, 4, and 5 hours) to the control at 0 hours, which had no visible expression bands. However, when comparing the 3, 4, and 5-hour induction times to each other, there was no substantial visual difference in the intensity of the bands, indicating that protein expression was already strong at the 3-hour mark and remained consistent through the later induction times.

sds-barbie1-mad10.png

Figure 3. Analysis of protein expression of Barbie1-RFP-3xMad10 over different induction times (0h, 3h, 4h, 5h) via SDS-PAGE. Lanes: 1- Ladder (molecular weight marker), 2- Barbie1-RFP-3xMad10 (0h), 3- Barbie1-RFP-3xMad10 (3h), 4- Barbie1-RFP-3xMad10 (4h), 5- Barbie1-RFP-3xMad10 (5h), 6- BLS1 (untransformed control).

Final expression and purification

During the sample preparation, we observed a high level of viscosity, which complicated the process. To address this, we centrifuged the sample multiple times and filtered it twice before loading it onto the equipment. Unfortunately, a significant portion of the sample was lost during this additional processing. Despite these challenges, the chromatogram still showed a clear peak corresponding to Barbie1-RFP-Mad10.

After purification, we sought to understand the cause of the excessive viscosity in our sample. To investigate, we performed a computational analysis using Aggrescan4D (A4D) to predict the aggregation propensities of the protein fold state. The results indicate that Barbie1-RFP-3xMad10 exhibits a higher aggregation propensity compared to BaCBM2-RFP-3xMad10. Figure 4 shows the scores for each residue in both proteins, where residues with negative values are considered "soluble residues", and those with positive values are classified as "aggregation-prone residues". The difference between the two proteins is particularly noticeable in the initial residues, corresponding to BaCBM2 and Barbie1.

imagem-2024-10-01-091300553.png

Figure 4. Scores for each residue in both proteins (BaCBM2-RFP-3xMad10 and Barbie1-RFP-3xMad10). Residues with negative values are identified as "soluble residues", while those with positive values are considered "aggregation-prone residues”.

Based on this finding, we decided to perform new rounds of expression and purification, this time incorporating detergents to improve protein solubilization. In subsequent cycles, we adapted our strategy to account for the aggregation tendency of Barbie1-RFP-Mad10, aiming to minimize this characteristic and improve the purification process.

chromatogram-barbie1-mad10.png

Figure 5. Chromatogram of Barbie1-RFP-Mad10 purification using IMAC (Immobilized Metal Affinity Chromatography) on a Ni-column. The peak corresponding to the BARBIE1-RFP-Mad10 protein is highlighted, indicating successful, although partial, purification.

Characterization

SEC-MALS

After purifying Barbie1-RFP-3xMad10, SDS-PAGE analysis revealed that the sample was eluted along with other E. coli proteins. To better characterize the sample and improve its purity, a 30 kDa concentrator was used to separate the larger proteins from the smaller ones. With the sample more concentrated, we proceeded with SEC-MALS analysis. Unlike the results observed for BaCBM2-RFP-3xMad10, the SEC-MALS outcome for Barbie1-RFP-3xMad10 displayed three distinct peaks in the curve. However, none of the three peaks had an estimated molecular mass consistent with the expected mass of Barbie1-RFP-3xMad10 (46.47 kDa).

As noted earlier, during sample preparation for purification, we detected the formation of large aggregates. Additionally, the aggregation score for Barbie1-RFP-3xMad10, calculated by Aggrescan4D, is significantly higher than that of BaCBM2-RFP-3xMad10. Considering the SEC-MALS results along with the SDS-PAGE findings, which confirmed the presence of Barbie1-RFP-3xMad10 in the collected sample, we believe that the peak with a molecular mass of 157.7 kDa corresponds to a protein aggregate.

sec-mals-barbie1.png

Figure 6. SEC-MALS result of 6His-BARBIE1-RFP-3xMAD10.

Circular Dichroism (CD)

CD can provide crucial insights due to the presence of beta sheets in their structures. Typically, β-sheets exhibit a characteristic negative band around 218 nanometers (nm) and a positive band near 195 nm, while α-helices show a negative band around 222 nm and a positive band around 190 nm. However, a limitation of the equipment used restricts the analyzed wavelength range to between 200 and 250 nm.

The pipeline acquisition of the spectrum was followed by a 2 seconds of integration for each point in each 6 spectrum calculated at room temperature (20 ºC). From the collected data, it is calculated its average, subtracted from the water absorbance, normalized, and smoothed using a Savitzky-Golay filter.

imagem-2024-09-23-150429757.png

Figure 7. CD acquisition of Barbie1-RFP-3xMad10 protein at 20 ºC.

Barbie1-RFP-3xMad10 Heating Ramp

An essential evaluation for the protein is the heating ramp, in which we can evaluate how the protein acts in different temperatures and their melting point. For doing so, the same procedure was followed, in which each protein was heated from 20 to 95ºC, as shown in Figure 8. The data acquisition was made in steps of 5ºC.

The CD absorbance is shown as the figure’s color map. It is evident that the protein show absorbance in the secondary structure regions, indicating a thermal resistivity.

imagem-2024-09-23-151406532.png

Figura 8. Heating ramp for Barbie1-RFP-3xMad10.

Barbie1-RFP-3xMad10 Pontual Spectra

With this result, it is interesting to observe that, even after the protein engineering, the protein could still maintain great results for tests like the described. To further explore it, we also ran a protein pontual spectra, in which we measured CD absorbance for only 215 nanometers, but each measure made at 1ºC increase.

Figure 9 shows the result of the experiment, which emphasizes the decrease in the CD absorbance for Barbie1-RFP-3xMad10 as already expected from Figure 8. Although Barbie1-RFP-3xMad10 loses part of its absorbance with the increase of temperature, it is not possible to confirm its denaturation. A viable hypothesis is some reduction in the beta-sheets structures.

imagem-2024-09-23-151807308.png

Figure 9. Pontual spectra for Barbie1-RFP-3xMad10 at 215 nm.

In the first instant, this behavior might show a bad characteristic of Barbie1-RFP-3xMad10. However, the loss of function at high temperatures can be used to properly design the water filter. Once saturated with plastics, the system can be heated to high temperatures making Barbie1-RFP-3xMad10 lose part of its structure. Thus, knowing the function is correlated to the form, the protein might lose its plastic affinity, allowing the removal of microplastics and further processing for reutilizing it.

Heating and Cooling Barbie1-RFP-3xMad10 Ramp

With this possibility in our sights, the following experiment is heating and cooling the protein, in a way we can observe how its structures will behave. Illustrated in Figure 10, the plot was rotated for better visualizing the CD absorbance in function of temperature variation. From left to right, it is possible to firstly see how Barbie1-RFP-3xMad10behaves when heated and then cooled.

As a result, protein’s CD absorbance in the cooling ramp still presented a signal, indicating it maintained its structure and a consequently renovelation of the possible denoveled parts.

heating-cooling-barbie.png

Figure 10. Heating and cooling ramp CD absorbance of Barbie1-RFP-3xMad10.

Specific Wavelength Analysis

If we deepen our analysis into 210 and 222 nm, which are in the region of secondary structure absorbance, we can better understand the variancy. In Figure 11, the graphs are related to the 210 nm wavelength, whereas the graphs from below to the 222 nm.

In a dashed gray line, an average initial absorbance for Barbie1-RFP-3xMad10 was added to highlight the loss of the original absorbance. Additionally, both wavelengths show an interesting behavior, that is an important trend change between 55 and 60 ºC, which may be a critical point for the protein. From the cooling graphs, it is notable how the protein loses part of the original signal, indicating a conformational alteration from the temperature change.

For future adaptation of Barbie1-RFP-3xMad10, it is fundamental to assess and design the protein in a way it can be further reused for multiple filtration.

imagem-2024-09-23-160301522.png

Figure 11. Two specific wavelengths, 210 nm above and the 222 nm below, plots for the heating and cooling ramp for Barbie1-RFP-3xMad10.

Barbie1-RFP-3xMad10 Melting Temperature

From the previous result of analyzing a specific wavelength, it is possible to note the curve’s behavior, which is usual for a protein denaturation. Thus, a viable way to uncover the protein’s melting temperature is by fitting the Boltzmann curve in the collected data, which is given by:

$$A2 + (A1 - A2) / (1 + np.exp((x - x0) / dx)),$$

Where:

  • A1 is the protein’s native state;
  • A2 is the completely unfolded protein’s state;
  • Tm is the protein’s melting temperature;
  • k is a constant that represents the curve inclination.

After fitting the curve in Figure 12, the melting point found for Barbie1-RFP-3xMad10 is approximately 58.81ºC. This information is fundamental for building the filter, since it represents stability in a higher temperature. In particular, it is known that in water purifiers the temperatures do not exceed 30 to 35ºC, making the protein appropriate to the system.

fitting-curve-barbie.png

Figure 12. Boltzmann fit in the absorbance values of Barbie1-RFP-3xMad10 in a heating ramp.

Plastic Interaction Analysis

Through secondary structures absorbance, it was possible to observe that Barbie1-RFP-3xMad10 structure do not change in the presence of plastic nanoparticles

Finally, the last experimental evaluation at CD was with the polystyrene nanoplastics. We used a 100 nanometer particle which was inserted into a solution with Barbie1-RFP-3xMad10.

Nevertheless, as reported from the computational simulations, we do not expect any conformational change. Since CD is only evaluating the absorbance in secondary structures, there may not happen major differences in the experiment.

Consequently, Figure 13 shows as projected. When the Barbie1-RFP-3xMad10a spectra is plotted with the nanoparticle solution spectra, minor changes can be observed. Specifically, this contrast might be a result of the polystyrene nanoparticle scattering, creating noise in the measurement.

barbie-np-cd.png

Figure 13. Plastic binding proteins Barbie1-RFP-3xMad10and BaCBM2-RFP-3xMad10 spectra acquisition with 100 nm polystyrene nanoparticles.

Considering we were not able to assess the plastic-protein binding, further studies might be done. One important form to evaluate this is the infrared spectroscopy, which allows studying the vibrational frequencies of molecules.

Sequence and Features


Assembly Compatibility:
  • 10
    COMPATIBLE WITH RFC[10]
  • 12
    COMPATIBLE WITH RFC[12]
  • 21
    INCOMPATIBLE WITH RFC[21]
    Illegal BglII site found at 594
  • 23
    COMPATIBLE WITH RFC[23]
  • 25
    INCOMPATIBLE WITH RFC[25]
    Illegal AgeI site found at 88
  • 1000
    COMPATIBLE WITH RFC[1000]