AND
OR
NOT
NAND
XOR
NOR
GENETIC LOGIC GATES · 6 OPERATORS
G = f(I₁, I₂)

Biochemical-Circuit-Cell-Programming

生化電路與細胞編程:基因邏輯閘、回饋控制網路與生物計算,將細胞重新編程為活的感測器與微型化工廠。

CIRCUIT PARTS~10⁴ BiobricksiGEM Registry中已表徵的標準化基因零件數量
RESPONSE TIME8–45 min典型轉錄調控電路的信號傳播延遲範圍
NOISE LEVELCV 0.15–0.6基因表達雜訊的變異係數,影響電路可靠性
LOGIC DEPTH≤ 7 layers目前可實現的最大基因邏輯閘串聯層數

將細胞重新編程:合成生物學的運算隱喻

如果細胞是硬體,DNA就是其作業系統——而合成生物學的目標是為這個作業系統撰寫新的應用程式。過去二十年,基因工程的典範從「複製貼上」(將一個基因從物種A移到物種B)進化到了「從頭編程」(de novo programming)——工程師開始像設計電子電路一樣設計基因電路:有邏輯閘、有回饋迴路、有振盪器與雙穩態開關。細胞不再只是生產重組蛋白的工廠,而是可編程的活的計算機

這一轉變的關鍵推動力來自三個方向的交匯:DNA合成成本每五年下降一個數量級(現已低於每鹼基對0.05美元)、CRISPR-Cas技術提供了可定址的基因組讀寫介面、以及標準化基因零件庫(如iGEM Registry)使電路設計從手工藝邁向工程化。當一個大學生團隊可以在一個暑假中設計並建造一個檢測飲用水砷含量的活細菌感測器時,我們正站在生物工程民主化的門檻上。

基因邏輯閘:DNA編碼的布林運算

電子電路使用電壓的高低來代表0與1;基因電路則使用轉錄因子(transcription factor)的濃度作為邏輯訊號。一個基本的AND閘由兩個串聯的啟動子(promoter)構成——只有當轉錄因子A與轉錄因子B同時存在時,下游的報導基因(如GFP螢光蛋白)才會被轉錄。NOR閘則利用兩個並聯的抑制子(repressor)——任一輸入訊號的存在都會關閉輸出。

2012年,MIT的Voigt實驗室在E. coli中實現了首個完整的4-bit基因邏輯運算單元,包含AND、OR、NOT、NAND、NOR、XOR六種基本閘門。更令人印象深刻的是,這些閘門的輸出可以被串聯——一個閘門的輸出蛋白質成為下一個閘門的輸入轉錄因子——形成多層級的邏輯運算網路,其理論深度僅受限於訊號雜訊的累積與蛋白質半衰期導致的延遲。

Genetic circuit design
Fig 1. 合成生物學中的基因電路設計:啟動子、核糖體結合位點與編碼序列的標準化組裝Source: Unsplash

回饋控制與動態行為

如果基因邏輯閘是靜態的運算單元,回饋控制則賦予基因電路動態行為。正回饋(positive feedback)可以創造雙穩態開關(toggle switch)——細胞在兩個穩定狀態之間切換,一如電子正反器(flip-flop)的記憶功能。負回饋(negative feedback)則可以抑制雜訊、加速響應時間與產生持續振盪——著名的Repressilator電路(Elowitz & Leibler, 2000)利用三個串聯的抑制子形成環形負回饋,在E. coli中產生了週期約160分鐘的穩定螢光振盪。

更複雜的控制架構正在被實現:PID-like控制回路可以將基因表達水平精確維持在設定點,抵禦環境波動;前饋回路(feed-forward loop)可以產生脈衝響應與加速關閉;而細胞間的通訊系統——如群體感應(quorum sensing)的醯基高絲胺酸內酯訊號——使得多細胞協同計算成為可能,一個菌落可以作為一個分散式計算網路來執行單一細胞無法完成的任務。

基因電路雜訊模擬

以下Python程式模擬了轉錄雜訊對基因邏輯閘輸出可靠性的影響,使用Gillespie演算法進行隨機模擬。

GeneCircuitNoiseSim.pyPYTHON 3.11
import numpy as np

class GeneCircuitNoiseSim:
    """使用Gillespie演算法模擬基因邏輯閘的轉錄雜訊"""

    def __init__(self, k_tx=0.15, k_deg=0.02, input_level=10):
        self.k_tx = k_tx    # 基礎轉錄速率 (s⁻¹)
        self.k_deg = k_deg  # mRNA降解速率 (s⁻¹)
        self.input_level = input_level

    def simulate(self, duration=3600, gate_type='AND'):
        """Gillespie模擬基因閘門輸出 (mRNA分子數軌跡)"""
        t, mrna = 0, 0
        history = [(t, mrna)]
        while t < duration:
            # 轉錄傾向 (依賴閘門類型與輸入)
            if gate_type == 'AND':
                tx_propensity = self.k_tx * (self.input_level / 20) ** 2
            elif gate_type == 'OR':
                tx_propensity = self.k_tx * (1 - np.exp(-self.input_level / 5))
            else:  # NOT gate
                tx_propensity = self.k_tx * np.exp(-self.input_level / 5)
            total_rate = tx_propensity + self.k_deg * mrna
            if total_rate <= 0: break
            tau = np.random.exponential(1 / total_rate)
            t += tau
            if np.random.rand() < tx_propensity / total_rate:
                mrna += 1
            else:
                mrna -= 1
            history.append((t, mrna))
        return history

sim = GeneCircuitNoiseSim(k_tx=0.15, input_level=10)
traj = sim.simulate(duration=1800, gate_type='AND')
print(f"最終mRNA數: {traj[-1][1]} | 平均: {np.mean([p[1] for p in traj]):.1f}")

從感測器到微型工廠:應用的邊界

基因電路的第一波商業應用聚焦於全細胞生物感測器——工程菌在檢測到特定分子(重金屬、病原體、爆炸物前體)後產生可量測的螢光或色素訊號。這些活感測器的優勢在於其自複製性(zero-power deployment)與對複雜樣品基質的耐受性。在環境監測中,工程菌被固定在微流體晶片上,連續監測河流或工業廢水中的特定汙染物。

更具顛覆性的應用在於生物製造——基因電路可以在發酵過程中動態調節代謝通量,在生長期最大化生物量積累、在生產期將全部代謝資源轉向目標產物的合成。這種動態調控可將發酵產率提升30%至200%。而長遠來看,基因電路的最終目標是實現真正的可編程細胞療法——CAR-T細胞在識別到癌細胞後,透過內建的邏輯電路決定是否啟動細胞毒性、是否分泌細胞激素招募更多免疫細胞、以及在清除癌細胞後自動進入休眠。

免責聲明:本文所述基因電路為研究階段概念,實際基因表達受啟動子強度、質體拷貝數、宿主代謝負荷等多因素影響。任何臨床應用需經嚴格的藥理學與倫理審查。