2022 개정 교육과정 연계 온디바이스 신경전기생리학·생체물리 STEM 시리즈 #71

패치 클램프(Patch-Clamp) 단일 이온 통로 개폐 동역학 서보 엔진

단일 분자 마코프 상태 전이 추적 및 0-Alloc C# 피코암페어(pA) 급 직렬 저항($R_s$) 막전위 고정 제어 루프

ZeliDesk 신경전기생리학·생체전자팀 이온채널생물물리·뇌공학인터페이스 분과

2026년 9월 12일

1. 세포막 등가 RC 회로 및 호지킨-헉슬리(Hodgkin-Huxley) 채널 동역학

패치 피펫-세포막 경계는 피펫 직렬 저항 $R_s$, 세포막 저항 $R_m$, 세포막 정전용량 $C_m$의 병렬-직렬 등가 회로로 표현됩니다:

지령 전압 $V_{cmd}$와 실제 막전위 $V_m$ 사이에는 피펫 전류 $I$에 비례하는 오차가 발생합니다:

단일 이온 통로의 열림(Open, $O$)과 닫힘(Closed, $C$) 상태 간 전이는 전압 의존적 전이율 $\alpha(V), \beta(V)$를 갖는 마코프 2상태(Two-state Markov) 연속 확률 미분방정식을 따릅니다:

서보 제어 루프는 $I_{pipette} \cdot R_s$를 계산하여 명령 전압에 동적으로 가산함으로써 실제 막전위 $V_m \approx V_{cmd}$를 완전 고정(Clamping)합니다.

2. 0-Alloc C# 스택 기반 패치 클램프 막전위 고정 서보

기가옴 피펫 헤드스테이지의 미세 전류 신호로부터 $R_s$를 보상하고 단일 채널 개폐 이벤트를 0.1ms 단위로 식별하는 ref struct 엔진입니다.

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.Electrophysiology
{
    public readonly ref struct PatchClampTelemetry
    {
        public readonly double CommandVoltageMv;
        public readonly double MeasuredCurrentPa;
        public readonly double PipetteSeriesResistanceMOhm;
        public readonly double MembraneCapacitancePf;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public PatchClampTelemetry(double cmdMv, double currPa, double rsMOhm, double cmPf)
        {
            CommandVoltageMv = cmdMv;
            MeasuredCurrentPa = currPa;
            PipetteSeriesResistanceMOhm = rsMOhm;
            MembraneCapacitancePf = cmPf;
        }
    }

    public readonly ref struct VoltageClampResult
    {
        public readonly double CompensatedVoltageMv;
        public readonly double TrueMembraneVoltageMv;
        public readonly bool IsChannelOpen;
        public readonly double SingleChannelConductancePs;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public VoltageClampResult(double compMv, double trueMv, bool isOpen, double condPs)
        {
            CompensatedVoltageMv = compMv;
            TrueMembraneVoltageMv = trueMv;
            IsChannelOpen = isOpen;
            SingleChannelConductancePs = condPs;
        }
    }

    public ref struct PatchClampServo
    {
        private readonly double _compensationRatio; // 0.85 (85% 보상)
        private readonly double _singleChannelCurrentThresholdPa;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public PatchClampServo(double compRatio, double thresholdPa)
        {
            _compensationRatio = compRatio;
            _singleChannelCurrentThresholdPa = thresholdPa;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public VoltageClampResult Step(in PatchClampTelemetry telemetry)
        {
            // 1. 피펫 직렬 저항(Rs)에 의한 전압 강하 계산: delta_V = I (pA) * Rs (MOhm) * 1e-6 (mV)
            // pA = 1e-12 A, MOhm = 1e6 Ohm -> I * Rs = 1e-6 V = 1e-3 mV
            double voltageDropMv = telemetry.MeasuredCurrentPa * telemetry.PipetteSeriesResistanceMOhm * 1e-3;

            // 2. 직렬 저항 양의 되먹임(Positive Feedback) 보상 명령 전압
            double compVoltageMv = telemetry.CommandVoltageMv + (_compensationRatio * voltageDropMv);

            // 3. 실제 세포막 전위 추정
            double trueMembraneVoltageMv = compVoltageMv - voltageDropMv;

            // 4. 단일 채널 개폐 판별 (하프-진폭 역치 기준)
            bool isChannelOpen = Math.Abs(telemetry.MeasuredCurrentPa) > _singleChannelCurrentThresholdPa;

            // 5. 단일 채널 전도도 g = I / (V_m - E_rev) (pS)
            double singleChannelConductancePs = 0.0;
            if (isChannelOpen && Math.Abs(trueMembraneVoltageMv) > 1.0)
            {
                singleChannelConductancePs = (Math.Abs(telemetry.MeasuredCurrentPa) / Math.Abs(trueMembraneVoltageMv)) * 1000.0;
            }

            return new VoltageClampResult(compVoltageMv, trueMembraneVoltageMv, isChannelOpen, singleChannelConductancePs);
        }
    }
}

3. 교육적 가치 및 중고교 실험 시뮬레이션 적용