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
| Aim | Title | Focus |
|---|---|---|
| 1 | Bioinformatics Discovery & Genetic Design | Metagenomics, structural prediction, circuit design |
| 2 | Wet Lab Validation Under Zambian Conditions | Transformation, metal assays, pH & stress testing |
| 3 | ZAMGEL Containment & Commercial Product Design | Hydrogel 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:
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.
| Database | Purpose |
|---|---|
| NCBI SRA | Primary source for Zambian mine metagenome FASTQ files |
| MG-RAST | Mine microbiome metagenomes with functional annotation |
| IMG/M | Integrated Microbial Genomes — metal resistance gene clusters |
| UniProt/SwissProt | Reference 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
.pdbfiles 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
- 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:
| Assay | Method | Pass Criterion |
|---|---|---|
| Colony PCR | MT-specific primers flanking insert; 30 cycles, 55°C annealing | Band at expected insert size |
| Sanger Sequencing | Sequence full insert with M13 forward/reverse primers | 100% identity to designed cassette |
| SDS-PAGE + Western Blot | Anti-His-tag antibody; 4h induction at 37°C | Band at ~6 kDa (49 AA protein) |
| GFP Fluorescence Microscopy | Image 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):
| Cu²⁺ Concentration | Environmental Context | Measurements |
|---|---|---|
| 0 mg/L | Negative control | GFP baseline, OD600, ICP-MS |
| 0.5 mg/L | WHO drinking water limit | GFP, OD600, ICP-MS |
| 5 mg/L | WHO industrial discharge limit | GFP, OD600, ICP-MS |
| 50 mg/L | Typical Konkola drainage concentration | GFP, OD600, ICP-MS |
| 500 mg/L | Peak Copperbelt leachate concentration | GFP, OD600, ICP-MS, survival rate |
| 1000 mg/L | Toxicity threshold — LD50 determination | Colony 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.
| pH | Environmental Context (Zambia) | Measurements |
|---|---|---|
| 2.5 | Active acid mine drainage leachate | GFP, OD600, ICP-MS |
| 3.5 | Tailing pond runoff | GFP, OD600, ICP-MS |
| 4.5 | Near-tailing agricultural soil leachate | GFP, OD600, ICP-MS |
| 5.5 | Mildly acidic Copperbelt soil | GFP, OD600, ICP-MS |
| 6.5 ★ | Neutral control (laboratory standard) | GFP, OD600, ICP-MS |
| 7.5 | Borehole drinking water (Kitwe) | GFP, OD600, ICP-MS |
| 8.5 | Alkaline mine neutralisation runoff | GFP, OD600, ICP-MS |
| 9.0 | Extreme 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.
| Stressor | Zambia-Specific Condition | Test Parameters | Output Measured |
|---|---|---|---|
| Temperature | Avg 24°C; dry season peak 38°C | 20, 28, 37, 42°C | OD600, GFP, %BSE |
| Co-metal toxicity | Cu²⁺ + Zn²⁺ + Co²⁺ + Pb²⁺ co-contamination | Single vs cocktail, 50 mg/L each | ICP-MS all ions, GFP |
| Desiccation | Dry season soil water activity < 0.85 | aw 0.85, 0.90, 0.95 via NaCl | OD600, colony viability |
| UV exposure | High solar UV at 12–15°S latitude | UV-C 254 nm: 0, 10, 30, 60 s pulse | Colony survival, DNA damage gel |
| Competing microbiome | Indigenous Copperbelt soil microbiome | 10% v/v heat-killed soil extract | GFP, 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:
| Layer | Composition | Function | Sourcing |
|---|---|---|---|
| Outer shell | Calcium alginate + CaCO₃ nanoparticles | pH buffering: neutralises acidic mine leachate to pH 5.5–6.5 before bacteria are exposed; structural integrity in soil | Food-grade alginate; CaCO₃ from local limestone |
| Middle membrane | Cellulose nanofibre + chitosan crosslink | Size-selective filter: 200 nm pores allow Cu²⁺ ions (0.73 Å) to enter freely; bacteria (1–2 µm) physically cannot escape | Local agricultural waste cellulose; chitosan import |
| Inner core | PVA + gelatin hydrogel + activated charcoal | Bacteria viability matrix at 10⁸ CFU/mL; activated charcoal provides passive metal co-adsorption during biological lag phase | Commercial PVA/gelatin; charcoal from local Copperbelt source |
Sub-aim 3b: Containment Validation & Kill-Switch Integration
Containment Validation
| Test | Protocol | Pass Threshold |
|---|---|---|
| Bacterial escape | Plate surrounding water on LB agar at 7, 14, 30 days | < 1 CFU/mL at 30 days |
| Ion permeability | ICP-MS of surrounding fluid vs bead interior after 24h Cu²⁺ exposure | Cu²⁺ enters freely; bacteria absent in external fluid |
| Mechanical durability | Compression to 50 kPa (equivalent to 30 cm soil overburden) | No structural failure; containment maintained |
| Biodegradation rate | Bury spent beads in Zambian soil analogue at 28°C; measure mass loss weekly | Full 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.
Sub-aim 3c: Commercial Product Formats & Digital Predictor App
| Format | Description | Use Case | Deployment |
|---|---|---|---|
| ZAMGEL Beads | 3–5 mm spheres, ~10⁸ CFU/bead | Mine water treatment ponds | Broadcast by hand or machine |
| ZAMGEL Sheets | 10×10 cm biodegradable mats | Soil surface tailing cap treatment | Lay directly on contaminated soil |
| ZAMGEL Cartridges | Inline filter column packed with beads | Borehole and drainage pipe treatment | Install 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
| Week | Aim | Activity |
|---|---|---|
| 1 | 1a | SRA/MG-RAST/IMG/M search for Konkola, Nchanga, Mufulira mine datasets; quality trim with fastp |
| 2 | 1a | MEGAHIT assembly → Prodigal ORF prediction → BLASTp + Prokka annotation of metal resistance genes |
| 3 | 1a | Cys-X-Cys motif filter → top 5 candidates selected; IQ-TREE 2 maximum-likelihood phylogenetic tree |
| 4 | 1b | AlphaFold3 structure prediction for all 5 candidates; retrieve .pdb files |
| 5 | 1b | PyMOL/ChimeraX binding pocket quantification: volume, Cys coordination geometry, pLDDT mapping |
| 6 | 1b | Benchling codon optimisation + Cello 2.0 logic gate design + Twist Bioscience DNA order |
| 7 | 2a | B. subtilis 168 transformation; colony PCR; Sanger sequencing verification |
| 8 | 2a | SDS-PAGE + western blot + GFP fluorescence microscopy to confirm MT expression |
| 9 | 2b | Cu²⁺ concentration gradient assays (0–1000 mg/L); ICP-MS; GFP plate reader; dose-response curve |
| 10 | 2c | pH stress assays (pH 2.5–9.0) at 50 mg/L Cu²⁺; identify operational pH window |
| 11 | 2d | Multi-stressor factorial experiment: temperature × co-metals × UV × desiccation × microbiome cocktail |
| 12 | 3a | ZAMGEL prototype fabrication: alginate outer shell + chitosan membrane + PVA/gelatin inner core |
| 13 | 3b | Containment validation: LB plating, ICP-MS permeability, compression testing, biodegradation assay |
| 14 | 3b | MazF/MazE kill-switch chromosomal integration + aTc withdrawal 48h death assay; CcdB/CcdA backup |
| 15 | 3c | Streamlit app prototype; ZEMA regulatory pathway draft; final in silico feasibility report |
Validation Criteria & Contingency Plans
| Experiment | Pass Threshold | If Fail — Contingency |
|---|---|---|
| AlphaFold3 pLDDT (binding domain) | > 85 on core domain; ipTM > 0.80 | Use SmtA (Synechococcus PCC 7942) as positive control scaffold; re-run with AlphaFold2 |
| GFP activation in Cu²⁺ media | > 5× fluorescence over background | Redesign 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 window | Active sequestration at pH 4.5–8.0 | Increase CaCO₃ loading in ZAMGEL outer shell; add internal carbonate buffer inside PVA core |
| ZAMGEL containment (30 days) | < 1 CFU/mL in surrounding medium | Increase chitosan crosslink density; reduce pore size to 100 nm |
| Kill-switch efficacy | 100% cell death within 48h of aTc removal | Switch 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