Molekularne determinanty specyficznego i niespecyficznego wiązania białek z DNA

Identyfikator grantu: PT00818

Kierownik projektu: Jacek Czub

Realizatorzy:

  • Cyprian Kleist
  • Miłosz Wieczór
  • Michał Badocha
  • Michał Jurkowski
  • Mateusz Kogut
  • Michał Kowalczyk
  • Hubert Biniek
  • Abolfazl Shiroudi
  • Mikołaj Ławicki
  • Maciej Łada

Politechnika Gdańska

Wydział Chemiczny

Gdańsk

Data otwarcia: 2020-04-14

Streszczenie projektu

Oddziaływania białek z DNA pełnią kluczową rolę w procesach powielania, utrzymywania oraz ekspresji informacji genetycznej. W zależności od swojej funkcji komórkowej białka wiążą się z DNA w sposób niezależny od sekwencji (wiązanie niespecyficzne) lub wykazują istotne powinowactwo do konkretnych motywów sekwencyjnych (wiązanie specyficzne). W obecnym projekcie zamierzamy zbadać wybrane molekularne mechanizmy leżące u podstaw specyficznego i niespecyficznego wiązania się białek z DNA. W szczególności, chcemy zbadać w jakim stopniu za rozpoznawanie konkretnych motywów odpowiadać może zależne od sekwencji zróżnicowanie elastyczności podwójnej helisy DNA. Na podstawie licznych eksperymentów twierdzi się, że białka, których związanie do DNA wymaga istotnego wygięcia podwójnej helisy mają wyższe powinowactwo do sekwencji DNA, dla których taka deformacja jest energetycznie korzystniejsza. Wykorzystując metody "alchemicznej" predykcji różnic energii swobodnej w symulacjach dynamiki molekularnej w połączeniu z odpowiednimi cyklami termodynamicznymi zamierzamy zdekomponować całkowite powinowactwo wybranych białek do DNA na wkłady odpowiadające i) oddziaływaniu z zasadami i szkieletem fosfocukrowym, ii) zależnej od sekwencji deformacji elastycznej DNA, iii) zależnej od sekwencji desolwatacji DNA. Taka dekompozycja pozwoliłaby po raz pierwszy ilościowo ustalić udział niebezpośrednich czynników w rozpoznaniu sekwencji DNA. Ponadto mamy zamiar przetestować postawioną przez nas poprzednio hipotezę, wedle której wszystkie reszty aminokwasowe zaangażowane w rozpoznanie sekwencji DNA można podzielić na dwie grupy. Pierwsza z nich zwiększa absolutne powinowactwo białka do danej sekwencji poprzez obniżenie energii swobodnej wiązania białka z tą konkretną sekwencją (odpowiada ona zatem za tzw. "pozytywną" selekcję). Druga grupa z kolei zwiększa relatywne powinowactwo białka do danej sekwencji poprzez zwiększenie energii swobodnej wiązania do konkurencyjnych miejsc wiązania wzdłuż podwójnej helisy DNA (odpowiada zatem za tzw. "negatywną" selekcję, redukując powinowactwo do sekwencji nietargetowych). W założonym celu zamierzamy wykorzystać m.in. szereg technik przyspieszonego próbkowania w celu wyznaczenia wkładów do energii swobodnej wiązania pochodzących od poszczególnych reszt aminokwasowych.Trzecim elementem planowowanych badań jest zrozumienie molekularnych podstaw rozpoznawania przez helikazy z rodziny DEAH/RHA niekanonicznych motywów DNA typu G-kwadrupleksów. Helikazy te wykazują niezwykle silne powinowactwo do G-kwadrupleksów zarówno DNA jak i RNA a wykazując zdolność do ich rozwijania biorą udział w mediowanej przez G-kwadrupleksy kontroli transkrypcyjnej i post-transkrypcyjnej. Aby zrozumieć mechanizm rozpoznania G-kwadrupleksów przez helikazy typu DEAH/RHA wyznaczymy profile energii swobodnej dla wiązania wybranego G-kwadrupleksu o topologii równoległej z helikazą DHX36 oraz jej poszczególnymi domenami zaangażowanymi w oddziaływanie z cząsteczką kwasu nukleinowego.
Ze względu na konieczność precyzyjnej analizy specyficznych oddziaływań białko-DNA, wszystkie planowane obliczenia wymagają wykorzystania pełnoatomowej dynamiki molekularnej. Badane układy (z wodą traktowaną eksplicite) złożone będą z 40--200 tys. atomów, przez co ich wydajne pełnoatomowe symulacje dynamiki molekularnej wymagają zrównoleglonych wersji silników MD oraz zasobów obliczeniowych dostępnych wyłącznie w centrach komputerowych. Ponieważ skale czasowe odpowiadające wiązaniu białek do DNA o kilka rzędów przewyższają te dostępne dla pełnoatomowej dynamiki molekularnej, zachodzi konieczność stosowania kosztownych obliczeniowo metod przyspieszających próbkowanie przestrzeni konformacyjnej, takich jak umbrella sampling czy wymiana replik. Wszystkie symulacje dynamiki molekularnej prowadzone będą przy pomocy pakietu Gromacs.

Publikacje

  1. M. Wieczór, J. Czub., Telomere uncapping by common oxidative guanine lesions: Insights from atomistic models. , Free Rad. Bio. Med. 148, (2020) 162-169
  2. B. Zalewska-Piątek, M. Olszewski, T. Lipniacki, S. Błoński, M. Wieczór, P. Bruździak, A. Skwarska, B. Nowicki, S. Nowicki, R. Piątek., A shear stress micromodel of urinary tract infection by the Escherichia coli producing Dr adhesin, PLoS Path. 16, (2020) e1008247
  3. K.A. Hossain, M. Jurkowski, J. Czub, M. Kogut, Mechanism of recognition of parallel G-quadruplexes by DEAH/RHAU helicase DHX36 explored by molecular dynamics simulations. , Comput. Struct. Biotechnol. J. 19, (2021) 2526-2536
  4. Michał Jurkowski, Dlaczego G-kwadrupleksy wykazują silną tendencję do przyjmowania konformacji prawoskrętnych?, praca magisterska -, (2021) -
  5. Kazi Amirul Hossain, Mateusz Kogut , Joanna Słabonska , Subrahmanyam Sappati, Miłosz Wieczór, Jacek Czub, How acidic amino acid residues facilitate DNA target site selection, Proc Natl Acad Sci USA 120, (2023) e2212501120


← Powrót do spisu projektów

CONTACT

Our consultants help future and novice users of specialized software installed on High Performance Computers (KDM) at the TASK IT Center.

Contact for High Performance Computers, software / licenses, computing grants, reports:

kdm@task.gda.pl

Administrators reply to e-mails on working days between 8:00 am – 3:00 pm.