Part:BBa_K5033001
OncoBiotica: mFadA[B]_GSLinker_CDA[Alpha]
If you are interested in an overview of the parts designed by the iGEM Team Aachen 2024, visit our Parts page.
This part, developed by iGEM Aachen 2024, is our 'Alpha' mutant of the basic part BBa_K5033000. Within the codA cytosine deaminase (CDA) domain, it contains a D314A mutation (the aspartic acid on position 314 has been exchanged with alanine). In the context of our fusionprotein it is a D352A mutation. This mutant has been selected for research in the lab with the help of literature research and modeling.
Its main use is to analyze the catalytic behaviour of our fusionprotein by replicating already described mutations of the CDA within our fusionprotein.
It is an optimized version of BBa_K5033000 which served as the foundation for exploring the concept of microbiota-directed cancer therapy. This part encodes a fusion protein designed to combine two functionalities. Binding specific bacteria and having an optimized enzymatic function. This part is to be cloned into a vector based on an inducable expression system. iGEM Aachen 2024 used a pET21b(+) vector.
iGEM Aachen 2024 successfully demonstrated that the enzymatic function is enhanced in comparison to the CDA-wild type fusionprotein. This is one of five mutants analyzed by iGEM Aachen 2024 in addition to the CDA-wild type fusionprotein.
See the other four variants: 'Beta' (BBa_K5033002), 'Gamma' (BBa_K5033003), 'Epsilon' (BBa_K5033004) and 'Theta' (BBa_K5033005).
Contents
Part Composition
The first protein domain is derived from the part BBa_K4990002 but has been codon optimized for expression in E. coli. It is the mFadA B-domain, found in various Fusobacterium strains. This part has already been well described by the iGEM23_CPU-CHINA team. This domain should be able to bind to FadA pili on Fusobacterium nucleatum and its former subspecies Fusobacterium nucleatum, F. polymorphum, F. vincentii, F. animalis via self assembly.
To further investigate the binding domain's functionality, iGEM Aachen 2024 created the basic part BBa_K5033006. This variant replaces the enzyme in our fusion protein with eGFP as a reporter protein.
The second functional protein domain is linked to the mFadA B-domain by a synthetic flexible linker consisting of Glycin and Serine in alternating order. This linker is eleven amino acids long.
This second functional domain of the fusionprotein is a mutant of the codA cytosine deaminase (CDA) that is native to E. coli. This mutant contains the D314A mutation of the enzyme. This Enzyme converts cytosine to uracil in its host organism but it is also able to convert 5-fluorocytosine (5-FC) to 5-fluorouracil (5-FU).[1]
In this case the Enzyme can be used for an enzyme directed prodrug therapy. To be precise, 5-FC is the non-toxic substrate and 5-FU is the active chemotherapeutic agent.
The fusionprotein encoded by this part also contains a downstream hexa-histidine tag for protein purification.
Protein Modeling
To find interesting mutations that shall be investigated in the lab, our team used a research based modeling approach.
Before transformation of this biological part (cloned into the pET21b(+) plasmid backbone), the structure of the expected fusionprotein was modeled.
Selection of the Mutant
During literature research we found the D314A mutation. This is a very promising mutant as it was reported to have a high specificity and efficiency compared to the wild type (WT) enzyme [2]. This mutant was the most stable variant (ΔΔGfold = -6.8 kcal/mol) calculated with the use of computer-assisted recombination (CompassR) [3].
Biochemical Properties
The fundamental biochemical properties like molecular mass and extinction coefficient are important for a lot of SynBio work done with proteins. To see an overview of these properties, have a look at figure 2.
Protein Structure Prediction
The tertiary structure has been predicted using AlphaFold2 by DeepMind. In this case it is especially important, that the binding domain and the His-Tag are freely available.
Modeling of Substrate and Active Site Interaction
Why RoseTTAFold All-Atom?
Proteins rarely act alone. Although substantial progress in the prediction of protein structures has been made, modeling of proteins and their ligands still remains challenging. The development of RoseTTAFold All-Atom (RFAA) aims to tackle this issue by building a neural network that is trained to accurately model general biomolecules containing a wide range of non-protein components. In contrast to other tools that only include sequence based modeling, RFAA incorporates a graphical representation that models non-protein molecules at the atomic level, capturing their chemical bonds and interactions. In combination with the training data set that also includes ligand-bound protein structures from the Protein Data Bank (pdb), it allows RFAA to predict protein structures, ions and non-protein ligands. Interestingly, during our project, DeepMind released a new AlphaFold version (v3) that includes selected ions and ligands. However, an earlier release would not have been advantageous for us, as 5-FC and cytosine are not among the selected ligands that AlphaFold3 includes. Nevertheless, this shows that the improvements made this year mark a significant step forward, paving the way for more refined and accurate modeling of proteins and ligands in the future.
We observed a good overall structural alignment of the wild type enzymes' crystal structure [4] to the RoseTTAFold All-Atom model. Upon closer inspection of the active site, we noticed small differences in torsion angles of the side chains which naturally led to slight differences in bond lengths between amino acids and the ligand. However, these differences are inherent to the modeling process and do not reflect significant deviations. Therefore, they do not compromise the reliability of our approach for predicting structural changes in the mutants. This allowed us to apply the approach to the generated mutants by CompassR.
Modeling Results
All critical amino acids (E217, H246, and D313) are predominantly in the correct conformation relative to the iron and substrate, allowing the reaction to occur. Aspartic acid at position 314 lies on a flexible loop that undergoes conformational changes during ligand binding [2]. This results in a interaction with the pyrimidine ring of 5-FC or cytosine [2]. The substitution changes how the side chain interacts with the ligand, particularly with the fluorine atom in 5-FC. Compared to the carboxyl group bearing side chain of aspartic acid, alanine has a hydrophobic side chain due to it's methyl group. The partially positively charged hydrogen atoms of the methyl group may form stabilizing interactions with the highly electronegative fluorine, in contrast to the potential destabilizing repulsion between the negatively charged carboxyl group and fluorine in aspartic acid. Additionally, this change of the functional group results in a favorable van der Waals interaction between the fluorine atom and the methyl group.
Looking at our model generated by RFAA, fluorine and the side chain at position 314 exhibit a distance in a range of possible interactions (3.446 Å). This could explain the reported slight improvement of the Km value (2.8 mM for D314A, 3.3 mM for the WT enzyme) [2]. It appears that the mutant results in a loss of interaction between Q156 and the ligand. However, this can be an artifact due to incorrect torsion angles predicted by RFAA as stated above.
Cloning of the Plasmid
To build the plasmid containing the gene for our Alpha variant. We used the plasmid we already had for our WT-Fusionprotein (BBa_K5033000; pET21b(+)_mFad[A]_GSLinker_CDA[WT]). The gene sequence for this part contains a BamHI restriction site between the linker and the enzyme. The backbone contains a XhoI restriction site at the end of the gene insert.
After modeling of the Alpha variant we ordered the gene fragment, encoding this variant. We made sure to include the correct restriction sites.
The backbone was prepared using the BamHI and XhoI restriction enzymes. After digestion, the cut backbone was cleaned up using an agarose gel and a gel extraction kit. The same was done for the Insert.
After gel cleanup the cut backbone and insert were ligated using the T4 Ligase.
To enhance the efficiency of the plasmid transformation into E. coli BL21 (DE3) the plasmid was first propagated via transformation in E. coli DH5α.
The propagated pET21b(+)_mFadA[B]_GSLinker_CDA[Alpha] plasmid could then be purified with a plasmid miniprep kit and used for transformation into the production organism E. coli BL21 (DE3).
Producing the Fusionprotein
After successful transformation of the pET21b(+)_mFadA[B]_GSLinker_CDA[Alpha] plasmid into the production organism E. coli BL21 (DE3) the protein could be expressed and purified. The pET21b(+) backbone has a lac operon (including the lacI repressor), which can be induced with IPTG (IUPAC: Propan-2-yl 1-thio-β-D-galactopyranoside).
Expression and Purification of the Fusionprotein
The fusionprotein was expressed by adding IPTG to the medium to a final concentration of 1mM.
The His-tagged protein was purified using a Protino Ni-IDA 2000 packed column by Macherey & Nagel®.
The fusionptrotein is expected to have a molecular weight of 52.08kDA (cf. Fig. 2). This corresponds to the big bands visible on the gel.
This gel shows that the E10 fraction still has a lot of impurities. The E50 fraction was desalted and stored in 50mM TRIS buffer, to use for the kinetic assays.
Kinetic Assays
If you are interested in the methods used, take a look at our Experiments page.
UV-Vis
An assay using a spectrometer was devised. This measures the peak of the absorbance spectrum, which shifts to lower wavelengths as the ratio of 5-FU to 5-FC increases. We established a linear relationship, which we calibrated using standards. More details can be found in the UV-Vis Peak Shift Assay protocol on our Experiments page.
Note that this method does not give a high resolution beyond 5% 5-FU/5-FC intervals.
High-Performance Liquid Chromatographie (HPLC)
We used Reverse Phase High Performance Liquid Chromatography for quantitative Analysis of 5-fluorocytosine and 5-fluorouracil in mutual solution. The results seen below were all measured with the same method. (see Experiments page)
Standards at between 10 µM and 500 µM were made to translate the peak area into compound concentration.
After measuring, the chromatograms were evaluated with “OpenChrom” by Lablicate. For this, a baseline subtraction filter was applied, after this the standard first derivitave peak detector and trapezoid peak integrator were run. We identified para-aminobenzoicacid as a potential internal standard, but no problems which would necessitate the use of an internal standard arose.
Enzyme Kinetics
Similar to the wild type enzyme, the mutant cytosine deaminase followed Michaelis-Menten kinetics. However, due to the limited solubility of 5-fluorocytosine (5-FC), achieving substrate concentrations sufficient for determining the mutants' maximum velocity (Vmax) was difficult, leading to variability in kinetic measurements. This made the direct determination of Km and Vmax unreliable. To assess the mutants' performance, three key metrics were employed: (1) the time required to reach half-maximal product formation t1/2RP, (2) the initial reaction rate (initial velocity) v0, and (3) the total product formed over time as a measure of relative efficiency RFA.
For each mutant and method (HPLC, NMR, and UV-Vis Spectroscopy), relative substrate and product concentrations were plotted against time, and these graphs formed the basis of the analysis. The time points at which half-maximal product formation was reached were graphically determined, and substrate and product concentrations for each variant are plotted against time in figure 7. As with the wild type, the relative efficiency of the mutant was calculated by determining the area under the graph within a predefined time interval, assuming linearity between data points. This area represents the total product formed over time and provides a valuable comparison of efficiency of the mutant.
The initial velocities were determined by linear approximation of relative product formation during the first 80 seconds for both the HPLC and Uv-Vis methods. Each mutant's performance was directly compared to the wild type by plotting against the time, with the slope of the regression line representing the relative velocity for each condition.
The relative activity of each mutant was calculated as the mean of the two key parameters: the time to half-maximal product formation and the initial velocity, with the standard deviation providing an estimate of error. This combined measure of relative activity, allowed for a direct comparison of each variant's performance (cf. fig. 9).
This detailed approach allowed us to quantify the relative activity and efficiency of each mutant, highlighting the differences in their performance compared to the wild type.
Conclusion
We were able to show, that the Alpha variant has approximately double catalytic performance in comparison to the wild type fusionprotein (cf. fig. 10).
Literature has already described this single point mutation within the codA CDA to improve its performance. By expressing and analyzing our Alpha variant (fusionprotein) we showed that this mutation has similar effects within our fusionprotein.
References
[1] Aučynaitė, A., Rutkienė, R., Tauraitė, D., Meškys, R., Urbonavičius, J., 2018. Discovery of Bacterial Deaminases That Convert 5-Fluoroisocytosine Into 5-Fluorouracil. Frontiers in Microbiology 9.. https://doi.org/10.3389/fmicb.2018.02375
[2] Sheri D. Mahan, Greg C. Ireton, Catherine Knoeber, Barry L. Stoddard, Margaret E. Black, Random mutagenesis and selection of Escherichia coli cytosine deaminase for cancer gene therapy, Protein Engineering, Design and Selection, Volume 17, Issue 8, August 2004, Pages 625–633, https://doi.org/10.1093/protein/gzh074
[3] Cui, H., Cao, H., Cai, H., Jaeger, K., Davari, M.D., Schwaneberg, U., 2020. Computer‐Assisted Recombination (CompassR) Teaches us How to Recombine Beneficial Substitutions from Directed Evolution Campaigns. Chemistry – A European Journal 26, 643–649.. https://doi.org/10.1002/chem.201903994
[4] PDB Entry - 1RA0. https://doi.org/10.2210/pdb1RA0/pdb
Sequence and Features
- 10COMPATIBLE WITH RFC[10]
- 12COMPATIBLE WITH RFC[12]
- 21INCOMPATIBLE WITH RFC[21]Illegal BamHI site found at 115
Illegal XhoI site found at 1396 - 23COMPATIBLE WITH RFC[23]
- 25INCOMPATIBLE WITH RFC[25]Illegal NgoMIV site found at 1254
Illegal NgoMIV site found at 1341 - 1000COMPATIBLE WITH RFC[1000]
None |