Spotlight #42 – Biologia komórki i biofizyka: błony, osmoza, podział komórki i maszyny molekularne

Biofizyka komórki leży na przecięciu fizyki, chemii i biologii. Te same równania ośrodków ciągłych, które opisują przepływ płynów (Naviera-Stokesa), dyfuzję (drugie prawo Ficka) i deformację sprężystą (prawo Hooke’a), rządzą transportem przez błony, mechaniką cytoszkieletu i propagacją sygnałów elektrycznych wzdłuż aksonów. Ten Spotlight rozwija ilościowe ramy dla czterech obszarów: struktury błony i osmozy, mechaniki cytoszkieletu, faz podziału komórki oraz motorów molekularnych.

1. Dwuwarstwa lipidowa i płynność błony

Błona komórkowa to dwuwarstwa w modelu płynnej mozaiki (Singer i Nicolson, 1972): dwa listki amfipatycznych cząsteczek fosfolipidów samoorganizują się tak, że hydrofobowe ogony kwasów tłuszczowych skierowane są do wnętrza błony, a hydrofilowe grupy fosforanowe — do wodnej cytoplazmy i przestrzeni zewnątrzkomórkowej. Błona ma grubość ok. 7–8 nm oraz moduł ściśliwości powierzchniowej K_A ≈ 0,24 N m−¹ i moduł zginania κ ≈ 20 k_BT.

Dyfuzja lateralna lipidów w obrębie listka podlega dwuwymiarowemu ruchowi Browna ze współczynnikiem dyfuzji D ≈ 1 µm² s−¹ — lipid przemieszcza się przez komórkę o średnicy 10 µm w około 100 sekund. Przejścia fazowe między fazą ciekłą uporządkowaną (L_o) a ciekłą nieuporządkowaną (L_d) leżą u podstaw tworzenia się tratw lipidowych; temperatura przejścia T_m rośnie wraz z długością łańcucha i stopniem nasycenia.

F_bend = ∫∫ [(κ/2)(c₁ + c₂ − c₀)² + κ̅ c₁c₂] dA

c₁, c₂  = krzywizny główne (m⁻¹)
c₀      = krzywizna spontaniczna
κ       ≈ 20 k₂T  (moduł zginania, typowa dwuwarstwa)
κ̅       ≈ −κ   (moduł gaussowski)

2. Ciśnienie osmotyczne i turgor

Osmoza to wypadkowy przepływ wody przez błonę półprzepuszczalną, napędzany gradientem potencjału chemicznego. Dla roztworów rozcieńczonych ciśnienie osmotyczne określa równanie van ’t Hoffa:

π = iMRT

π  = ciśnienie osmotyczne  [Pa]
i   = współczynnik van ’t Hoffa  (1 dla glukozy, 2 dla NaCl przy idealnym rozcieńczeniu)
M   = stężenie molowe [mol L⁻¹]
R   = 8,314 J mol⁻¹ K⁻¹
T   = temperatura bezwzględna [K]

Przy 37°C, 300 mOsm/kg daje π ≈ 7,7 atm ≈ 780 kPa

Dla komórki roślinnej wypadkową siłą napędową pobierania wody jest potencjał wody Ψ:

Ψ = Ψₛ + Ψ𝗣        (potencjał wody = składowa rozpuszczonych substancji + ciśnienie)

Ψₛ = −π    (składowa osmotyczna/rozpuszczonych substancji, zawsze ujemna)
Ψ𝗣 = P     (składowa ciśnienia/turgoru, dodatnia w komórce turgorowej)

Przepływ netto: Jᵨ = L𝗣 × (Ψ₀ᵩᵴ − Ψ𝚢ᵦᵬ)
L𝗣 = przewodność hydrauliczna błony

Komórka zwiotczała: Ψ𝗣 = 0, Jᵨ > 0 (napływ netto)
Komórka turgorowa:  Ψ = 0, Jᵨ = 0 (równowaga)
Plazmoliza:  Ψ𝗣 = 0, Ψₛ𝚢ᵦᵬ < Ψₛ₀ᵩᵴ (odpływ netto)

Interaktywna symulacja osmozy modeluje ten układ za pomocą dwóch równoczesnych wizualizacji: płótna komórkowego pokazującego 50 animowanych cząsteczek wody wykonujących błądzenie Browna, przekraczających błonę probabilistycznie na podstawie gradientu potencjału chemicznego i zmieniających stężenie wewnątrzkomórkowe oraz objętość wakuoli w czasie rzeczywistym; oraz towarzyszącego wykresu liniowego objętości w funkcji czasu. Pięć presetów obejmuje scenariusze fizjologiczne spotykane na kursach fizjologii roślin.

3. Mechanika cytoszkieletu

Cytoszkielet to sieć trzech typów filamentów polimerowych, które nadają komórkom sztywność mechaniczną i pełnią funkcję torów dla białek motorycznych. Ich właściwości fizyczne obejmują cztery rzędy wielkości sztywności na zginanie:

• Filamenty aktynowe (F-aktyna): średnica ≈ 7 nm, długość persystencji L_p ≈ 17 µm. Powstają przez bieżnikowanie (ATP-aktyna dołącza się na końcu kolczastym, ADP-aktyna dysocjuje na końcu spiczastym). Moduł Younga wzdłuż filamentu E ≈ 2 GPa.
• Mikrotubule (MT): średnica ≈ 25 nm, L_p ≈ 5,2 mm (sztywne pręty w skali komórkowej). Wykazują dynamiczną niestabilność: tubulina-GTP polimeryzuje na końcu plus; katastrofa następuje, gdy czapeczka GTP zostaje utracona. Tworzą wrzeciono mitotyczne.
• Filamenty pośrednie (FP): średnica ≈ 10 nm, L_p ≈ 1 µm (przypominające liny). Wimentyna, keratyna, lamina. Opierają się dużym deformacjom; lamina jądrowa (lamina A/C/B) używa filamentów pośrednich do mechanicznej ochrony chromatyny.

Efektywny moduł Younga całej komórki, mierzony metodą wgłębnikowania mikroskopem sił atomowych (AFM) lub aspiracją mikropipetą, waha się od ~0,1 kPa (leukocyty) przez ~10 kPa (komórki nabłonkowe) po ~40 kPa (chondrocyty) — o rzędy wielkości miększe niż większość materiałów inżynierskich.

4. Mitoza i mejoza

Podział komórki eukariotycznej obejmuje ściśle regulowany rozpad i ponowne złożenie otoczki jądrowej, dynamiczne wychwytywanie chromosomów przez mikrotubule wrzeciona mitotycznego oraz fizyczne rozdzielenie chromatyd siostrzanych (mitoza) lub par chromosomów homologicznych (mejoza I).

Fazy mitozy

Siedem etapów mitozy animowanych przez symulację:

1. Interfaza (G2): chromosomy zreplikowane, wciąż w postaci rozciągniętej chromatyny. Centrosomy zdublowane i zaczynają się rozdzielać.
2. Profaza: chromatyna kondensuje się w widoczne chromosomy (chromatydy siostrzane połączone w centromerze). Otoczka jądrowa zaczyna się rozpadać.
3. Prometafaza: otoczka jądrowa całkowicie rozpuszczona. Mikrotubule kinetochorowe (kMT) z obu biegunów wychwytują kinetochory. Rozpoczyna się kongresja chromosomów.
4. Metafaza: chromosomy ustawiają się na płytce metafazalnej (płaszczyzna równikowa). Punkt kontrolny składania wrzeciona (SAC) zapewnia, że wszystkie kinetochory znajdują się pod odpowiednim napięciem, zanim dozwolona zostanie anafaza.
5. Anafaza A+B: kohezyna jest przecinana przez separazę (aktywowaną degradacją sekuryny za pośrednictwem APC/C). Chromatydy siostrzane segregują do biegunów (anafaza A, skracanie kMT). Wrzeciono wydłuża się dzięki antyrównoległemu nakładaniu mikrotubul (anafaza B).
6. Telofaza: otoczki jądrowe formują się ponownie wokół zestawu chromosomów każdego bieguna. Chromosomy zaczynają się dekondensować.
7. Cytokineza: aktomiozynowy pierścień kurczliwy (w komórkach zwierzęcych) lub tworzenie płytki komórkowej (w komórkach roślinnych) fizycznie dzieli cytoplazmę na dwie komórki potomne.

5. Białka motoryczne i mechanochemia

Białka motoryczne cytoszkieletu przekształcają energię swobodną hydrolizy ATP (ΔG ≈ −54 kJ mol−¹ w warunkach fizjologicznych) w ukierunkowany ruch mechaniczny. Istnieją trzy nadrodziny, każda wyspecjalizowana w innym rodzaju toru filamentowego:

Kinezyna-1 (anterogradowa, skierowana ku końcowi plus):
  - Rozmiar kroku: 8 nm na cykl hydrolizy ATP
  - Siła zatrzymania: Fₛ ≈ 7 pN (mierzona pęsetą optyczną)
  - Prędkość: v ≈ 800 nm s⁻¹ przy nasyceniu ATP
  - Wydajność mechanochemiczna: η = Fₛ d / ΔGₐᵩᵰ ≈ 7,0 pN x 8 nm / 90 zJ ≈ 62%

Miozyna II (mięsień):
  - Skok roboczy: d ≈ 5–10 nm na cykl ATPazy
  - Siła na mostek poprzeczny: F ≈ 1–5 pN
  - Skok siłowy zależy od obrotu ramienia dźwigni: θₛ ≈ 70°

Dyneina (retrogradowa, skierowana ku końcowi minus):
  - Rozmiar kroku: zmienny 8/16/24/32 nm
  - Siła zatrzymania: Fₛ ≈ 6–7 pN
  - Niezbędna dla transportu wewnątrzwiciowego (IFT) oraz kMT kinetochorów

Ruchliwość pojedynczych cząsteczek kinezyny po raz pierwszy rozwiązali Howard, Hudspeth i Vale (1989) za pomocą testu ślizgu in vitro. Dziś eksperymenty z pęsetą optyczną bezpośrednio wizualizują pojedyncze kroki o długości 8 nm, potwierdzając mechanizm „ręka za rękę” (naprzemienne głowy) z termicznie napędzanym dyfuzyjnym poszukiwaniem między cyklami katalitycznymi.

6. Kanały jonowe i bioelektryczność

Zachowanie elektryczne błon pobudliwych (neuronów, mięśni) wynika z bramkowanych napięciem kanałów jonowych — transbłonowych porów białkowych, które przełączają się między stanem otwartym a zamkniętym w odpowiedzi na napięcie błonowe. Model Hodgkina-Huxleya (1952) pozostaje ilościowym fundamentem:

C𝑥 dV/dt = I𝉥ᵴ − 𝑔ⱺ𝑴 m³h(V − Eⱺ𝑴) − 𝑔𝑘 n⁴(V − E𝑘) − 𝑔ℒ (V − Eℒ)

m, h, n  = zmienne bramkowania (0–1)
dm/dt   = α𝔦(V)(1−m) − β𝔦(V) m   (aktywacja Na⁺, szybka)
dh/dt   = αₕ(V)(1−h) − βₕ(V) h   (inaktywacja Na⁺, wolna)
dn/dt   = αₙ(V)(1−n) − βₙ(V) n   (aktywacja K⁺)

Potencjał spoczynkowy:  Vₐ ≈ −65 mV  (olbrzymi akson kałamarnicy)
Potencjał czynnościowy:   V𝔪𝔭𝔩 ≈ +40 mV, czas trwania ≈ 1 ms

Schemat bramkowania m³h dla kanału Na⁺ odzwierciedla trzy niezależne bramki aktywacyjne (każda aktywuje się szybko przy depolaryzacji) oraz jedną bramkę inaktywacyjną (zamyka się wolniej po aktywacji). Iloczyn m³h daje dzwonowaty kształt przewodnictwa Na⁺, generujący narastanie potencjału czynnościowego. Człon n⁴ dla K⁺ generuje wolniejsze przewodnictwo repolaryzujące.

Łącząc skale

Biofizyka komórki to fizyka przekraczająca skale. Prawa, które rządzą osmotycznym transportem wody, są tymi samymi prawami termodynamicznymi, które rządzą zachowaniem dowolnej błony półprzepuszczalnej w równowadze. Moduły sprężystości cytoszkieletu podlegają tej samej fizyce polimerów w modelu łańcucha robaczkowego, który stosuje się do DNA lub syntetycznych polimerów w roztworze. Równania Hodgkina-Huxleya to układ równań różniczkowych zwyczajnych o tej samej strukturze matematycznej, co kinetyka chemiczna czy modele dynamiki populacji.

To, co czyni komórkę niezwykłą, to nie fakt, że podlega innym prawom, lecz to, że wykorzystuje nierównowagową termodynamikę — nieustanną hydrolizę ATP, utrzymywane gradienty jonowe, bieżnikujące polimery — do tworzenia struktur i funkcji, które są o rzędy wielkości bardziej złożone niż jakikolwiek pojedynczy proces fizyczny mógłby wytworzyć w izolacji.