Analiza kontroli metabolicznej

Analiza kontroli metabolicznej ( MCA ) to matematyczne ramy opisu szlaków metabolicznych , sygnalizacyjnych i genetycznych . MCA określa ilościowo, w jaki sposób zmienne, elastyczne, takie jak strumienie i stężenia gatunków , zależą od parametrów sieci . W szczególności jest w stanie opisać, w jaki sposób właściwości zależne od sieci, zwane współczynnikami kontrolnymi , zależą od lokalnych właściwości zwanych elastycznościami lub współczynnikami elastyczności .

MCA został pierwotnie opracowany w celu opisania kontroli w szlakach metabolicznych, ale został później rozszerzony w celu opisania sieci sygnalizacyjnych i genetycznych . MCA była czasami nazywana Teorią Kontroli Metabolicznej, ale tej terminologii zdecydowanie sprzeciwiał się Henrik Kacser , jeden z założycieli [ potrzebne źródło ] .

Nowsze prace wykazały, że MCA można bezpośrednio odwzorować na klasyczną teorię sterowania i jako takie są one równoważne.

Teoria systemów biochemicznych jest podobnym formalizmem , choć ma raczej inne cele. Obie są ewolucją wcześniejszej analizy teoretycznej Josepha Higginsa.

Współczynniki kontrolne

Współczynnik kontrolny mierzy względną zmianę stanu ustalonego zmiennej systemowej, np. przepływu szlaku (J) lub stężenia metabolitu (S), w odpowiedzi na względną zmianę parametru , np. aktywności enzymu lub szybkości stanu ustalonego ( ) kroku . Dwa główne współczynniki kontrolne to współczynniki kontrolne strumienia i stężenia. Współczynniki kontroli strumienia są określone przez

i współczynniki kontroli stężenia wg

Współczynniki kontrolne mierzą wpływ zaburzeń w enzymie na obserwowalny stan stacjonarny, taki jak strumień lub stężenie metabolitu. Należy zauważyć, że efekt perturbacji może być pozytywny lub negatywny w zależności od kontekstu. Na rysunku zakłada się, że zaburzenie występuje w kroku trzecim. Rysunek zmodyfikowany i przerysowany z

.

Twierdzenia o sumowaniu

sumowaniu kontroli strumienia zostało odkryte niezależnie przez grupę Kacsera/Burnsa i grupę Heinricha/Rapoporta na początku lat 70. i pod koniec lat 60. XX wieku. Twierdzenie o sumowaniu kontroli strumieni sugeruje , że strumienie metaboliczne są właściwościami systemowymi i że ich kontrola jest wspólna dla wszystkich reakcji w systemie. Kiedy pojedyncza reakcja zmienia swoją kontrolę nad strumieniem, jest to kompensowane przez zmiany w kontroli tego samego strumienia przez wszystkie inne reakcje.

Współczynniki sprężystości

Współczynnik elastyczności mierzy lokalną reakcję enzymu lub innej reakcji chemicznej na zmiany w jego środowisku. Takie zmiany obejmują czynniki, takie jak substraty, produkty lub stężenia efektorów. Więcej informacji można znaleźć na dedykowanej stronie dotyczącej współczynników elastyczności .

Elastyczności to lokalne wielkości, które mierzą wpływ substratów, produktów i efektorów na daną reakcję. Pionowe czarne strzałki reprezentują katalizę enzymatyczną. Rysunek przerysowany i zmodyfikowany z

.

Twierdzenia o łączności

Twierdzenia o łączności to specyficzne zależności między elastycznościami a współczynnikami kontrolnymi. Są przydatne, ponieważ podkreślają ścisły związek między kinetycznymi poszczególnych reakcji a właściwościami systemowymi szlaku. Istnieją dwa podstawowe zestawy twierdzeń, jeden dotyczący strumienia, a drugi dotyczący koncentracji. { \ .

Równania kontrolne

Możliwe jest połączenie sumowania z twierdzeniami o łączności w celu uzyskania wyrażeń zamkniętych , które wiążą współczynniki kontrolne ze współczynnikami elastyczności. Rozważmy na przykład najprostszą nietrywialną ścieżkę:

Zakładamy, że są gatunkami ustalonej granicy ścieżka może osiągnąć stan współczynnik, a drugi krok . Koncentrując się na współczynnikach kontroli strumienia, możemy napisać jedno sumowanie i jedno twierdzenie o łączności dla tej prostej ścieżki:

Korzystając z tych dwóch równań, możemy rozwiązać współczynniki kontroli strumienia

Korzystając z tych równań, możemy przyjrzeć się prostym ekstremalnym zachowaniom. Załóżmy na przykład, że pierwszy krok jest całkowicie niewrażliwy na swój produkt (tj. Nie reaguje z nim), S, wtedy . W takim przypadku współczynniki kontrolne zmniejszają się do

To wszystko kontrola (lub czułość) jest na pierwszym stopniu. Ta sytuacja reprezentuje klasyczny etap ograniczania szybkości , o którym często wspomina się w podręcznikach. Przepływ przez ścieżkę jest całkowicie zależny od pierwszego kroku. W tych warunkach żaden inny krok na ścieżce nie może wpływać na strumień. Efekt jest jednak zależny od całkowitej niewrażliwości pierwszego kroku na jego produkt. Taka sytuacja prawdopodobnie będzie rzadkością w rzeczywistych ścieżkach. W rzeczywistości klasyczny etap ograniczania szybkości prawie nigdy nie był obserwowany eksperymentalnie. Zamiast tego obserwuje się zakres ograniczeń, przy czym niektóre kroki mają więcej ograniczeń (kontroli) niż inne.

Możemy również wyprowadzić współczynniki kontroli stężenia dla prostej dwuetapowej ścieżki:

Ścieżka w trzech krokach

Rozważ prostą trzyetapową ścieżkę:

gdzie są gatunkami o ustalonej granicy, równania kontrolne dla tej ścieżki można wyprowadzić w podobny sposób jak przypadku prostej dwuetapowej ścieżki, bardziej nudne.

gdzie D mianownik jest dany przez

Zauważ, że każdy wyraz w liczniku pojawia się w mianowniku, co zapewnia spełnienie twierdzenia o sumowaniu współczynnika kontroli strumienia.

Podobnie można wyprowadzić współczynniki kontroli stężenia, dla

I dla

Zauważ, że mianowniki pozostają takie same jak poprzednio i zachowują się jak czynnik normalizujący .

Wyprowadzenie za pomocą perturbacji

Równania kontrolne można również wyprowadzić, biorąc pod uwagę wpływ zakłóceń na system. Weź pod , że szybkości reakcji i są określane przez dwa enzymy i odpowiednio. Zmiana któregokolwiek z enzymów spowoduje zmianę poziomu stanu ustalonego szybkości reakcji w stanie ustalonym . niewielką zmianę wielkości } Będzie to miało wiele skutków, zwiększy się, się co wzrośnie W końcu system osiągnie nowy stan ustalony. Możemy opisać te zmiany, skupiając się na zmianie w i . Zmiana w wyznaczamy , w wyniku zmiany \ rozważamy tylko małe zmiany, możemy wyrazić zmianę kategoriach korzystając z relacji

pochodna jak reaguje na Pochodną można obliczyć, jeśli znamy prawo kursu dla . Na przykład, jeśli założymy, że prawo stopy wynosi to pochodna wynosi . Możemy również użyć podobnej w wyniku Tym razem zmiana w wynikiem dwóch zmian, zmiany w zmiany w . Możemy wyrazić te zmiany, sumując dwa indywidualne wkłady:

{ \ . Ponieważ pozwoliliśmy systemowi ustabilizować się do nowego stanu ustalonego, możemy również stwierdzić, że zmiana szybkości reakcji musi być taka sama (inaczej nie byłby w stanie ustalonym). Oznacza to, że możemy stwierdzić, że . Mając to na uwadze, zrównujemy oba równania i piszemy

Rozwiązując stosunek otrzymujemy:

W granicy, gdy zmniejszamy i zmniejszamy zmianę zbiega się do pochodnej :

Możemy pójść o krok dalej i skalować pochodne, aby wyeliminować jednostki. Mnożenie obu stron przez dzielenie obu stron przez daje skalowane pochodne: mi

po prawej stronie to elastyczności, a skalowany składnik po lewej stronie to skalowany współczynnik czułości lub współczynnik

Możemy jeszcze bardziej uprościć to wyrażenie. Szybkość reakcji funkcją liniową . Na przykład w równaniu Briggsa-Haldane'a szybkość reakcji jest określona wzorem . Różniczkując to prawo szybkości ze względu na i skalowanie daje .

Korzystanie z tego wyniku daje:

analizę można przeprowadzić, gdy . W tym przypadku uzyskujemy czułość odniesieniu do :

Powyższe wyrażenia mierzą, jak bardzo enzymy \ kontrolują stężenie półproduktu w . Możemy również rozważyć, w jaki sposób na szybkości reakcji w stanie ustalonym v wpływają perturbacje w i . Jest to często ważne dla inżynierów metabolicznych, którzy są zainteresowani zwiększeniem tempa produkcji. W stanie ustalonym szybkości reakcji są często nazywane strumieniami i określane skrótem i . W przypadku ścieżki liniowej, takiej jak ten przykład, oba strumienie są równe w stanie ustalonym, tak że strumień przez prostu określany jako . Wyrażając zmianę strumienia w wyniku zaburzenia mi i biorąc granicę jak przed otrzymaniem

Powyższe wyrażenia mówią nam, jak bardzo enzymy mi kontrolują strumień stanu ustalonego Kluczową kwestią jest tutaj to, że zmiany stężenia enzymu lub równoważnie aktywności enzymu muszą być spowodowane działaniem zewnętrznym.

Właściwości ścieżki liniowej

Liniowy łańcuch etapów reakcji katalizowanych przez enzymy bez pętli ujemnego sprzężenia zwrotnego jest najprostszą ścieżką do rozważenia. Poniższy rysunek przedstawia trzystopniowy łańcuch liniowy.

Linear chain of four reactions catalyzed by enzymes e1 to e4

jest odwracalna i że gatunki graniczne są ustalone, tak że ustalony

Rozwiązania analityczne współczynników kontrolnych można uzyskać, jeśli przyjmiemy prostą kinetykę masowego działania na każdym etapie reakcji:

gdzie i to odpowiednio stałe szybkości do przodu i do tyłu. i . Jeśli przypomnimy sobie, że stała równowagi dla tej prostej reakcji wynosi:

możemy zmodyfikować równanie kinetyczne działania masy w następujący sposób:

Mając dane szybkości reakcji, można opisać równania różniczkowe opisujące szybkości zmian gatunku. Na przykład tempo zmian będzie równe:

Ustawiając równania różniczkowe na zero, można wyprowadzić stężenie w stanie ustalonym dla gatunku, z którego można również wyznaczyć równanie strumienia szlaku . W trzyetapowej ścieżki stężenia w stanie stacjonarnym i są określone wzorem:

Wstawienie albo lub do jednego z praw szybkości da stały strumień ścieżki, :

W tym równaniu można zobaczyć wzór taki, że ogólnie dla ścieżki liniowej złożonej z kroków strumień ścieżki w stanie ustalonym jest określony wzorem:

Należy zauważyć, że strumień ścieżki jest funkcją wszystkich parametrów kinetycznych i termodynamicznych. Oznacza to, że nie ma jednego parametru, który całkowicie określałby strumień. Jeśli , to każdy enzym na szlaku ma pewien wpływ na przepływ Biorąc pod uwagę wyrażenie strumienia, możliwe jest wyprowadzenie współczynników sterowania strumieniem przez różniczkowanie i skalowanie wyrażenia strumienia. Można to zrobić w ogólnym przypadku :

Wynik ten daje dwa wnioski:

  • Suma współczynników kontroli strumienia wynosi jeden. Potwierdza to twierdzenie o sumowaniu.
  • Wartość indywidualnego współczynnika kontroli strumienia w liniowym łańcuchu reakcji jest większa lub mniejsza niż jeden:

Dla trzystopniowego łańcucha liniowego współczynniki kontroli strumienia są określone wzorem:

gdzie jest podane przez:

Biorąc pod uwagę te wyniki, nasuwają się natychmiastowe spostrzeżenia:

  • duże stałe równowagi, to znaczy , to , a pozostałe współczynniki do dążyć do zera.
  • Jeśli stałe równowagi są mniejsze, kontrola ma tendencję do rozkładania się na wszystkie trzy etapy.

Powodem, dla którego kontrola staje się bardziej rozproszona, jest to, że przy bardziej umiarkowanych stałych równowadze perturbacje mogą łatwiej przemieszczać się zarówno w górę, jak iw dół rzeki. Na przykład zaburzenie na ostatnim etapie ma lepszy wpływ na szybkość reakcji w górę, co skutkuje zmianą strumienia w

Ważny wynik można uzyskać, jeśli ustawimy wszystkie na równe sobie. W tych warunkach współczynnik kontroli strumienia jest proporcjonalny do licznika. To jest:

Jeśli założymy, że wszystkie stałe równowagi są większe niż 1,0, to ponieważ wcześniejsze kroki mają więcej że wcześniejsze kroki będą miały na ogół wysokie, większe współczynniki kontroli strumienia. W liniowym łańcuchu etapów reakcji kontrola strumienia będzie miała tendencję do przesunięcia w kierunku przodu ścieżki. Z punktu widzenia inżynierii metabolicznej lub celowania leków, należy preferować ukierunkowanie na wcześniejsze etapy szlaku, ponieważ mają one największy wpływ na przepływ szlaku. Należy zauważyć, że ta zasada dotyczy tylko ścieżek bez pętli ujemnego sprzężenia zwrotnego.

Oprogramowanie do analizy kontroli metabolicznej

Istnieje wiele narzędzi programowych, które mogą bezpośrednio obliczać elastyczność i kontrolować współczynniki:

Związek z klasyczną teorią sterowania

Klasyczna teoria sterowania to dziedzina matematyki zajmująca się sterowaniem systemami dynamicznymi w inżynieryjnych procesach i maszynach. W 2004 roku Brian Ingalls opublikował artykuł, w którym wykazał, że klasyczna teoria kontroli i analiza kontroli metabolicznej są identyczne. Jedyna różnica polegała na tym, że analiza kontroli metabolicznej była ograniczona do zerowych odpowiedzi częstotliwościowych, gdy rzutowano ją w domenie częstotliwości, podczas gdy klasyczna teoria kontroli nie nakłada takiego ograniczenia. Inną istotną różnicą jest to, że klasyczna teoria sterowania nie ma pojęcia stechiometrii i zachowania masy, co czyni ją bardziej kłopotliwą w użyciu, ale oznacza również, że nie rozpoznaje właściwości strukturalnych właściwych sieciom stechiometrycznym, które dostarczają użytecznych informacji biologicznych.

Zobacz też

Linki zewnętrzne