Coding

Part:BBa_K5033001

Designed by: Jonas Martin Westphal   Group: iGEM24_Aachen   (2024-09-16)

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).

GOI
Figure 1: Schematic view of the fusion protein's coding sequence.

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.

properties
Figure 2: Biochemical properties of mFadA[B]_GSLinker_CDA[Alpha] as described by Benchling.


Protein Structure Prediction

3D
Figure 3: Tertiary structure of mFadA[B]_GSLinker_CDA[Alpha] as predicted by AlphaFold2. From left to right: mFadA[B] as an alpha helix (blue), the flexible linker, the enzyme-mutant and the freely accessible His-Tag (orange).

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

5-FC
Figure 4: Model of interactions of the amino acids in the active site of the fusionprotein with 5-fluorocytosine. As predicted using RoseTTAFold All-Atom.

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®.

purification
Figure 5: SDS page showing the proteins in the elution fractions. The number corresponds to the imidazole concentration (in mM) in the elution buffer respectively. Example: E10 is an elution buffer with 10mM imidazole.
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.

example HPLC
Figure 6: Example HPLC Chromatogram obtained from measurement using 400 µM 5-FC and 200 µM 5-FU; the first peak being the 5-FC, the second 5-FU.


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.


HPLC results
Figure 7: Representation of the relative amount of substance nrel for the Alpha mutant against the time t obtained from HPLC measurements.


uvvis results
Figure 8: Representation of the relative amount of substance nrel for the Alpha mutant obtained from Uv-Vis measurement against the time t.

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).

NMR results
Figure 9: Summary of the obtained RPA, t(1/2)RP and v0 values from Uv-Vis measurements.


results overwiev
Figure 10: Comparison of the relative enyzme performance for the wild type and Alpha variant; data obtained from HPLC.

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


Assembly Compatibility:
  • 10
    COMPATIBLE WITH RFC[10]
  • 12
    COMPATIBLE WITH RFC[12]
  • 21
    INCOMPATIBLE WITH RFC[21]
    Illegal BamHI site found at 115
    Illegal XhoI site found at 1396
  • 23
    COMPATIBLE WITH RFC[23]
  • 25
    INCOMPATIBLE WITH RFC[25]
    Illegal NgoMIV site found at 1254
    Illegal NgoMIV site found at 1341
  • 1000
    COMPATIBLE WITH RFC[1000]
[edit]
Categories
Parameters
None