Individual Final Project

cover image cover image

Subsections of Individual Final Project

Week 1 HW: Principles and Practices

Zambia Mineral-Waste Bioremediation Predictor

From Metagenome to Marketable Bioremediation Product

HTGAA 2026 Final Project · Elsa Muleya · SynBio USFQ Node


Project Rationale

Zambia’s Copperbelt Province faces severe heavy metal contamination from decades of copper mining at Konkola, Nchanga, Mufulira, and Chingola. Cu²⁺, Zn²⁺, Co²⁺, and Pb²⁺ leach from mine tailings into groundwater and agricultural soils at concentrations far exceeding WHO limits, with no affordable or accessible remediation solution for affected communities.

This project designs, validates, and packages a living biological solution: engineered Bacillus subtilis carrying a novel metallothionein (MT) gene discovered from Zambian mine-associated bacterial genomes, encapsulated in a field-deployable dual-layer hydrogel biocontainment system — ZAMGEL — that can be commercially produced and applied without specialist equipment or laboratory infrastructure.


Three-Aim Project Structure

AimTitleFocus
1Bioinformatics Discovery & Genetic DesignMetagenomics, structural prediction, circuit design
2Wet Lab Validation Under Zambian ConditionsTransformation, metal assays, pH & stress testing
3ZAMGEL Containment & Commercial Product DesignHydrogel bioencapsulation, kill-switch, market pathway

Aim 1: Bioinformatics Discovery & Genetic Design

Goal: Identify and structurally validate novel metallothioneins from Zambian mine-associated bacterial genomes, and design a complete synthetic expression cassette ready for wet lab transformation.

Sub-aim 1a: Metagenomic Mining of Zambian Copperbelt Sequences

Mine publicly available sequencing datasets from NCBI SRA, MG-RAST, and IMG/M targeting the Konkola, Nchanga, and Mufulira mine regions. The full computational pipeline:

FASTQ → fastp (QC trim) → MEGAHIT (assembly) → Prodigal (ORF prediction) → BLASTp + Prokka (annotation)

Filter candidates by the presence of the Cys-X-Cys motif — the canonical Cu/Zn coordination fingerprint in prokaryotic metallothioneins — and cross-reference against known prokaryotic MT families (SmtA-like, BmtA-like, CzcA operons, CopA ATPases). Build a maximum-likelihood phylogenetic tree using IQ-TREE 2 to confirm novelty.

DatabasePurpose
NCBI SRAPrimary source for Zambian mine metagenome FASTQ files
MG-RASTMine microbiome metagenomes with functional annotation
IMG/MIntegrated Microbial Genomes — metal resistance gene clusters
UniProt/SwissProtReference MT homology and Cys-X-Cys motif validation

Sub-aim 1b: Structural Validation & Synthetic Expression Cassette Design

For the top 5 MT candidates from Sub-aim 1a, simultaneously validate 3D structural integrity and design the full synthetic genetic system.

Structural Validation

  • Submit top candidate sequences to AlphaFold3 to generate .pdb files and visualise cysteine-rich metal-binding pockets
  • Pass threshold: pLDDT > 85 across the metal-binding domain; ipTM > 0.80 for confident fold prediction
  • Quantify binding pocket geometry in PyMOL / ChimeraX: pocket volume (ų), solvent accessibility, Cys coordination angle, and closest Cys–Cys distance (target < 6 Å for effective Cu²⁺ coordination)
  • Calculate predicted dissociation constant: Kd = e^(ΔG/RT) at T = 310 K (37°C); expected range 10⁻¹³ to 10⁻¹⁵ M for high-performance prokaryotic MTs
  • Compare all candidates against reference proteins (SmtA from Synechococcus PCC 7942; BmtA from Pseudomonas) on Kd, Cys count, and pLDDT

Expression Cassette Design (Benchling)

  • Codon-optimise the best-scoring MT sequence for B. subtilis 168 using Benchling’s built-in optimiser
  • Design a metal-responsive synthetic circuit in Cello 2.0: Cu²⁺ sensor (PcorA or PmtA promoter) → NOT gate logic → MT expressed only when Cu²⁺ exceeds threshold
  • Include eGFP fluorescent reporter downstream of MT as a real-time visual proxy for circuit activation
5'─[PcopA/PmtA]─[RBS B0034]─[MT_Bsubtilis_optimised]─[eGFP]─[T_B0015]─3'
    Cu²⁺ sensor   strong RBS    codon-optimised         reporter  terminator
  • Verify BioBrick RFC10 compatibility in Benchling
  • Submit all sequences through Twist Bioscience biosecurity screening (“Green” classification required before synthesis order)

Aim 2: Wet Lab Validation Under Zambian Environmental Conditions

Goal: Transform the computationally designed system into a living, functional biosensor-remediator and rigorously stress-test it against the real environmental conditions of the Zambian Copperbelt.

Sub-aim 2a: Chassis Construction & Verification

Transform B. subtilis 168 with the assembled MT expression plasmid and confirm successful integration using three independent assays before proceeding to metal exposure experiments:

AssayMethodPass Criterion
Colony PCRMT-specific primers flanking insert; 30 cycles, 55°C annealingBand at expected insert size
Sanger SequencingSequence full insert with M13 forward/reverse primers100% identity to designed cassette
SDS-PAGE + Western BlotAnti-His-tag antibody; 4h induction at 37°CBand at ~6 kDa (49 AA protein)
GFP Fluorescence MicroscopyImage colonies in Cu²⁺-spiked media at Ex 488 / Em 510 nm> 5× fluorescence over water control

Sub-aim 2b: Metal Ion Concentration Response Assays

Expose the engineered B. subtilis to a full Cu²⁺ concentration gradient spanning real Copperbelt mine drainage (reported range: 0.5–500 mg/L). Measure metal removal using ICP-MS on growth media supernatant and calculate Bio-Sequestration Efficiency (%BSE):

%BSE = ([Metal]₀ − [Metal]f) ÷ [Metal]₀ × 100
Cu²⁺ ConcentrationEnvironmental ContextMeasurements
0 mg/LNegative controlGFP baseline, OD600, ICP-MS
0.5 mg/LWHO drinking water limitGFP, OD600, ICP-MS
5 mg/LWHO industrial discharge limitGFP, OD600, ICP-MS
50 mg/LTypical Konkola drainage concentrationGFP, OD600, ICP-MS
500 mg/LPeak Copperbelt leachate concentrationGFP, OD600, ICP-MS, survival rate
1000 mg/LToxicity threshold — LD50 determinationColony viability, LD50 endpoint

Sub-aim 2c: pH Stress Testing

Zambian mine tailings range from pH 2.5–4.5 (active acid mine drainage) to pH 8–9 (alkaline neutralisation runoff). Test bacteria across this full range at fixed 50 mg/L Cu²⁺ to define the operational pH window and inform ZAMGEL outer shell buffer design.

pHEnvironmental Context (Zambia)Measurements
2.5Active acid mine drainage leachateGFP, OD600, ICP-MS
3.5Tailing pond runoffGFP, OD600, ICP-MS
4.5Near-tailing agricultural soil leachateGFP, OD600, ICP-MS
5.5Mildly acidic Copperbelt soilGFP, OD600, ICP-MS
6.5 ★Neutral control (laboratory standard)GFP, OD600, ICP-MS
7.5Borehole drinking water (Kitwe)GFP, OD600, ICP-MS
8.5Alkaline mine neutralisation runoffGFP, OD600, ICP-MS
9.0Extreme alkaline drainage (worst case)GFP, OD600, ICP-MS

Sub-aim 2d: Multi-Stressor Environmental Simulation

Real Copperbelt soil presents multiple co-occurring stresses. Bacteria must survive all of these simultaneously to be field-deployable. Each stressor is tested at fixed Cu²⁺ = 50 mg/L and pH 6.5 to isolate the effect; a final cocktail experiment combines all worst-case stressors simultaneously.

StressorZambia-Specific ConditionTest ParametersOutput Measured
TemperatureAvg 24°C; dry season peak 38°C20, 28, 37, 42°COD600, GFP, %BSE
Co-metal toxicityCu²⁺ + Zn²⁺ + Co²⁺ + Pb²⁺ co-contaminationSingle vs cocktail, 50 mg/L eachICP-MS all ions, GFP
DesiccationDry season soil water activity < 0.85aw 0.85, 0.90, 0.95 via NaClOD600, colony viability
UV exposureHigh solar UV at 12–15°S latitudeUV-C 254 nm: 0, 10, 30, 60 s pulseColony survival, DNA damage gel
Competing microbiomeIndigenous Copperbelt soil microbiome10% v/v heat-killed soil extractGFP, OD600, ICP-MS

Aim 3: ZAMGEL Containment System & Commercial Product Design

Goal: Design a biomaterial containment system that physically and genetically contains the engineered bacteria inside a field-deployable carrier, preventing environmental escape while maintaining full metal-sequestration function — creating a product that can be commercially sold and applied without ecological risk.

Sub-aim 3a: ZAMGEL Dual-Layer Hydrogel Bioencapsulation

The ZAMGEL biocapsule is a three-layer biomaterial architecture. Each layer performs a distinct function, together creating a self-contained living bioreactor deployable directly onto mine tailings:

LayerCompositionFunctionSourcing
Outer shellCalcium alginate + CaCO₃ nanoparticlespH buffering: neutralises acidic mine leachate to pH 5.5–6.5 before bacteria are exposed; structural integrity in soilFood-grade alginate; CaCO₃ from local limestone
Middle membraneCellulose nanofibre + chitosan crosslinkSize-selective filter: 200 nm pores allow Cu²⁺ ions (0.73 Å) to enter freely; bacteria (1–2 µm) physically cannot escapeLocal agricultural waste cellulose; chitosan import
Inner corePVA + gelatin hydrogel + activated charcoalBacteria viability matrix at 10⁸ CFU/mL; activated charcoal provides passive metal co-adsorption during biological lag phaseCommercial PVA/gelatin; charcoal from local Copperbelt source

Sub-aim 3b: Containment Validation & Kill-Switch Integration

Containment Validation

TestProtocolPass Threshold
Bacterial escapePlate surrounding water on LB agar at 7, 14, 30 days< 1 CFU/mL at 30 days
Ion permeabilityICP-MS of surrounding fluid vs bead interior after 24h Cu²⁺ exposureCu²⁺ enters freely; bacteria absent in external fluid
Mechanical durabilityCompression to 50 kPa (equivalent to 30 cm soil overburden)No structural failure; containment maintained
Biodegradation rateBury spent beads in Zambian soil analogue at 28°C; measure mass loss weeklyFull degradation in 90–180 days; no persistent residue

Genetic Kill-Switch (MazF/MazE Toxin-Antitoxin)

A MazF/MazE kill-switch is integrated into the B. subtilis chromosome (not plasmid, to prevent loss). MazE antitoxin is expressed under a Ptet promoter requiring anhydrotetracycline (aTc) to remain active. When aTc is withdrawn (ZAMGEL retrieved or degraded at end of life), MazE degrades, MazF mRNA interferase cleaves all mRNA, and all bacteria die within 48 hours. A secondary CcdB/CcdA kill-switch on the plasmid backbone provides an orthogonal safety layer.

aTc present → MazE expressed → MazF neutralised → Bacteria LIVE
aTc absent  → MazE degraded  → MazF active      → Bacteria DEAD within 48h

Sub-aim 3c: Commercial Product Formats & Digital Predictor App

FormatDescriptionUse CaseDeployment
ZAMGEL Beads3–5 mm spheres, ~10⁸ CFU/beadMine water treatment pondsBroadcast by hand or machine
ZAMGEL Sheets10×10 cm biodegradable matsSoil surface tailing cap treatmentLay directly on contaminated soil
ZAMGEL CartridgesInline filter column packed with beadsBorehole and drainage pipe treatmentInstall in drainage infrastructure

A Streamlit-based mobile web app (offline-capable PWA) allows community members and mine site managers to input local soil Cu²⁺ concentration, pH, temperature, and treatment area, and receive a data-driven treatment recipe — number of ZAMGEL beads, predicted %BSE, and estimated remediation timeline — based on dose-response curves generated in Aim 2. No laboratory equipment required.

Regulatory pathway: Zambia Environmental Management Agency (ZEMA) contained-use application under Biosafety Act No. 10 of 2007; Nagoya Protocol compliance for use of indigenous Zambian microbial genetic resources; community consent framework with Copperbelt mining communities. Primary commercial client: ZCCM-IH.


15-Week Project Timeline

WeekAimActivity
11aSRA/MG-RAST/IMG/M search for Konkola, Nchanga, Mufulira mine datasets; quality trim with fastp
21aMEGAHIT assembly → Prodigal ORF prediction → BLASTp + Prokka annotation of metal resistance genes
31aCys-X-Cys motif filter → top 5 candidates selected; IQ-TREE 2 maximum-likelihood phylogenetic tree
41bAlphaFold3 structure prediction for all 5 candidates; retrieve .pdb files
51bPyMOL/ChimeraX binding pocket quantification: volume, Cys coordination geometry, pLDDT mapping
61bBenchling codon optimisation + Cello 2.0 logic gate design + Twist Bioscience DNA order
72aB. subtilis 168 transformation; colony PCR; Sanger sequencing verification
82aSDS-PAGE + western blot + GFP fluorescence microscopy to confirm MT expression
92bCu²⁺ concentration gradient assays (0–1000 mg/L); ICP-MS; GFP plate reader; dose-response curve
102cpH stress assays (pH 2.5–9.0) at 50 mg/L Cu²⁺; identify operational pH window
112dMulti-stressor factorial experiment: temperature × co-metals × UV × desiccation × microbiome cocktail
123aZAMGEL prototype fabrication: alginate outer shell + chitosan membrane + PVA/gelatin inner core
133bContainment validation: LB plating, ICP-MS permeability, compression testing, biodegradation assay
143bMazF/MazE kill-switch chromosomal integration + aTc withdrawal 48h death assay; CcdB/CcdA backup
153cStreamlit app prototype; ZEMA regulatory pathway draft; final in silico feasibility report

Validation Criteria & Contingency Plans

ExperimentPass ThresholdIf Fail — Contingency
AlphaFold3 pLDDT (binding domain)> 85 on core domain; ipTM > 0.80Use SmtA (Synechococcus PCC 7942) as positive control scaffold; re-run with AlphaFold2
GFP activation in Cu²⁺ media> 5× fluorescence over backgroundRedesign Cello promoter with stronger RBS; increase plasmid copy number
ICP-MS metal removal (%BSE)> 60% BSE at 50 mg/L Cu²⁺Increase MT copy number via multi-copy plasmid (pHT01); co-express CopA copper ATPase
pH operational windowActive sequestration at pH 4.5–8.0Increase CaCO₃ loading in ZAMGEL outer shell; add internal carbonate buffer inside PVA core
ZAMGEL containment (30 days)< 1 CFU/mL in surrounding mediumIncrease chitosan crosslink density; reduce pore size to 100 nm
Kill-switch efficacy100% cell death within 48h of aTc removalSwitch to CcdB/CcdA system; add second orthogonal kill-switch on separate chromosome locus

Why This Project Matters

Existing Copperbelt remediation approaches — lime neutralisation, chemical precipitation, pump-and-treat — are capital-intensive, infrastructure-dependent, and inaccessible to subsistence communities adjacent to mine tailings. The ZAMGEL system offers:

  • No electricity or specialist infrastructure required — scatter-and-forget deployment
  • Zero environmental release — physically contained by 200 nm membrane; genetically contained by dual kill-switch
  • Self-regulating — MT only expressed when Cu²⁺ exceeds threshold; GFP reporter confirms activity in real time
  • Locally grounded — MT gene discovered from Zambian mine-associated bacterial genomes
  • Commercially viable — manufacturable from locally sourced materials; approvable under existing Zambian biosafety law
  • Community-facing — Streamlit app enables treatment planning without laboratory equipment or expertise