PL 
EnglishDomyślny🧬 Mechanika DNA: długość perzystencji i łańcuch robaczkowy
Interaktywna symulacja modelu łańcucha robaczkowego DNA pod wpływem wirtualnej pęsety optycznej. Zobacz, jak termiczne wichrowanie zapada się w rozciągniętą cząsteczkę, i porównaj bieżącą krzywą siła–wydłużenie ze wzorem Marko-Siggi.
O tej symulacji
Ta symulacja traktuje DNA jako półelastyczny polimer opisany modelem łańcucha robaczkowego (WLC). Łańcuch koralikowo-prętowy złożony z 32 segmentów jest propagowany za pomocą przetłumionej dynamiki Langevina: każde złącze niesie rzeczywistą karę zginania Kratky'ego-Poroda proporcjonalną do długości perzystencji Lp podzielonej przez długość segmentu, sąsiednie koraliki utrzymuje w odległości sztywna sprężyna rozciągająca, a każdy koralik otrzymuje termiczne szturchnięcie zgodne z twierdzeniem fluktuacyjno-dysypacyjnym w wybranej temperaturze. W efekcie powstaje łańcuch, który wichruje się zgodnie z poprawną statystyką skorelowanej krzywizny prawdziwego półelastycznego polimeru, a nie sztywnego pręta czy w pełni losowego łańcucha o swobodnych złączach.
Jeden koniec cząsteczki jest zakotwiczony; drugi trzyma wirtualna pęseta optyczna, wywierająca stałą siłę ciągnącą. Wraz ze wzrostem siły bieżący punkt siła–wydłużenie jest nanoszony na wykres względem wzoru interpolacyjnego Marko-Siggi, F(x) = (kT/Lp)[¼(1−x/L)⁻² − ¼ + x/L] — podręcznikowego wyniku zweryfikowanego bezpośrednio przez eksperymenty z pęsetami optycznymi i magnetycznymi rozciągającymi pojedyncze cząsteczki DNA. Przy małej sile łańcuch opiera się rozciąganiu entropowo — to fluktuacje termiczne, a nie wiązania chemiczne, przeciwdziałają wyprostowaniu. Blisko pełnej długości konturu krzywa gwałtownie się stroma, gdy cząsteczka zbliża się do reżimu entalpowego, w którym zaczyna się rzeczywiste rozciąganie wiązań szkieletu.
Najczęściej zadawane pytania
Czym jest długość perzystencji?
Długość perzystencji Lp to skala długości, na której kierunek polimeru „zapomina” swoją początkową orientację na skutek termicznego zginania. Dla dwuniciowego DNA wynosi ona około 50 nm, czyli mniej więcej 150 par zasad — znacznie sztywniej niż wiotki łańcuch, ale znacznie bardziej elastycznie niż sztywny pręt w skali całych chromosomów.
Czym jest model łańcucha robaczkowego (WLC)?
Model WLC, zwany też modelem Kratky'ego-Poroda, traktuje polimer jako w sposób ciągły elastyczny pręt z energią zginania, która przeciwstawia się krzywiźnie. Jest to standardowy statystyczno-mechaniczny opis półelastycznych biopolimerów, takich jak DNA, poprawnie przewidujący zarówno równowagowe zwijanie nierozciągniętej cząsteczki, jak i jej reakcję na przyłożoną siłę.
Czym jest wzór Marko-Siggi?
To wzór interpolacyjny, F(x) = (kT/Lp)[¼(1−x/L)⁻² − ¼ + x/L], który przybliża siłę potrzebną do utrzymania łańcucha robaczkowego przy danym ułamkowym wydłużeniu x/L. Wyprowadzili go John Marko i Eric Siggia, a jego słuszność potwierdzono bezpośrednio w eksperymentach rozciągania pojedynczych cząsteczek DNA.
Dlaczego DNA zachowuje się jak entropowa guma przy małej sile?
Przy małym wydłużeniu rozciąganie dwóch końców zmniejsza liczbę dostępnych termicznie zgiętych konfiguracji, jakie może przyjąć łańcuch. Siła przywracająca wynika z tej utraty entropii konfiguracyjnej, dokładnie jak przy rozciąganiu gumki, a nie z rozciągania jakiegokolwiek wiązania chemicznego. Dlatego dolny zakres siły na krzywej nazywa się reżimem entropowym.
Co dzieje się blisko pełnej długości konturu?
Gdy wydłużenie zbliża się do długości konturu L, niemal cały termiczny luz znika, a siła potrzebna do dalszego rozciągania gwałtownie rośnie — to człon (1−x/L)⁻² we wzorze Marko-Siggi. Bardzo blisko pełnego wydłużenia prawdziwe DNA przechodzi w reżim entalpowy, w którym zaczynają się rozciągać same wiązania szkieletu cukrowo-fosforanowego.
Jak pęsety optyczne i magnetyczne mierzą to w laboratorium?
Pęsety optyczne wykorzystują skupioną wiązkę lasera do uwięzienia mikroskopijnego koralika przyczepionego do jednego końca pojedynczej cząsteczki DNA, podczas gdy drugi koniec jest zamocowany do powierzchni lub drugiego koralika; przesuwanie pułapki rozciąga cząsteczkę, a przekaz pędu lasera wskazuje siłę. Pęsety magnetyczne zamiast tego wykorzystują magnetyczny koralik i gradient pola, aby przyłożyć kontrolowaną, stałą siłę ciągnącą. Obie techniki pozwalają zmierzyć siły rzędu pikoniutonów i wydłużenia rzędu nanometrów na pojedynczej cząsteczce.
Jak temperatura zmienia symulację?
Podniesienie temperatury zwiększa termiczne szturchnięcia działające na każdy koralik (twierdzenie fluktuacyjno-dysypacyjne wiąże siłę szumu bezpośrednio z kT), podczas gdy wewnętrzna sztywność zginania cząsteczki pozostaje stała. Efektem netto jest krótsza efektywna długość perzystencji w wyższej temperaturze — łańcuch wichruje się bardziej przy danej sile ciągnącej, dokładnie tak jak prawdziwe DNA staje się względnie bardziej wiotkie po podgrzaniu.
Dlaczego energia zginania jest proporcjonalna do Lp/l?
W zdyskretyzowanym modelu Kratky'ego-Poroda kara za zginanie między dwoma sąsiednimi segmentami o długości l jest tak dobrana, by uśredniona po zespole korelacja między kierunkami segmentów zanikała jak exp(−l/Lp). Wymaga to sztywności rzędu kT·Lp/l na złącze — dokładnie tej zależności użyto do napędzania segmentów w tej symulacji.
Czy ta symulacja jest ilościowo dokładna?
Potencjał zginania, termostat Langevina i prawo siły Marko-Siggi to prawdziwa, publikowana fizyka stosowana w nauce o polimerach i biofizyce pojedynczych cząsteczek. Rozstaw koralików i dynamika rozciągania są uproszczone na potrzeby wizualizacji w czasie rzeczywistym, ale zgodność między symulowanym punktem siła–wydłużenie a krzywą teoretyczną jest prawdziwa, a nie zaskryptowana.