Consultas y Reporte de Experiencia
This paper extends the PIEZO1/PIEZO2 Markov kinetic model established in Caballero (2026) by formalizing the biophysical translation bridge between peripheral mechanotransduction (Hardware) and central neurochemical modulation (Software). Replacing unstable geometric approximations, we introduce a compartmental diffusive nanodomain model governed by empirical time constants for local Calcium exchange ([Ca2+]nano). We demonstrate that local calcium accumulation coupled with a fourth-order (n=4) allosteric Hill kinetics triggers a massive quantum exocytosis of GABA. By integrating this flux into the Goldman-Hodgkin-Katz (GHK) electrodiffusion equation, we report a critical hyperpolarization of the membrane potential (Vm ≈ -84 mV). This electrical shunt block provides a rigorous mechanistic basis for amygdalar decoupling, explaining the clinical emergence of long-term memories devoid of stress-induced somatic markers.
Introducción Matemática de la Sección de Goldman (GHK) Para mapear la transición desde la señalización química extracelular hasta el silenciamiento eléctrico celular, el sistema debe abandonar las aproximaciones de potencial fijo y adoptar la Teoría de Campo Eléctrico Constante de Goldman-Hodgkin-Katz (GHK).
El potencial transmembrana dinámico (Vm) se formula a partir del equilibrio de flujo difusivo y electroforético de los iones dominantes (K+, Na+, Cl-), expresado como:
Vm(t) = (RT / F) * ln( (P_K * [K+]out + P_Na * [Na+]out + P_Cl(t) * [Cl-]in) / (P_K * [K+]in + P_Na * [Na+]in + P_Cl(t) * [Cl-]out) )
Donde:
RT / F ≈ 26.72 mV a temperatura corporal (310.15 K).
P_K, P_Na representan las permeabilidades relativas basales de la membrana en reposo.
[Cl-]out es condicionado por la Capa 0 del protocolo (ingesta de solución salina acida), incrementando la disponibilidad de Cloro en el fluido intersticial a ~ 125 mEq/L.
Modulación Alostérica de la Permeabilidad al Cloro (P_Cl) En este marco expandido, la permeabilidad al Cloro no es una constante estática, sino una función acoplada a la densidad de ocupación del receptor GABA_A en la hendidura sináptica. El receptor ionotrópico se modela mediante una cinética de saturación no lineal derivada de la disponibilidad transitoria de GABA (GABA):
P_Cl(t) = P_Cl,base + Delta P_Cl,max * ( GABA / (GABA + K_m) )
Donde P_Cl,base = 0.10, Delta P_Cl,max = 1.5, y K_m = 5.0 uM representa la constante de activación media del canal conductor de Cloro.
Cuando la estimulación mecánica sostenida (Fase PIEZO1) eleva GABA por encima del umbral de saturación, el término de permeabilidad iónica al Cloro domina el denominador de la ecuación GHK debido al elevado gradiente extracelular de [Cl-]out. El flujo neto de cargas negativas hacia el interior celular fuerza matemáticamente a Vm a colapsar desde su potencial de reposo de -70 mV hasta el potencial de equilibrio de inversión del Cloro (-84 mV a -87 mV), induciendo un estado de shunting inhibitorio completo.
El presente trabajo constituye la expansión formal del modelo cuantitativo de transducción de señales de la Ingeniería Cognitiva Humana (IC-H). En esta etapa se formaliza la transición biofísica entre el hardware periférico mecanosensible y la modulación del software neuroquímico central.
El marco conceptual se estructura en un sistema de cuatro capas analíticas acopladas que rigen la transferencia de información biológica:
Capa 0 - Entorno Mecánico y Precondicionamiento: Incorpora la dinámica del estímulo físico entrante. Se define por un régimen de percusión rítmica (Fase I a 2.1 Hz) seguido de una pausa y una fase de presión estática sostenida (Fase III). Esta entrada modula la tensión elástica basal de la membrana y prepara el andamiaje para la apertura de los canales PIEZO.
Capa 1 - Dinámica Nanodominio y Citosólica de Calcio (Ca2+): Modela la entrada de calcio provocada por la deformación mecánica. En lugar de asumir un volumen de difusión homogéneo, se establece un nanodominio local ([Ca2+]nano) acoplado al citosol ([Ca2+]bulk) con constantes de tiempo empíricas de intercambio y aclaramiento.
Capa 2 - Transducción Quimio-Acústica y Liberación de GABA: El pico de concentración local de calcio en el nanodominio actúa como regulador alostérico directo de la exocitosis cuántica de GABA. Esta relación se gobierna mediante una cinética cooperativa de Hill de orden cuatro (n = 4), reflejando la ocupación de los sitios de unión del sensor presináptico.
Capa 3 - Ecuación de Goldman-Hodgkin-Katz (GHK) y Potencial Transmembrana (Vm): La acumulación de GABA en la hendidura sináptica incrementa dinámicamente la permeabilidad relativa al cloro (P_Cl). Al integrar esta función en la ecuación de electrodifusión de GHK, la entrada masiva de cargas negativas fuerza un colapso del potencial de membrana hacia valores de hiperpolarización (~ -84 mV a -87 mV). Este fenómeno genera un desacoplamiento o bloqueo eléctrico que desarticula la actividad redundante de las redes centrales de hipervigilancia.
A continuación se adjunta el script computacional en Python que simula de forma explícita la integración temporal de estas cuatro capas y grafica la respuesta dinámica del sistema:
import numpy as np
import matplotlib.pyplot as plt
# =============================================================================
# 1. PARÁMETROS GENERALES
# =============================================================================
dt = 0.01
T_total = 60.0
time = np.arange(0, T_total, dt)
N = len(time)
# Constantes físicas
R = 8.314
T_K = 310.15
F = 96485.0
RT_F = (R * T_K / F) * 1000.0 # mV
# =============================================================================
# 2. CAPA 0: ENTRADA MECÁNICA DEL PROTOCOLO
# =============================================================================
stim_mec = np.zeros(N) # estímulo mecánico normalizado (0 a ~1)
for i, t in enumerate(time):
if 0.0 <= t < 35.0:
# Activación pulsátil a 2.1 Hz
stim_mec[i] = max(0.0, np.sin(2.0 * np.pi * 2.1 * t))
elif 35.0 <= t < 45.0:
# Pausa respiratoria / descarga mecánica
stim_mec[i] = 0.0
elif 45.0 <= t < 55.0:
# Presión sostenida
stim_mec[i] = 1.0
else:
stim_mec[i] = 0.0
# Ganancia mecánica por “conductividad” tisular relativa
G_tisular = 1.0
stim_mec = G_tisular * stim_mec
# =============================================================================
# 3. CAPA 1: DINÁMICA DE CALCIO EN NANODOMINIO Y CITOSOL
# =============================================================================
Ca_bulk = np.full(N, 0.1) # uM
Ca_nano = np.full(N, 0.1) # uM
f_Ca = 0.02 # fracción efectiva de corriente que entra como Ca2+
tau_bulk = 1.2 # s, aclaramiento citosólico
tau_nano = 0.08 # s, intercambio local nanodominio-citosol
Ca_basal = 0.1 # uM
# Flujo mecano-dependiente de Ca2+ en unidades arbitrarias escaladas a uM/s
Ca_drive_max = 0.8
# =============================================================================
# 4. CAPA 2: SEÑAL QUÍMICA DEPENDIENTE DE CALCIO
# =============================================================================
GABA = np.zeros(N) # uM
GABA_max = 120.0 # uM/s equivalente fenomenológico
K_gaba = 1.5 # uM
n_hill = 4
tau_GABA = 0.25 # s
tau_clear = 4.0 # s
# =============================================================================
# 5. CAPA 3: POTENCIAL DE MEMBRANA POR GHK
# =============================================================================
V_membrana = np.zeros(N)
P_K = 1.0
P_Na = 0.04
P_Cl_base = 0.10
K_out = 4.5
K_in = 140.0
Na_out = 142.0
Na_in = 14.0
Cl_out = 125.0
Cl_in = 6.0
V_rest = -70.0
# =============================================================================
# 6. INTEGRACIÓN TEMPORAL
# =============================================================================
for i in range(1, N):
# --- Capa 1: Ca2+ ---
Ca_influx = Ca_drive_max * stim_mec[i-1] * f_Ca
dCa_nano = (Ca_influx - (Ca_nano[i-1] - Ca_bulk[i-1]) / tau_nano) * dt
dCa_bulk = ((Ca_nano[i-1] - Ca_bulk[i-1]) / tau_nano - (Ca_bulk[i-1] - Ca_basal) / tau_bulk) * dt
Ca_nano[i] = max(Ca_basal, Ca_nano[i-1] + dCa_nano)
Ca_bulk[i] = max(Ca_basal, Ca_bulk[i-1] + dCa_bulk)
# --- Capa 2: GABA ---
gaba_release = GABA_max * (Ca_nano[i-1]**n_hill) / (K_gaba**n_hill + Ca_nano[i-1]**n_hill)
dGABA = (gaba_release - GABA[i-1] / tau_clear) * dt
GABA[i] = max(0.0, GABA[i-1] + dGABA)
# --- Capa 3: GHK y conductancia de Cl- modulada ---
P_Cl_dyn = P_Cl_base + 1.5 * (GABA[i-1] / (GABA[i-1] + 5.0))
num = P_K * K_out + P_Na * Na_out + P_Cl_dyn * Cl_in
den = P_K * K_in + P_Na * Na_in + P_Cl_dyn * Cl_out
V_ghk = RT_F * np.log(num / den)
# Mezcla fenomenológica con reposo
V_membrana[i] = 0.7 * V_rest + 0.3 * V_ghk
# =============================================================================
# 7. GRÁFICOS
# =============================================================================
fig, axes = plt.subplots(4, 1, figsize=(11, 9), sharex=True)
axes[0].plot(time, stim_mec, color='blue', label='Estímulo mecánico')
axes[0].set_ylabel('u.a.')
axes[0].grid(True)
axes[0].legend()
axes[1].plot(time, Ca_nano, color='crimson', label='Ca2+ nanodominio')
axes[1].plot(time, Ca_bulk, color='orange', label='Ca2+ citosólico')
axes[1].set_ylabel('[Ca2+] (uM)')
axes[1].grid(True)
axes[1].legend()
axes[2].plot(time, GABA, color='purple', label='GABA')
axes[2].set_ylabel('[GABA] (uM)')
axes[2].grid(True)
axes[2].legend()
axes[3].plot(time, V_membrana, color='black', lw=2, label='Vm')
axes[3].axhline(V_rest, color='gray', linestyle='--', label='Reposo')
axes[3].set_ylabel('Vm (mV)')
axes[3].set_xlabel('Tiempo (s)')
axes[3].grid(True)
axes[3].legend()
plt.suptitle('Simulación IC-H: mecánica, Ca2+, señal química y Vm', fontsize=12, fontweight='bold')
plt.tight_layout()
plt.show()