Week 1 HW: Principles and Practices

I am interested in developing an engineered probiotic system designed to release the lactase enzyme on demand in the human gut for individuals with lactose intolerance. This project is entirely in silico, combining concepts from synthetic biology, microbiome modeling, and systems biology without any wet-lab implementation.

The system would simulate a probiotic chassis such as Lactobacillus or Bifidobacterium, equipped with virtual genetic circuits inspired by lactose metabolism. These circuits would model regulatory control of lactase expression based on local lactose concentration, using logic-gate–like behavior and feedback mechanisms. Enzyme production would increase when lactose is present and decrease once lactose is depleted, allowing adaptive and resource-efficient regulation.

Because this project represents an early, in-silico design phase, its governance goals focus on the responsible framing, communication, and interpretation of computational results rather than regulation of a finalized biological product.

From my perspective, scenario-based risk modeling can be prioritized over the other governance options, because all three approaches address public health and safety either directly or indirectly. Scenario-based analysis explicitly explores what could go wrong if an in-silico model were physically implemented, making it the most direct mechanism for anticipating risks to gut microbiome balance or unintended metabolic effects. However, maintaining scientific integrity also plays a critical indirect role in protecting public health: by avoiding overclaiming the safety or effectiveness of a purely conceptual model, the transition from simulation to real-world application becomes more cautious, accurate, and oriented toward appropriate experimental validation, thereby reducing the likelihood of harmful misinterpretations. Similarly, ensuring ethical transparency through clear and accurate documentation of modeling assumptions, parameters, and limitations improves how the model is interpreted and reused by others, helping prevent incorrect applications that could ultimately pose health risks.

in this homework, AI ChatGPT assisted me in organizing and clearly articulating my answers and descriptions, ensuring that the content is well-structured and easy to understand.

Sources:

• Emerging biotechnologies: Technology, choice and the public good. (n.d.). Nuffield Council on Bioethics. Retrieved February 9, 2026, from https://www.nuffieldbioethics.org/publication/emerging-biotechnologies-technology-choice-and-the-public-good/
• Gingold-Belfer, R., Levy, S., Layfer, O., Pakanaev, L., Niv, Y., Dickman, R., & Perets, T. T. (2020). Use of a Novel Probiotic Formulation to Alleviate Lactose Intolerance Symptoms-a Pilot Study. Probiotics and Antimicrobial Proteins, 12(1), 112–118. https://doi.org/10.1007/s12602-018-9507-7
• iGEM Responsibility. (n.d.). Retrieved February 9, 2026, from https://responsibility.igem.org/
• Khalil, A. S., & Collins, J. J. (2010). Synthetic biology: Applications come of age. Nature Reviews Genetics, 11(5), 367–379. https://doi.org/10.1038/nrg2775
• Lactose Intolerance—NIDDK. (2024, January 30). National Institute of Diabetes and Digestive and Kidney Diseases. https://www.niddk.nih.gov/health-information/digestive-diseases/lactose-intolerance

Assignment (Week 2 Lecture Prep):

Homework Questions from Professor Jacobson:

1. Error rate and genome context • From the slide N°= 8 , DNA polymerase has an error rate of ~1 in 10⁶ bases. • With the human genome of ~3 × 10⁹ bp, this would result in ~3,000 errors per replication without repair. • Biology reduces this discrepancy with proofreading activity of DNA polymerase (3′→5′ exonuclease) and post-replication mismatch repair like MutS, NER, BER…, which collectively reduce the final error rate to ~1 in 10⁹–10¹⁰.
2. Human protein: ~1036 bp (~345 amino acids), With ~3 codons per amino acid on average, the number of possible DNA sequences for an average human protein is ~3³⁴⁵ (~10¹⁶⁴ possible sequences). Not all sequences work in practice because of Mutations: Insertions, deletions, transitions, and transversions that can introduce frameshifts or premature stop codons, making the protein non-functional. Also, there are some mechanism of regulations that make some Sequences creating unwanted secondary structures in mRNA, affect splicing, or introduce cryptic signals that disrupt translation.

Homework Questions from Dr. LeProust:

1. Most commonly used method for oligo synthesis Today, almost all synthetic DNA is made using phosphoramidite solid-phase synthesis. This method adds one nucleotide at a time on a solid support and is reliable, efficient, and easy to automate, which is why it became the standard for modern DNA synthesizers. https://biolabmix.ru/en/info/detail/oligonucleotide-synthesis/#:~:text=The%20most%20common%20approach%20to,for%20example%2C%20by%20attaching%20fluorophores.

2. Why it’s hard to make oligos longer than ~200 nt Each step in chemical DNA synthesis is very efficient but not perfect, so small errors happen every time a base is added. As the oligo gets longer, these errors pile up, and beyond about 200 nucleotides it becomes very difficult to get a clean, full-length sequence. https://pubs.rsc.org/en/content/articlepdf/2025/sc/d4sc06958g

3. Why you can’t directly synthesize a 2000 bp gene Making a 2000-base gene in one piece would accumulate too many chemical errors and damaged bases to be useful. Instead, companies synthesize short oligos and then assemble them enzymatically, followed by cloning and sequence checking to make sure the gene is correct. https://www.pnas.org/doi/10.1073/pnas.2237126100#:~:text=The%20broader%20implications%20of%20the,without%20multiple%20repair/selection%20steps.

Homework Question from George Church:

All animals require the same 10 essential amino acids because they cannot synthesize them and must obtain them from their diet. These are: histidine, isoleucine, leucine, lysine, methionine, phenylalanine, threonine, tryptophan, valine, and arginine (arginine is essential for all animals and conditionally essential in adult humans). The “lysine contingency” refers to the idea that lysine is often the limiting essential amino acid in plant-based diets, especially those dominated by cereals like wheat, rice, or maize. Since animals cannot make lysine, their growth and health are directly constrained by how much lysine is available in their food. So knowing that all animals share the same essential amino acid requirements makes lysine’s importance stand out even more. It shows that lysine is not just nutritionally important but evolutionarily critical.

https://www.kemin.com/ap/en/blog/animal/amino-acids-for-animal-health#:~:text=Essential%20amino%20acids:%20These%20are,essential)%2C%20leucine%20and%20lysine

Week 2 HW: DNA Read, Write, and Edit

• In this part, I imported The complete 48,502 bp linear genome of bacteriophage lambda from NCBI GenBank into Benchling. This sequence corresponds to the Lambda DNA sold by NEB (N3011) and will be used for in-silico restriction digestion.

• Then simulated restriction enzyme digestion using EcoRI, HindIII, BamHI, KpnI, EcoRV, SacI, and SalI. By running in-silico gel electrophoresis . The resulting virtual gel shows discrete bands corresponding to these fragments, which demostrates how sequence information maps to physical separation in gel electrophoresis.

• To create a pattern in the style of Paul Vanouse’s work, I experimented with different combinations of restriction enzymes to control the gel band patterns. By adjusting the number and length of the resulting DNA fragments, I explored how these parameters influence the final visual outcome. Through this process, I ultimately obtained a gel pattern resembling a butterfly shape.

• This helped me understand how restriction digests and gels work before doing any real lab experiment. I treated this as both a technical exercise and a creative exploration, inspired by DNA gel art concepts.

• For the DNA design challenge, I chose a protein related to my project interest in engineered probiotics and conditional enzyme release in the gut.The enzyme β-galactosidase is well-characterized and commonly expressed in Escherichia coli, making it an ideal candidate for computational DNA design and expression modeling.
• I first searched online database UniProt to obtain the amino acid sequence of the protein.

• the amino acid equence was as follow:

>sp|P00722|BGAL_ECOLI Beta-galactosidase OS=Escherichia coli (strain K12) OX=83333 GN=lacZ PE=1 SV=2
MTMITDSLAVVLQRRDWENPGVTQLNRLAAHPPFASWRNSEEARTDRPSQQLRSLNGEWR
FAWFPAPEAVPESWLECDLPEADTVVVPSNWQMHGYDAPIYTNVTYPITVNPPFVPTENP
TGCYSLTFNVDESWLQEGQTRIIFDGVNSAFHLWCNGRWVGYGQDSRLPSEFDLSAFLRA
GENRLAVMVLRWSDGSYLEDQDMWRMSGIFRDVSLLHKPTTQISDFHVATRFNDDFSRAV
LEAEVQMCGELRDYLRVTVSLWQGETQVASGTAPFGGEIIDERGGYADRVTLRLNVENPK
LWSAEIPNLYRAVVELHTADGTLIEAEACDVGFREVRIENGLLLLNGKPLLIRGVNRHEH
HPLHGQVMDEQTMVQDILLMKQNNFNAVRCSHYPNHPLWYTLCDRYGLYVVDEANIETHG
MVPMNRLTDDPRWLPAMSERVTRMVQRDRNHPSVIIWSLGNESGHGANHDALYRWIKSVD
PSRPVQYEGGGADTTATDIICPMYARVDEDQPFPAVPKWSIKKWLSLPGETRPLILCEYA
HAMGNSLGGFAKYWQAFRQYPRLQGGFVWDWVDQSLIKYDENGNPWSAYGGDFGDTPNDR
QFCMNGLVFADRTPHPALTEAKHQQQFFQFRLSGQTIEVTSEYLFRHSDNELLHWMVALD
GKPLASGEVPLDVAPQGKQLIELPELPQPESAGQLWLTVRVVQPNATAWSEAGHISAWQQ
WRLAENLSVTLPAASHAIPHLTTSEMDFCIELGNKRWQFNRQSGFLSQMWIGDKKQLLTP
LRDQFTRAPLDNDIGVSEATRIDPNAWVERWKAAGHYQAEAALLQCTADTLADAVLITTA
HAWQHQGKTLFISRKTYRIDGSGQMAITVDVEVASDTPHPARIGLNCQLAQVAERVNWLG
LGPQENYPDRLTAACFDRWDLPLSDMYTPYVFPSENGLRCGTRELNYGPHQWRGDFQFNI
SRYSQQQLMETSHRHLLHAEEGTWLNIDGFHMGIGGDDSWSPSVSAEFQLSAGRYHYQLV
WCQK

• After selecting the protein, I converted the amino acid sequence of β-galactosidase (1024 residues) into the corresponding DNA sequence using the Sequence Manipulation Suite Reverse Translate tool. Because the genetic code is degenerate, multiple codons can encode the same amino acid. The resulting 3072 bp DNA sequence represents one valid nucleotide sequence capable of encoding the β-galactosidase protein.

• the resulted DNA sequence was as follow:

>reverse translation of sp|P00722|BGAL_ECOLI Beta-galactosidase OS=Escherichia coli (strain K12) OX=83333 GN=lacZ PE=1 SV=2 to a 3072 base sequence of most likely codons.
atgaccatgattaccgatagcctggcggtggtgctgcagcgccgcgattgggaaaacccg
ggcgtgacccagctgaaccgcctggcggcgcatccgccgtttgcgagctggcgcaacagc
gaagaagcgcgcaccgatcgcccgagccagcagctgcgcagcctgaacggcgaatggcgc
tttgcgtggtttccggcgccggaagcggtgccggaaagctggctggaatgcgatctgccg
gaagcggataccgtggtggtgccgagcaactggcagatgcatggctatgatgcgccgatt
tataccaacgtgacctatccgattaccgtgaacccgccgtttgtgccgaccgaaaacccg
accggctgctatagcctgacctttaacgtggatgaaagctggctgcaggaaggccagacc
cgcattatttttgatggcgtgaacagcgcgtttcatctgtggtgcaacggccgctgggtg
ggctatggccaggatagccgcctgccgagcgaatttgatctgagcgcgtttctgcgcgcg
ggcgaaaaccgcctggcggtgatggtgctgcgctggagcgatggcagctatctggaagat
caggatatgtggcgcatgagcggcatttttcgcgatgtgagcctgctgcataaaccgacc
acccagattagcgattttcatgtggcgacccgctttaacgatgattttagccgcgcggtg
ctggaagcggaagtgcagatgtgcggcgaactgcgcgattatctgcgcgtgaccgtgagc
ctgtggcagggcgaaacccaggtggcgagcggcaccgcgccgtttggcggcgaaattatt
gatgaacgcggcggctatgcggatcgcgtgaccctgcgcctgaacgtggaaaacccgaaa
ctgtggagcgcggaaattccgaacctgtatcgcgcggtggtggaactgcataccgcggat
ggcaccctgattgaagcggaagcgtgcgatgtgggctttcgcgaagtgcgcattgaaaac
ggcctgctgctgctgaacggcaaaccgctgctgattcgcggcgtgaaccgccatgaacat
catccgctgcatggccaggtgatggatgaacagaccatggtgcaggatattctgctgatg
aaacagaacaactttaacgcggtgcgctgcagccattatccgaaccatccgctgtggtat
accctgtgcgatcgctatggcctgtatgtggtggatgaagcgaacattgaaacccatggc
atggtgccgatgaaccgcctgaccgatgatccgcgctggctgccggcgatgagcgaacgc
gtgacccgcatggtgcagcgcgatcgcaaccatccgagcgtgattatttggagcctgggc
aacgaaagcggccatggcgcgaaccatgatgcgctgtatcgctggattaaaagcgtggat
ccgagccgcccggtgcagtatgaaggcggcggcgcggataccaccgcgaccgatattatt
tgcccgatgtatgcgcgcgtggatgaagatcagccgtttccggcggtgccgaaatggagc
attaaaaaatggctgagcctgccgggcgaaacccgcccgctgattctgtgcgaatatgcg
catgcgatgggcaacagcctgggcggctttgcgaaatattggcaggcgtttcgccagtat
ccgcgcctgcagggcggctttgtgtgggattgggtggatcagagcctgattaaatatgat
gaaaacggcaacccgtggagcgcgtatggcggcgattttggcgataccccgaacgatcgc
cagttttgcatgaacggcctggtgtttgcggatcgcaccccgcatccggcgctgaccgaa
gcgaaacatcagcagcagttttttcagtttcgcctgagcggccagaccattgaagtgacc
agcgaatatctgtttcgccatagcgataacgaactgctgcattggatggtggcgctggat
ggcaaaccgctggcgagcggcgaagtgccgctggatgtggcgccgcagggcaaacagctg
attgaactgccggaactgccgcagccggaaagcgcgggccagctgtggctgaccgtgcgc
gtggtgcagccgaacgcgaccgcgtggagcgaagcgggccatattagcgcgtggcagcag
tggcgcctggcggaaaacctgagcgtgaccctgccggcggcgagccatgcgattccgcat
ctgaccaccagcgaaatggatttttgcattgaactgggcaacaaacgctggcagtttaac
cgccagagcggctttctgagccagatgtggattggcgataaaaaacagctgctgaccccg
ctgcgcgatcagtttacccgcgcgccgctggataacgatattggcgtgagcgaagcgacc
cgcattgatccgaacgcgtgggtggaacgctggaaagcggcgggccattatcaggcggaa
gcggcgctgctgcagtgcaccgcggataccctggcggatgcggtgctgattaccaccgcg
catgcgtggcagcatcagggcaaaaccctgtttattagccgcaaaacctatcgcattgat
ggcagcggccagatggcgattaccgtggatgtggaagtggcgagcgataccccgcatccg
gcgcgcattggcctgaactgccagctggcgcaggtggcggaacgcgtgaactggctgggc
ctgggcccgcaggaaaactatccggatcgcctgaccgcggcgtgctttgatcgctgggat
ctgccgctgagcgatatgtataccccgtatgtgtttccgagcgaaaacggcctgcgctgc
ggcacccgcgaactgaactatggcccgcatcagtggcgcggcgattttcagtttaacatt
agccgctatagccagcagcagctgatggaaaccagccatcgccatctgctgcatgcggaa
gaaggcacctggctgaacattgatggctttcatatgggcattggcggcgatgatagctgg
agcccgagcgtgagcgcggaatttcagctgagcgcgggccgctatcattatcagctggtg
tggtgccagaaa

• After reverse translation, I verified the identity of the resulting nucleotide sequence by performing a BLASTn search against the reference lacZ gene from Escherichia coli K-12. The alignment showed 100% query coverage with an E-value of 0.0, confirming a highly significant match. The percent identity was ~84%, which is expected because reverse translation produces a synonymous DNA sequence that differs at the codon level while still encoding the same β-galactosidase protein. This result confirmed that the reverse-translated sequence correctly corresponds to the lacZ gene.

Next, I performed codon optimization of the sequence originates from E. coli K-12 to improve expression efficiency in a Lactobacillus probiotic strain (delbrueckii subsp. Bulgaricus), as this organism is the intended chassis for conditional lactase expression in the human gut, to ensure efficient translation in the final probiotic host organism. Codon optimization was performed using a host-specific algorithm using the Vector Builder codon orimisation tool that adjusts synonymous codon usage to match the preferred codons of L. delbrueckii while preserving the original amino acid sequence.

• the resulted optimised sequence is as following:

Host organism: Lactobacillus delbruekii susbsp. Bulgaricus ATCC 11842 = JCM 1002 Original Sequence: GC=59.80%, CAI=0.72 Optimized Sequence: GC=60.16%, CAI=0.89

Improved DNA[1]: GC=60.16%, CAI=0.89
ATGACTATGATCACCGACAGCCTGGCAGTTGTTTTGCAACGGCGGGACTGGGAAAACCCGGGCGTCACTCAGTTGAACCGGCTGGCCGCCCACCCACCATTTGCCAGCTGGCGCAACTCCGAAGAAGCCCGGACCGACCGGCCGAGCCAGCAACTGAGAAGCTTGAACGGCGAATGGCGTTTCGCCTGGTTTCCGGCCCCGGAAGCCGTCCCAGAAAGCTGGTTGGAATGCGACCTCCCGGAAGCCGATACCGTCGTGGTGCCGAGCAACTGGCAAATGCACGGCTATGACGCCCCCATCTACACCAATGTTACCTACCCAATTACCGTCAACCCGCCATTTGTCCCGACCGAAAACCCGACTGGTTGCTATAGCTTGACCTTCAACGTTGACGAAAGCTGGCTGCAAGAAGGCCAGACCCGCATTATTTTTGACGGCGTTAACAGCGCCTTCCACTTGTGGTGCAACGGCCGCTGGGTCGGCTACGGCCAGGACAGCCGCTTGCCATCCGAATTTGACCTGAGTGCTTTCTTGCGGGCCGGCGAAAACCGTCTGGCCGTCATGGTCCTGCGCTGGAGCGACGGCAGCTACCTGGAAGACCAAGACATGTGGCGGATGTCCGGCATTTTCCGGGACGTCAGCCTGCTGCACAAGCCGACCACCCAGATTTCCGACTTTCACGTTGCAACCCGGTTCAACGACGACTTCTCTCGGGCTGTGCTGGAAGCTGAAGTCCAGATGTGCGGCGAATTGCGGGACTACCTGCGGGTTACTGTTTCATTGTGGCAGGGCGAAACCCAGGTTGCCTCAGGCACCGCCCCGTTTGGCGGTGAAATTATCGACGAACGCGGCGGGTACGCCGACCGGGTTACCTTGAGACTGAACGTGGAAAACCCGAAGTTGTGGAGCGCCGAAATCCCAAATCTGTACCGCGCCGTCGTCGAATTGCACACCGCTGACGGCACCCTGATCGAAGCCGAAGCCTGCGACGTTGGCTTCCGGGAAGTCCGCATCGAAAACGGCTTGCTGCTCCTGAACGGCAAGCCACTGCTGATCCGGGGCGTTAACCGGCACGAACACCACCCATTGCACGGCCAAGTCATGGACGAACAGACTATGGTCCAGGACATCCTGCTGATGAAGCAGAACAACTTCAACGCTGTTCGTTGCTCACACTATCCAAACCATCCACTGTGGTACACTCTGTGCGACCGGTACGGCCTGTACGTTGTGGACGAAGCCAACATCGAAACTCACGGCATGGTTCCGATGAACCGGCTGACCGACGACCCGAGATGGCTGCCAGCCATGAGCGAACGGGTTACTCGCATGGTTCAACGCGACCGGAACCACCCATCCGTTATTATCTGGAGCCTGGGGAACGAAAGCGGCCACGGCGCCAATCACGACGCTCTGTACCGGTGGATCAAGTCCGTCGACCCATCCCGCCCTGTTCAGTACGAAGGCGGCGGCGCCGATACGACCGCCACCGACATCATCTGCCCAATGTACGCCCGGGTTGATGAAGACCAGCCGTTTCCGGCTGTCCCAAAGTGGAGCATCAAGAAGTGGCTGAGCCTGCCAGGCGAAACTCGGCCGCTGATCCTGTGCGAATACGCCCACGCCATGGGCAACTCCCTGGGCGGCTTTGCCAAGTACTGGCAGGCTTTTCGCCAGTATCCACGGTTGCAGGGCGGCTTTGTTTGGGACTGGGTCGACCAAAGCCTGATCAAGTACGACGAAAACGGCAACCCGTGGAGCGCCTACGGCGGCGACTTTGGCGACACCCCGAACGACCGCCAGTTTTGCATGAACGGTCTGGTTTTCGCTGACCGGACGCCACACCCGGCCCTGACCGAAGCCAAGCACCAGCAGCAGTTCTTCCAGTTCCGGCTGTCAGGCCAGACCATCGAAGTGACTAGCGAATACCTGTTTCGCCACTCCGACAACGAATTGTTGCACTGGATGGTCGCCCTGGACGGCAAGCCACTGGCCAGCGGCGAAGTTCCGCTGGACGTTGCCCCACAGGGCAAGCAGCTGATCGAATTGCCGGAACTGCCGCAGCCGGAAAGCGCCGGCCAACTGTGGCTGACTGTTCGGGTCGTTCAGCCGAACGCCACTGCCTGGTCTGAAGCCGGGCACATCTCAGCCTGGCAGCAGTGGCGCCTGGCCGAAAACTTGAGCGTTACGCTGCCGGCCGCCAGCCACGCCATCCCACACCTGACTACTAGCGAAATGGACTTTTGCATCGAATTGGGCAACAAGCGGTGGCAATTCAACCGGCAGAGCGGCTTTCTGAGCCAGATGTGGATCGGCGACAAGAAGCAGTTGCTGACCCCACTGCGGGATCAGTTCACCCGGGCCCCGCTGGACAACGACATCGGCGTCAGCGAAGCCACTCGGATCGACCCAAACGCCTGGGTCGAACGCTGGAAGGCCGCCGGCCACTACCAGGCCGAAGCCGCTCTGCTGCAATGTACCGCTGATACGCTGGCTGACGCCGTCTTGATTACTACCGCTCACGCCTGGCAACACCAGGGCAAGACTTTGTTTATCAGCCGGAAGACCTACCGGATTGACGGCAGCGGTCAGATGGCCATCACAGTCGATGTCGAAGTTGCCAGCGACACCCCGCACCCGGCACGGATCGGCCTGAACTGCCAGCTGGCCCAGGTTGCCGAACGGGTTAACTGGCTGGGCCTGGGCCCTCAGGAAAACTACCCAGACCGTTTGACGGCTGCCTGCTTTGACCGGTGGGACTTACCGTTGAGCGATATGTACACTCCATACGTCTTTCCGTCCGAAAACGGCCTGCGGTGCGGCACCAGAGAACTGAACTATGGCCCGCACCAGTGGCGCGGTGACTTTCAATTCAACATCAGCCGGTACTCCCAGCAGCAGTTGATGGAAACCAGCCACCGCCACCTGCTGCACGCCGAAGAAGGGACGTGGTTGAACATCGACGGCTTTCACATGGGCATCGGCGGCGACGACTCATGGAGCCCGAGCGTTAGCGCTGAATTCCAGTTGAGCGCCGGCCGGTACCACTACCAGTTGGTTTGGTGCCAGAAG

• To produce the protein from this DNA sequence, I would use a cell-dependent expression system based on bacterial transformation and expression. In this approach, This gene is then placed into an expression cassette with the necessary regulatory elements so it can be used by a biological system.
• To produce the protein, I would use a cell-dependent expression system through bacterial cloning. The designed DNA sequence is inserted into a plasmid and introduced into a bacterial host by transformation. Inside the cell, the gene is transcribed into mRNA under the control of the selected promoter. The mRNA is then translated by ribosomes, which read the codons starting at the start codon and assemble the corresponding amino acids into the lactase protein. This approach follows the natural flow of genetic information (DNA to RNA to protein) and allows controlled production of the enzyme in living cells.

• A lactose-inducible promoter was selected to enable conditional expression of lactase in response to lactose availability in the gut. The PlacA promoter region was extracted from the Lactococcus lactis lac operon upstream of the native ribosome binding site, with preserving lactose-responsive regulation.

AATCGTCGTTTTTTGTTCATATGAAGACTTTCTTTCATAAAGTAATTTTTTTCCAAAGATAATTCTCTTT
TAATTGTATCATAAAAGATAATATTTTCAAGGTAAAACAAACAATTTCAAACAAAAACAAACGTTAGATG
ATGAAATAAGAACAGAGGATTGACGTATATTAGCTTAGGTCAGATTTTGTATAAGACGAAAATAAAGTAG
GACCTCTTAATCAGTAAGTTATAGAAAGTAAAAGACTTTTGTAATACCTGAATAGATATTTCACGTCCAT
TTTGTGATGGATTAAATGAACAAAAATGAACAATAATTTAACGGTGTTATCTATTTTTTAAAAAAACAAA
TAAAAAAAAACAAAAAATTAACAAAAATAGTTGCGTTTTGTTTGAATGTTTGATATCATATAAACAAAGA
AATGATGAAAACGTTATCTTGAACATTTTGCAAAATATTTTCTACTTCTACGTAGCATTTCTTTTTAAAA
TTTAGGAGGTAGTCCAA

• For the RBS, I chose to keep the native Lactococcus lactis ribosome binding site (RBS) derived from the lacA operon which is the region immediately upstream of the coding sequence (CDS) and preserved its original spacer length to ensure efficient translation initiation in the probiotic host. Maintaining native RBS spacing is critical in Gram-positive bacteria, as ribosome binding and translation efficiency are highly sensitive to the distance between the Shine–Dalgarno sequence and the start codon.

• the RBS sequence is as follow:

AGGAGGTAGTCCAA

• I selected the transcription terminator from the tpi gene of Lactococcus lactis, a highly expressed native housekeeping gene, to ensure efficient and reliable transcription termination in the probiotic host. While two related annotations are present in GenBank for this region, both correspond to the same rho-independent transcription terminator. Therefore, I chose the complete annotated terminator region (positions 958–988), which includes both the inverted repeat and the downstream poly-T tract, to ensure proper formation of the termination hairpin and robust termination of transcription.

• A transcription terminator was included downstream of the lactase coding sequence to ensure proper termination of transcription. This prevents transcriptional read-through into adjacent sequences and improves the stability and predictability of gene expression, independent of promoter regulation.

• ATG used as start codon and AAG as stop codon

• From the selected elements, I built a linear expression cassette in Benchling containing a lactose-regulated promoter, native LAB ribosome binding site, codon-optimized lacZ, and a native transcription terminator. I exported this sequence as a FASTA file. Cassette_link_to_Benchling

• When I first uploaded my expression cassette FASTA file to Twist Bioscience, I encountered an initial error related to the FASTA header name. The header exceeded the maximum allowed length (32 characters), which caused the sequence to be rejected. I fixed this issue by shortening the header name and re-uploading the file. After this correction, the sequence was accepted for further analysis.

However, after re-uploading the corrected file, additional synthesis warnings appeared. These warnings were related to large GC content variation, repetitive regions, and overall sequence complexity. These issues are mainly due to the codon-optimized lacZ gene and the presence of multiple regulatory elements such as the ribosome binding site and transcription terminator. Twist flagged these features as potential manufacturability risks. Unfortunately, I was not able to resolve these additional issues at this stage. Fixing them would have required re-optimizing the enzyme sequence, possibly changing the host organism for codon optimization, and redesigning the regulatory architecture of the cassette. Due to time constraints and because this assignment focuses on learning the design and ordering workflow rather than producing a synthesis-ready construct, I chose not to redesign the sequence further.

For this exercise, I proceeded by selecting a Twist clonal vector (pTwist Amp High Copy) to complete the plasmid design. Although the insert sequence still contained manufacturability warnings. However, In a real DNA synthesis order, additional sequence optimization would be required to reduce GC content extremes and repetitive regions to meet synthesis constraints.

DNA read:

• I would want to sequence DNA used for digital data storage. In my knowledge, this technology enables the storage of digital information such as text, images, or files by encoding them into DNA sequences instead of being stored on hard drives. DNA is extremely stable and can store a huge amount of information in a very small space, which makes it interesting for long-term data storage. Reading this DNA by sequencing is necessary to retrieve the stored information and check that the data has not been damaged or changed over time.
• For this porpose, I would use Illumina sequencing because it is very accurate and well suited for reading short DNA fragments, which is how DNA data storage is usually organized. this strategy can be performed following 4 crusial steps: Image_adress

1. Generation This method is a second-generation sequencing technology. It sequences millions of short DNA fragments in parallel, which makes it fast and reliable, but it cannot read very long DNA molecules in one piece.

2. Input and preparation The input is DNA that contains the encoded digital data. To prepare it: The DNA is fragmented into short pieces, Adapters are added to both ends of the fragments, The fragments are amplified using PCR, The prepared DNA is loaded onto a flow cell

3. How the technology reads DNA (base calling) Each DNA fragment is copied one base at a time using fluorescently labeled nucleotides. A camera records the color added at each step, and the machine translates these signals into DNA letters (A, T, C, G).

4. Output The output is a large number of short DNA sequence reads saved as digital files. These reads are then assembled and decoded to recover the original stored data.

DNA write:

• I am particularly interested in the genes in human genomic DNA related to pharmacogenomics and pharmacogenetics. These fields study how genetic variation affects how people respond to drugs. So, I would want to synthesize genes encoding drug-metabolizing enzymes, like human cytochrome P450 enzymes. Since, these genes are central to pharmacogenetics as variations in them strongly influence how drugs are processed in the body. Synthesizing these genes allows them to be studied, expressed, and tested in controlled systems.

• So in order to synthetizing them , I would use chemical DNA synthesis combined with gene assembly, which is the standard approach used by commercial DNA synthesis companies.

1. Essential steps

• DNA synthesis starts with the digital design of the DNA sequence. This is followed by the chemical synthesis of short oligonucleotides, which are then assembled into full-length genes (for example, using Gibson Assembly). The synthesized genes are cloned into plasmids and finally sequence-verified to confirm their accuracy before use.

• This DNA synthesis method is easy to use and works well for many projects.However, it can sometimes make mistakes during the process. Parts of DNA that have lots of G and C letters or repeated sequences are harder to make. Very long DNA pieces also need to be built from many shorter fragments, which can be tricky and may cause errors.

DNA Edit:

• I would want to edit DNA in human cell lines used for drug testing, focusing on genes that affect how drugs work. Changing these genes helps researchers see how different genetic variants influence drug effects and side effects, which is useful in pharmacogenomics.

• The modification can be realised by CRISPR for editing because it allows precise and programmable changes to DNA. this stratigy works by using a guide RNA to find a specific DNA sequence. The Cas enzyme then makes a cut or nick, and the cell repairs it, introducing the change we want.

• To use CRISPR, you need to design guide RNAs, prepare the CRISPR components (DNA, RNA, or protein), deliver them into cells, and then check which cells were correctly edited.

• However, there are some limitations, like different editing efficiencies depending on cell type, and ethical or regulatory concerns when working with human cells.

in this homework, AI ChatGPT assisted me in organizing and clearly articulating my answers and descriptions, ensuring that the content is well-structured and easy to understand.

Sources:

• Ahmad, E., Mahapatra, V., M, V. V., & Nagaraja, V. (2022). Intrinsic and Rho-dependent termination cooperate for efficient transcription termination at 3’ untranslated regions (p. 2022.07.21.500918). bioRxiv. https://doi.org/10.1101/2022.07.21.500918
• Amin, A. A., Olama, Z. A., & Ali, S. M. (2023). Characterization of an isolated lactase enzyme produced by Bacillus licheniformis ALSZ2 as a potential pharmaceutical supplement for lactose intolerance. Frontiers in Microbiology, 14, 1180463. https://doi.org/10.3389/fmicb.2023.1180463 Bioinformatic Tools | VectorBuilder. (n.d.). Retrieved February 17, 2026, from https://en.vectorbuilder.com/tool/overview.html
• Coenen, T. M. M., Bertens, A. M. C., de Hoog, S. C. M., & Verspeek-Rip, C. M. (2000). Safety evaluation of a lactase enzyme preparation derived from Kluyveromyces lactis. Food and Chemical Toxicology, 38(8), 671–677. https://doi.org/10.1016/S0278-6915(00)00053-3
• De Jesus, L. C. L., Aburjaile, F. F., Sousa, T. D. J., Felice, A. G., Soares, S. D. C., Alcantara, L. C. J., & Azevedo, V. A. D. C. (2022). Genomic Characterization of Lactobacillus delbrueckii Strains with Probiotics Properties. Frontiers in Bioinformatics, 2, 912795. https://doi.org/10.3389/fbinf.2022.912795
• de Vrese, M., Stegelmann, A., Richter, B., Fenselau, S., Laue, C., & Schrezenmeir, J. (2001). Probiotics—Compensation for lactase insufficiency123. The American Journal of Clinical Nutrition, 73(2), 421s–429s. https://doi.org/10.1093/ajcn/73.2.421s
• How can I find the promoter sequence of a gene on NCBI? (n.d.). ResearchGate. Retrieved February 17, 2026, from https://www.researchgate.net/post/How_can_I_find_the_promoter_sequence_of_a_gene_on_NCBI
• Lactase—An overview | ScienceDirect Topics. (n.d.). Retrieved February 17, 2026, from https://www.sciencedirect.com/topics/biochemistry-genetics-and-molecular-biology/lactase
• Reverse Translate. (n.d.). Retrieved February 17, 2026, from https://www.bioinformatics.org/sms2/rev_trans.html
• Saqib, S., Akram, A., Halim, S. A., & Tassaduq, R. (2017). Sources of β-galactosidase and its applications in food industry. 3 Biotech, 7(1), 79. https://doi.org/10.1007/s13205-017-0645-5

Week 3 HW: Lab Automation

I created two different agar art designs using two Arabic calligraphy styles. For the first design, I used a simple calligraphy style and created it directly using Python scripting in a Google Colab notebook. For the second design, I used the Opentrons Automation Art interface to design the calligraphy and obtain the coordinates.

I used the Google Gemini AI tool in Colab to understand the logic of the example Opentrons scripts provided in the lab. It helped me understand how coordinates, loops, and pipetting commands work. I also used Gemini AI to help identify and correct mistakes in my Python script, such as indentation errors. I reviewed the suggestions and edited the final code myself.

1. Featured Article: Automated Assembly of Programmable RNA-Based Sensors

The research aimed to solve the challenge of rapidly designing and building large libraries of RNA sensors that can “sense” specific viral RNA signatures. These sensors are crucial for diagnostic applications and understanding RNA-protein interactions. The authors focused on the biological validation of these sensors in both in vivo (bacteria) and cell-free systems.

They used the following lab automation:

• Hardware: Hamilton Microlab STAR liquid-handling workstation.
• Software: Custom Python scripts integrated with the liquid handler’s control software to manage complex plate layouts and reaction conditions.

The researchers used the automated system as a tool to facilitate:

• High-Throughput Plasmid Assembly: The authors needed to construct 144 unique plasmids encoding different riboregulator designs. Doing this manually would be prone to pipetting errors and extremely time-consuming.
• Library Preparation: Automation was used to prepare DNA libraries and reaction mixes for cell-free protein synthesis assays, ensuring consistent reagent volumes across hundreds of samples.
• Normalization and Dilution: The Hamilton system handled the precise normalization of DNA concentrations across plates, which is critical for accurate comparative screening of sensor performance.

The study successfully identified several high-performing RNA sensors capable of detecting viral targets. The use of automation allowed the team to scale their construction phase by nearly 10-fold compared to manual workflows, enabling them to test a much wider range of biological designs than previously possible. For understanding the content of this artical and which type of Lab automation the authors used in their research , i used the AI tool “SCISPACE”.

2. Final project Lab Automation:

My final project focuses on developing an in silico model of a lactose-responsive probiotic that produces lactase only when lactose is present. The physical implementation of this model would allow laboratory automation to verify its predicted results through experimental tests. A liquid-handling robot such as Opentrons could be used to prepare a multi-well plate containing a gradient of lactose concentrations. The robot would then inoculate each well with the engineered probiotic strain, and perform timed sampling to measure lactase activity or reporter output. The automated workflow enables scientists to perform systematic and repeatable tests on lactose responses of the genetic circuit. This helps them match their experimental results with their computer-based model. The project currently exists as a computational project which will use automation as a future extension of the project which does not require automation for its current research activities.

1st Idea: In-Silico Model of an Engineered Probiotic Producing Lactase in Response to Lactose

• Problematic:

Many people cannot digest lactose because they lack enough lactase in their intestine. A possible solution is to use engineered probiotics that produce lactase only when lactose is present. Before building such probiotics in the laboratory, it is important to understand how the genetic system would behave. So, without computational modeling, designing these systems requires trial-and-error experiments that are slow and expensive. There is a need for a simple computational model that can predict how a lactose-responsive genetic circuit would control lactase production over time. Image_ref

• Objectives:

The project is based on a lactose-responsive genetic cassette, who’s dynamic behavior is modeled as a genetic circuit in silico. The objectives of this project are:

1. –> To build an in-silico model of a lactose-responsive genetic circuit.
2. –> To simulate how lactose stimulate lactase production.
3. –> To study how changing key parameters affects lactose degradation.
4. –> To explore system behavior completely in silico.

• Project Description:

The project develops a purely computational model of an engineered probiotic strain. The model is based on a lactose-responsive genetic cassette, whose dynamic behavior is represented in silico as a genetic circuit: – a lacA promoter, operator, and native RBS from Lactococcus lactis for lactose sensing and regulation, – the lacZ gene from Escherichia coli K-12, encoding β-galactosidase (lactase).

I used an AI tool (ChatGPT) to guide me about the repression mechanism I should use, and its response was as follows: To ensure realistic behavior in the model, the lacA promoter includes a native operator, normally repressed by a LacR-like protein in Lactococcus lactis. In the simulation, a repression term is included to prevent unnecessary accumulation of lacZ (lactase) when lactose is absent. The model simulates how the presence of lactose activates the promoter, leading to lactase production, and how this enzyme then degrades lactose over time. No DNA construction or wet-lab experiments are performed. All behavior is represented mathematically and simulated using a computer.

• Steps to Achieve the Project:

1. Define simplified biological assumptions (single strain, constant environment).
2. Represent lactose as the input signal.
3. Model promoter activation based on lactose concentration.
4. Model lactase production and degradation over time.
5. Model lactose degradation by lactase.
6. Run simulations to observe system behavior.
7. Change parameters to study different scenarios.

• Limitations:

1. The model does not include other gut microbes.
2. The gut environment is assumed constant.
3. Results are predictive, not experimentally validated.

2nd Idea:Engineering an E. coli Reporter Strain to Monitor Protein Aging During Heterologous Expression Using a Fluorescent Timer Protein

• Problematic:

Escherichia coli BL21(DE3) is one of the most widely used hosts for heterologous protein expression in research and biotechnology. Although protein expression levels can be easily measured, there are very limited tools to determine how long the expressed protein molecules have persisted inside the cell. During prolonged induction, proteins may accumulate, age, misfold, or lose functionality, even when expression appears successful. Most current methods detect protein quality only after purification, making optimization of expression conditions slow and inefficient. So, there is a need for a genetically encoded reporter system that can estimate protein aging in living cells during expression. Image_ref

• Objectives:

This project is based on a fluorescent timer protein–based reporter system integrated into a heterologous protein expression strain. The objectives are:

1. –> To engineer a reporter strain capable of estimating protein age in vivo.
2. –> To use a fluorescent timer protein to distinguish newly synthesized and older proteins.
3. –> To monitor protein aging during prolonged heterologous expression.
4. –> To provide a practical tool for optimizing protein expression conditions.

• Project Description:

The project focuses on the genetic engineering of a protein expression strain of E. coli BL21(DE3). The reporter system is based on a genetic fusion between:

• a protein of interest (POI) expressed under the T7 promoter, and
• a fluorescent timer protein whose emission spectrum changes over time after synthesis. The genetic construct consists of:
• a T7 promoter and ribosome binding site,
• the gene encoding the protein of interest,
• a flexible linker sequence,
• the fluorescent timer protein gene,
• a transcriptional terminator. Image_ref

After induction, newly synthesized POI–timer fusion proteins initially emit one fluorescent signal. As time progresses, the timer protein matures and shifts to a second fluorescent signal. The ratio of the two fluorescence signals provides an estimate of the age distribution of the expressed protein population.

I used AI tool (ChatGPT) version to refine questions related to the necessary genetic elements required for T7-based heterologous expression in Escherichia coli BL21(DE3) and to determine the appropriate placement of a fluorescent timer gene for monitoring the age of the expressed protein.

• Steps to Achieve the Project:

1. Select a heterologous protein suitable for expression in E. coli.
2. Design a genetic fusion between the protein of interest and a fluorescent timer protein.
3. Clone the fusion construct under a T7 promoter into an expression plasmid.
4. Transform the plasmid into E. coli BL21(DE3).
5. Induce protein expression using IPTG.
6. Monitor fluorescence signals over time using appropriate excitation/emission settings.
7. Calculate fluorescence signal ratios to estimate protein aging.
8. Compare protein aging under different induction times and expression conditions.

• Limitations:

1. Fusion of the timer protein may affect protein folding or function.
2. Protein damage mechanisms are not directly measured.

3rd Idea:Engineering Houseplants for Atmospheric Carbon Monoxide (CO) Capture

• Problematic:

Carbon monoxide (CO) is a toxic gas produced by cars, heaters, and incomplete combustion. It is dangerous for humans, especially in indoor environments. Current solutions such as CO detectors can detect the gas but cannot remove it. Some bacteria naturally use CO as an energy source and convert it into carbon dioxide (CO₂). However, common houseplants cannot metabolize CO. If plants could be engineered to convert CO into CO₂, they could act as natural biological air filters. Image_ref

• Objectives:

The objectives of this project are:

1. –> To engineer a houseplant capable of converting carbon monoxide into carbon dioxide.
2. –> To use microbial genes that naturally perform CO oxidation.
3. –> To ensure the system works safely in oxygen-rich (indoor) environments.
4. –> To allow the produced CO₂ to be reused by the plant’s normal photosynthesis.
5. –> To design a genetically stable and safe indoor plant system.

• Project Description:

This project engineers a plant to express a bacterial enzyme called carbon monoxide dehydrogenase (CODH). This enzyme converts carbon monoxide (CO) into carbon dioxide (CO₂). The CO₂ produced by this reaction is not wasted. Instead, it enters the plant’s natural photosynthetic pathway (Calvin cycle), where it can be fixed into sugars. The plant therefore detoxifies CO while continuing its normal metabolism. The system is designed to work only when CO is present, to avoid unnecessary energy use.

• Genetic Elements for construct design:

1. CO Oxidation Enzymes

The core of the system is the carbon monoxide dehydrogenase (CODH) enzyme, which is responsible for converting carbon monoxide (CO) into carbon dioxide (CO₂). This enzyme is composed of three subunits encoded by the genes coxL, coxM, and coxS. The coxL gene encodes the large catalytic subunit, coxM encodes a subunit involved in electron transfer, and coxS encodes a structural subunit that stabilizes the enzyme complex. These genes originate from Oligotropha carboxidovorans, a bacterium that can oxidize CO in the presence of oxygen, making it suitable for expression in plant cells.

2. Promoter (Gene Expression Control)

To drive the expression of the CODH genes in plant cells, the CaMV 35S promoter is used. This promoter originates from the Cauliflower mosaic virus and is one of the most widely used promoters in plant biotechnology. It enables strong and constitutive gene expression across many plant tissues and is well characterized, making it a reliable choice for this project.

3. Subcellular Targeting Signal

A chloroplast transit peptide is included to ensure that the CODH proteins are transported into the chloroplast after synthesis. This targeting signal is derived from the small subunit of the plant enzyme Rubisco, which naturally localizes to the chloroplast. By directing the CODH enzymes to the chloroplast, the CO₂ produced from CO oxidation is generated close to the photosynthetic machinery, allowing it to be efficiently reused by the plant during photosynthesis.

5. Transcription Terminator

The NOS terminator is used to ensure proper termination of transcription and stable gene expression. This terminator originates from Agrobacterium tumefaciens and is commonly used in plant genetic constructs. Its function is to signal the end of transcription, improving mRNA stability and ensuring reliable expression of the introduced genes.

• Steps to Achieve the Project:

1. Select CO-oxidation genes from aerobic bacteria.
2. Adapt bacterial gene sequences for plant expression (codon optimization).
3. Design a plant expression construct containing:

• Plant promoter
• CODH genes
• Chloroplast targeting signal
• Transcription terminator Image_ref

4. Introduce the construct into the plant genome.
5. Confirm expression of CODH proteins in plant cells.
6. Evaluate CO removal and plant health in controlled conditions.
7. Assess whether produced CO₂ supports normal photosynthesis.

• Limitations:

1. Plant genetic engineering is slow and complex.
2. CO uptake by plants may be limited.
3. CO metabolism efficiency may be low.

in this homework, AI ChatGPT also assisted me in organizing and clearly articulating my answers and descriptions, ensuring that the content is well-structured and easy to understand.

Sources:

• Bidwell, R., & Bebee, G. (2011). Carbon monoxide fixation by plants. Canadian Journal of Botany, 52, 1841–1847. https://doi.org/10.1139/b74-236 Buckley, S. (2024, September 12). Bioengineered Plants Offer Superior Indoor Air Purification. Sustainable Brands. https://sustainablebrands.com/read/bioengineered-plants-indoor-air-purification
• DesignWorksGarage. (n.d.). NCHH. NCHH. Retrieved February 24, 2026, from https://nchh.org/information-and-evidence/learn-about-healthy-housing/health-hazards-prevention-and-solutions/carbon-monoxide/
• Hartl, F. U., Bracher, A., & Hayer-Hartl, M. (2011). Molecular chaperones in protein folding and proteostasis. Nature, 475(7356), 324–332. https://doi.org/10.1038/nature10317
• Heiss, S., Hörmann, A., Tauer, C., Sonnleitner, M., Egger, E., Grabherr, R., & Heinl, S. (2016). Evaluation of novel inducible promoter/repressor systems for recombinant protein expression in Lactobacillus plantarum. Microbial Cell Factories, 15(1), 50. https://doi.org/10.1186/s12934-016-0448-0
• Orina, F., Amukoye, E., Bowyer, C., Chakaya, J., Das, D., Devereux, G., Dobson, R., Dragosits, U., Gray, C., Kiplimo, R., Lesosky, M., Loh, M., Meme, H., Mortimer, K., Ndombi, A., Pearson, C., Price, H., Twigg, M., West, S., & Semple, S. (2024). Household carbon monoxide (CO) concentrations in a large African city: An unquantified public health burden? Environmental Pollution, 351, 124054. https://doi.org/10.1016/j.envpol.2024.124054
• Promoter and terminator considerations for gene integration into the genome of E.coli? (n.d.).ResearchGate. Retrieved February 24, 2026, from https://www.researchgate.net/post/Promoter_and_terminator_considerations_for_gene_integration_into_the_genome_of_Ecoli
• Robson, J. M., Arevalos, N. R., & Green, A. A. (2025). Automated Assembly of Programmable RNA-Based Sensors. bioRxiv, 2025.08.12.669972. https://doi.org/10.1101/2025.08.12.669972

week 04 HW: protein design-part-I

1. Amino Acid Count in 500g Meat: Meat is roughly 20% protein by mass. (Human Nutrition - Protein, Vitamins, Minerals | Britannica, n.d.)
   • 500g meat x 0.20 = 100g protein.
   • Using an average mass of 100 Daltons (Da) per amino acid: 100g / 100 Daltons (or g/mol) = 1 moles of amino acids
   • 1x 6.022 x 10^{23 = 6.022 x 10}23 molecules /1 mole.
2. Why we don’t become cows: When we eat protein, our digestive system breaks it down into individual amino acids. Our body then uses its own DNA information to reassemble those amino acids into human proteins. The information which is coded by the sequence of AA is destroyed, but the building blocks or AA are reused.
3. Why only 20 amino acids: In nature, the use of 20 amino acids is often explained as a “frozen accident” that originated in the early RNA World. This set worked well very early in Earth’s history and then became fixed. These 20 amino acids were good enough to build strong and functional proteins. Even though many other amino acids exist, this small group provides enough variety to perform many functions while remaining simple, stable, and efficient for cells to use. (Doig, 2017)
4. Non-natural amino acids: Yes, scientists can make non-natural (unnatural) amino acids. They do this using chemical methods and special genetic tools that allow new amino acids to be added to proteins. These new amino acids can give proteins new properties that natural amino acids do not have. (Young & Schultz, 2010) For example, A new amino acid could be made by taking a normal amino acid, like alanine, and adding a fluorine atom to its side chain. This fluorinated amino acid would make proteins more stable and less likely to break down, which is useful for drug design. (Adhikari et al., n.d.)
5. Pre-life origins of amino acids: According to Gutiérrez-Preciado, Romero, and Peimbert (2010) Before enzymes and living organisms existed, amino acids were probably formed naturally on early Earth. Energy from lightning, UV light, and volcanic heat helped simple gases react to make amino acids. Some amino acids were also brought to Earth by meteorites and comets. Together, these processes created a “primordial soup” of basic organic molecules. (Amino Acids, Evolution | Learn Science at Scitable, n.d.)
6. D-amino acid α-helix: In nature, L-amino acids form right-handed helices. If you used only D-amino acids, the stereochemistry would be mirrored, resulting in a left-handed $\alpha$-helix. (Zotti et al., n.d.)
7. Additional helices: Yes, additional helical structures besides the standard α-helix can be found in proteins. Studies show that other types of helices occur in many proteins, but they are often overlooked or mistaken for small distortions in α-helices. These helices are especially common in membrane proteins and are found in a significant number of known protein structures.(Vieira-Pires & Morais-Cabral, 2010)
8. Why right-handed helices: because this shape is the most stable for the natural building blocks of life. L-amino acids and D-sugars fit together best in a right-handed twist, which allows strong hydrogen bonds and reduces crowding between atoms. Left-handed helices are usually less stable or hard to form. (Right-Handed Alpha-Helix - an Overview | ScienceDirect Topics, n.d.)
9. β -sheet aggregation: β-sheets tend to aggregate because their edges have exposed hydrogen-bonding groups that easily stick to other β-strands. The main driving forces are hydrogen bonding between strands and the hydrophobic effect, which together make the stacked β-sheet structure very stable and allow fibrils to form.(Gsponer & Vendruscolo, 2006)

1. I chose the protein mCherry because it is a small red fluorescent protein that is easy to visualize and analyze using 3D protein visualization software. Its structure is well known and has a clear β-barrel shape, which makes it easy to study secondary structure, amino acid distribution, and surface features. This makes mCherry a good example protein for learning basic protein sequence and structure analysis.

• The mCherry protein analyzed here is the standard red fluorescent protein and does not function as a fluorescent timer. However, according to the fluorescent protein database (FPbase), mCherry is the parent fluorescent protein for several timer-based reporters, including the medium fluorescent timer planned for my final project. Therefore, mCherry is used in this assignment as a reference protein to understand the structure and sequence properties of fluorescent proteins before working with fluorescent timer variants.

2. I obtained the amino acid sequence of mCherry from the FPbase, which links laboratory fluorescent protein names to biological databases. FPbase provided the UniProt identifier X5DSL3, which is now stored in UniParc (UPI000046F63B) because the UniProtKB entry was removed. And also, the same database provided the genebank identifier for this protein AAV52164, which from where I got the sequence in fasta format.

• This is the obtained sequence :

>AAV52164.1 monomeric red fluorescent protein [synthetic construct]
MVSKGEEDNMAIIKEFMRFKVHMEGSVNGHEFEIEGEGEGRPYEGTQTAKLKVTKGGPLPFAWDILSPQF
MYGSKAYVKHPADIPDYLKLSFPEGFKWERVMNFEDGGVVTVTQDSSLQDGEFIYKVKLRGTNFPSDGPV
MQKKTMGWEASSERMYPEDGALKGEIKQRLKLKDGGHYDAEVKTTYKAKKPVQLPGAYNVNIKLDITSHN
EDYTIVEQYERAEGRHSTGGMDELYK

• The protein sequence is 236 amino acids long and a molecular mass of approximately 26.7 kDa. It has been confirmed at the protein level, although the UniProt entry is currently unreviewed (TrEMBL). Using the provided Colab notebook, I analyzed the amino acid composition of the sequence and found that glycine (G) is the most frequent amino acid, appearing 25 times.

• To identify protein sequence homologs of mCherry, I used the BLAST tool available on UniProt.

• Using the BLAST tool in UniProt, a total of 227 homologous protein sequences were identified for mCherry in the UniProtKB database. Among these results, 13 sequences are reviewed (Swiss-Prot) and 214 are unreviewed (TrEMBL). The homologs show a wide range of sequence identities, from about 23.6% up to 100%, with very low E-values (as low as 4.4 × 10⁻¹⁷⁵), indicating strong evolutionary relatedness.

• Most homologous proteins have sequence lengths between 200 and 400 amino acids, which is similar to mCherry (236 amino acids). Many homologs originate from marine organisms, especially corals and sea anemones such as Porites lobata, Pocillopora meandrina, and Discosoma species, which are known natural sources of GFP-like fluorescent proteins. Some homologs also appear in bacteria and other organisms, reflecting that mCherry is an engineered protein that has been widely introduced into different hosts for research purposes. Overall, these results confirm that mCherry belongs to a well-conserved GFP-like fluorescent protein family with broad biological and biotechnological use.

• The mCherry protein belongs to a known protein family. According to UniProt family and domain analysis, mCherry is part of the green fluorescent protein (GFP)-like family , even though it emits red light. This classification is supported by several databases, including InterPro, Pfam, Gene3D, and PRINTS, all of which identify mCherry as a GFP or GFP-related protein. Proteins in this family share a conserved structure and chromophore-forming mechanism.

3. The structure of the selected protein mCherry is available in the RCSB Protein Data Bank under the PDB ID 2H5Q, titled “Crystal structure of mCherry.” This structure represents the red fluorescent protein mCherry derived from Discosoma species and expressed in Escherichia coli. The structure was solved using X-ray diffraction and was deposited in May 2006 and released in August 2006.

• The quality of this structure is very high. It was solved at a resolution of 1.36 Å, which is much better than the 2.70 Å threshold typically used to define a good-quality structure. Lower resolution values indicate more detailed and accurate atomic positions, so a resolution of 1.36 Å means the structure is very reliable. In addition, the reported R-values (R-work ≈ 0.15 and R-free ≈ 0.19) further support that this is a well-refined and high-quality crystal structure.

• Besides the protein itself, the solved structure contains a modified residue that corresponds to the mature chromophore of mCherry. This chromophore is formed from amino acids within the protein chain and is responsible for fluorescence. No additional ligands, cofactors, or external small molecules are present. The biological assembly is a single monomer, which means that the protein functions as one chain and does not require binding to other protein subunits.

• According to SCOP (Structural Classification of Proteins), mCherry is classified within the fluorescent protein family and the GFP-like superfamily. SCOP groups proteins based on their three-dimensional structure rather than their biological function or expression host. In this classification, mCherry contains a single domain (residues 6–224) that forms the characteristic β-barrel fold shared by GFP-like proteins. This confirms that mCherry belongs to the same structural superfamily as other green and red fluorescent proteins that use a similar fold to support fluorescence.

4. The protein was visualized using cartoon, ribbon, and ball-and-stick representations to examine overall fold and atomic details.

• Coloring by secondary structure shows that mCherry contains many β-sheets about 11 β-sheets and very few α-helices (only 3 helices) . The protein is dominated by a β-barrel fold, which is typical for GFP-like fluorescent proteins.

• Using the PyMOL command line, I colored the hydrophobic residues yellow and the hydrophilic residues red. The resulting structure shows a clear alternating pattern along the β-strands, where hydrophilic side chains face the exterior to interact with the aqueous environment (supported by the presence of surrounding water molecules), while hydrophobic side chains face the interior. This internal hydrophobic core effectively shields the chromophore from the solvent, which is essential for its fluorescence.

• Based on the surface visualization of the mCherry protein (PDB: 2H5Q), the protein does not show any clear holes or binding pockets. The surface is compact and smooth, forming a closed β-barrel structure that surrounds the chromophore inside the protein. Although small bumps and grooves are visible on the surface due to amino acid side chains, there are no deep openings that lead into the protein core. This sealed structure is important for mCherry’s function, because it protects the internal chromophore from water or oxygen that could interfere with fluorescence. The closed surface therefore supports the role of mCherry as a stable fluorescent protein.

C1 Protein Language Modeling

1. Deep Mutational Scans

To analyze how different mutations affect my protein, I used the ESM-2 protein language model to generate a deep mutational scan. The output is shown as a heatmap, where each color represents how favorable or unfavorable a specific mutation is. The score (z-value) reflects how likely the mutation is to be stable: positive values mean the mutation is well tolerated, while negative values suggest the mutation may damage the protein.

To Understand the Heatmap Colors, these are some exmples:

• The darkest color (black) represents the most harmful mutations. For example, the mutation at position 92 to Cysteine (C) has a very low score (z = −5.01). This position is buried deep inside the protein. Changing it to cysteine is predicted to strongly disrupt the protein, likely causing misfolding or aggregation.
• The dark blue color represents very risky mutations. An example is position 180 mutated to Proline (P) with a score of z = −3.08. This residue lies in a β-strand. Proline is known to break regular protein structures, so inserting it here would likely distort or break the β-barrel.
• The green color indicates neutral mutations. For example, position 183 mutated to Threonine (T) has a score of z = 0, meaning the model predicts little to no effect on protein stability.
• The yellow color represents favorable mutations. At position 45 mutated to Valine (V), the score is z = 3.04, suggesting this mutation may slightly improve protein stability compared to the original amino acid.

When looking at the entire heatmap, many positions appear as vertical dark bands. These positions do not tolerate most mutations and are therefore highly conserved. These residues usually form the hydrophobic core of the protein and point inward to build the β-barrel structure. Because mCherry has a tightly sealed β-barrel, mutations in these regions can disrupt proper folding or destabilize the barrel. If the β-barrel is damaged or becomes leaky, the chromophore inside can no longer be protected, which would stop the protein from fluorescing. So, this explains why mutations in these regions are strongly disfavored by the model.

My protein of interest is the Medium-FT variant, which is related to my final project and works as a protein “aging timer.” This behavior is controlled by specific mutations that change the chromophore chemistry without breaking the overall protein structure. To explore the functional mutations in the parent protein mCherry (PDB: 2H5Q, the one I used to represent the heatmap), I focused on two important mutations: K69R, and A224S (F. V. Subach et al., 2009; O. M. Subach et al., 2022). So as indicated in the heatmap, they showed positive scores (z = 0.75 ; z = 1.08) respectively.

Both mutations appear as light green to yellow on the heatmap, meaning they are well tolerated. This confirms that these changes do not disrupt the β-barrel or overall stability. they adjust the protein’s function by slowing down fluorescence maturation while keeping the main structure intact.

2. Latent Space Analysis

To perform latent space analysis, I used the provided dataset of protein sequences from the SCOP database and generated numerical embeddings for each sequence using the ESM-2 protein language model, which results in a three-dimensional map where each point represents one protein.

When analyzing the resulting map, proteins do not appear randomly distributed. Instead, they form local neighborhoods where nearby points correspond to proteins with similar structural properties. These neighborhoods approximate similarities in protein fold and secondary structure rather than biological function. This shows that the language model organizes proteins based on shared “structural rules,” such as how alpha helices and beta sheets are arranged, even when the proteins come from different organisms or have different functions.

For example, the protein d2cw3a1 a.2.11.0 (A:4–90) from Perkinsus marinus has three closest neighbors that come from very different organisms, including Escherichia coli and cow. These neighboring proteins also have very different biological functions.

My protein of interest, mCherry (PDB: 2H5Q) which is represented by bleu dot, is located in a neighborhood dominated by proteins rich in β-sheet structures. Its closest neighbors include proteins such as the β-propeller domain of the enzyme PepX, the β-barrel domain of the chaperone protein Sis1, and other β-sheet–containing domains like transferrin-binding protein and latexin. Although these proteins perform very different biological roles, they share similar β-sheet-based structural architectures. The close proximity of mCherry to these proteins confirms that the ESM-2 model groups proteins based on structural similarity, correctly placing mCherry among other β-sheet and β-barrel-like proteins in the latent space.

C2. Protein Folding

The predicted coordinates matched the original structure very well. The overall shape, especially the β-barrel structure, was preserved, and the folding pattern looked almost identical. This shows that ESMFold can accurately predict the structure of mCherry from its amino acid sequence.

Next, I changed the protein sequence by introducing several mutations, including small amino acid changes and changes spread across the sequence. After folding the mutated sequence with ESMFold, the structure showed noticeable changes compared to the original protein. While the general β-barrel shape was still present, some regions were slightly distorted. This indicates that mCherry is partly resilient to mutations, but too many or poorly placed mutations can affect proper folding and reduce structural stability.

C3. Protein Generation

I used ProteinMPNN to do inverse folding on the mCherry protein (PDB: 2H5Q). I used the default settings and turned off the homomer option because this protein has only one chain. ProteinMPNN uses the 3D shape of the protein and suggests new amino acid sequences that can keep the same shape.

The output includes a probability heatmap, which shows the model’s confidence for each amino acid at every position in the sequence. In the heatmap, bright colors (yellow/green) indicate amino acids that are highly preferred at a specific position, while dark colors (blue/purple) indicate unlikely choices. Some positions show a strong preference for one amino acid, meaning they are important for maintaining the protein structure. Other positions show more flexibility, suggesting they can tolerate different amino acids without disrupting the fold.

ProteinMPNN generated a new sequence candidate with a sequence recovery of about 47.93 %, meaning nearly half of the amino acids are identical to the original mCherry sequence. The designed sequence received a lower score (0.8107) compared to the native sequence score (1.3913). Because lower scores indicate a better statistical fit to the backbone, this suggests that the designed sequence is predicted to be highly compatible and stable for the 11-stranded β-barrel structure of mCherry.

• The native protein sequence and its score are shown below:

>2H5Q, score=1.3913, fixed_chains=[], designed_chains=['A'], model_name=v_48_020
NMAIIKEFMRFKVHMEGSVNGHEFEIEGEGEGRPYEGTQTAKLKVTKGGPLPFAWDILSPQFXXXSKAYVKHPADIPDYLKLSFPEGFKWERVMNFEDGGVVTVTQDSSLQDGEFIYKVKLRGTNFPSDGPVMQKKTMGWEASSERMYPEDGALKGEIKQRLKLKDGGHYDAEVKTTYKAKKPVQLPGAYNVNIKLDITSHNEDYTIVEQYERAEGRHST

• The newly generated protein sequence and its evaluation metrics are shown below:

>T=0.1, sample=0, score=0.8107, seq_recovery=0.4793
VDDIVKPVQKYTVNLDGSVNGHKFKIKGEGIGTPYEGKYEVDLEVTEGGPLPFSFDILAPLFXXXAQQFTKYPADIPDYVKQAFPEGYTEERTADYEDGGKLKSTKTVTLKDGVVVQEIEADGSNFPADGPVMTKKTAGWEPVVWHCYPKDGALYCEADAALKLKDGGTYKAKVTAKIKPNHPVPLPGPFEIDEKLTVTDHNADETKVKLSKEAVARRAS

I attempted to refold the newly designed sequence using ESMFold in order to compare the predicted structure with the original mCherry structure. However, ESMFold requires GPU resources, and GPU access was not available at the time of execution. As a result, a direct structural comparison could not be performed. Despite this limitation, the strong sequence score and conserved structural regions indicate that the designed sequence would likely fold into a structure very similar to the original β-barrel if GPU resources were available.

Gemini AI tools integrated with Google Colab were used to help explain code errors, interpret the generated outputs such as heatmaps, and analyze the latent space by identifying the closest neighboring proteins through distance calculations between my protein and other sequences.

One-Page Proposal

in this homework, AI ChatGPT assisted me in organizing and clearly articulating my answers and descriptions, ensuring that the content is well-structured and easy to understand.

Sources:

• Adhikari, A., Bhattarai, B. R., Aryal, A., Thapa, N., KC, P., Adhikari, A., Maharjan, S., Chanda, P. B., Regmi, B. P., & Parajuli, N. (n.d.). Reprogramming natural proteins using unnatural amino acids. RSC Advances, 11(60), 38126–38145. https://doi.org/10.1039/d1ra07028b
• Chamakura, K. R., Tran, J. S., & Young, R. (2017). MS2 Lysis of Escherichia coli Depends on Host Chaperone DnaJ. Journal of Bacteriology, 199(12), e00058-17. https://doi.org/10.1128/JB.00058-17
• Doig, A. J. (2017). Frozen, but no accident – why the 20 standard amino acids were selected. The FEBS Journal, 284(9), 1296–1305. https://doi.org/10.1111/febs.13982
• Gsponer, J., & Vendruscolo, M. (2006). Theoretical Approaches to Protein Aggregation. Protein & Peptide Letters, 13(3), 287–293. https://doi.org/10.2174/092986606775338407
• Human nutrition—Protein, Vitamins, Minerals | Britannica. (n.d.). Retrieved February 27, 2026, from https://www.britannica.com/science/human-nutrition/Meat-fish-and-eggs
• Mezhyrova, J., Martin, J., Börnsen, C., Dötsch, V., Frangakis, A. S., Morgner, N., & Bernhard, F. (2023). In vitro characterization of the phage lysis protein MS2-L. Microbiome Research Reports, 2(4), 28. https://doi.org/10.20517/mrr.2023.28
• Phage Therapy: Past, Present and Future. (n.d.). ASM.Org. Retrieved March 2, 2026, from https://asm.org:443/articles/2022/august/phage-therapy-past,-present-and-future
• Right-Handed Alpha-Helix—An overview | ScienceDirect Topics. (n.d.). Retrieved February 27, 2026, from https://www.sciencedirect.com/topics/chemistry/right-handed-alpha-helix
• Strathdee, S. A., Hatfull, G. F., Mutalik, V. K., & Schooley, R. T. (2023). Phage therapy: From biological mechanisms to future directions. Cell, 186(1), 17–31. https://doi.org/10.1016/j.cell.2022.11.017
• Subach, F. V., Subach, O. M., Gundorov, I. S., Morozova, K. S., Piatkevich, K. D., Cuervo, A. M., & Verkhusha, V. V. (2009). Monomeric fluorescent timers that change color from blue to red report on cellular trafficking. Nature Chemical Biology, 5(2), 118–126. https://doi.org/10.1038/nchembio.138
• Subach, O. M., Tashkeev, A., Vlaskina, A. V., Petrenko, D. E., Gaivoronskii, F. A., Nikolaeva, A. Y., Ivashkina, O. I., Anokhin, K. V., Popov, V. O., Boyko, K. M., & Subach, F. V. (2022). The mRubyFT Protein, Genetically Encoded Blue-to-Red Fluorescent Timer. International Journal of Molecular Sciences, 23(6), 3208. https://doi.org/10.3390/ijms23063208
• Vieira-Pires, R. S., & Morais-Cabral, J. H. (2010). 310 helices in channels and other membrane proteins. The Journal of General Physiology, 136(6), 585–592. https://doi.org/10.1085/jgp.201010508
• Young, T. S., & Schultz, P. G. (2010). Beyond the Canonical 20 Amino Acids: Expanding the Genetic Lexicon. The Journal of Biological Chemistry, 285(15), 11039–11044. https://doi.org/10.1074/jbc.R109.091306

week-05-HW-protein-design-part-II

Part 1: Generate Binders with PepMLM

For the PepMLM analysis, the amino acid sequence of the normal Superoxide Dismutase 1 protein was obtained from UniProt using the accession number P00441. To simulate the disease-associated variant, the A4V mutation was then manually introduced into the sequence to generate the mutant form of the protein used for the peptide design experiments. This mutation corresponds to the substitution of alanine by valine at position 4 of the protein sequence.

• Original Superoxide dismutase 1 (SOD1) sequence from Uniprot:

sp|P00441|SODC_HUMAN Superoxide dismutase [Cu-Zn] OS=Homo sapiens OX=9606 GN=SOD1 PE=1 SV=2

MATKAVCVLKGDGPVQGIINFEQKESNGPVKVWGSIKGLTEGLHGFHVHEFGDNTAGCTS AGPHFNPLSRKHGGPKDEERHVGDLGNVTADKDGVADVSIEDSVISLSGDHCIIGRTLVV HEKADDLGKGGNEESTKTGNAGSRLACGVIGIAQ

• Mutant A4V sequence (A4V means Alanine → Valine at position 4):

sp|P00441|SODC_HUMAN Superoxide dismutase [Cu-Zn] OS=Homo sapiens OX=9606 GN=SOD1 PE=1 SV=2

MATKVVCVLKGDGPVQGIINFEQKESNGPVKVWGSIKGLTEGLHGFHVHEFGDNTAGCTS AGPHFNPLSRKHGGPKDEERHVGDLGNVTADKDGVADVSIEDSVISLSGDHCIIGRTLVV HEKADDLGKGGNEESTKTGNAGSRLACGVIGIAQ

After running the PepMLM model using mutant A4V sequence, I obtained four candidate peptide binder sequences. In the generated results, each sequence ended with the amino acid symbol ‘X’, which represents an undefined residue predicted by the model.

To proceed with the structural analysis and fold the predicted binders with the mutated SOD1 protein, I needed to assign a specific amino acid at this position. For this reason, I replaced the ‘X’ residue with alanine (A) in each sequence. I chose alanine because it is a small and neutral amino acid that generally has minimal effects on protein structure and interactions. This allowed me to obtain complete peptide sequences that could be used for the subsequent folding and interaction prediction with the A4V mutant SOD1 protein.

Part 2: Evaluating Binders with AlphaFold3

The predicted peptide binders were evaluated using structural modeling. Each peptide was folded together with the A4V mutant of SOD1 to evaluate the potential protein–peptide interactions.

• The known binder (FLYRWLPSRRGG) showed a relatively low binding confidence, with an ipTM score of 0.28. The peptide binds mainly on the surface of the SOD1 β-barrel, close to the electrostatic loop and the zinc-binding loop. It does not bind near the N-terminus, where the A4V mutation is located, and it also does not interact with the dimer interface. The peptide remains mostly surface-bound rather than buried inside the protein structure. Several residues help stabilize the interaction. For example, Trp5 and Tyr3 can form aromatic contacts with the protein surface, while Arg8 and Arg9 may form hydrogen bonds with nearby residues of SOD1. However, the peptide does not form a strong or compact binding interface, which suggests that the interaction may be weak or transient.

• Binder 1 (WLSPAAGVEWKA) showed a clear improvement compared with the known binder, with an ipTM score of 0.39. The peptide binds on a hydrophobic groove on the surface of the SOD1 β-barrel. In this interaction, Trp1 acts as an important anchoring residue, helping the peptide attach to a hydrophobic pocket on the protein surface. Other residues such as Ser3 and Pro4 help position the peptide backbone against the protein surface. In addition, Glu9 forms stabilizing hydrogen bonds with nearby residues on SOD1. Because of these interactions, the peptide forms a more compact and organized binding conformation than the known binder. Although the peptide still binds away from the A4V mutation site, the higher ipTM score and the stronger interaction network suggest that Binder 1 may represent a more promising peptide candidate.

• Binder 2 (WHYYVVVVRHWA) showed a moderate interaction with SOD1, with an ipTM score of 0.33, which is higher than the known binder but lower than Binder 1. The peptide binds on a surface patch of the SOD1 β-barrel region. Several residues appear to contribute to this interaction. Trp1 participates in both hydrogen bonding and aromatic interactions with the protein surface, helping to anchor the peptide. Tyr3 and Arg9 also participate in hydrogen bonding that stabilizes the peptide orientation. In addition, the terminal residue Ala12 contributes to stabilizing the peptide backbone through hydrogen bonding with the protein surface. Compared with the known binder, Binder 2 shows a more localized and organized binding mode, although the peptide still binds mainly on the surface of the protein rather than deeply inside the structure.

• Binder 3 (WRSYVVVVELKA) showed the lowest binding confidence among the designed peptides, with an ipTM score of 0.20, which is even lower than the known binder. The peptide still localizes on the surface of the SOD1 β-barrel, but the interaction appears weak and poorly defined. The interaction is mainly supported by Arg2 and Lys11, which can form hydrogen bonds with residues on the SOD1 surface. In addition, Tyr4 may contribute through aromatic interactions with the protein surface. However, the peptide forms only a limited number of stabilizing contacts, and the interaction appears less stable compared with Binder 1 and Binder 2. These results suggest that Binder 3 may not be a strong candidate for stable binding to the SOD1 mutant.

• Binder 4 (WRYPAVVAAHGA) showed a moderate structural confidence, with an ipTM score of 0.33, similar to Binder 2 and higher than the known binder. The peptide binds on the surface of the SOD1 β-barrel region. Several residues contribute to this interaction. Trp3, Val6, and Gly11 appear to form hydrogen bonds with residues on the SOD1 surface, helping stabilize the interaction. In addition, an internal hydrogen bond between Val6 and His10 helps stabilize the peptide backbone and maintain its conformation. Compared with Binder 3, this peptide forms more defined interactions with the protein surface, which explains its higher predicted binding confidence. Although the peptide still binds away from the A4V mutation site, the interaction appears more organized and stable than the known binder.

To further explore whether peptide length influences binding stability, the same structural analysis was also performed using 11-residue versions of the peptides obtained by removing the final alanine that replaced the unknown residue X. For Binder 1, the ipTM score decreased from 0.39 (12 aa) to 0.27 (11 aa), indicating that the twelfth residue likely helps stabilize the interaction with the SOD1 surface. In contrast, Binder 2 showed a small increase in structural confidence, where the score changed from 0.33 (12 aa) to 0.35 (11 aa), suggesting that the slightly shorter peptide may adopt a somewhat better orientation on the protein surface. Binder 3 showed the strongest negative effect of shortening the peptide, with the score decreasing from 0.20 (12 aa) to 0.13 (11 aa), confirming that this peptide already forms weak interactions and becomes even less stable when shortened. Interestingly, Binder 4 showed the opposite trend, where the 11-residue version reached the highest score of all tested peptides (0.44) compared with 0.33 for the 12-residue version, suggesting that removing the last residue may allow the peptide to adopt a more favorable binding conformation. Overall, these exploratory results suggest that peptide length can influence binding stability, but the effect is sequence-dependent, since shortening the peptide reduced stability for some binders (Binder 1 and Binder 3) while improving it for others (Binder 2 and Binder 4).

Part 3: Evaluate Properties of Generated Peptides in the PeptiVerse

I evaluated the therapeutic properties of all 12-residue peptide binders using PeptiVerse. The known binder FLYRWLPSRRGG showed weak binding affinity (pKd 5.97), good solubility, very low hemolysis probability (0.047), and a positive net charge of 2.76. Among the PepMLM-generated peptides, Binder 1 (WLSPAAGVEWKA) had weak binding affinity (pKd 5.61), excellent solubility, very low hemolysis (0.037), and a near-neutral net charge (-0.24). Binder 2 (WHYYVVVVRHWA) exhibited medium binding affinity (pKd 7.12), fully soluble, non-hemolytic (0.115), and slightly positive net charge (0.93). Binder 3 (WRSYVVVVELKA) showed weak binding (pKd 6.28), soluble, non-hemolytic (0.115), and net charge 0.76. Binder 4 (WRYPAVVAAHGA) had weak binding (pKd 5.22), soluble, non-hemolytic (0.037), and net charge 0.85.

Comparing these properties with the structural confidence from AlphaFold, we see that higher ipTM scores do not always directly match stronger predicted binding. For example, Binder 1 had the highest ipTM (0.39) but only weak predicted binding, while Binder 2 had slightly lower ipTM (0.33) but showed medium predicted binding. All generated peptides are soluble and non-hemolytic, which is favorable for therapeutic use. Considering both structural confidence and predicted properties, Binder 1 (WLSPAAGVEWKA) is the most promising overall: it has the highest structural stability on SOD1, is non-hemolytic, fully soluble, and has a near-neutral charge that may support safe and effective binding in a biological context.

Part 4: Generating Optimized Peptides with moPPIt

Mutations such as the A4V variant can destabilize the structure of Superoxide Dismutase 1, increasing the probability of protein misfolding, dissociation of the dimer, and toxic aggregation, which are processes associated with Amyotrophic Lateral Sclerosis (ALS). For this reason, the design strategy of this step focused on generating short peptides that can bind simultaneously to both monomers at the dimer interface, effectively acting as a molecular bridge that reconnects and stabilizes the two subunits. By reinforcing the interaction between the chains, these peptides may help restore a conformation closer to the native functional state of the SOD1 complex, while reducing the structural instability caused by the mutation.

To do so, several design parameters were selected before generating peptides. The peptide length was fixed at 12 amino acids. The motif position was focused on residues 3–10, meaning the central region of the peptide was encouraged to interact with the target protein. In addition, affinity guidance and solubility optimization were enabled, and hemolysis prediction was considered to reduce potential toxicity. These settings allow the model to design peptides that not only bind the protein but also have better therapeutic properties.

After generating the following peptides, their structures were evaluated using AlphaFold to predict how they interact with SOD1.

• The generated sequences:

• Binder 1 – KRQCEIFNQFMA (ipTM: 0.91) This peptide localizes at the dimer interface of the SOD1 homodimer and acts as a molecular bridge between the two monomers. Several residues of the peptide participate in stabilizing the interaction. Gln3 forms a hydrogen bond with residues on the first monomer, while Cys4 interacts with a cysteine residue on the second monomer. In addition, Asn8 forms multiple hydrogen bonds with residues on Chain A. These multiple contacts allow the peptide to connect both monomers simultaneously, which could help stabilize the dimer structure of SOD1.
• Binder 2 – EKDNKWVITSQF (ipTM: 0.86) This peptide also binds at the dimer interface and connects the two monomers. The interaction is mainly driven by the N-terminal region of the peptide. Glu1 forms several hydrogen bonds with residues on Chain B, creating a strong anchoring point. In addition, Ser10 interacts with residues on Chain A. Through these interactions with both monomers, the peptide may help maintain the stability of the SOD1 dimer.
• Binder 3 – VCQFDYKTLFKK (ipTM: 0.87) This peptide spans the interface between the two monomers, forming stabilizing contacts with both chains. Val1 forms a hydrogen bond with residues on Chain B, while Phe4 interacts with residues on Chain A. These interactions allow the peptide to bridge the two monomers and stabilize the interface region.
• Binder 4 – GQQSLFKTKTLD (ipTM: 0.89) Unlike the previous peptides, this binder attaches to the outer surface of a single SOD1 monomer, particularly near the β-barrel structure. The interaction is mainly driven by residues near the C-terminus of the peptide. Thr10 forms a hydrogen bond with the monomer, while Asp12 forms two hydrogen bonds with residues on Chain A. Lys9 also contributes to stabilization by forming an additional hydrogen bond. This peptide does not bridge the dimer but instead stabilizes the surface structure of the monomer. The four peptides show two different binding strategies:

Three peptides (KRQCEIFNQFMA, EKDNKWVITSQF, and VCQFDYKTLFKK) bind at the dimer interface, where they interact with residues from both monomers. These peptides may help stabilize the SOD1 dimer by acting as a bridge between the two chains. In contrast, GQQSLFKTKTLD binds only to one monomer, specifically on the β-barrel surface. Instead of bridging the two chains, this peptide may stabilize the structure of the individual monomer.

Among the peptides, KRQCEIFNQFMA shows the highest ipTM score (0.91), suggesting the strongest predicted interaction with the protein complex. When comparing the peptides generated by PepMLM and moPPIt, the main difference lies in the design strategy. PepMLM mainly samples possible peptide sequences that could bind to the target protein based on patterns learned from protein sequence data. However, it does not allow the user to control exactly where the peptide should bind on the protein. As a result, the generated peptides are plausible binders, but their binding location and biochemical properties are not specifically optimized.

In contrast, moPPIt enables guided peptide design. In this approach, the user can select specific residues or regions on the protein where the peptide should bind, such as the dimer interface of Superoxide Dismutase 1 or regions near the A4V mutation. The model also optimizes several properties simultaneously, including binding affinity, solubility, hemolysis risk, and motif placement. Because of this multi-objective optimization, moPPIt peptides are designed to better satisfy several therapeutic requirements at the same time.

For this part, unfortunately, I was unable to access the BRD4 Drug Discovery Platform, as the access was not granted to me despite my request.

Option 1 : Improve autofolding and lysis efficiency

The goal of this part was to design mutations in the L-protein in order to improve its function. Two main objectives were considered. The first objective was to improve the autofolding ability of the L-protein so that it can fold correctly without strong dependence on host chaperones. The second objective was to improve the lysis efficiency of the protein by enhancing its ability to form pores in the E. coli membrane and promote faster or more efficient bacterial lysis.

To identify possible mutations, the provided mutation scoring notebook was used. This notebook evaluates all possible amino-acid substitutions in the L protein and assigns a score to each mutation.

After running the notebook, the resulted mutation predictions are presented in the following dataset:

The predicted mutations were compared with the experimental dataset of L-protein mutants provided in the course material. This dataset contains mutations that were experimentally tested and their effect on lysis activity.

The goal of this comparison was to determine whether mutations with high prediction scores correspond to mutations that show improved lysis in experimental studies. This step helps evaluate the reliability of the prediction model. The results of this comparison revealed a limited overlap between predicted beneficial mutations and experimentally tested mutations. Two mutations, C29R and K50I, appeared in both datasets. However, experimental data indicated that these substitutions did not improve lysis activity. This suggests that, while the protein language model captures sequence compatibility, it does not fully predict functional outcomes such as lysis efficiency. For this raison, experimental validation remains essential to confirm computational predictions.

To avoid mutations that could disrupt essential protein functions, sequence conservation analysis was performed. Multiple sequences related to the MS2 L protein obtained from the BLAST results provided in the course folder were uploaded to Clustal Omega and aligned. The conserved regions of the L protein were identified after analyzing the multiple sequence alignment results. Highly conserved residues, which are the same across all sequences, were marked with stars (*) in the alignment output. while colon (:) indicates residues with strongly similar chemical properties. These positions were considered critical for protein function, so mutations at these residues were avoided. The remaining positions, which showed variability among sequences, were classified as non-conserved and were selected as potential sites for mutation. This approach ensured that the chosen mutations would minimize disruption of essential protein structure and function. Mutations were selected using the resulted mutation scoring predictions and evolutionary conservation analysis. Only residues located in non-conserved positions were chosen in order to reduce the risk of disrupting essential protein functions. The selected mutations (F5Q, S9Q, F22S, Y27L, and A45L) -as represented in the following table- are distributed between the N-terminal region, the central region, and the transmembrane domain of the L protein. This distribution allows the exploration of potential effects on protein autofolding and membrane activity, while maintaining the overall structural integrity of the protein.

Because the L gene overlaps with other genes in the MS2 genome, the nucleotide changes corresponding to the selected mutations were checked to ensure that they do not introduce stop codons in the overlapping reading frames. The mutations F5Q and S9Q are located in the region overlapping with the coat protein (CP) gene, near its C-terminal end, while the mutation A45L is located in the region overlapping with the replicase (Rep) gene, near its N-terminal region. For each mutation, the possible codon substitutions were examined and confirmed not to generate stop codons in the overlapping genes. Therefore, these mutations are considered compatible with the genome organization of MS2.

Option 2: Achieve DnaJ independence

Here the goal was to reduce or eliminate the dependence of the L-protein on the host chaperone DnaJ. By designing mutations in the soluble N-terminal domain of the L-protein, i aimed to weaken its interaction with DnaJ while maintaining proper folding. This approach could potentially allow the phage to function even if DnaJ is mutated or absent in the host.

To study the interaction, i used the AlphaFold2-Multimer notebook in ColabFold to co-fold the soluble domain of the L-protein with the full sequence of E. coli DnaJ. The sequences used were:

• DnaJ sequence:

MAKQDYYEILGVSKTAEEREIRKAYKRLAMKYHPDRNQGDKEAEAKFKEIKEAYEVLTDSQKRAAYDQYGHAAFEQGGMGGGGFGGGADFSDIFGDVFGDIFGGGRGRQRAARGADLRYNMELTLEEAVRGVTKEIRIPTLEECDVCHGSGAKPGTQPQTCPTCHGSGQVQMRQGFFAVQQTCPHCQGRGTLIKDPCNKCHGHGRVERSKTLSVKIPAGVDTGDRIRLAGEGEAGEHGAPAGDLYVQVQVKQHPIFEREGNNLYCEVPINFAMAALGGEIEVPTLDGRVKLKVPGETQTGKLFRMRGKGVKSVRGGAQGDLLCRVVVETPVGLNERQKQLLQELQESFGGPTGEHNSPRSKSFFDGVKKFFDDLTR

• The soluble domain of Lysis protein (N terminal Domain):

METRFPQQSQQTPASTNRRRPFKHEDYPCRRQQRSS

The first co-folding run generated five ranked models with the following parameters. Here, pLDDT reflects the confidence in the predicted structure of the L–DnaJ complex, pTM indicates the overall predicted quality of the complex, and ipTM estimates the predicted strength of the interaction between L-protein and DnaJ:

Rank_004_alphafold2_multimer_v3_model_1: pLDDT=79.3 pTM=0.567 ipTM=0.198 Rank_005_alphafold2_multimer_v3_model_2 : pLDDT=78.9 pTM=0.573 ipTM=0.161 Rank_003_alphafold2_multimer_v3_model_3: pLDDT=78.4 pTM=0.572 ipTM=0.265
Rank_001_alphafold2_multimer_v3_model_4: pLDDT=78.9 pTM=0.583 ipTM=0.373
Rank_002_alphafold2_multimer_v3_model_5: pLDDT=77.9 pTM=0.572 ipTM=0.287

• Plots for L_DnaJ_complex

By comparing the different predicted models, the fourth model was identified as the best-ranked model because it showed the highest ipTM score, indicating the strongest predicted interaction between the L-protein and DnaJ. Using PyMol I Analysed the predicted L–DnaJ complex of the best predicted model (rank 4) and the results revealed multiple interaction residues located in the N-terminal region of the L protein as summarized in the following table:

Key residues such as Arg4, Thr3, Pro6, Phe5, Arg18, Lys23, His24, and Tyr27 were found to interact with several residues of DnaJ, including Asp116, Leu117, Arg118, Glu122, and Glu266. These interactions involve a combination of electrostatic, hydrophobic, and hydrogen-bond contacts. Residues forming multiple contacts were considered potential targets for mutagenesis aimed at reducing the dependence of the L protein on the DnaJ chaperone.

Two hydrophobic residues (Pro6 and Phe22), two positively charged residues (Arg4 and Arg18), and two negatively charged residues (Glu2 and Asp26) were selected for mutational analysis. These residues participate in multiple contacts with DnaJ and represent different physicochemical interaction types involved in stabilizing the L–DnaJ interface.

To evaluate the contribution of different interaction types at the L–DnaJ interface, selected residues were substituted with alanine using an alanine-scanning approach in order to remove their side-chain interactions while minimizing structural perturbation.

The resulting N-terminal sequence of the lysis protein was used to re-predict the interaction with the DnaJ protein in order to evaluate whether the introduced mutations could reduce the dependence of the lysis protein on the host chaperone:

MATAFAQQSQQTPASTNARRPAKHEAYPCRRQQRSS

The mutated L-protein was co-folded again with DnaJ using AlphaFold2-Multimer. The five ranked models obtained were:

When we compared the new models with the wild-type complex, we can see clearly that the ipTM values were slightly lower. In the wild-type prediction, the best model showed an ipTM value of 0.373, while after mutation the highest ipTM value decreased to 0.291. Since ipTM reflects the predicted strength of interaction between two proteins, this decrease suggests that the interaction between the L-protein and DnaJ became weaker after the mutations were introduced. This reduction is consistent with the mutation strategy, where several key residues involved in hydrophobic and electrostatic contacts were replaced with alanine in order to remove their side-chain interactions.

Despite these changes, the overall structural confidence of the models (pLDDT values were 78.9 to 78.6) remained similar to the wild-type predictions, indicating that the L-protein is still likely to fold correctly. Therefore, these results suggest that the designed mutations may reduce the dependence of the L-protein on the DnaJ chaperone while maintaining a stable protein structure. This computational approach demonstrates how targeted mutagenesis combined with AlphaFold2-Multimer predictions can be used to design L-protein variants with potentially lower chaperone dependency.

in this homework, AI ChatGPT assisted me in organizing and clearly articulating my answers and descriptions, ensuring that the content is well-structured and easy to understand.

Week 06 HW: genetic circuits part-I

Answer 01: image_ref The Phusion High-Fidelity PCR Master Mix is a ready-to-use solution used for PCR amplification with high accuracy. It already contains the main components needed for DNA amplification. image_ref

• Phusion DNA Polymerase

Phusion DNA polymerase is the enzyme that copies the DNA during PCR. It synthesizes new DNA strands using the template DNA. This polymerase has proofreading activity, which helps detect and correct errors during DNA synthesis, making the amplification very accurate.

• Reaction Buffer

The reaction buffer provides the optimal chemical environment for the polymerase to function properly. It maintains the correct pH and salt conditions needed for efficient DNA amplification. Different buffers can be used depending on the DNA template, such as HF buffer for standard templates or GC buffer for GC-rich DNA.

• MgCl₂ (Magnesium Chloride)

Magnesium ions are essential cofactors for DNA polymerase activity. They help stabilize the interaction between the enzyme, the primers, and the DNA template during DNA synthesis.

• dNTPs (Deoxynucleotide Triphosphates)

dNTPs are the building blocks used to synthesize new DNA strands. They include dATP, dTTP, dCTP, and dGTP. During PCR, the polymerase adds these nucleotides to the growing DNA strand according to complementary base-pairing rules.

• Additional additives (e.g., DMSO)

Some reactions may include additives such as DMSO, which helps improve amplification of GC-rich DNA by reducing secondary structures and improving primer binding.

Answer 02:

The primer annealing temperature (Ta) is the temperature at which primers bind to the DNA template during PCR. It mainly depends on the primer melting temperature (Tm). In practice, the annealing temperature is usually set about 3–5 °C lower than the lowest primer Tm so that the primers can bind correctly to the DNA template.

Important factors (from the lab and lecture)

-Primer melting temperature (Tm): image_ref The melting temperature is the temperature at which 50 % of the primer–DNA duplex separates into single strands. It is the main factor used to determine the annealing temperature.

-Primer length (18–22 nucleotides): image_ref Primers are usually designed with a length of 18–22 bases. This length provides good specificity and stable binding to the template DNA.

-GC content (40–60 %)

Primers should contain about 40–60 % GC bases. GC pairs form stronger bonds than AT pairs, which increases the stability of primer binding.

-GC clamp (≤3 GC bases at the 3′ end) image_ref A small number of GC bases at the 3′ end of the primer (called a GC clamp) helps the primer bind more strongly to the template DNA and improves PCR efficiency.

-Primer secondary structures image_ref Primers should avoid forming hairpins, self-dimers, or cross-dimers. These structures prevent primers from binding properly to the template DNA. Recommended primer Tm range (52–58 °C).

Primers are usually designed to have a melting temperature between 52 °C and 58 °C, which allows efficient and specific amplification.

-GC sequence composition

Primers with higher GC content bind more strongly because GC base pairs form three hydrogen bonds, while AT pairs form only two.

Additional factors

-Ionic environment (Mg²⁺ and salt concentration)

Ions such as Mg²⁺ and other salts stabilize the DNA double strand and influence primer binding. Changes in these concentrations can affect the optimal annealing temperature.

-Primer concentration

Higher primer concentrations increase the probability that primers will bind to the DNA template, which can influence the optimal annealing temperature.

-Optimization using gradient PCR

In many experiments, scientists perform gradient PCR to test different annealing temperatures and find the best one for efficient and specific amplification.

Answer 03:

PCR (Polymerase Chain Reaction) and restriction enzyme digestion are two common molecular biology methods used to produce linear DNA fragments, but they work in very different ways. PCR works by amplifying a specific DNA sequence, while restriction enzymes cut existing DNA at specific recognition sites. Both techniques are widely used in cloning and DNA assembly experiments, including methods such as Gibson Assembly. image_ref PCR is a technique used to copy a specific region of DNA many times. It starts with a small amount of template DNA and uses specific primers that bind to the target sequence. A heat-stable DNA polymerase enzyme (such as Taq or Phusion polymerase) then synthesizes new DNA strands. The reaction takes place in a thermal cycler, which repeatedly changes the temperature through three main steps: denaturation, where the DNA strands separate; annealing, where primers bind to the template DNA; and extension, where the polymerase enzyme copies the DNA. After many cycles, PCR produces large amounts of a specific linear DNA fragment. One advantage of PCR is that researchers can design primers to add new sequences to the ends of the fragment, such as restriction sites or overlapping regions for Gibson Assembly. image_ref Restriction enzyme digestion works differently. Instead of amplifying DNA, it cuts existing DNA molecules at specific short sequences called recognition sites. Restriction enzymes recognize these sequences and cut the DNA at or near those locations. In a typical protocol, the DNA (for example a plasmid) is mixed with the restriction enzyme and a specific buffer, and the reaction is incubated at a constant temperature, usually around 37 °C, for about one hour. The enzyme then cuts the DNA to produce linear fragments. Depending on the enzyme, the cut DNA can produce sticky ends (overhangs) or blunt ends, which can be used for cloning.

These two methods are used in different situations depending on the goal of the experiment. PCR is preferred when the DNA is present in low concentration, because it can amplify a very small amount of template into large quantities. PCR is also useful when researchers want to introduce new sequences, mutations, or overlaps into the DNA fragment. For example, primers can be designed to add restriction sites, promoter sequences, or homologous overlaps needed for Gibson Assembly. PCR is also commonly used when scientists want to isolate a specific gene or region from genomic DNA.

Restriction enzyme digestion is more suitable when the DNA is already available in large quantities, such as a purified plasmid. It is commonly used when researchers want to cut DNA at precise and known locations to isolate fragments or prepare a plasmid for cloning. Restriction enzymes are also often used for diagnostic analysis, such as verifying plasmid identity or checking the size of DNA fragments through restriction mapping.

Answer 04:

To ensure that DNA fragments produced by PCR or restriction digestion are suitable for Gibson Assembly, several preparation and verification steps must be followed. Gibson Assembly joins DNA fragments that contain overlapping homologous sequences, so the fragments must be designed carefully and purified before the assembly reaction.

The first and most important step is primer design. Primers used in PCR should include overlapping sequences of about 20–40 base pairs that match the ends of the neighboring DNA fragment. These overlaps allow the fragments to align and assemble correctly during the Gibson reaction. The overlapping regions should have similar melting temperatures (Tm) to allow stable annealing during the isothermal reaction. It is also important to design overlaps with a balanced GC content and to avoid strong secondary structures such as hairpins, because these structures can reduce assembly efficiency.

Another important step is using a high-fidelity DNA polymerase, such as Phusion or Q5 polymerase, during PCR amplification. These enzymes have proofreading activity and reduce the number of mutations introduced during amplification. This is important because Gibson Assembly is often used to construct precise DNA sequences or multi-fragment plasmids.

After PCR amplification, the DNA fragments should be verified using agarose gel electrophoresis to confirm that the fragments have the expected size. The correct DNA bands are then purified from the gel to remove primers, nucleotides, enzymes, and non-specific products that might interfere with the assembly reaction.

To reduce background contamination from the original template plasmid, PCR products can be treated with the restriction enzyme DpnI, which digests methylated template DNA but does not affect the newly synthesized PCR fragments.

If a plasmid backbone is used, the vector must be completely linearized before Gibson Assembly. This can be done by restriction enzyme digestion or PCR. When restriction enzymes are used, it is important to ensure that the digestion is complete so that no circular plasmid remains, because this could produce unwanted background colonies during transformation.

Another important step is DNA quantification. The concentration of each DNA fragment should be carefully measured using methods such as fluorometric quantification (for example Qubit) or gel analysis. The correct molar ratio of vector to insert fragments, often about 1:2 or 1:3, helps improve assembly efficiency.

Finally, after Gibson Assembly and bacterial transformation, the resulting plasmid constructs are usually verified by DNA sequencing to confirm that the fragments assembled correctly and that no mutations were introduced during PCR.

Answer 05:

Plasmid DNA enters Escherichia coli cells during a process called bacterial transformation. In this process, the bacterial cells must first be made competent, meaning their membranes become temporarily able to allow DNA molecules to enter. image_ref In the most common method, called chemical transformation, the cells are treated with a solution containing calcium chloride (CaCl₂). The calcium ions (Ca²⁺) play an important role because they neutralize the negative charges on both the plasmid DNA and the phospholipids of the bacterial membrane. Normally, DNA and the membrane repel each other because they are both negatively charged. The calcium ions reduce this repulsion and allow the plasmid DNA to attach to the surface of the bacterial cell.

After mixing the plasmid DNA with the competent cells, the mixture is kept on ice (around 0 °C) for a short time. The cells are then exposed to a brief heat shock, usually at about 42 °C for 30–60 seconds. This sudden temperature change creates a strong thermal gradient between the cold cells and the warm environment. As a result, the bacterial membrane becomes temporarily destabilized and small pores form, allowing the plasmid DNA to pass into the cell.

Immediately after the heat shock, the cells are placed back on ice. This rapid cooling helps close the pores and stabilize the membrane again. The cells are then transferred into a nutrient recovery medium and incubated for a short period. During this recovery step, the cells repair their membranes and begin expressing the antibiotic resistance gene carried by the plasmid.

Finally, the bacteria are plated on agar plates containing the appropriate antibiotic. Only the cells that successfully received the plasmid DNA will survive and form colonies.

Another alternative method used to introduce plasmid DNA into E. coli is electroporation. In this method, competent bacterial cells are mixed with plasmid DNA and placed in a special electroporation cuvette. A short electrical pulse is then applied using an electroporator. The electrical pulse temporarily creates small pores in the bacterial cell membrane, allowing the plasmid DNA to pass directly into the cell.

After the pulse, the membrane quickly reseals and the cells recover in a nutrient medium. Electroporation is often more efficient than chemical transformation and is commonly used when transforming difficult DNA constructs or when very high transformation efficiency is required.

Answer 06:

Another DNA assembly method is Golden Gate Assembly, which allows several DNA fragments to be joined together in a single reaction. This technique uses special restriction enzymes called Type IIS restriction enzymes, such as BsaI or BsmBI, together with T4 DNA Ligase. Unlike traditional restriction enzymes, Type IIS enzymes cut outside their recognition sequence, which allows scientists to design custom 4-base pair overhangs at the ends of DNA fragments. These overhangs are designed so that fragments can only join with the correct neighboring fragment, ensuring the correct order and orientation of the assembled DNA. During the reaction, the restriction enzyme cuts the DNA fragments and creates the overhangs, and the DNA ligase joins the fragments together. The recognition sites of the restriction enzyme are removed during assembly, which means the final DNA construct cannot be cut again by the same enzyme. The digestion and ligation steps occur in the same tube using alternating temperatures, making Golden Gate Assembly a very efficient method for assembling multiple DNA fragments, especially in synthetic biology and modular cloning experiments. image_ref

This diagram is a clear example of Golden Gate Assembly, a cloning method that joins several DNA fragments in one reaction. In the example, three DNA parts — Promoter (Fragment A), ORF (Fragment B), and Terminator (Fragment C) — are assembled into a final plasmid called the destination vector. The process uses the Type IIS restriction enzyme BsaI together with T4 DNA Ligase.

In the first step, each fragment is present in an entry vector that contains the BsaI recognition site (GGTCTC). Unlike classical restriction enzymes, BsaI cuts outside of its recognition site, generating specific 4-base pair sticky ends (overhangs). Because the cut occurs outside the recognition sequence, the recognition site is removed during assembly and does not remain in the final DNA construct.

The fragments are designed with specific overhangs so they connect in the correct order. For example, Fragment A ends with the overhang CCAC, which matches the beginning of Fragment B. Fragment B ends with CGAT, which matches the start of Fragment C. These complementary overhangs act like puzzle pieces, ensuring that the fragments assemble correctly and in the proper orientation.

All fragments, the destination vector, BsaI, and T4 DNA Ligase are mixed in a single tube. During the reaction, BsaI cuts the DNA fragments to create sticky ends, and T4 DNA ligase joins fragments with matching overhangs. The reaction cycles between temperatures that allow DNA digestion and ligation, gradually assembling the correct construct.

Once fragments are ligated together, the BsaI recognition sites are no longer present, so the final product cannot be cut again by the enzyme. This makes the process efficient and irreversible, allowing the formation of a seamless DNA construct containing Fragment A + Fragment B + Fragment C in the destination plasmid. image_ref

• Modeling Golden Gate Assembly in Benchling

In this part, I modeled a Golden Gate Assembly to construct a genetic circuit for my second project, which is the engineering of an Escherichia coli reporter strain to monitor protein aging using a fluorescent timer protein.

First, I selected all the genetic elements needed for my construct. The backbone plasmid was obtained from Addgene, and it already contains a T7 promoter, a ribosome binding site (RBS), and a T7 terminator, which are very suitable for strong expression of the inserted gene. This vector also includes the GST (Glutathione S-Transferase from Schistosoma japonicum), which I used as the protein of interest because it has stable folding and is suitable for initial testing of my genetic system.

Then, I designed two additional fragments: a flexible linker (Gly₄Ser)₃ and a fluorescent timer (FT) protein (Medium FT). Their sequences were also obtained from Addgene. The linker allows proper folding between the GST protein and the fluorescent timer, while the FT protein provides a signal that changes over time, allowing estimation of protein age inside the cell.

• The full sequence of the pET28a-GST-P(11)4:

GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCTGTCTTCGGTATCGTCGTATCCCACTACCGAGATATCCGCACCAACGCGCAGCCCGGACTCGGTAATGGCGCGCATTGCGCCCAGCGCCATCTGATCGTTGGCAACCAGCATCGCAGTGGGAACGATGCCCTCATTCAGCATTTGCATGGTTTGTTGAAAACCGGACATGGCACTCCAGTCGCCTTCCCGTTCCGCTATCGGCTGAATTTGATTGCGAGTGAGATATTTATGCCAGCCAGCCAGACGCAGACGCGCCGAGACAGAACTTAATGGGCCCGCTAACAGCGCGATTTGCTGGTGACCCAATGCGACCAGATGCTCCACGCCCAGTCGCGTACCGTCTTCATGGGAGAAAATAATACTGTTGATGGGTGTCTGGTCAGAGACATCAAGAAATAACGCCGGAACATTAGTGCAGGCAGCTTCCACAGCAATGGCATCCTGGTCATCCAGCGGATAGTTAATGATCAGCCCACTGACGCGTTGCGCGAGAAGATTGTGCACCGCCGCTTTACAGGCTTCGACGCCGCTTCGTTCTACCATCGACACCACCACGCTGGCACCCAGTTGATCGGCGCGAGATTTAATCGCCGCGACAATTTGCGACGGCGCGTGCAGGGCCAGACTGGAGGTGGCAACGCCAATCAGCAACGACTGTTTGCCCGCCAGTTGTTGTGCCACGCGGTTGGGAATGTAATTCAGCTCCGCCATCGCCGCTTCCACTTTTTCCCGCGTTTTCGCAGAAACGTGGCTGGCCTGGTTCACCACGCGGGAAACGGTCTGATAAGAGACACCGGCATACTCTGCGACATCGTATAACGTTACTGGTTTCACATTCACCACCCTGAATTGACTCTCTTCCGGGCGCTATCATGCCATACCGCGAAAGGTTTTGCGCCATTCGATGGTGTCCGGGATCTCGACGCTCTCCCTTATGCGACTCCTGCATTAGGAAGCAGCCCAGTAGTAGGTTGAGGCCGTTGAGCACCGCCGCCGCAAGGAATGGTGCATGCAAGGAGATGGCGCCCAACAGTCCCCCGGCCACGGGGCCTGCCACCATACCCACGCCGAAACAAGCGCTCATGAGCCCGAAGTGGCGAGCCCGATCTTCCCCATCGGTGATGTCGGCGATATAGGCGCCAGCAACCGCACCTGTGGCGCCGGTGATGCCGGCCACGATGCGTCCGGCGTAGAGGATCGAGATCTCGATCCCGCGAAATTAATACGACTCACTATAGGGGAATTGTGAGCGGATAACAATTCCCCTCTAGAAATAATTTTGTTTAACTTTAAGAAGGAGATATACCATGGGCCATCATCATCATCATCATAGCCCGATCCTGGGTTACTGGAAAATCAAGGGCCTGGTGCAACCGACCCGCCTGCTGCTGGAATACCTGGAGGAAAAATACGAGGAACACCTGTATGAGCGTGACGAAGGCGATAAGTGGCGTAACAAGAAATTCGAGCTGGGTCTGGAATTTCCGAACCTGCCGTACTATATTGACGGCGATGTGAAACTGACCCAGAGCATGGCGATCATTCGTTACATCGCGGACAAACACAACATGCTGGGTGGCTGCCCGAAGGAGCGTGCGGAAATTAGCATGCTGGAGGGCGCGGTGCTGGATATTCGTTACGGTGTTAGCCGTATCGCGTATAGCAAAGACTTCGAAACCCTGAAGGTGGATTTTCTGAGCAAACTGCCGGAGATGCTGAAGATGTTCGAGGACCGTCTGTGCCACAAAACCTATCTGAACGGTGACCACGTTACCCACCCGGATTTTATGCTGTACGACGCGCTGGATGTGGTTCTGTATATGGACCCGATGTGCCTGGATGCGTTCCCGAAGCTGGTTTGCTTTAAGAAACGTATCGAGGCGATTCCGCAAATCGACAAGTACCTGAAAAGCAGCAAGTATATTGCGTGGCCGCTGCAAGGTTGGCAAGCGACCTTTGGTGGCGGTGATCACCCGCCGAAGGGTGGCGGTGGTAGCGGCGGCGGCGGCAGCCAACAGCGTTTTGAATGGGAATTTGAACAGCAGTAATAACTCGAGCACCACCACCACCACCACTGAGATCCGGCTGCTAACAAAGCCCGAAAGGAAGCTGAGTTGGCTGCTGCCACCGCTGAGCAATAACTAGCATAACCCCTTGGGGCCTCTAAACGGGTCTTGAGGGGTTTTTTGCTGAAAGGAGGAACTATATCCGGATTGGCGAATGGGACGCGCCCTGTAGCGGCGCATTAAGCGCGGCGGGTGTGGTGGTTACGCGCAGCGTGACCGCTACACTTGCCAGCGCCCTAGCGCCCGCTCCTTTCGCTTTCTTCCCTTCCTTTCTCGCCACGTTCGCCGGCTTTCCCCGTCAAGCTCTAAATCGGGGGCTCCCTTTAGGGTTCCGATTTAGTGCTTTACGGCACCTCGACCCCAAAAAACTTGATTAGGGTGATGGTTCACGTAGTGGGCCATCGCCCTGATAGACGGTTTTTCGCCCTTTGACGTTGGAGTCCACGTTCTTTAATAGTGGACTCTTGTTCCAAACTGGAACAACACTCAACCCTATCTCGGTCTATTCTTTTGATTTATAAGGGATTTTGCCGATTTCGGCCTATTGGTTAAAAAATGAGCTGATTTAACAAAAATTTAACGCGAATTTTAACAAAATATTAACGCTTACAATTTAGGTGGCACTTTTCGGGGAAATGTGCGCGGAACCCCTATTTGTTTATTTTTCTAAATACATTCAAATATGTATCCGCTCATGAATTAATTCTTAGAAAAACTCATCGAGCATCAAATGAAACTGCAATTTATTCATATCAGGATTATCAATACCATATTTTTGAAAAAGCCGTTTCTGTAATGAAGGAGAAAACTCACCGAGGCAGTTCCATAGGATGGCAAGATCCTGGTATCGGTCTGCGATTCCGACTCGTCCAACATCAATACAACCTATTAATTTCCCCTCGTCAAAAATAAGGTTATCAAGTGAGAAATCACCATGAGTGACGACTGAATCCGGTGAGAATGGCAAAAGTTTATGCATTTCTTTCCAGACTTGTTCAACAGGCCAGCCATTACGCTCGTCATCAAAATCACTCGCATCAACCAAACCGTTATTCATTCGTGATTGCGCCTGAGCGAGACGAAATACGCGATCGCTGTTAAAAGGACAATTACAAACAGGAATCGAATGCAACCGGCGCAGGAACACTGCCAGCGCATCAACAATATTTTCACCTGAATCAGGATATTCTTCTAATACCTGGAATGCTGTTTTCCCGGGGATCGCAGTGGTGAGTAACCATGCATCATCAGGAGTACGGATAAAATGCTTGATGGTCGGAAGAGGCATAAATTCCGTCAGCCAGTTTAGTCTGACCATCTCATCTGTAACATCATTGGCAACGCTACCTTTGCCATGTTTCAGAAACAACTCTGGCGCATCGGGCTTCCCATACAATCGATAGATTGTCGCACCTGATTGCCCGACATTATCGCGAGCCCATTTATACCCATATAAATCAGCATCCATGTTGGAATTTAATCGCGGCCTAGAGCAAGACGTTTCCCGTTGAATATGGCTCATAACACCCCTTGTATTACTGTTTATGTAAGCAGACAGTTTTATTGTTCATGACCAAAATCCCTTAACGTGAGTTTTCGTTCCACTGAGCGTCAGACCCCGTAGAAAAGATCAAAGGATCTTCTTGAGATCCTTTTTTTCTGCGCGTAATCTGCTGCTTGCAAACAAAAAAACCACCGCTACCAGCGGTGGTTTGTTTGCCGGATCAAGAGCTACCAACTCTTTTTCCGAAGGTAACTGGCTTCAGCAGAGCGCAGATACCAAATACTGTCCTTCTAGTGTAGCCGTAGTTAGGCCACCACTTCAAGAACTCTGTAGCACCGCCTACATACCTCGCTCTGCTAATCCTGTTACCAGTGGCTGCTGCCAGTGGCGATAAGTCGTGTCTTACCGGGTTGGACTCAAGACGATAGTTACCGGATAAGGCGCAGCGGTCGGGCTGAACGGGGGGTTCGTGCACACAGCCCAGCTTGGAGCGAACGACCTACACCGAACTGAGATACCTACAGCGTGAGCTATGAGAAAGCGCCACGCTTCCCGAAGGGAGAAAGGCGGACAGGTATCCGGTAAGCGGCAGGGTCGGAACAGGAGAGCGCACGAGGGAGCTTCCAGGGGGAAACGCCTGGTATCTTTATAGTCCTGTCGGGTTTCGCCACCTCTGACTTGAGCGTCGATTTTTGTGATGCTCGTCAGGGGGGCGGAGCCTATGGAAAAACGCCAGCAACGCGGCCTTTTTACGGTTCCTGGCCTTTTGCTGGCCTTTTGCTCACATGTTCTTTCCTGCGTTATCCCCTGATTCTGTGGATAACCGTATTACCGCCTTTGAGTGAGCTGATACCGCTCGCCGCAGCCGAACGACCGAGCGCAGCGAGTCAGTGAGCGAGGAAGCGGAAGAGCGCCTGATGCGGTATTTTCTCCTTACGCATCTGTGCGGTATTTCACACCGCAATGGTGCACTCTCAGTACAATCTGCTCTGATGCCGCATAGTTAAGCCAGTATACACTCCGCTATCGCTACGTGACTGGGTCATGGCTGCGCCCCGACACCCGCCAACACCCGCTGACGCGCCCTGACGGGCTTGTCTGCTCCCGGCATCCGCTTACAGACAAGCTGTGACCGTCTCCGGGAGCTGCATGTGTCAGAGGTTTTCACCGTCATCACCGAAACGCGCGAGGCAGCTGCGGTAAAGCTCATCAGCGTGGTCGTGAAGCGATTCACAGATGTCTGCCTGTTCATCCGCGTCCAGCTCGTTGAGTTTCTCCAGAAGCGTTAATGTCTGGCTTCTGATAAAGCGGGCCATGTTAAGGGCGGTTTTTTCCTGTTTGGTCACTGATGCCTCCGTGTAAGGGGGATTTCTGTTCATGGGGGTAATGATACCGATGAAACGAGAGAGGATGCTCACGATACGGGTTACTGATGATGAACATGCCCGGTTACTGGAACGTTGTGAGGGTAAACAACTGGCGGTATGGATGCGGCGGGACCAGAGAAAAATCACTCAGGGTCAATGCCAGCGCTTCGTTAATACAGATGTAGGTGTTCCACAGGGTAGCCAGCAGCATCCTGCGATGCAGATCCGGAACATAATGGTGCAGGGCGCTGACTTCCGCGTTTCCAGACTTTACGAAACACGGAAACCGAAGACCATTCATGTTGTTGCTCAGGTCGCAGACGTTTTGCAGCAGCAGTCGCTTCACGTTCGCTCGCGTATCGGTGATTCATTCTGCTAACCAGTAAGGCAACCCCGCCAGCCTAGCCGGGTCCTCAACGACAGGAGCACGATCATGCGCACCCGTGGGGCCGCCATGCCGGCGATAATGGCCTGCTTCTCGCCGAAACGTTTGGTGGCGGGACCAGTGACGAAGGCTTGAGCGAGGGCGTGCAAGATTCCGAATACCGCAAGCGACAGGCCGATCATCGTCGCGCTCCAGCGAAAGCGGTCCTCGCCGAAAATGACCCAGAGCGCTGCCGGCACCTGTCCTACGAGTTGCATGATAAAGAAGACAGTCATAAGTGCGGCGACGATAGTCATGCCCCGCGCCCACCGGAAGGAGCTGACTGGGTTGAAGGCTCTCAAGGGCATCGGTCGAGATCCCGGTGCCTAATGAGTGAGCTAACTTACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC

• Linker sequence (Gly₄Ser)₃ :

GGTGGCGGTGGCTCGGGCGGTGGTGGGTCGGGTGGCGGCGGATCT

• Medium-FT sequence:

atggtgagcaagggcgaggaggataacatggccatcatcaaggaattcatgcgtttcaaggtgcacctggagggctccgtggacggccacgagttcgagatcgagggcgagggcgagggccgcccctacgagggcacccagagcgccaagctgaaggtgaccaagggtggccccctgcccttcgcctgggacatcctgtcccctcagttcatgtacggctccagggcctacgtgaagcaccccgccgacatccccgactactggaagctgtccttccccgagggcttcaagtgggagcgcgtgatgaacttcgaggatggcggcgtggtgaccgtgacccaggactcctccctgcaggacggcgagttcatctacaaggtgaagctgcgcggcaccaacttcccttccgacggccccgtaatgcagaagaagaccatgggctgggaggcctcctccgagcggatataccccgaggacggcgccctgaagggcgagatcaagcagaggctgaagctgaaggacggcggccactacgacgctgaggtcaagaccacctacaaggccaagaagcccgtgctgctgcccggcgcctacaacgtcaacatcaagatggacatcacctcccacaacgaggactacaccatcgttgaacagtgcgaacgcgccgagggccaccattccaccggcggcatggacgagctgtacaagtaa

At the beginning, I manually designed the overhangs based on the coding sequence. I assumed that the last four nucleotides of the GST sequence (GAAG) would serve as the correct overhang to connect with the next fragment. Based on this assumption, I designed the linker fragment to have a compatible overhang (GAAG, GGTA). Similarly, I defined the overhangs between the linker and the fluorescent timer protein (GGTA) in order to maintain a continuous reading frame. During this step, I also verified that no frameshift was introduced at the junctions and that the coding sequence remained in frame across all fragments as indicated in the following table:

In addition, I checked that all BsaI restriction enzyme recognition sites were positioned outside of the fragments that would be recovered after digestion, ensuring that the internal sequences of the inserts would not be disrupted during the assembly process. the designed overhangs are as the following:

The designed overhangs are supposed to orient the assembly in the following order: the linker is placed immediately after the GST sequence, and the Medium FT is positioned just before the C-terminal His tag, as indicated in the following diagram:

After preparing the vector and all fragments, the designed vector digestion cuts were defined as follows: The designed linker fragment sticky ends were defined as follows: The designed Medium FT sticky ends were defined as follows:

Be careful !! A critical point to consider during the design is the correct placement of BsaI restriction enzyme recognition sites. For the inserted fragments, the BsaI sites must be located outside of the sequences of interest so that they are removed during digestion and do not remain in the final construct. In contrast, for the backbone vector, the BsaI sites must be positioned within the region to be replaced, so that digestion removes this segment and allows the insertion of the designed fragments.

It is also essential to ensure that the BsaI recognition sites are oriented correctly (inverted orientation) to generate the desired overhangs and to cut the backbone precisely at the intended insertion site. Any incorrect placement or orientation of these sites can lead to incompatible sticky ends and result in assembly failure.

I imported all sequences into Benchling and created a new assembly using the Golden Gate cloning option. I selected the pET-28 plasmid as the backbone and added the designed fragments, including the linker and the fluorescent timer protein, as inserts. I specified the use of the BsaI restriction enzyme and defined the final construct as circular. Since all sequences were already designed with appropriate BsaI recognition sites, I selected the option to use existing restriction sites for fragment generation. I then attempted to run the assembly.

However, the assembly failed, and Benchling returned an error indicating that the sticky ends were incompatible. Specifically, the system showed a mismatch between the overhangs “AAGC” (from the vector) and “GAAG” (from the insert). This result indicated that the fragments could not ligate properly.

After analyzing this issue, I realized that the mistake came from misunderstanding how Golden Gate Assembly works. I initially assumed that the overhang corresponds directly to visible nucleotides in the sequence. In reality, the overhang is determined by the position of the BsaI cutting site, not simply by the sequence at the end of the gene. Since BsaI cuts outside of its recognition site, the actual generated overhang in the vector was “AAGC” and not “GAAG” as I had expected.

This mismatch between expected and real overhangs caused the failure of the assembly. Additionally, the cloning workflow in Benchling does not automatically correct or reinterpret overhangs; it strictly checks for compatibility. Therefore, any small design error leads to a complete assembly failure.

In order to overcome the limitations encountered in the first approach, I tried another method available in Benchling by using the Assembly tool dedicated to multi-fragment cloning. This method is specifically designed to simulate Golden Gate Assembly in a more automated and flexible way, allowing better handling of fragment compatibility and overhang generation.

First, I opened the Assembly tool from the bottom toolbar and created a new assembly. I then added all the required DNA sequences, including the pET-28 plasmid as the backbone and the designed fragments (linker and Medium FT) as inserts. After that, I selected the BsaI restriction enzyme as the Type IIS restriction enzyme used for the assembly.

Unlike the previous method, this approach automatically analyzed the positions of the BsaI recognition sites and simulated the digestion process. It generated the correct sticky ends based on the actual cutting positions of the enzyme and evaluated the compatibility between fragments. This allowed the system to correctly align and assemble the different parts according to their matching overhangs.

After running the assembly, the construct was successfully generated as a circular plasmid. I carefully verified that all fragments were assembled in the correct order and orientation. I also confirmed that no frameshift was introduced across the junctions and that the reading frame was maintained from the GST sequence through the linker and into the fluorescent timer protein. In addition, I checked that no unwanted BsaI sites remained inside the final construct and that all restriction sites had been properly removed during the assembly process.

this is the direct link to benchling for this assembly: using the assembly tool

in this homework, AI ChatGPT assisted me in organizing and clearly articulating my answers and descriptions, ensuring that the content is well-structured and easy to understand.

Sources:

• Activity 3-2—Primer Design and Barcoding via Bioinformatics. (2025, March 8). Biology LibreTexts. https://bio.libretexts.org/Courses/Irvine_Valley_College/Lab_manual%3A_Molecular_biology_and_Genetic_Engineering_(Biot_275)/11%3A_Activity_3-2_-_Primer_Design_and_Sequence_Verification_Using_Bioinformatics
• Alberts, B., Johnson, A., Lewis, J., Raff, M., Roberts, K., & Walter, P. (2002). Isolating, Cloning, and Sequencing DNA. In Molecular Biology of the Cell. 4th edition. Garland Science. https://www.ncbi.nlm.nih.gov/books/NBK26837/
• Biocat.com/bc/files/Gibson_Guide_V2_101417_web_version_8.5_x_11_FINAL.pdf. (n.d.). Retrieved March 30, 2026, from https://www.biocat.com/bc/files/Gibson_Guide_V2_101417_web_version_8.5_x_11_FINAL.pdf
• Bird, J. E., Marles-Wright, J., & Giachino, A. (2022). A User’s Guide to Golden Gate Cloning Methods and Standards. ACS Synthetic Biology, 11(11), 3551–3563. https://doi.org/10.1021/acssynbio.2c00355
• Competent Cells: Principle, Methods & Functions Explained. (n.d.). VEDANTU. Retrieved March 30, 2026, from https://www.vedantu.com/biology/competent-cells
• DNA Quantification with a Plate Reader | BMG LABTECH. (n.d.). Retrieved March 30, 2026, from https://www.bmglabtech.com/en/blog/dna-quantification/
• Erjavec, M. S. (2019). Annealing Temperature of 55°C and Specificity of Primer Binding in PCR Reactions. In Synthetic Biology—New Interdisciplinary Science. IntechOpen. https://doi.org/10.5772/intechopen.85164
• Froger, A., & Hall, J. E. (2007). Transformation of Plasmid DNA into E. coli Using the Heat Shock Method. Journal of Visualized Experiments : JoVE, (6), 253. https://doi.org/10.3791/253
• Golden Gate Assembly—Snapgene. (n.d.). Retrieved March 30, 2026, from https://www.snapgene.com/guides/golden-gate-assembly
• Hoseini, S. S., & Sauer, M. G. (2015). Molecular cloning using polymerase chain reaction, an educational guide for cellular engineering. Journal of Biological Engineering, 9, 2. https://doi.org/10.1186/1754-1611-9-2
• How are competent bacterial cells transformed with a plasmid? (n.d.). Retrieved March 30, 2026, from https://worldwide.promega.com/resources/pubhub/enotes/how-are-competent-bacterial-cells-transformed-with-a-plasmid/
• Jacob Elmer. (2019, August 29). Chemical Transformation of E. coli [Video recording]. https://www.youtube.com/watch?v=xcc2ywDASag Med.unc.edu/pharm/sondeklab/wp-content/uploads/sites/868/2019/10/gibson-cloning.pdf. (n.d.). Retrieved March 30, 2026, from https://www.med.unc.edu/pharm/sondeklab/wp-content/uploads/sites/868/2019/10/gibson-cloning.pdf
• Molecular Cloning Methods – Benchling. (n.d.). Retrieved March 30, 2026, from https://help.benchling.com/hc/en-us/articles/9684255457805-Molecular-Cloning-Methods
• Molecular cloning using Gibson assembly. (n.d.). EPFL. Retrieved March 30, 2026, from https://www.epfl.ch/labs/lpbs/internal/general/molecular-cloning-using-gibson-assembly/
• PCR conditions | Primer annealing specificity | PCR buffers. (n.d.). Retrieved March 30, 2026, from https://www.qiagen.com/us/knowledge-and-support/knowledge-hub/bench-guide/pcr/introduction/pcr-conditions
• PCR Using Phusion® High-Fidelity PCR Master Mix with HF Buffer (NEB #M0531) | NEB. (n.d.). Retrieved March 30, 2026, from https://www.neb.com/en/protocols/protocol-phusion-high-fidelity-pcr-master-mix-with-hf-buffer-m0531
• (PDF) Study of Transformation in Escherichia coli. (n.d.). ResearchGate. https://doi.org/10.13140/RG.2.2.22413.45285 Polymerase Chain Reaction—An overview | ScienceDirect Topics. (n.d.). Retrieved March 30, 2026, from https://www.sciencedirect.com/topics/food-science/polymerase-chain-reaction
• Protocol for PhusionTM High-Fidelity PCR Master Mix with GC Buffer | NEB. (n.d.). Retrieved March 30, 2026, from https://www.neb.com/en/protocols/protocol-for-phusion-high-fidelity-pcr-master-mix-with-gc-buffer-m0532
• Restriction Enzyme Digestion | NEB. (n.d.). Retrieved March 30, 2026, from https://www.neb.com/en/applications/cloning-and-synthetic-biology/dna-preparation/restriction-enzyme-digestion
• SnapGene. (2022, June 30). Introduction à l’assemblée de Golden Gate [Video recording]. https://www.youtube.com/watch?v=aBcqev1NMMo Sorida, M., & Bonasio, R. (2023). An efficient cloning method to expand vector and restriction site compatibility of Golden Gate Assembly. Cell Reports Methods, 3(8), 100564. https://doi.org/10.1016/j.crmeth.2023.100564
• Staff, L. in the L. (2025, May 29). Gibson Assembly 101: Expert Cloning Tips You Need to Know. Life in the Lab. https://www.thermofisher.com/blog/life-in-the-lab/gibson-assembly-101-expert-cloning-tips-you-need-to-know/
• The Different Types of PCR Methods | Pipette.com. (n.d.). Retrieved March 30, 2026, from https://pipette.com/blog/types-of-pcr Universal Annealing Temperature in PCR and its Impact on Amplification Results. (n.d.).

Week 07 HW: genetic circuits part-II

Traditional genetic circuits work using Boolean logic, where the output is binary (either ON or OFF). The output depends on whether the input signals pass a fixed threshold. For example, in a genetic AND gate, a protein is only produced when both transcription factors A and B are present above a certain level. If one of them is missing or below the threshold, the output is zero. This type of system is useful for simple decisions, but it has important limitations because real biological signals are usually continuous, variable, and noisy, not strictly ON or OFF.

Intracellular Artificial Neural Networks (IANNs) solve these limitations by mimicking artificial neural networks inside the cell. Instead of treating inputs as binary, IANNs assign a continuous weight to each input. These weighted inputs are then summed, and the result is passed through a biological activation mechanism (such as a riboswitch or a protease-regulated system) to generate a graded output.

IANN approach provides several important advantages:

• Continuous output resolution: Unlike Boolean circuits that only produce ON or OFF outputs, IANNs generate different levels of expression depending on the strength of the inputs. This allows cells to respond in a more precise and dose-dependent way, which is important for applications like metabolic regulation or controlled therapeutic delivery.

• Weighted signal integration: Each input does not contribute equally. Instead, every signal has a specific weight that determines how much it influences the final output. This allows the system to prioritize certain signals over others, which is not possible in traditional AND/OR gates where all inputs are treated equally.

• Robustness to biological noise: Cellular environments are naturally noisy, and signals can vary between cells. Because IANNs work with continuous values rather than strict thresholds, they are more tolerant to noise and variability, making them more reliable in real biological conditions.

• Greater computational power: A multilayer IANN can act as a universal function approximator, meaning it can represent very complex relationships between inputs and outputs. In contrast, Boolean circuits are limited to simple logical combinations, which restricts the complexity of decisions they can perform.

• Rational tunability: The weights and biases in an IANN can be adjusted through DNA design (for example, by modifying promoters or regulatory elements) or improved through directed evolution. This makes it possible to “train” the system to recognize complex patterns, such as a specific combination of biomarkers, with much higher precision than traditional Boolean circuits.

• Application of an IANN: Smart Lactase-Producing Probiotic System

As a highly relevant and practical application of an Intracellular Artificial Neural Network (IANN) I chose to apply it in the engineering of a probiotic bacterium capable of context-aware lactase production for the management of lactose intolerance. Unlike conventional synthetic circuits that respond to a single input in a binary manner, this system integrates multiple physiological signals from the gastrointestinal environment to produce a graded and condition-dependent enzymatic response.

-Input Layer: Multidimensional Environmental Sensing

The system incorporates multiple biologically relevant input signals, each representing a distinct physiological parameter of the gastrointestinal environment:

X1: Lactose concentration This serves as the primary input signal, directly reflecting the presence and abundance of the substrate requiring enzymatic degradation.

X2: pH level This input provides spatial context by distinguishing between different regions of the gastrointestinal tract. The acidic pH of the stomach versus the near-neutral pH of the intestine allows the system to restrict activation to physiologically appropriate locations, thereby preventing premature or energetically wasteful enzyme production.

X3: Inflammatory biomarkers Molecules such as nitric oxide, reactive oxygen species or cytokine-associated metabolites act as indicators of intestinal stress or dysbiosis. This input enables modulation of the system’s response based on host physiological state, allowing adaptive tuning of output under pathological conditions.

Lactose sensitivity can be increased using a strong promoter or high-affinity regulator, corresponding to a positive weight. pH sensitivity may be implemented through a regulatory element that suppresses output under acidic conditions, corresponding to a negative or inhibitory weight. Inflammatory signals could be integrated via modulatory promoters or regulatory RNAs that amplify output under stress conditions, acting as an adjustable positive or negative weight depending on the desired response.

At the molecular level, each input is transduced into regulatory signals (e.g., transcription factors, small RNAs, or protease-mediated regulators). These signals are then integrated through combinatorial gene regulation, where promoter strengths, ribosome binding site efficiencies, and degradation dynamics collectively encode the effective weights.

The aggregated signal undergoes a transformation through a biological activation function, which may be implemented via nonlinear regulatory elements such as riboswitches, cooperative transcriptional regulators, or proteolytic cascades. This step introduces thresholding and saturation effects analogous to activation functions in artificial neural networks, thereby enabling continuous and nonlinear input–output relationships.

Output Layer: Graded Lactase Expression

The final output of the system is the expression of the lactase enzyme, with expression levels determined by the integrated and nonlinearly transformed input signal

This enables a spectrum of responses:

• Sub-threshold activation: (e.g., low lactose concentration or inhibitory pH conditions) result in negligible or no enzyme production.
• Intermediate activation: moderate enzyme expression
• High activation: (e.g., high lactose concentration under optimal pH conditions, potentially combined with inflammatory signals) drive maximal enzyme production.

Functional Behavior and Decision-Making Capability

The system effectively implements a context-dependent decision-making process, wherein output is not determined by a single condition but by the weighted combination of multiple environmental cues. For example:

• The presence of lactose alone is insufficient to trigger activation under acidic conditions, thereby preventing inappropriate expression in the stomach.
• Under intestinal pH, lactose induces activation in a concentration-dependent manner.
• In the presence of both high lactose and inflammatory signals, the system can upregulate lactase production, potentially enhancing digestive efficiency under stress conditions.

Limitations and Practical Constraints

Despite its conceptual advantages, the implementation of such an IANN-based system faces several challenges:

• Stochastic gene expression: Intrinsic and extrinsic noise can introduce variability in circuit performance across individual cells. Parameter tuning complexity: Precise calibration of weights and activation thresholds through genetic elements (e.g., promoters, RBSs) remains experimentally demanding.
• Kinetic limitations: Transcriptional and translational processes impose temporal delays, limiting the speed of system response.
• Regulatory crosstalk: Interactions between synthetic and endogenous pathways may lead to unintended behaviors.
• Metabolic burden: The expression of complex regulatory networks can reduce host fitness and stability.
• Environmental variability: Dynamic and heterogeneous gut conditions may challenge the robustness and predictability of the system.

Implementation of a Multilayer Perceptron Using Endoribonucleases

To implement a multilayer perceptron in a biological system, the output of one computational layer must regulate the activity of the next. This can be achieved using a cascade of endoribonucleases, where each layer processes inputs and produces a regulatory molecule that serves as the input for the subsequent layer.

• Input Representation

The system integrates multiple biological inputs represented as molecular signals:

X1: Csy4 endoribonuclease (constitutively or inducibly expressed)

X2: an additional regulatory signal (e.g., inducible promoter or transcriptional activator)

X3: environmental or metabolic signal (e.g., pH, or inflammatory markers such as nitric oxide)

These inputs are converted into regulatory effects at the gene expression level, analogous to numerical inputs in an artificial neural network.

• Layer 1: Intermediate Processing In the first layer, the inputs jointly regulate the expression of an intermediate endoribonuclease (e.g., Cas6a). The mRNA encoding this enzyme is engineered to contain specific recognition sites for Csy4. As a result:

–> The presence of Csy4 (X1) induces cleavage of the mRNA, leading to repression of Cas6a expression

–> The second input (X2) can act as an activator, promoting transcription of the Cas6a gene

Thus, Layer 1 integrates activating and inhibitory signals. The resulting expression level of Cas6a reflects a balance between these opposing regulatory effects, analogous to a weighted sum followed by a nonlinear activation function in a perceptron.

• Layer 2: Output Generation The output of Layer 1 (Cas6a protein) serves as the regulatory input for Layer 2.

The mRNA encoding a reporter protein (e.g., GFP) is engineered to contain Cas6a recognition sites. Consequently:

–> High levels of Cas6a lead to cleavage of GFP mRNA and repression of fluorescence

–> Low levels of Cas6a allow GFP expression

This establishes a second computational layer in which the input is not external, but derived from the processed output of the first layer.

• System-Level Behavior

This cascading architecture enables hierarchical signal processing within the cell. –> When Csy4 levels are high, Cas6a production is suppressed, allowing GFP expression

–> When Csy4 levels are low and activation dominates, Cas6a is produced and represses GFP

Therefore, the final output depends on both the original inputs and the intermediate computation performed in Layer 1.

Fungal Materials and Their Uses

Fungi, especially their root-like networks called mycelium, are becoming a surprisingly powerful source of sustainable materials. Unlike plastics or leather, which require heavy manufacturing and chemicals, mycelium-based materials are grown. The fungi take agricultural waste—things like sawdust, rice husks, or hemp—and weave it into solid, structured materials. It’s almost like nature is doing the 3D printing for us.

One of the most familiar uses today is in packaging. Mycelium can form lightweight, shock-absorbing foams that replace Styrofoam or plastic inserts. Fragile items like electronics, furniture, or delicate goods can be safely packed in these eco-friendly alternatives. Companies like Ecovative and IKEA have already begun experimenting with this approach, showing that sustainable materials don’t have to compromise practicality. image_ref

Fungi are also stepping into the fashion world through myco-leather. By processing mats of mycelium into flexible sheets, it’s possible to make shoes, bags, and even clothing without harming animals. Myco-leather is fully biodegradable, reduces chemical waste, and offers a much lower environmental footprint than traditional leather production. It’s a great example of how biology can meet design. image_ref

In construction, fungal materials are finding their place as well. Mycelium boards can provide thermal and acoustic insulation, or serve as lightweight panels for ceilings and walls. They are naturally fire-resistant, resistant to pests like termites, and completely biodegradable. This means that even in building applications, mycelium offers both functionality and sustainability.

The versatility doesn’t stop there. Designers and researchers are exploring fungal foams, textiles, and even furniture, taking advantage of mycelium’s ability to grow into complex shapes. People are also experimenting with specialty applications, like wearable electronics, wound dressings, filters, and acoustic panels. Fungi aren’t just materials—they’re living factories that can be shaped, molded, and sometimes even programmed to do more. image_ref

Advantages and Disadvantages of Fungal Materials

Genetic Engineering and Synthetic Biology in Fungi

Fungi are not just fascinating organisms—they are also incredibly versatile tools for engineering. If I were to genetically engineer fungi, I would aim to enhance the properties that currently limit their use while maximizing their natural strengths. For instance, one limitation of mycelium-based materials is their mechanical strength, which can make them less competitive compared to plastics or synthetic foams. I would focus on modifying the cell wall composition or growth patterns to produce stronger, more durable materials, making fungi a realistic alternative for packaging, textiles, and construction.

Another area I would target is environmental resilience. Fungal materials are naturally biodegradable, which is a huge advantage, but they can degrade too quickly in humid or wet environments. By engineering fungi to better tolerate moisture or extreme temperatures, it would be possible to create materials that maintain their structure and functionality in a wider range of conditions, expanding their practical applications.

Beyond materials, fungi can also be engineered for functional enhancements. I would consider adding traits like pigmentation for natural coloring, antimicrobial properties to extend shelf life, or even self-healing abilities so that minor damage doesn’t ruin the product. These modifications could transform mycelium into “smart materials” that are not only sustainable but also highly functional.

Why Use Fungi Instead of Bacteria?

Fungi offer several important advantages over bacteria when it comes to synthetic biology. First, as eukaryotic organisms, they have more advanced cellular machinery. This allows them to properly fold and modify complex proteins through processes like glycosylation, which is essential for many pharmaceuticals and functional biomolecules.

Another major advantage is their filamentous growth. Many fungi grow as long branching structures (hyphae), which makes them very efficient at secreting enzymes and other products into their environment. This simplifies downstream processing because the desired product is often already outside the cell.

Fungi also have a much richer and more diverse metabolism compared to most bacteria. They naturally produce a wide range of secondary metabolites, which means they can be engineered to generate a broader variety of useful compounds, from drugs to pigments to bioactive molecules. In addition, fungi are generally more robust in industrial settings. They can grow on cheap, low-quality substrates like agricultural waste and tolerate harsher conditions than many bacteria, making them more practical for large-scale, sustainable production.

That said, working with fungi can be more complex. They tend to grow more slowly than bacteria, and genetic engineering tools are less standardized. However, despite these challenges, their unique capabilities make them extremely valuable for applications where bacteria fall short.

for this part, I have developed three potential ideas for my final project and would greatly appreciate your feedback to help refine my direction. While I am still open to suggestions, I currently find myself most aligned with my second idea, as it feels both biologically intuitive and well-matched to the techniques we have learned throughout the course.

The idea I am leaning toward is focused on engineering an E. coli reporter system to monitor protein aging during heterologous expression using a fluorescent timer protein. I am particularly drawn to this concept because it allows me to integrate multiple core synthetic biology tools, including DNA construct design, protein engineering, and computational structure prediction, while also remaining experimentally feasible within the scope of the course. In addition, the system is mechanistically clear, which makes it easier to design, test, and interpret.

I have further refined this idea into a more specific and functional design: a time-dependent protein quality control system in which a fluorescent timer regulates the exposure of a degron, leading to the selective degradation of aged proteins. In this system, a protein of interest is fused to a fluorescent timer and a C-terminal degron. As the protein matures and the timer shifts from its “young” to “old” fluorescence state, conformational or structural changes are expected to increase the accessibility of the degron. This, in turn, allows recognition by the host proteolytic machinery, enabling targeted degradation of older protein populations. The key modification from the original idea is the addition of a functional outcome—degradation—rather than only monitoring protein age. This transforms the system from a passive reporter into an active quality control mechanism. The purpose of this change is to address a limitation in current heterologous expression systems, where proteins can accumulate in misfolded or non-functional states over time. By selectively degrading older or potentially damaged proteins, this system could improve overall protein quality and stability.

The broader gap I am attempting to address is the lack of dynamic, time-resolved control over protein lifespan in bacterial systems. Most current approaches either measure protein expression statically or rely on constitutive degradation signals that do not account for protein age. This project introduces a strategy to link protein function, age, and degradation in a single genetically encoded system.

At this stage, I would greatly value any feedback on the conceptual design, feasibility, or potential improvements. In particular, I would appreciate input on whether the proposed mechanism for degron exposure is realistic, and whether there are alternative design strategies that could strengthen the system. Any suggestions on experimental design, protein choice, or construct optimization would also be extremely helpful.

Please feel free to share feedback through any preferred channel, including email or whatsApp. Thank you for your time and guidance.

Designing the isert sequence in Benchling:

for this idea i designed the genetic construct in Benchling that encodes a fusion protein consisting of GST as the protein of interest, followed by a flexible linker, a fluorescent timer protein, a second short linker, and a C-terminal ssrA degron whose sequences are represented in the following table:

ATG

AGCCCGATCCTGGGTTACTGGAAAATCAAGGGCCTGGTGCAACCGACCCGCCTGCTGCTGGAATACCTGGAGGAAAAATACGAGGAACACCTGTATGAGCGTGACGAAGGCGATAAGTGGCGTAACAAGAAATTCGAGCTGGGTCTGGAATTTCCGAACCTGCCGTACTATATTGACGGCGATGTGAAACTGACCCAGAGCATGGCGATCATTCGTTACATCGCGGACAAACACAACATGCTGGGTGGCTGCCCGAAGGAGCGTGCGGAAATTAGCATGCTGGAGGGCGCGGTGCTGGATATTCGTTACGGTGTTAGCCGTATCGCGTATAGCAAAGACTTCGAAACCCTGAAGGTGGATTTTCTGAGCAAACTGCCGGAGATGCTGAAGATGTTCGAGGACCGTCTGTGCCACAAAACCTATCTGAACGGTGACCACGTTACCCACCCGGATTTTATGCTGTACGACGCGCTGGATGTGGTTCTGTATATGGACCCGATGTGCCTGGATGCGTTCCCGAAGCTGGTTTGCTTTAAGAAACGTATCGAGGCGATTCCGCAAATCGACAAGTACCTGAAAAGCAGCAAGTATATTGCGTGGCCGCTGCAAGGTTGGCAAGCGACCTTTGGTGGCGGTGATCACCCGCCGAAG

GGTGGCGGTGGCTCGGGCGGTGGTGGGTCGGGTGGCGGCGGATCT

ATGGTGAGCAAGGGCGAGGAGGATAACATGGCCATCATCAAGGAATTCATGCGTTTCAAGGTGCACCTGGAGGGCTCCGTGGACGGCCACGAGTTCGAGATCGAGGGCGAGGGCGAGGGCCGCCCCTACGAGGGCACCCAGAGCGCCAAGCTGAAGGTGACCAAGGGTGGCCCCCTGCCCTTCGCCTGGGACATCCTGTCCCCTCAGTTCATGTACGGCTCCAGGGCCTACGTGAAGCACCCCGCCGACATCCCCGACTACTGGAAGCTGTCCTTCCCCGAGGGCTTCAAGTGGGAGCGCGTGATGAACTTCGAGGATGGCGGCGTGGTGACCGTGACCCAGGACTCCTCCCTGCAGGACGGCGAGTTCATCTACAAGGTGAAGCTGCGCGGCACCAACTTCCCTTCCGACGGCCCCGTAATGCAGAAGAAGACCATGGGCTGGGAGGCCTCCTCCGAGCGGATATACCCCGAGGACGGCGCCCTGAAGGGCGAGATCAAGCAGAGGCTGAAGCTGAAGGACGGCGGCCACTACGACGCTGAGGTCAAGACCACCTACAAGGCCAAGAAGCCCGTGCTGCTGCCCGGCGCCTACAACGTCAACATCAAGATGGACATCACCTCCCACAACGAGGACTACACCATCGTTGAACAGTGCGAACGCGCCGAGGGCCACCATTCCACCGGCGGCATGGACGAGCTGTACAAGTAA

GGTGGTGGTGGTAGC

GCTGCTAACGACGAAAACTACGCTCTGGCTGCT

TAA

This design enables time-dependent exposure of the degron, allowing selective degradation of aged proteins by the host proteolytic system.

The designed insert will be cloned into a pET28 expression vector for protein expression in Escherichia coli BL21(DE3). This vector provides a T7 promoter, ribosome binding site, transcription terminator, and an N-terminal His₆ tag for protein purification. Therefore, only the coding sequence of the fusion protein was designed in Benchling.

the direct link to Benchling: GST_Timer_Degron_Insert

In this homework, ChatGPT helped me structure and write the answers and descriptions clearly, while Cloud AI generated the diagrams comparing Boolean genetic circuits and INNAs, the example illustrating a multilayer perceptron application, and the diagram describing my final project idea proposal.

Sources:

• Brophy, J. A. N., & Voigt, C. A. (2014). Principles of Genetic Circuit Design. Nature Methods, 11(5), 508–520. https://doi.org/10.1038/nmeth.2926 Differences Between SLP and MLP | PDF | Theoretical Computer Science | Machine Learning. (n.d.). Retrieved March 30, 2026, from https://fr.scribd.com/document/858039220/Single-Layer-Perceptron-and-Multilayer-Perceptron
• Gandia, A., van den Brandhof, J. G., Appels, F. V. W., & Jones, M. P. (2021). Flexible Fungal Materials: Shaping the Future. Trends in Biotechnology, 39(12), 1321–1331. https://doi.org/10.1016/j.tibtech.2021.03.002
• Halužan Vasle, A., & Moškon, M. (2024). Synthetic biological neural networks: From current implementations to future perspectives. BioSystems, 237, 105164. https://doi.org/10.1016/j.biosystems.2024.105164
• Hinneburg, H., Gu, S., & Naseri, G. (2025). Fungal Innovations—Advancing Sustainable Materials, Genetics, and Applications for Industry. Journal of Fungi, 11(10). https://doi.org/10.3390/jof11100721
• Lim, H. G., Jang, S., Jang, S., Seo, S. W., & Jung, G. Y. (2018). Design and optimization of genetically encoded biosensors for high-throughput screening of chemicals. Current Opinion in Biotechnology, Analytical Biotechnology, 54, 18–25. https://doi.org/10.1016/j.copbio.2018.01.011
• Mattern, D. J., Valiante, V., Unkles, S. E., & Brakhage, A. A. (2015). Synthetic biology of fungal natural products. Frontiers in Microbiology, 6, 775. https://doi.org/10.3389/fmicb.2015.00775
• Moorman, A., Samaniego, C. C., Maley, C., & Weiss, R. (2019). A Dynamical Biomolecular Neural Network. 2019 IEEE 58th Conference on Decision and Control (CDC), 1797–1802. https://doi.org/10.1109/CDC40024.2019.9030122
• Parhizi, Z., Dearnaley, J., Kauter, K., Mikkelsen, D., Pal, P., Shelley, T., & Burey, P. (Polly). (2025). The Fungus Among Us: Innovations and Applications of Mycelium-Based Composites. Journal of Fungi, 11(8), 549. https://doi.org/10.3390/jof11080549
• Seak, L. C. U., Lo, O. L. I., Suen, W. C.-W., & Wu, M.-T. (2021). Next-generation biocomputing: Mimicking artificial neural network with genetic circuits (p. 2021.03.12.435120). bioRxiv. https://doi.org/10.1101/2021.03.12.435120
• Secret fungi in everyday life | Kew. (n.d.). Retrieved March 30, 2026, from https://www.kew.org/read-and-watch/everyday-fungi-food-medicine Stock, C. H., Harvey, S. E., Ocko, S. A., & Ganguli, S. (2022). Synaptic balancing: A biologically plausible local learning rule that provably increases neural network noise robustness without sacrificing task performance. PLoS Computational Biology, 18(9), e1010418. https://doi.org/10.1371/journal.pcbi.1010418
• Undecided with Matt Ferrell. (2021, June 22). Is Mycelium Fungus the Plastic of the Future? [Video recording]. https://www.youtube.com/watch?v=cApVVuuqLFY
• van der Linden, A. J., Pieters, P. A., Bartelds, M. W., Nathalia, B. L., Yin, P., Huck, W. T. S., Kim, J., & de Greef, T. F. A. (2022). DNA Input Classification by a Riboregulator-Based Cell-Free Perceptron. ACS Synthetic Biology, 11(4), 1510–1520. https://doi.org/10.1021/acssynbio.1c00596
• Wang, X., Chen, Y.-Z., Qiu, X.-D., Chen, L., Teng, Y.-M., Ding, C., Huang, Y.-T., Wang, S.-Y., Liu, S.-Y., Ding, B., Laborda, P., & Zhu, S.-Q. (2026). Bioactivity and mechanisms of Ewingella americana for the control of Alternaria leaf spot on peanut. Physiological and Molecular Plant Pathology, 142, 103088. https://doi.org/10.1016/j.pmpp.2025.103088
• Yang, P., Condrich, A., Lu, L., Scranton, S., Hebner, C., Sheykhhasan, M., & Ali, M. A. (2024). Genetic Engineering in Bacteria, Fungi, and Oomycetes, Taking Advantage of CRISPR. DNA, 4(4), 427–454. https://doi.org/10.3390/dna4040030

Week 09 HW: Cell Free Systems

Homework Part A: General and Lecturer-Specific Questions

1. Main Advantages of Cell-Free Protein Synthesis (CFPS) Over Traditional In Vivo Methods

Cell-free protein synthesis removes the constraints of using living cells. You are working in a test tube, which gives you direct control over the reaction environment without worrying about cell viability. image ref

Flexibility and Control:

• Direct manipulation: You can easily change pH, salt concentration, redox potential, or add detergents, chaperones, or unnatural amino acids at any time. In living cells, these changes would kill the cells or fail to enter.
• No cell walls or membranes: You add DNA directly to the extract. There is no need for transformation, selection, or cell lysis steps. This saves hours or days.
• Toxic protein production: You can synthesize proteins that would kill living cells (e.g., membrane proteins, proteases, toxins).
• Speed and efficiency: Protein production takes 2–4 hours instead of days. All energy goes into making your target protein, not cell growth. Two cases where CFPS is more beneficial than cell production:

1. High-throughput screening of enzyme variants or genetic circuits – Because reactions are fast and can be done in 96- or 384-well plates, you can test hundreds of conditions or mutants in a single afternoon.

2. Production of toxic membrane proteins (e.g., GPCRs, viral ion channels) – These proteins kill E. coli or insect cells when produced in vivo. In CFPS, you can add detergents or nanodiscs directly to the reaction to keep the protein soluble and stable.

3. Main Components of a Cell-Free Expression System and Their Roles

A cell-free system combines cellular machinery with necessary nutrients and energy. Below are the key components and what each does.

image ref

These components are combined either as a crude extract (fast and cheap) or a PURE system (reconstituted from purified components, cleaner but more expensive).

3. Why Energy Provision Regeneration Is Critical and a Method to Ensure Continuous ATP Supply

• Why it is critical?

Cell-free systems lack the metabolic networks of living cells that continuously generate ATP. Translation consumes ATP rapidly – each peptide bond uses 2 ATP equivalents. Without regeneration, ATP drops to zero within 10–20 minutes, and protein synthesis stops. To produce protein for 2–6 hours, you need a way to keep making ATP from ADP.

image ref

Method for continuous ATP supply: Phosphoenolpyruvate (PEP) / Pyruvate Kinase system

• What you add: Phosphoenolpyruvate (PEP) and the enzyme pyruvate kinase.
• How it works: Pyruvate kinase transfers a high-energy phosphate group from PEP to ADP, regenerating ATP and producing pyruvate as a byproduct.
• Why it works well: PEP has a higher phosphate transfer potential than ATP, so the reaction favors ATP formation. It is reliable and commonly used in E. coli systems.

Alternative methods (if PEP causes problems):

If the PEP system presents limitations, other options can be used:

• Creatine phosphate / creatine kinase: Converts ADP + creatine phosphate → ATP + creatine. Very common and stable.
• Glucose / hexokinase or maltodextrin – cheaper but can cause pH drops.

4. Comparison of Prokaryotic vs. Eukaryotic Cell-Free Expression Systems

Cell-free expression systems can be broadly divided into prokaryotic and eukaryotic platforms, and the choice between them mainly depends on the complexity of the target protein.

• Prokaryotic system (e.g., Escherichia coli)

These systems are typically derived from E. coli and are widely used because they are fast, cost-effective, and produce high protein yields in a short time. However, they lack the machinery needed for post-translational modifications such as glycosylation, and they often have difficulty forming correct disulfide bonds and folding complex proteins properly.

• Eukaryotic system (e.g., rabbit reticulocyte lysate, wheat germ extract)

These systems provide a more suitable environment for protein folding. They contain molecular chaperones and can support disulfide bond formation and, in some cases, post-translational modifications. However, they are generally more expensive, slower, and produce lower yields compared to prokaryotic systems.

• Choosing proteins for each system

–> For prokaryotic systems:

The general rule is to choose proteins that are simple, relatively small, and do not require post-translational modifications or complex folding. These proteins should be able to fold easily in the cytoplasm. Based on these criteria, bacterial luciferase is a suitable choice. This enzyme produces a measurable light signal, making it very useful as a reporter protein. It does not require glycosylation and can be efficiently expressed and folded in E. coli, allowing easy detection through luminescence assays.

–> For eukaryotic systems:

The selection criteria are different. Proteins are usually more complex, may contain multiple domains, require disulfide bonds, or need chaperones for correct folding. Some are also membrane proteins and need a suitable environment to function. Membrane proteins, such as G protein-coupled receptors (GPCRs), are good examples. These proteins have complex structures with multiple transmembrane domains and require proper folding machinery and membrane-like conditions. Such requirements cannot be met by prokaryotic systems, while eukaryotic systems can support their correct folding and functionality

5. Designing a Cell-Free Experiment to Optimize Membrane Protein Expression

Optimizing the expression of a membrane protein in a cell-free system requires careful consideration of the protein’s complexity, folding requirements, and membrane integration. Membrane proteins are challenging to produce because of their hydrophobic transmembrane domains, tendency to aggregate, and need for a membrane-like environment and proper chaperones.

• Choosing the right Expression System

The choice of a cell-free system depends on the nature and complexity of the membrane protein:

• Prokaryotic system (e.g., Escherichia coli): Suitable for simpler membrane proteins with few transmembrane domains that do not require complex folding or post-translational modifications. Advantages include fast expression, high yield, and low cost. However, proper folding must be supported using membrane mimics such as liposomes, nanodiscs, or mild detergents.

• Eukaryotic system (e.g., rabbit reticulocyte lysate, wheat germ extract): Preferable for complex membrane proteins with multiple transmembrane domains or disulfide bonds. These systems contain molecular chaperones and provide a more natural folding environment, reducing aggregation and increasing the chance of functional protein production. Limitations include higher cost, slower expression, and lower yields.

• Providing a Membrane-Like Environment

Membrane proteins require an environment that mimics a lipid bilayer. In both prokaryotic and eukaryotic systems, this can be achieved by:

• Adding liposomes or nanodiscs
• Using mild detergents carefully optimized to prevent aggregation

This ensures proper insertion of the protein into a membrane-like environment, which is critical for correct folding and functionality.

• Optimizing Folding and Expression

To further improve expression and functionality:

• Add chaperones if the protein is prone to misfolding
• Adjust reaction conditions such as temperature, Mg²⁺ concentration, and DNA template concentration
• Use a continuous ATP regeneration system (e.g., PEP/pyruvate kinase) to sustain protein synthesis
• Employ a Continuous Exchange Cell-Free (CECF) setup to extend reaction time up to 24 hours. This setup constantly provides fresh energy (ATP/GTP) and removes inhibitory byproducts, which significantly improves protein yield and folding efficiency

Challenges and how to address them:

Optimization checklist:

• Titrate magnesium (8–16 mM) – critical for ribosome function.
• Test temperatures (20°C, 25°C, 30°C, 37°C).
• Try 2–3 different detergents or lipid preparations.
• Run a small-scale (10 µL) screening reaction before scaling up.

6. Low Yield of Target Protein – Three Possible Reasons and Troubleshooting

If your cell-free reaction produces very low yield protein, check these common issues:

Reason 1: Low quality of DNA template

The DNA may contain inhibitors (salts, ethanol, phenol, agarose) or be degraded by nucleases. Without a good template, no mRNA is made.

• Troubleshooting:

✅ Purify DNA using a spin column kit (not just alcohol precipitation).

✅ Avoid using DNA cut from agarose gels – re-extract if necessary.

✅ Check DNA concentration and run an agarose gel to see if it is intact.

✅ Use 10–20 µg of plasmid or 5–10 µg of linear PCR product per 1 mL reaction.

Reason 2: Codon bias (rare codons in the target gene)

If your gene contains many codons that are rare in the host (e.g., human gene expressed in E. coli extract), ribosomes stall or terminate early. This produces truncated or no protein.

• Troubleshooting:

✅ Re-synthesize the gene with codons optimized for your extract (E. coli or wheat germ). Many online tools and services do this.

✅ Use an extract from a strain that supplies extra rare tRNAs (e.g., E. coli Rosetta or BL21 CodonPlus).

✅ Switch to a PURE system, which is less sensitive to codon bias.

Reason 3: Rapid energy depletion

ATP runs out after 30–60 minutes because the energy regeneration system is weak or missing. The reaction stops while plenty of template and amino acids remain.

• Troubleshooting:

✅ Switch to a Continuous Exchange Cell-Free (CECF) format (dialysis membrane or two-chamber system). This constantly supplies fresh energy and removes waste.

✅ Increase the concentration of your energy regeneration components (e.g., double creatine phosphate from 50 mM to 100 mM).

✅ Use a more efficient energy source: PEP/pyruvate kinase or a maltodextrin-based system.

✅ Check the pH after the reaction – if it dropped below 7.0, your energy system may be producing acid. Switch to creatine phosphate (less pH drop).

Additional common reasons (if the above don’t help):

• Protein aggregation: Lower temperature to 20–25°C. Add 0.5% detergent or 1 mM DTT.
• RNase contamination: Use nuclease-free tubes, add RNase inhibitor (e.g., murine RNase inhibitor at 1 U/µL), and wear gloves.
• Wrong magnesium concentration: Test a range from 8 to 16 mM Mg²⁺. Too low and ribosomes dissociate; too high and they lock up.

1. Pick a function and describe it.

a. What would your synthetic cell do? What is the input and what is the output?

My synthetic minimal cell (SMC) is a “killer biosensor” that detects the presence of Staphylococcus aureus and responds by producing and secreting lysostaphin, a specific anti-staphylococcal enzyme.

• Input: AIP-1 (autoinducing peptide-1), a quorum sensing molecule secreted by S. aureus (Group I strains) when it reaches high cell density.
• Output: Lysostaphin (27 kDa zinc metalloprotease from Simulans staphylolyticus), which specifically cleaves the pentaglycine cross-bridges in the S. aureus cell wall, causing bacterial lysis.

Overall function: The SMC acts as a sentinel that detects S. aureus quorum signaling and releases a targeted killer, preventing infection, biofilm formation, and the spread of antibiotic-resistant strains.

b. Could this function be realized by cell-free Tx/Tl alone, without encapsulation?

No. Without encapsulation, the cell-free reaction would produce lysostaphin immediately and continuously, regardless of whether AIP-1 is present. The SMC would release its output constitutively, wasting the enzyme and providing no sensing function. Encapsulation creates a barrier that allows the system to wait for the input signal before producing the output. Additionally, without a membrane: The membrane-bound receptor AgrC could not be properly inserted and oriented and Lysostaphin would diffuse away uncontrollably instead of being released only after detecting S. aureus.

c. Could this function be realized by genetically modified natural cell?

Yes, in principle, but with significant drawbacks compared to a synthetic minimal cell (SMC). Natural GMOs can grow, divide, and potentially spread in the environment, and they may transfer genes to other bacteria through horizontal gene transfer. They can also mutate over time and lose their function, and the produced antibacterial molecule (e.g., lysostaphin) might harm the host cell itself. In contrast, SMCs do not replicate, cannot transfer genes, and do not evolve, making them safer and more stable. Additionally, their activity is more controlled, since the toxic compound is produced only when needed and released outward, which makes SMCs more suitable for applications such as medical treatments or topical use.

d. Describe the desired outcome of your synthetic cell operation.

In the presence of S. aureus (which secretes AIP-1), the synthetic cell detects AIP-1 via the membrane-bound AgrC receptor. This triggers a phosphorylation cascade that activates AgrA, which then binds the P2 promoter and drives transcription of the lysostaphin gene. Lysostaphin is produced inside the vesicle and secreted into the environment. The released lysostaphin specifically cleaves the pentaglycine bridges in the S. aureus cell wall, causing bacterial lysis and death.

In the absence of AIP-1 (no S. aureus), the synthetic cell remains inactive. The P2 promoter is “off” (no leak), and no lysostaphin is produced. This ensures the toxin is only made when and where it is needed.

2. Design all components that would need to be part of your synthetic cell.

a. What would the membrane be made of?

The membrane needs to be stable but also allow the AgrC histidine kinase (a transmembrane protein) to insert properly. A suitable choice is liposomes composed of:

• POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) – 60 mol% → Main structural lipid of the membrane
• Cholesterol – 30 mol% → Increases membrane stability and reduces leakage
• DOPG (1,2-dioleoyl-sn-glycero-3-phospho-(1’-rac-glycerol)) – 10 mol% → Adds negative charge, which helps the insertion and function of membrane proteins like AgrC

b. What would you encapsulate inside? Enzymes, small molecules.

Inside the synthetic cell, i would encapsulate the basic components needed for protein production and function. First, a cell-free transcription–translation system from Escherichia coli is included, which contains all the machinery such as ribosomes, tRNAs, enzymes, and T7 RNA polymerase to make proteins.

I also add the DNA templates: agrC and agrA genes (from Staphylococcus aureus) under a constitutive promoter to sense the signal and activate the response, and the lysostaphin (lys) gene (from Staphylococcus simulans) under the P2 promoter to produce the antibacterial protein. A secretion signal is fused to the lys gene so the protein can be exported outside the cell.

In addition, small molecules like:

• ATP, GTP, CTP, UTP (nucleotide triphosphates for transcription)
• 20 amino acids (building blocks for protein synthesis)
• Creatine phosphate + creatine kinase (energy regeneration system)
• Magnesium acetate (10–14 mM) – critical for ribosome function
• Potassium glutamate (100–150 mM) – maintains ionic strength
• DTT (1–2 mM) – maintains reducing environment
• RNase inhibitors – protect mRNA from degradation

c. Which organism will your Tx/Tl system come from?

The Tx/Tl system will come from a bacterial source, specifically an Escherichia coli extract. This is because the AgrC/AgrA system is naturally bacterial and works well in an E. coli cell-free system, where AgrC can insert into liposomes properly. In addition, lysostaphin is a bacterial enzyme that does not require complex modifications, so it can be produced efficiently in this system. Finally, using a bacterial extract is simpler, faster, and cheaper than using a mammalian system, which is not needed in this case.

d. How will your synthetic cell communicate with the environment?

This synthetic cell communicates with its environment in a simple and efficient way using natural bacterial mechanisms:

For input, the signaling molecule AIP-1 does not need to enter the cell; instead, it binds directly to the external part of AgrC, a membrane protein embedded in the liposome. This means the sensor is already on the surface, so no channels are needed.

For output, lysostaphin (a relatively large protein, about 27 kDa) cannot pass through the membrane by diffusion. To solve this, a secretion signal peptide is added to lysostaphin, which directs it to the membrane during its synthesis. The protein is then transported across the membrane through the SecYEG translocon, a natural protein channel present in the Escherichia coli extract. This allows the protein to be released outside the synthetic cell in a controlled and efficient way, without needing artificial pores.

3. Experimental details

a. List all lipids and genes (specific names).

• Lipids:

• Genes:

Cell-free Tx/Tl system: All machinery for transcription and translation from E. coli extract (ribosomes, tRNAs, aminoacyl-tRNA synthetases, initiation/elongation/termination factors, T7 RNA polymerase).

b. How will you measure the function of your system?

• Measurement 1: AIP-1 sensing (dose–response)

The synthetic cells can be exposed to different concentrations of AIP-1. After a few hours, lysostaphin production is measured using methods like ELISA, Western blot, or an enzyme activity assay. If the system works properly, higher AIP-1 levels should lead to higher lysostaphin production.

• Measurement 2: Lysostaphin production (fluorescent reporter)

The lysostaphin gene can be replaced with GFP (green fluorescent protein) under the same promoter. The synthetic cells can then be monitored over time using a plate reader to measure fluorescence. Higher fluorescence indicates stronger gene expression.

• Measurement 3: Killing of Staphylococcus aureus (functional assay)

The synthetic cells can be incubated with live bacteria in culture medium. After several hours, bacterial growth can be measured using OD600, colony counting (CFU), or live/dead staining. Reduced growth shows that the system is effective.

• Measurement 4: Secretion efficiency

The synthetic cells can be centrifuged to separate them from the surrounding liquid. Lysostaphin activity is then measured both in the supernatant (outside) and inside the cells. A good system will show most of the protein in the supernatant.

• Measurement 5: Promoter leakiness (control test)

The synthetic cells can be tested without adding AIP-1 to check background expression. Ideally, very little lysostaphin should be produced. If significant production is observed, the promoter may be leaky and require optimization.

1. One-sentence pitch

A wall paint containing synthetic minimal cells that detect toxic mold signals in damp walls and produce enzymes to neutralize mycotoxins and inhibit mold growth.

image ref

2. How will the idea work?

The paint is embedded with microcapsules containing freeze-dried synthetic minimal cells (SMCs). When the wall becomes damp, the SMCs are activated by chemical signals released by mold, such as those from Stachybotrys chartarum. Once triggered, the SMCs produce enzymes or antimicrobial proteins that either degrade mycotoxins or prevent further mold growth. This creates a self-protecting coating that actively reduces mold and mycotoxin levels in real-time, improving indoor air safety.

3. What societal challenge or market need does this address?

Toxic wall moisture is a serious indoor environmental problem. Persistent dampness encourages growth of black mold, which releases mycotoxins harmful to human health, causing respiratory issues, chronic fatigue, and neurological problems. Current paints only act as passive barriers and do not remove toxins. This smart paint provides active protection, reducing health risks and the need for costly remediation.

4. How will you address limitations of cell-free systems?

The SMCs are freeze-dried within protective microcapsules, remaining inactive until moisture activates them. Microcapsules shield the system during storage and paint application. Activation only occurs when mold is present, ensuring efficient use. The one-time-use limitation is addressed by applying fresh paint layers during regular maintenance, keeping the wall continuously protected.

1. One-sentence summary pitc

We will use a freeze-dried cell-free system to test how microgravity affects protein production using a GFP reporter, providing insight into reduced collagen synthesis in space.

2. How the idea works

Freeze-dried cell-free reactions containing a GFP reporter gene will be prepared in sealed chambers. In space, they will be rehydrated and incubated under microgravity conditions. GFP fluorescence will act as a direct indicator of protein synthesis efficiency. By comparing fluorescence levels between space and Earth conditions, we can determine whether microgravity directly affects the molecular machinery responsible for producing proteins such as collagen.

4. Societal challenge / market need

Long-duration space missions lead to bone loss and tissue weakening in astronauts, partly due to reduced production of structural proteins like collagen. Understanding whether this reduction is caused by fundamental limits in protein synthesis will help develop countermeasures for bone loss, injury prevention, and tissue regeneration, improving astronaut health during missions to Mars and beyond.

6. Limitation of cell-free reactions and how to address them

Cell-free reactions are single-use and require activation by water. To overcome this, we will freeze-dry the reactions in sealed chambers, ensuring long-term stability. The experiment will be activated by rehydration in space, allowing controlled and efficient protein production measurements under microgravity conditions.

8. Molecular or genetic target

Green Fluorescent Protein (GFP) gene used as a reporter to measure protein synthesis efficiency linked to collagen-related biological processes.

10. How the target relates to the space biology challenge

Collagen is essential for maintaining bone and tissue structure, but its production decreases in microgravity. Instead of directly expressing collagen, which is complex, GFP is used as a reporter to measure overall protein synthesis efficiency. If microgravity reduces GFP production, it suggests that the basic machinery needed to produce proteins like collagen is affected. This helps determine whether tissue weakening in space is caused by direct physical effects on protein production or by cellular regulation, providing clearer insight into astronaut health challenges.

12. Hypothesis or research goal

We hypothesize that microgravity reduces protein synthesis efficiency, which contributes to decreased production of structural proteins such as collagen in astronauts. The goal is to measure GFP production in a cell-free system under microgravity and Earth conditions. Since cell-free systems isolate transcription and translation from cellular signaling, any observed decrease in GFP fluorescence would indicate that physical factors—such as altered diffusion, molecular interactions, or protein folding—directly impact protein synthesis. This would suggest that microgravity imposes fundamental constraints on biological processes, helping explain tissue weakening. The results could guide the development of targeted countermeasures to maintain astronaut health during long-duration missions.

14. Experimental plan

Freeze-dried BioBits® reactions containing GFP DNA will be used. Samples include: (1) microgravity test reactions, (2) Earth-based positive controls, and (3) negative controls without DNA. Reactions will be rehydrated and incubated using the miniPCR®. GFP fluorescence will be measured with the P51 Molecular Fluorescence Viewer. Fluorescence intensity will be compared between conditions to determine whether microgravity reduces protein synthesis efficiency.

• For this homework, I used DeepSeek and Google as sources of information. ChatGPT was used to improve the structure and clarity of the writing, while Cloud AI was used to generate the illustration of the synthetic minimal cell function.

Sources:

• A Comparative Guide: Prokaryotic vs. Eukaryotic Cell-Free Expression Systems for Eukaryotic Proteins—CD Biosynsis. (n.d.). Retrieved April 5, 2026, from https://www.biosynsis.com/a-comparative-guide-prokaryotic-vs-eukaryotic-cell-free-expression-systems-for-eukaryotic-proteins.html
• Advantages of Cell-Free Protein Expression. (n.d.). Retrieved April 5, 2026, from https://www.biocompare.com/Editorial-Articles/594727-Advantages-of-Cell-Free-Protein-Expression/
• Biology Goal. (2023, November 29). Autoinducing Peptides (AIP) in Staphylococcus Quorum Sensing [Video recording]. https://www.youtube.com/watch?v=nXXVLU7Dh-w
• Black Mold: Health Risks, Costs & Identification Guide. (n.d.). Mold Busters. Retrieved April 5, 2026, from https://www.bustmold.com/resources/about-mold/types-of-mold/black-mold/
• Brookwell, A., Oza, J. P., & Caschera, F. (2021). Biotechnology Applications of Cell-Free Expression Systems. Life, 11(12), 1367. https://doi.org/10.3390/life11121367
• Cell-Free Expression Support—Troubleshooting—DZ. (n.d.). Retrieved April 5, 2026, from https://www.thermofisher.com/ng/en/home/technical-resources/technical-reference-library/protein-expression-support-center/cell-free-expression-support/cell-free-expression-support-troubleshooting.html
• Cell-Free Protein Expression—DZ. (n.d.). Retrieved April 5, 2026, from https://www.thermofisher.com/ng/en/home/life-science/protein-biology/protein-biology-learning-center/protein-biology-resource-library/pierce-protein-methods/cell-free-protein-expression.html
• Cell-Free Systems for Protein Production: Advantages Over Living Cells. (n.d.). Cytion. Retrieved April 5, 2026, from https://www.cytion.com/ca/About-Cytion/Knowledge-Hub/Blog/Cell-Free-Systems-for-Protein-Production-Advantages-Over-Living-Cells/
• Chong, S. (2014). Overview of Cell-Free Protein Synthesis: Historic Landmarks, Commercial Systems, and Expanding Applications. Current Protocols in Molecular Biology / Edited by Frederick M. Ausubel … [et Al.], 108, 16.30.1-16.30.11. https://doi.org/10.1002/0471142727.mb1630s108
• Contact us via LiveChat! (n.d.). Retrieved April 5, 2026, from https://secure.livechatinc.com/customer/action/open_chat?license_id=7496091&group=0&embedded=1&widget_version=3&unique_groups=0&organization_id=eef41f83-49a8-4805-8843-da337d811599&use_parent_storage=1&x-region=us-south1
• Don’t Ignore Mould In The Home | EnviroVent. (n.d.). Retrieved April 5, 2026, from https://www.envirovent.com/blog/dont-ignore-mould-in-the-home/ Jiang, S., Çelen, G., Glatter, T., Niederholtmeyer, H., & Yuan, J. (2024). A cell-free system for functional studies of small membrane proteins. The Journal of Biological Chemistry, 300(11), 107850. https://doi.org/10.1016/j.jbc.2024.107850
• Jiang, Y., Geng, M., & Bai, L. (2020). Targeting Biofilms Therapy: Current Research Strategies and Development Hurdles. Microorganisms, 8(8), 1222. https://doi.org/10.3390/microorganisms8081222
• Khambhati, K., Bhattacharjee, G., Gohil, N., Braddick, D., Kulkarni, V., & Singh, V. (2019). Exploring the Potential of Cell-Free Protein Synthesis for Extending the Abilities of Biological Systems. Frontiers in Bioengineering and Biotechnology, 7, 248. https://doi.org/10.3389/fbioe.2019.00248
• Lysostaphin—An overview | ScienceDirect Topics. (n.d.). Retrieved April 5, 2026, from https://www.sciencedirect.com/topics/neuroscience/lysostaphin Meyer, C., Arizzi, A., Henson, T., Aviran, S., Longo, M. L., Wang, A., & Tan, C. (2025). Designer artificial environments for membrane protein synthesis. Nature Communications, 16(1), 4363. https://doi.org/10.1038/s41467-025-59471-1
• NASA OSDR: Open Science for Life in Space. (n.d.). Retrieved April 6, 2026, from https://osdr.nasa.gov/bio/repo/search?q=&data_source=cgene,alsda,esa&data_type=study
• Optimization of Energy Regeneration in Cell-free Protein Synthesis. (n.d.). PatSnap Eureka. Retrieved April 5, 2026, from https://eureka.patsnap.com/report-optimization-of-energy-regeneration-in-cell-free-protein-synthesis
• Seki, E., Matsuda, N., Yokoyama, S., & Kigawa, T. (2008). Cell-free protein synthesis system from Escherichia coli cells cultured at decreased temperatures improves productivity by decreasing DNA template degradation. Analytical Biochemistry, 377(2), 156–161. https://doi.org/10.1016/j.ab.2008.03.001
• Sitaraman, K., Esposito, D., Klarmann, G., Le Grice, S. F., Hartley, J. L., & Chatterjee, D. K. (2004). A novel cell-free protein synthesis system. Journal of Biotechnology, 110(3), 257–263. https://doi.org/10.1016/j.jbiotec.2004.02.014
• Szaflarski, W., & Nierhaus, K. H. (2007). Question 7: Optimized Energy Consumption for Protein Synthesis. Origins of Life and Evolution of Biospheres, 37(4), 423–428. https://doi.org/10.1007/s11084-007-9091-4
• Troubleshooting Guide for NEBExpressTM Cell-free E. coli Protein Synthesis System (NEB #E5360) | NEB. (n.d.). Retrieved April 5, 2026, from https://www.neb.com/en/tools-and-resources/troubleshooting-guides/troubleshooting-guide-for-nebexpress-cell-free-e-coli-protein-synthesis-system-neb-e5360
• Williams, P., Hill, P., Bonev, B., & Chan, W. C. (2023). Quorum-sensing, intra- and inter-species competition in the staphylococci. Microbiology, 169(8), 001381. https://doi.org/10.1099/mic.0.001381
• Wu, J. A., Kusuma, C., Mond, J. J., & Kokai-Kun, J. F. (2003). Lysostaphin disrupts Staphylococcus aureus and Staphylococcus epidermidis biofilms on artificial surfaces. Antimicrobial Agents and Chemotherapy, 47(11), 3407–3414. https://doi.org/10.1128/AAC.47.11.3407-3414.2003

Week 10 HW: Imaging And Measurement

My project aims to express the carbon monoxide dehydrogenase (CODH) pathway from Oligotropha carboxidovorans in Nicotiana tabacum (tobacco) using a two-plasmid system. I need to measure whether the system works at every level — from DNA integration to enzyme function to plant health. Below i included what I will measure, how I will measure it, and the technologies I will use:

1. Confirming DNA Integration and Sequence

What I measure: Whether the seven CODH genes are present in the tobacco genome and whether their sequences are correct.

How I measure it:

• Genomic PCR: Extract DNA from leaves, design primers specific to each of my seven codon-optimized genes, run PCR, and look for bands on an agarose gel.
• Border-specific PCR: Use one primer in the T-DNA border (LB or RB) and one primer in my gene to confirm the entire T-DNA integrated.
• Sanger sequencing: Send PCR products to a sequencing facility, align the returned sequences against my Benchling design using SnapGene. image ref

Technologies: PCR thermocycler, agarose gel electrophoresis, UV transilluminator, Sanger sequencing service, sequence alignment software.

2. Confirming mRNA Transcription

What I measure: Whether the seven genes are being transcribed into mRNA, and whether the three structural subunits (CoxL, CoxM, CoxS) are expressed at balanced levels.

How I measure it:

• Extract total RNA from leaves using an RNA extraction kit.
• Treat with DNase to remove genomic DNA.
• Convert mRNA to cDNA using reverse transcriptase.
• Run qPCR with gene-specific primers and SYBR Green.
• Include reference genes for normalization.
• Compare Ct values across the three structural subunits.

Technologies: RNA extraction kit, DNase, reverse transcriptase, qPCR machine, SYBR Green. image ref

3. Confirming Protein Presence and Assembly

What I measure: Whether CoxL, CoxM, and CoxS are present, whether the chloroplast transit peptide was cleaved, and whether the three subunits assemble into the complex.

How I measure it:

• Isolate intact chloroplasts using Percoll gradient centrifugation.
• Lyse chloroplasts gently and perform Co-IP using anti-FLAG magnetic beads (FLAG is on CoxS).
• Elute with FLAG peptide.
• Split eluate: run on Tricine-SDS-PAGE (silver stain) to see individual subunits at 88 kDa (CoxL), 32 kDa (CoxM), and 18 kDa (CoxS).
• Run on Blue Native PAGE (Coomassie stain) to see the assembled complex at ~280 kDa.
• For maturation proteins: run anti-FLAG Western (detects CoxD) and anti-His Western (detects His-tagged CoxF) on total chloroplast extract.

Technologies: Ultracentrifuge, anti-FLAG magnetic beads, PAGE equipment, silver stain, Coomassie stain, Western blot transfer system, chemiluminescence imager.

image ref

image ref

image ref

4. Confirming Chloroplast Targeting

What I measure: Whether the used chloroplast transit peptides direct proteins to the chloroplast.

How I measure it:

• Build a separate reporter construct: promoter + CTP + GFP + terminator.
• Transform into tobacco, select on hygromycin.
• Take fresh leaf samples, mount on slides with water.
• Observe under confocal microscope: GFP channel (green) and chlorophyll autofluorescence (red).
• Calculate Pearson’s correlation coefficient using ImageJ (target >0.7).

Technologies: Confocal laser scanning microscope, ImageJ software. image ref

5. Confirming CO Oxidation Activity

What I measure: Whether the assembled CODH enzyme can oxidize CO to CO₂.

How I measure it:

• Gas phase (whole plant): Place transformed plant in sealed transparent chamber, inject CO gas, record CO concentration in separate timelines using electrochemical CO sensor.
• Methylene blue (purified enzyme): Purify CODH complex via anti-FLAG Co-IP, add to reaction with methylene blue and CO in anaerobic cuvette, measure absorbance at 600 nm at different timelines. Calculate specific activity (μmol CO/min/mg protein).

Technologies: Sealed gas chamber, electrochemical CO sensor, spectrophotometer, anaerobic cuvettes. image ref

6. Confirming Cofactor Incorporation

What I measure: Whether the CODH complex contains molybdenum, copper, and iron-sulfur clusters.

How I measure it:

• ICP-MS: Send purified CODH complex to core facility. Measure Mo, Cu, and Fe content. Calculate metal-to-protein stoichiometry.
• UV-Vis spectroscopy: Measure absorbance spectrum of purified complex from 300-700 nm. Look for peak at 420 nm (Fe-S clusters). Technologies: ICP-MS instrument, UV-Vis spectrophotometer.

7. Confirming Electron Transfer Compatibility

What I measure: Whether electrons from CODH go to the photosynthetic electron transport chain or leak to oxygen.

How I measure it:

• Compare CO oxidation rate in light vs. dark using the gas chamber setup.
• Calculate light:dark ratio. Ratio >2 indicates electrons go to photosynthetic chain (requires light). Ratio ~1 indicates electrons go directly to oxygen (oxidative stress risk).

Technologies: Sealed gas chamber, electrochemical CO sensor, light source, dark cover.

8. Monitoring Plant Health

What I measure: Whether expressing CODH causes stress or benefits photosynthesis.

How I measure it:

• Chlorophyll fluorescence (Fv/Fm): Dark-adapt leaf for 20 minutes, measure with PAM fluorometer. Healthy plant = 0.80-0.83.
• CO₂ assimilation: Use infrared gas analyzer (IRGA) to measure net CO₂ uptake by leaf. Compare transformed vs. wild-type.
• Biomass: Dry plants at 70°C for 48 hours, weigh shoot and root. Compare transformed vs. wild-type.
• ROS detection: Stain leaf discs with NBT (detects superoxide, turns blue) and DAB (detects H₂O₂, turns brown). Photograph and quantify staining.

Technologies: PAM fluorometer, LI-COR IRGA, analytical balance, NBT/DAB staining, light microscope, ImageJ. image ref

image ref Histochemical detection of H2O2 by DAB staining (a), superoxide radical by NBT staining (b)

9. Monitoring Silencing Over Time

What I measure: Whether expression remains stable across generations (T0 → T1 → T2).

How I measure it:

• Grow T0 plants (primary transformants), measure mRNA by RT-qPCR.
• Self-pollinate T0 to obtain T1 seeds.
• Grow T1 plants, repeat RT-qPCR.
• Grow T2 plants, repeat RT-qPCR.
• Calculate silencing index = Expression(T1)/Expression(T0). Index >0.8 = stable.

Technologies: RT-qPCR, plant growth facilities.

1. Theoretical Molecular Weight Calculation

The theoretical molecular weight of eGFP was calculated using the online tool ExPASy Compute pI/Mw tool (Swiss Institute of Bioinformatics). The full amino acid sequence of eGFP, including the C-terminal His-tag (HHHHHH) and linker (LE), was entered into the calculator. The computed molecular weight obtained from this tool was: 28006.60 Da This value was used as the reference theoretical mass for comparison with the experimentally determined molecular weight obtained from LC-MS analysis.

2. Calculating the Experimental Molecular Weight (MW)

   2.1. Identification of Adjacent Charge States

Step 1: Identifing Two Adjacent Peaks from Figure 01

let’s use the following values from this figure:

m/z(n) = 903.7148

m/z(n+1) = 875.4421

Step 2: Solve for the Charge State (n)

The relationship between the two peaks is:

n = (m/z(n+1)−1)/(m/z(n)−m/z(n+1))

Let’s plug in our example numbers:

n = 875.4421 – 1 / (903.7148 – 875.4421)

n = 874.4421/ (28.2727)

n = 30.93

Since the charge state must be a whole integer, we round this to the nearest whole number. Therefore, n = 31. This means the peak at m/z 903.7148 is the +31 charge state. From this value, we can extract the charge state for the second adjacent peak: n+1 = 32, which means the peak at m/z 875.4421 is the +32 charge state.

2.2. Calculating (MW)

Now that we know n, we can calculate M using the following formula, which accounts for the mass of the protons that are adding the charge:

m/z = MW of protein + mass of all added protons / total number of charges (n)

MW of protein = (m/z x total number of charges (n)) – (mass of all added protons)

Note: mass of all added protons is: the total number of charges (n) x the mass of a proton (approximately 1.0078 Da) (H)

Using the charge state of the first peak:

MW = (m/z(n) x zn) – (zn x H)

MW = (875.4421 x 32) - (32 x 1.0078)

MW = 28014.1472 – 32.2496

MW = 27981.8976 Da

Using the second peak, I found: MW = 27983.917 Da, so the average experimental molecular weight of this protein is ≈ 27982.9073 Da By comparing the experimental result, we just calculated to the theoretical weight from Step 1, the resulted experimental molecular weight is approximate to the theoretical value calculated 28006.60 Da.

2.3. Calculating the Measurement Accuracy

The formula for accuracy is:

Accuracy = (|Calculated MW – Theoretical MW|) / (Theoretical MW) x 1,000,000

Accuracy = (|27982.9073 – 28006.60 |) / (28006.60) x 1,000,000

Accuracy = (9.56) / (28006.60) x 1,000,000

Accuracy = 845.96 ppm > 50 ppm

The measured accuracy (~846 ppm) is significantly higher than the acceptable threshold of 50 ppm.

This deviation is most likely due to instrumental factors, such as imperfect calibration of the mass spectrometer, which can lead to slight inaccuracies in measured m/z values. Since the theoretical mass was calculated directly from the provided amino acid sequence, it is unlikely that the discrepancy arises from errors in the protein sequence or its expression.

3. Charge State Determination (Zoomed Peak)

No, we cannot. The inability to determine the charge state from the zoomed-in peak is mainly due to the relationship between isotope spacing and instrument resolution. Proteins are made of atoms that exist in different isotopic forms, such as 12C and 13C, which create small differences in mass. In their neutral state, these isotopes are separated by about 1 Da. However, in mass spectrometry, we measure the mass-to-charge ratio (m/z), so the space between isotopic peaks becomes (1/z), where (z) is the charge. This means that as the charge increases, the spacing between peaks becomes smaller.

For large proteins like eGFP (approximately 28 kDa), the charge state is relatively high. As a result, the spacing between isotopic peaks becomes extremely small. For example, if the charge is around (z ≈ 19), the spacing between peaks is only about 0.05 (m/z). These very small differences are difficult for the instrument to detect.

The limitation comes from the resolution of the mass spectrometer. Resolution refers to the ability of the instrument to distinguish between two very close peaks. In this case, the required spacing (around 0.05 (m/z)) is smaller than what the instrument can clearly resolve. Instead of observing distinct isotopic peaks, the signals merge together and appear as a single broad and jagged peak.

Because the individual isotope peaks are not visible, it is not possible to measure their spacing and determine the charge state directly. Therefore, an alternative approach, such as the adjacent charge state method, must be used to calculate the charge and molecular weight.

1. the difference between native and denatured protein conformations

What happens when a protein unfolds?

In its native state, a protein such as eGFP is folded into a compact three-dimensional structure (often described as a beta-barrel). In this conformation, many basic amino acid residues (such as lysine and arginine) are buried inside the protein and are not easily accessible. When the protein becomes denatured, typically due to acidic or organic solvents, it loses this structure and unfolds into a more extended chain. This unfolding exposes a larger surface area and reveals previously hidden basic sites.

How is this determined with a Mass Spectrometer?

Mass spectrometry detects the charge-to-mass ratio (m/z). Because an unfolded protein has more surface area and more exposed basic sites, it can pick up a much higher number of protons (H+) during Electrospray Ionization (ESI). So, in simple way:

• Native (folded) protein: Compacted structure → Fewer exposed basic sites → Binds fewer protons → low charge state (low z)
• Denatured (unfolded) protein: Extended, flexible structure → More exposed basic sites → Binds more protons → high charge state (high z)

Changes Observed in the Mass Spectrum (Figure 2)

These differences in charge directly affect the mass-to-charge ratio (m/z): Since m/z= m x 1/z, a higher charge (z) results in a lower m/z

• Denatured (in Green): The peaks are shifted to the left (lower m/z). This is because the charge (z) is high. Since z is the denominator in m/z, a higher charge results in a lower m/z value. The distribution is also very broad, indicating many different charge states are possible for a flexible, unfolded chain.
• Native (in Red): The peaks are shifted to the right (higher m/z). A folded protein is “shielded,” so it can only pick up a few protons. Fewer protons mean a lower z, which results in a much higher m/z value.

2. When analyzing Figure 3 of the native mass spectrum of eGFP, I initially noticed a possible confusion in the question. The prompt refers to a zoomed-in region around m/z ~2800, however, the zoomed image shown in the figure is actually centered on the peak at m/z ~2545, not 2800. Because of this mismatch, I decided to carefully analyze the figure in two complementary ways to ensure a complete and correct interpretation.

Case 1: Analysis of the zoomed-in region (m/z ~2545)

Although the question mentions ~2800, the zoomed panel clearly shows the peak at m/z ≈ 2545. In this zoomed region, individual isotopic peaks are visible. This is important because isotopic resolution allows us to determine the charge state using peak spacing.

→ Method used: isotopic spacing

• In mass spectrometry, isotopic peaks of a given charge state are separated by: Δ(m/z) = 1/z
• From the zoomed spectrum, the spacing between adjacent isotopic peaks: Looking at the labeled values around ~2544–2545:

2544.8552 → 2544.7637 ≈ 0.0915 m/z

2544.7637 → 2544.6719 ≈ 0.0918 m/z

Average spacing ≈ 0.092 m/z

Calculation: z = 1/ 0.092 ≈ 10.86

Considering the measured values shown in the figure (around 2545.03–2545.22), the spacing is most consistent with: +11

Case 2: Interpretation of the peak at m/z ~2800 (main spectrum)

In the full (non-zoomed) spectrum, there is also a broader peak around m/z ~2800, but:

It is not zoomed in and the isotopic pattern is not resolved, Therefore, charge state cannot be directly read from spacing in this region What I did to solve this

Since isotopic resolution is not available at ~2800, I used the adjacent peak relationship between charge states in native mass spectrometry:

• Neighboring charge states follow predictable shifts in m/z
• Using the relationship between the 2545 peak and the 2799 peak:

n = (m/z(n+1)−1)/(m/z(n)−m/z(n+1))

n = 2545 – 1 / (2799 – 2545)

n = 2544 / (254)

n = 10.01

This indicates that the peak at ~2800 corresponds to the next charge state after +10.

1. Identification of Cleavage Sites (K and R residues)

To predict the tryptic digestion pattern of eGFP, I first analyzed the amino acid sequence and counted the number of lysine (K) and arginine (R) residues, since trypsin cleaves specifically after these amino acids. From the sequence analysis:

• Number of Lysine (K): 20
• Number of Arginine (R): 6

2. Prediction of Tryptic Peptides

To determine the number of peptides generated after digestion, I used the ExPASy PeptideMass tool by inputting the full eGFP sequence and applying trypsin cleavage conditions. The tool predicted a total of: 19 peptides

The theoretical molecular weight of eGFP used for reference was: Mw (average mass): 28006.60 Da

3. Chromatographic Peak Analysis

From the total ion chromatogram (TIC) shown in Figure 5a, I counted the number of peaks between 0.5 and 6 minutes, considering only peaks with a relative intensity greater than 10%. The number of observed peaks was: 18

4. Comparison Between Predicted Peptides and Observed Peaks

The theoretical digestion predicted 19 peptides, while the chromatogram shows 18 peaks. There is slight difference between the theoretical digestion and the chromatogram, but overall, the numbers are very close, indicating good agreement between theoretical prediction and experimental data.

5. Peptide Mass and Charge Determination

From Figure 5b, the most abundant peak was observed at: m/z = 525.76

By analyzing the isotope spacing:

• 526.25918 – 525.76712 = 0.49
• 526.76845 - 526.25918 = 0.50

Δm/z ≈ 0.5 → z = 1/ Δm/z = 1/ 0.5 = 2

Thus, the peptide is doubly charged (z=2).

The molecular weight was calculated using:

MW = (m/z x z) – (z x H)

MW = (525.76 x 2) – (2 x 1.0078)

MW = 1049.5044 ≈ 1050

6. Peptide Identification

Using the predicted peptide list from the ExPASy tool, I compared the calculated experimental mass (1049.5044 Da) with theoretical peptide masses. The closest match was:

• Peptide sequence: FEGDTLVNR with Theoretical mass: 1050.5214 Da

This confirms that the detected peptide corresponds to this sequence.

Then the mass accuracy was calculated using:

Accuracy = (|Calculated MW – Theoretical MW|) / (Theoretical MW) x 1,000,000

Accuracy = (|1049.5044 – 1050.5214|) / (1050.5214) x 1,000,000

Accuracy = (1.017) / (1050.5214) * 1,000,000

Accuracy = 968.09 ppm > 10

7. Sequence Coverage (Figure 6) From the coverage map shown in Figure 6, approximately: 88% of the eGFP sequence was identified This high coverage indicates that most of the protein sequence was successfully confirmed through peptide mapping.

Bonus part:

8. Peptide Sequence Confirmation Using Fragmentation

To confirm the identity of the peptide, I used the mass obtained from the LC-MS analysis and matched it with the predicted tryptic peptides. The peptide with the closest theoretical mass was identified as FEGDTLVNR, with a theoretical mass of 1050.52149 Da. To validate this identification, I used a fragmentation prediction tool to generate the expected b- and y-ion fragments of this peptide.
the resulted fragments are as following:
I then compared these predicted fragments with the experimental MS/MS spectrum shown in Figure 5c. Several peaks in the spectrum matched the predicted fragments, especially the y-ions, like :1050.52149; 903.45308; 774.41049; 602.36208, which confirms that the sequence FEGDTLVNR is correct. The experimental mass of the peptide was 1050.52438 Da, which is very close to the theoretical value. I calculated the mass accuracy using the ppm formula and obtained: accuracy ≈2.75 ppm
This very low error (well below 10 ppm) indicates high measurement accuracy and strong agreement between experimental and theoretical data.

9. Sequence Coverage and Protein Confirmation

To evaluate whether the results confirm the identity of the protein, I analyzed the sequence coverage shown in Figure 6. The coverage percentage was approximately 88%, indicating that a large portion of the eGFP sequence was successfully identified. Additionally, the identified peptide FEGDTLVNR (positions 115–123) is located within the covered regions of the sequence, confirming that this peptide contributes to the overall sequence identification. This high sequence coverage, along with the accurate peptide identification and fragmentation matching, confirms that the analyzed protein is indeed eGFP. Although some regions are not covered (likely due to peptides that are too small or poorly ionized), the overall results provide strong confidence in the protein identification.

Oligomer Identification Using CDMS

To determine the oligomeric states of Keyhole Limpet Hemocyanin (KLH), I used the subunit masses provided in Table 1 and calculated the expected total mass for each oligomeric form. The given subunit masses are:

• 7FU = 340 kDa
• 8FU = 400 kDa

Mass Calculations

For each oligomer, the total mass was calculated by multiplying the subunit mass by the number of subunits:

• 7FU Decamer (10 subunits): 10×340 = 3400kDa = 3.4MDa
• 8FU Didecamer (20 subunits): 20×400 = 8000kDa = 8 MDa
• 8FU 3-Decamer (30 subunits): 30×400 = 12000kDa = 12 MDa
• 8FU 4-Decamer (40 subunits): 40×400 = 16000kDa = 16 MDa

Note: While assigning the oligomeric peaks in the CDMS spectrum (Figure 7), I noticed that for the first three oligomers there are clear red peaks, but for the fourth one (~16 MDa), there is only a small blue signal without a corresponding red peak. This made me question why there are two different colors in the spectrum and why the fourth oligomer does not have a red peak.

After looking into this, I understood that the two colors represent different types of data:

• The blue line corresponds to the raw signal detected by the instrument. It includes all detected ions and therefore appears noisy and irregular.
• The red peaks correspond to a fitted model (Gaussian fit) generated by the software. This fit is applied to the raw data to determine the most accurate position (center) of each mass peak.

This means that the red peaks represent the most reliable mass values, while the blue signal shows all detected data, including weaker or less clear signals.

Using this understanding, I assigned the oligomers as follows:

• The peak at 3.4 MDa (red) corresponds to the 7FU decamer
• The peak at 8.33 MDa (red) corresponds to the 8FU didecamer
• The peak at 12.67 MDa (red) corresponds to the 8FU 3-decamer
• For the fourth oligomer (~16.0 MDa), I observed only a small blue “hump” in the region between 16–17 MDa, without any red fitted peak.

This can be explained by the fact that:

• The signal for this oligomer is much weaker compared to the others
• There may be fewer particles detected at this mass
• The signal may be too noisy or not well-defined

Because of this, the software was not able to confidently fit a Gaussian curve, and therefore no red peak was generated. Despite this, the presence of the blue signal at the expected mass range (~16 MDa) still indicates the existence of the 8FU 4-decamer, even if it is less abundant or less stable.

For this homework, I used AI tools such as ChatGPT and DeepSeek to help structure my ideas and improve the clarity of my writing. I also used NotebookLM to better understand the provided resources and supporting materials. For the final project measurements, DeepSeek suggested including the last four key measurements, which I integrated into my analysis.

Sources:

• EC:1.2.2.4—FACTA Search. (n.d.). Retrieved April 14, 2026, from https://www.nactem.ac.uk/facta/cgi-bin/facta3.cgi?query=EC%3A1.2.2.4%7C111111%7C0%7C0%7C233944%7C0%7C10
• Herminghaus, S., Schreier, P. H., McCarthy, J. E. G., Landsmann, J., Botterman, J., & Berlin, J. (1991). Expression of a bacterial lysine decarboxylase gene and transport of the protein into chloroplasts of transgenic tobacco. Plant Molecular Biology. https://agris.fao.org/search/en/records/65de1eb24c5aef494fd9fee5
• Kim, Y. M., & Hegeman, G. D. (1981a). Purification and some properties of carbon monoxide dehydrogenase from Pseudomonas carboxydohydrogena. Journal of Bacteriology, 148(3), 904–911. https://doi.org/10.1128/jb.148.3.904-911.1981
• Kim, Y. M., & Hegeman, G. D. (1981b). Purification and some properties of carbon monoxide dehydrogenase from Pseudomonas carboxydohydrogena. Journal of Bacteriology, 148(3), 904–911. https://doi.org/10.1128/jb.148.3.904-911.1981
• Matzke, M. A., & Matzke, A. J. (1998). Epigenetic silencing of plant transgenes as a consequence of diverse cellular defence responses. Cellular and Molecular Life Sciences: CMLS, 54(1), 94–103. https://doi.org/10.1007/s000180050128
• Maxwell, K., & Johnson, G. N. (2000). Chlorophyll fluorescence—A practical guide. Journal of Experimental Botany, 51(345), 659–668. https://doi.org/10.1093/jxb/51.345.659
• Pahlow, S., Ostendorp, A., Kr&#252, L., &#223, el, & Kehr, J. (2018). Phloem Sap Sampling from Brassica napus for 3D-PAGE of Protein and Ribonucleoprotein Complexes. JoVE (Journal of Visualized Experiments), (131), e57097. https://doi.org/10.3791/57097
• Recent Advances and Emerging Trends in Chlorophyll Fluorescence Parameter Fv/Fm. (2025). Phyton-International Journal of Experimental Botany, 94(9), 2615–2630. https://doi.org/10.32604/phyton.2025.069246
• Remelli, W., Villafiorita, F., Casazza, A. P., & Santabarbara, S. (2018). Comparative excitation-emission dependence of the FV/FM ratio in model green algae and cyanobacterial strains. https://iris.cnr.it/handle/20.500.14243/365902
• Schägger, H., & von Jagow, G. (1991). Blue native electrophoresis for isolation of membrane protein complexes in enzymatically active form. Analytical Biochemistry, 199(2), 223–231. https://doi.org/10.1016/0003-2697(91)90094-a
• Smith, P. K., Krohn, R. I., Hermanson, G. T., Mallia, A. K., Gartner, F. H., Provenzano, M. D., Fujimoto, E. K., Goeke, N. M., Olson, B. J., & Klenk, D. C. (1985). Measurement of protein using bicinchoninic acid. Analytical Biochemistry, 150(1), 76–85. https://doi.org/10.1016/0003-2697(85)90442-7
• Woo, J.-K., Hong, C. B., & Lee, J.-S. (1991). Chloroplast Targeting of Bacterial β-Glucuronidase with a Pea Transit Peptide in Transgenic Tobacco Plants. Molecules and Cells, 1(4), 451–457. https://doi.org/10.1016/S1016-8478(23)13893-3

Week 11 HW: Bioproduction & Cloud Labs

Contribution to the Collective Bioart Project

I contributed to several designs during the experiment. My final contribution was trying to create a geometric pattern inspired by Islamic geometric art in the bottom-right corner of the pixel canvas. The design did not stay until the end because other participants kept modifying it, but it was interesting to see how the artwork kept changing with everyone’s input. What I Liked About the Project

I really liked the collaborative aspect of this project. It was fun to work with others at the same time, contribute to different designs, and watch them change in real time. The canvas was dynamic and creative, and it encouraged experimentation and shared participation.

Suggestions for Improvement

One improvement could be to limit each participant to only one or two pixels. This would encourage more collaboration, because people would need to work together to create designs instead of working alone on bigger parts. It could make the final artwork more coordinated and truly collaborative.

1. Component Roles (20-Hour NMP–Ribose–Glucose System)

image ref

• E. coli Lysate

BL21 (DE3) Star Lysate: Provides the core cellular machinery required for gene expression, including ribosomes, tRNAs, aminoacyl-tRNA synthetases, and metabolic enzymes. The BL21 (DE3) strain also supplies T7 RNA polymerase, enabling strong transcription from T7 promoters.

• Salts / Buffer

Potassium Glutamate: Maintains proper ionic strength and mimics the natural intracellular environment, helping stabilize proteins and support enzyme activity.

HEPES-KOH (pH 7.5): Acts as a buffering agent to keep the pH stable, which is essential for maintaining enzyme function during long incubations.

Magnesium Glutamate: Provides Mg²⁺ ions, which are essential cofactors for ribosome stability, RNA polymerase activity, and interactions with nucleic acids.

Potassium Phosphate (monobasic/dibasic): Serves as a secondary buffer and provides inorganic phosphate needed for ATP regeneration and nucleotide metabolism.

• Energy / Nucleotide System

Ribose: Feeds into the pentose phosphate pathway to generate precursors (like PRPP) required for nucleotide synthesis.

Glucose: Acts as the main energy source, supporting ATP production through metabolic pathways such as glycolysis.

AMP, CMP, UMP: These nucleoside monophosphates (NMPs) are low-cost precursors that are enzymatically converted into NTPs (ATP, CTP, UTP) for RNA synthesis.

Guanine: Supplied as a nucleobase that is converted into GMP through salvage pathways, then further phosphorylated into GTP for transcription.

• Translation Mix (Amino Acids)

17 Amino Acid Mix + Tyrosine + Cysteine: Provide all amino acids required for protein synthesis. Tyrosine and cysteine are added separately because they are less stable or less soluble in standard mixtures.

• Additives / Backfill

Nicotinamide: Acts as a precursor for NAD⁺, an important cofactor in metabolic reactions that support long-term energy regeneration.

Nuclease-Free Water: Serves as the solvent to adjust final concentrations while preventing degradation of DNA or RNA by nucleases.

2. Differences Between 1-Hour PEP-NTP and 20-Hour NMP–Ribose–Glucose Systems

The 1-hour PEP-NTP system is designed for rapid protein production by providing ready-to-use NTPs and a high-energy phosphate donor (PEP), allowing fast transcription and translation but for a short duration due to quick depletion of resources. In contrast, the 20-hour NMP–Ribose–Glucose system uses cheaper precursors (NMPs, ribose, glucose) and relies on the lysate’s metabolic pathways to gradually regenerate NTPs and energy, enabling longer and more sustained protein production.

3. Bonus Question

Transcription can still occur without GMP because the system includes guanine, which is converted into GMP through the salvage pathway. In this process, guanine is combined with PRPP (derived from ribose metabolism) to form GMP, which is then phosphorylated into GDP and GTP. The produced GTP is then used by T7 RNA polymerase for RNA synthesis.

1. Fluorescent Protein Properties Affecting Cell-Free Expression

The biophysical properties of fluorescent proteins (FPs), including folding efficiency, maturation time, pH sensitivity, oxygen dependence, and structural stability, play a critical role in determining their fluorescence output in cell-free systems, especially during extended incubations such as 36 hours.

sfGFP (Superfolder GFP)

This protein exhibits very fast folding and high structural stability, with efficient chromophore maturation that is oxygen-dependent. Its resistance to misfolding and aggregation allows it to maintain strong and consistent fluorescence over long incubation periods, making it a reliable reference protein in cell-free systems.

mRFP1 (Monomeric Red Fluorescent Protein 1)

It is characterized by slow maturation kinetics and incomplete chromophore formation, which delays the appearance of fluorescence. Additionally, it may form non-fluorescent intermediates, leading to lower overall signal intensity compared to more advanced red fluorescent proteins.

mKO2 (Monomeric Kusabira Orange 2)

This protein shows relatively fast maturation and high brightness, but its chromophore formation is strongly dependent on oxygen availability and can also be influenced by temperature. In conditions with limited oxygen or suboptimal temperature, its fluorescence intensity may be reduced.

mTurquoise2 It has a complex maturation mechanism and high quantum yield but is sensitive to environmental conditions such as pH and oxygen levels. Acidic conditions can reduce fluorescence, while insufficient oxygen may limit proper chromophore formation.

mScarlet-I

This protein is known for its very high brightness due to an excellent extinction coefficient and quantum yield. However, its performance depends on proper folding, and it can be sensitive to temperature or conditions that promote misfolding, which may reduce fluorescence output.

Electra2

Electra2 is engineered for rapid maturation and improved performance under reducing conditions commonly found in cell-free systems. Its stability in such environments allows it to maintain fluorescence where other proteins may struggle, although its long-term stability or photostability may vary depending on conditions.

2. Hypothesis (Electra2 Optimization)

For Electra2, fluorescence output over a 36-hour incubation may be limited by the availability of nucleotides and the sustainability of transcription in the cell-free system. I hypothesize that increasing the concentrations of ribose and nucleoside monophosphates (AMP, CMP, UMP, and guanine) will enhance the regeneration of nucleoside triphosphates (NTPs) through the lysate’s metabolic pathways. Ribose can be converted into phosphoribosyl pyrophosphate (PRPP), which is required for nucleotide synthesis, while NMPs and guanine serve as precursors that are enzymatically converted into NTPs. By increasing these components, the system should maintain a continuous supply of NTPs, thereby sustaining transcription by T7 RNA polymerase and increasing mRNA production over time. As a result, this enhanced transcriptional activity is expected to support prolonged translation and lead to higher cumulative protein production and fluorescence intensity over the 36-hour period. This strategy is particularly suitable for Electra2, which is designed for rapid maturation and can efficiently convert increased protein synthesis into measurable fluorescence.

image ref

3. Master Mix Design: The Three-Well Strategy

To test this hypothesis, I designed three distinct reagent compositions to identify the “sweet spot” between fuel availability and metabolic stability.

• Mix 1: The “Maximized” Fuel Mix (Well Q4-H20)

Goal: To test the absolute capacity of the system by pushing precursors to the high end.

Key Adjustments: Ribose was increased to 19.0 g/L (+63.4 %) and NMPs (AMP/CMP/UMP) were increased by 60-100 %. Guanine was doubled to 0.313 mM to provide a surplus of base molecules for the salvage pathway.

• Mix 2: The “Intermediate” Mix (Well Q4-G21)

Goal: To establish a bridge between the standard mix and the maximum boost.

Key Adjustments: Ribose was set at 15.0 g/L (+29%) and NMPs/Guanine were increased by 20-33 %. This well helps determine if the “Max” mix is overkill or if a moderate increase is sufficient.

• Mix 3: The “Direct Supply” Mix (Well Q4-I21)

Goal: To test if bypassing the enzymatic salvage of Guanine improves initial speed.

Key Adjustments: While maintaining the Intermediate fuel levels, I added 0.500 mM of pure GMP. This tests whether providing a direct nucleotide (GMP) is more efficient for Electra2 than relying solely on Guanine-to-GMP conversion.

Final Concentration Comparison Table

4. Data Analysis Strategy

Once the 36-hour fluorescence data is returned, I will compare the slopes and peak intensities of these three wells.

• Validation: If Mix 1 > Mix 2 > Standard, the limiting factor was raw fuel.
• Metabolic Insights: If Mix 3 reaches a plateau faster than Mix 2, it proves the enzymatic conversion of Guanine was a kinetic bottleneck for Electra2 production.

Sources:

• Banks, A. M., Whitfield, C. J., Brown, S. R., Fulton, D. A., Goodchild, S. A., Grant, C., Love, J., Lendrem, D. W., Fieldsend, J. E., & Howard, T. P. (2022). Key reaction components affect the kinetics and performance robustness of cell-free protein synthesis reactions. Computational and Structural Biotechnology Journal, 20, 218–229. https://doi.org/10.1016/j.csbj.2021.12.013
• Burrington, L. R., Watts, K. R., & Oza, J. P. (2021). Characterizing and Improving Reaction Times for E. coli-Based Cell-Free Protein Synthesis. ACS Synthetic Biology, 10(8), 1821–1829. https://doi.org/10.1021/acssynbio.1c00195
• Deng, H., Callender, R., Schramm, V. L., & Grubmeyer, C. (2010). Pyrophosphate Activation in Hypoxanthine-Guanine Phosphoribosyltransferase with Transition State Analogue. Biochemistry, 49(12), 2705–2714. https://doi.org/10.1021/bi100012u
• Dopp, B. J. L., Tamiev, D. D., & Reuel, N. F. (2019). Cell-free supplement mixtures: Elucidating the history and biochemical utility of additives used to support in vitro protein synthesis in E. coli extract. Biotechnology Advances, 37(1), 246–258. https://doi.org/10.1016/j.biotechadv.2018.12.006
• Dudzinska, W., Lubkowska, A., Dolegowska, B., Safranow, K., & Jakubowska, K. (2010). Adenine, guanine and pyridine nucleotides in blood during physical exercise and restitution in healthy subjects. European Journal of Applied Physiology, 110(6), 1155–1162. https://doi.org/10.1007/s00421-010-1611-7
• Gregorio, N. E., Levine, M. Z., & Oza, J. P. (2019). A User’s Guide to Cell-Free Protein Synthesis. Methods and Protocols, 2(1), 24. https://doi.org/10.3390/mps2010024
• Hashimura, H., Nakagawa, H., & Sawai, S. (2025). Use of blue fluorescent protein Electra2 for live-cell imaging in Dictyostelium discoideum. microPublication Biology. https://doi.org/10.17912/micropub.biology.001774
• Hove-Jensen, B., Andersen, K. R., Kilstrup, M., Martinussen, J., Switzer, R. L., & Willemoës, M. (2016). Phosphoribosyl Diphosphate (PRPP): Biosynthesis, Enzymology, Utilization, and Metabolic Significance. Microbiology and Molecular Biology Reviews : MMBR, 81(1), e00040-16. https://doi.org/10.1128/MMBR.00040-16
• Jiang, L., Zhao, J., Lian, J., & Xu, Z. (2018). Cell-free protein synthesis enabled rapid prototyping for metabolic engineering and synthetic biology. Synthetic and Systems Biotechnology, 3(2), 90–96. https://doi.org/10.1016/j.synbio.2018.02.003
• Jiang, N., Ding, X., & Lu, Y. (2021). Development of a robust Escherichia coli-based cell-free protein synthesis application platform. Biochemical Engineering Journal, 165, 107830. https://doi.org/10.1016/j.bej.2020.107830
• Krinsky, N., Kaduri, M., Shainsky-Roitman, J., Goldfeder, M., Ivanir, E., Benhar, I., Shoham, Y., & Schroeder, A. (2016). A Simple and Rapid Method for Preparing a Cell-Free Bacterial Lysate for Protein Synthesis. PLOS ONE, 11(10), e0165137. https://doi.org/10.1371/journal.pone.0165137
• Vengut-Climent, E., Peñalver, P., Lucas, R., Gómez-Pinto, I., Aviñó, A., Muro-Pastor, A. M., Galbis, E., de Paz, M. V., Fonseca Guerra, C., Bickelhaupt, F. M., Eritja, R., González, C., & Morales, J. C. (2018). Glucose-nucleobase pairs within DNA: Impact of hydrophobicity, alternative linking unit and DNA polymerase nucleotide insertion studies †Electronic supplementary information (ESI) available. See DOI: 10.1039/c7sc04850e. Chemical Science, 9(14), 3544–3554. https://doi.org/10.1039/c7sc04850e
• Zhang, Y., Huang, Q., Deng, Z., Xu, Y., & Liu, T. (2018). Enhancing the efficiency of cell-free protein synthesis system by systematic titration of transcription and translation components. Biochemical Engineering Journal, 138, 47–53. https://doi.org/10.1016/j.bej.2018.07.001

Week 1 Lab: Pipetting

Individual Final Project

Project Title: Engineering Houseplants for Atmospheric Carbon Monoxide Capture: Chloroplast-Targeted Expression of the Bacterial CODH Enzyme Complex in Nicotiana tabacum

The Problem This Project Addresses

Carbon monoxide (CO) is a colorless, odorless, tasteless toxic gas that cannot be detected by human senses. It is produced whenever something burns incompletely — gas heaters, stoves, car engines, fireplaces, and wood-burning appliances all release CO. Indoors, CO accumulates silently and can reach dangerous or fatal concentrations before anyone notices. The current standard of protection is a battery-powered electrochemical CO detector. These devices are excellent at detecting CO and sounding an alarm , but they cannot remove the gas from the air. Once the alarm sounds, the occupants must evacuate and ventilate the space manually. Furthermore, CO detectors require regular battery replacement and eventually need to be replaced entirely. In low-income households worldwide, detectors are frequently absent, have dead batteries, or are past their useful lifespan.

–> This project proposes a fundamentally different approach: instead of detecting CO, make the plant remove it.

The Core Idea

Certain bacteria ,particularly Oligotropha carboxidovorans, have evolved the ability to use CO as a food source. They do this using an enzyme called Carbon Monoxide Dehydrogenase (CODH), which converts CO into CO₂ according to this reaction:

CO + H₂O → CO₂ + 2 electrons + 2 protons

The CO₂ produced by this reaction is not harmful at the quantities involved and supposed to be reused by a plant’s own photosynthesis through the Calvin cycle.

This project proposes to take the bacterial CODH system out of the bacterium and introduce it into a plant, specifically targeting it to the chloroplast (the organelle where photosynthesis happens). By placing CODH inside the chloroplast, two elegant outcomes occur simultaneously:

1. The plant actively breaks down CO from the surrounding air
2. The CO₂ produced by CODH is immediately captured by Rubisco and enters the Calvin cycle, making the plant slightly more productive

The scientific foundation for this idea is already established in the literature. Duffus et al. (2018) demonstrated that the complete CODH complex can be functionally expressed in Escherichia coli –> proving heterologous expression is achievable. South et al. (2019) demonstrated in Science that bacterial enzymes introduced into tobacco chloroplasts producing CO₂ directly in the stroma increased plant biomass by up to 40% –> proving that chloroplast-produced CO₂ is efficiently captured by photosynthesis. This project extends this logic to a new substrate: atmospheric CO.

The Complete Genetic System Required

The CODH enzyme from O. carboxidovorans is not a single protein. It is a complex system requiring seven genes organized into two functional groups:

1. Group 1 — Structural subunits (the enzyme itself):

coxL –> the large catalytic subunit (~88 kDa) where CO is actually oxidized. Contains the unique [CuSMoO₂] active site coxM –> the medium subunit (~30 kDa) containing FAD, responsible for electron transfer coxS –> the small subunit (~18 kDa) containing [2Fe-2S] iron-sulfur clusters, part of the electron relay chain

These three proteins assemble into a (CoxL·CoxM·CoxS)₂ heterohexamer — a complex of six protein subunits working together.

2. Group 2 — Maturation proteins (the assembly machinery):

coxD –> an AAA+ ATPase chaperone that acts as a “maturation protein,” responsible for the post-translational insertion of copper and the essential bridging sulfur into the apo-enzyme, converting it to active holo-enzyme. coxE, coxF and coxG –> “final processing” and “sulfur addition” are part of a complex pathway. According to research, coxF plays a role in copper acquisition/mobilization, and coxE and coxG are involved in the maturation pathway that leads to the properly sulfurated and copper-inserted active site. The exact individual functions of coxE and coxG are still being elucidated, though their role in the maturation complex is essential.

Overview of the Three Aims

AIM 1 — Computational Design and Validation of the Complete Genetic System

In simple terms: Design the complete genetic blueprint for the CO-capturing plant system on a computer, verify every element computationally, and produce a synthesis-ready design.

The seven bacterial genes cannot simply be pasted into a plant. They need to be comprehensively redesigned for plant expression:

• Their DNA sequences must be rewritten in “plant language” through codon optimization
• Each protein needs a molecular address label (chloroplast transit peptide) added to its beginning so it is directed to the correct location inside the plant cell
• The address labels must be verified to ensure the plant’s processing machinery will correctly remove them after the protein arrives
• Each gene needs its own promoter (an on-switch for gene expression) and terminator (an off-switch), carefully chosen to prevent the plant from silencing all the genes simultaneously
• Translation enhancer sequences must be added to maximize protein production
• Spacer sequences must be placed between genes to prevent one gene’s transcription from accidentally running into the next
• The complete system must be distributed across two separate transformation vectors

All of this is done computationally using Benchling, A codon optimization tool, ChloroP 1.1, Boltz, and the Asimov Kernel –> producing a complete verified design ready for DNA synthesis through Twist Biosciences.

AIM 2 — Wet Lab Transformation and Functional Validation (The next step — beyond this course)

In simple terms: Actually build the constructs in the lab, put them into tobacco plants, and prove the enzyme works. Aim 2 begins where Aim 1 ends. The Twist-synthesized multicassettes fragments are assembled into the pCAMBIA vectors using Gibson Assembly. The constructs are introduced into Nicotiana tabacum via Agrobacterium tumefaciens-mediated leaf disc transformation , the standard method for introducing genes into tobacco. Transgenic plants are selected on dual antibiotic medium (hygromycin + kanamycin, confirming both constructs integrated).

The experimental progression follows strict logic — each step must succeed before the next begins:

• –> Step 1 — Chloroplast targeting validation
• –> Step 2 — Gene integration and transcription
• –> Step 3 — Protein expression and CTP cleavage
• –> Step 4 — Complex assembly
• –> Step 5 — CO oxidation activity
• –> Step 6 — Plant health and photosynthesis

for more details, please take a look on part I of week 10 homework.

AIM 3 — Optimization, Transfer to Houseplants, and Real-World Deployment(The long-term vision)

In simple terms: Assuming Aim 2 succeeds, optimize the system, transfer it to real houseplants, and develop it toward real-world deployment. If Aim 2 demonstrates functional CO oxidation in tobacco, Aim 3 pursues three parallel directions:

Direction 1 — Transfer to real houseplants: The validated genetic architecture from tobacco is adapted for transformation into Epipremnum aureum (Pothos) and Spathiphyllum wallisii (Peace Lily) — widely kept, hardy, aesthetically acceptable houseplants. Agrobacterium-mediated transformation protocols established for tobacco are adapted for these species.

Direction 2 — System optimization: Several improvements are pursued to increase CO removal efficiency and operational range:

A CO-responsive inducible promoter system replaces constitutive promoters, activating CODH expression only when CO is present and saving plant energy otherwise Constitutively open stomata engineering to maintain CO uptake during nighttime hours when CO poisoning risk is highest Expression levels are optimized based on the quantitative CO removal model to increase per-plant removal capacity

Direction 3 — Safety, containment, and deployment:

Genetic Use Restriction Technology (GURT): To prevent seed viability and uncontrolled environmental spread, I will implement Genetic Use Restriction Technology (GURT). This ensures that any engineered plants cannot reproduce outside controlled environments. Additional containment strategy — chloroplast genome integration:

As an alternative or complement to GURT, I can integrate the transgenes into the chloroplast genome instead of the nuclear genome. Chloroplast DNA is maternally inherited in most flowering plants, including tobacco (Nicotiana tabacum). This means the transgenes are not transmitted via pollen, virtually eliminating the risk of gene flow to wild relatives. This is a well-established biosafety strategy for plant synthetic biology.

Regulatory pathway planning begins under USDA APHIS (Regulation of genetically engineered plantsand) EPA (Regulation of plants producing pesticidal substances (if applicable))frameworks.

The deployment target is refined based on the quantitative CO removal analysis: rather than acute emergency protection in homes (which requires too many plants), the primary application is chronic CO reduction in high-exposure industrial and semi-industrial environments like workshops, garages, underground parking facilities, and developing-world indoor cooking spaces where CO concentrations are higher and more sustained.

The ethical framework for commercial deployment ,including informed consent, false assurance prevention, equity of access, and environmental risk, is fully developed and integrated into regulatory submissions.

Sources:

• Bährle, R., Böhnke, S., Englhard, J., Bachmann, J., & Perner, M. (2023). Current status of carbon monoxide dehydrogenases (CODH) and their potential for electrochemical applications. Bioresources and Bioprocessing, 10(1), 84. https://doi.org/10.1186/s40643-023-00705-9
• Dent, M. R., Weaver, B. R., Roberts, M. G., & Burstyn, J. N. (2023). Carbon Monoxide-Sensing Transcription Factors: Regulators of Microbial Carbon Monoxide Oxidation Pathway Gene Expression. Journal of Bacteriology, 205(5), e00332-22. https://doi.org/10.1128/jb.00332-22
• Erb, T. J. (2024). Photosynthesis 2.0: Realizing New-to-Nature CO2-Fixation to Overcome the Limits of Natural Metabolism. Cold Spring Harbor Perspectives in Biology, 16(2), a041669. https://doi.org/10.1101/cshperspect.a041669
• Kaufmann, P., Duffus, B. R., Teutloff, C., & Leimkühler, S. (2018). Functional Studies on Oligotropha carboxidovorans Molybdenum–Copper CO Dehydrogenase Produced in Escherichia coli. Biochemistry, 57(19), 2889–2901. https://doi.org/10.1021/acs.biochem.8b00128
• Liu, C., Zhang, N., Sun, L., Gao, W., Zang, Q., & Wang, X. (2022). Potted plants and ventilation effectively remove pollutants from tobacco smoke. International Journal of Low-Carbon Technologies, 17, 1052–1060. https://doi.org/10.1093/ijlct/ctac081
• Park, S., Mani, V., Kim, J. A., Lee, S. I., & Lee, K. (2022). Combinatorial transient gene expression strategies to enhance terpenoid production in plants. Frontiers in Plant Science, 13, 1034893. https://doi.org/10.3389/fpls.2022.1034893
• Qin, S., Liu, Y., Yan, J., Lin, S., Zhang, W., & Wang, B. (2022). An Optimized Tobacco Hairy Root Induction System for Functional Analysis of Nicotine Biosynthesis-Related Genes. Agronomy, 12(2), 348. https://doi.org/10.3390/agronomy12020348
• Schübel, U., Kraut, M., Mörsdorf, G., & Meyer, O. (1995). Molecular characterization of the gene cluster coxMSL encoding the molybdenum-containing carbon monoxide dehydrogenase of Oligotropha carboxidovorans. Journal of Bacteriology, 177(8), 2197–2203. https://doi.org/10.1128/jb.177.8.2197-2203.1995
• Siebert, D., Busche, T., Metz, A. Y., Smaili, M., Queck, B. A. W., Kalinowski, J., & Eikmanns, B. J. (2020). Genetic Engineering of Oligotropha carboxidovorans Strain OM5—A Promising Candidate for the Aerobic Utilization of Synthesis Gas. ACS Synthetic Biology, 9(6), 1426–1440. https://doi.org/10.1021/acssynbio.0c00098
• Tao, Y., Chiu, L.-W., Hoyle, J. W., Dewhirst, R. A., Richey, C., Rasmussen, K., Du, J., Mellor, P., Kuiper, J., Tucker, D., Crites, A., Orr, G. A., Heckert, M. J., Godinez-Vidal, D., Orozco-Cardenas, M. L., & Hall, M. E. (2023). Enhanced Photosynthetic Efficiency for Increased Carbon Assimilation and Woody Biomass Production in Engineered Hybrid Poplar. Forests, 14(4), 827. https://doi.org/10.3390/f14040827
• Thagun, C., Odahara, M., Kodama, Y., & Numata, K. (2024). Identification of a highly efficient chloroplast-targeting peptide for plastid engineering. PLOS Biology, 22(9), e3002785. https://doi.org/10.1371/journal.pbio.3002785

Group Final Project