1 METODY MODELOWANIA CZĄSTECZKOWEGOCel nauczania: Zapoznanie studentów z wykorzystaniem metod modelowania cząsteczkowego w projektowaniu nowych leków. Omówienie wykorzystania metod komputerowych w wyznaczaniu zależności między budową leku a jego działaniem biologicznym. Program nauczania: 1. Wprowadzenie do modelowania cząsteczkowego poprzez omówienie metod mechaniki molekularnej, metod półempirycznych i metod ab initio. 2. Wykorzystanie metod modelowania cząsteczkowego do wyznaczania właściwości fizykochemicznych związków biologicznie aktywnych i określania zależności między budową leku a jego działaniem biologicznym (QSAR). 3. Określanie aktywnych konformacji leków przy wykorzystaniu analizy konformacyjnej i dynamiki molekularnej. 4. Poszukiwanie receptorów biologicznych, wyznaczanie budowy trójwymiarowego farmakoforu i struktur wiodących nowych leków w elektronicznych bazach danych. 5. Projektowanie leków z uwzględnieniem budowy miejsca wiążącego receptora biologicznego.

2

3

4

5

6

7

8

9 Bioinformatyka Nauka zajmująca się przepływem informacji w układach biologicznych Bioinformatyka zajmuje się tworzeniem i wykorzystaniem biologicznych baz danych. Łączy osiągnięcia biologii molekularnej i genetyki z metodami statystycznymi i osiągnięciami technologii informatycznej.

10 Bioinformatyka strukturalnaDział bioinformatyki zajmujacy się organizacją, przechowywaniem oraz analizą informacji strukturalnych białek, kwasów nukleinowych oraz ich kompleksów z różnymi ligandami w celu zrozumienia układów biologicznych na poziomie atomowym.

11 Wprowadzenie: 1954: wyznaczenie struktury białka (insulina - zespół Sangera) 1958: pierwsza 3D structura (X-ray) białka (mioglobina - Kendrew) 1972: pierwsze sekwencje DNA 1977: opracowanie szybkiej metody sekwencjowania DNA i białek (Gilbert i Sanger) 1986: PCR (szybka metoda kopiowania sekwencji DNA) 1992: wyznaczenie sekwencji chromosomu III drożdży (3*105 bp) 1995: wyznaczenie sekwencji bakterii Haemophilus influenzae (2*106 bp) 1998: odkrycie zjawiska interferencji mRNA (zjawisko wyciszania albo wyłączenia ekspresji genu przez dwuniciowy RNA) 1999: wyznaczenie sekwencji genomu wielokomórkowego organizmu (Caenorhabditis elegans) (108 bp) 2001: 92% wyznaczenie sekwencji ludzkiego genomu (3*109 bp) 2003: wyznaczenie sekwencji ludzkiego genomu (3*109 bp)

12 Przepływ informacji w układach biologicznychDNA 5’ 3’ Transkrypcja interferencja mRNA mRNA Translacja Poli-peptyd Fałdowanie białko Funkcja

13 Od genu poprzez strukturę do funkcjibiałko struktura funkcja

14

15 Bazy sekwencji DNA EMBL (Europe) http://www.ebi.ac.uk/embl/GenBank (USA) EMBL (Europe) DDBJ (Japan)

16

17

18

19 Mikromacierze DNA Metoda wyznaczania ekspresji mRNA w komórceKonstrukcja mikromacierzy: Przyłączanie krótkich sekwencji nukleotydowych do szkła w ściśle określonym położeniu Umożliwia jednoczesny pomiar aktywności około genów Nie wymaga dużej ilości mRNA Technika: precyzyjne roboty fotolitografia

20 Mikromacierz DNA - fotolitografia

21 Tworzenie mikromacierzy DNA1- próby (komórki) 2- ekstrakcja mRNA 3- przepisanie mRNA na oznakowany fluoroscencyjnie cDNA 4- hybrydyzacja 5- skanowanie 6- wizualizacja

22 Wyniki Affymetrix Genechip Mikromacierz DNA

23 Analiza wyników Korelacja pomiędzy ekspresją genu i jego funkcjąNajczęściej poszukujemy genów, których ekspresja może korelować z typem komórki (np. poszukujemy genowego wzorca komórek nowotworowych raka piersi MCF-7). Ten genowy „odcisk palca” powinien w najbliższym czasie stać się ważną metodą w diagnostyce chorób nowotworowych.

24

25

26

27

28 Topoisomerase II

29 Piękno nauki

30

31 Otrzymywania struktur trójwymiarowych kwasów nukleinowych i białek- Krystalografia rentgenowska Badania NMR mikroskopia elektronowa modelowanie cząsteczkowe (?)

32 Krystalografia rentgenowska

33 Poszukiwanie leków in silico

34 Wybór receptora biologicznegoStruktury trójwymiarowe ligand – receptor biologiczny Budowa modelu receptor – ligand Określenie kluczowych oddziaływań ligand - receptor projektowanie nowych struktur ligandów Analiza konformacyjna Solwatacja Określenie energii wiązań ligandów Synteza chemiczna Testy biologiczne

35 Hyperpowierzchnia energii potencjalnejOptymalizacja geometrii = minimalizacja energii Potrzebujemy zależności funkcyjnej wiążącej położenie atomów z energią: E = f(Q)

36 Hyperpowierzchnia energii potencjalnej

37 Hyperpowierzchnia energii potencjalnej

38 Cel: wyznaczenie stabilnej struktury cząsteczki lub układu cząsteczekWyznaczenie hyperpowierzchni energii potencjalnej układu molekularnego E = f(Q) minimum lokalne minimum globalne struktura przejściowa Cel: wyznaczenie stabilnej struktury cząsteczki lub układu cząsteczek

39 Metody obliczeniowe Teoria orbitali molekularnychMetody ab initio w przybliżeniu Hartree-Focka GAMESS, Gaussian Teoria funkcjonału gęstości (DFT) Metody półempiryczne CNDO, ZINDO, AM1, PM3, PM5 Mechanika molekularna MM, AMBER, CHARMM, SYBYL

40

41 Metody ab initio w przybliżeniu Hartree-FockaPrzybliżenie Borna-Oppenheimera ruch elektronów jest niezależny od ruchu jądra Przybliżenie Hartree-Focka elektrony znajdują się w uśrednionym polu innych elektronów Przybliżenie LCAO (orbitale cząsteczkowe są liniową kombinacją orbitali atomowych)

42 Metody ab initio w przybliżeniu Hartree-FockaDokładność wyników zależy od stopnia korelacji elektronowej i wielkości bazy funkcyjnej Koszt obliczeń wzrasta szybko jak wzrasta wielkość bazy funkcyjnej i ilość elektronów Czas CPU jest proporcjonalny do n4, gdzie n odpowiada liczbie funkcji bazowych

43 Mechanika kwantowa: wynikiGeometria cząsteczki: wyniki bardziej dokładne niż w innych metodach możliwa do wyznaczenia dla szerokiej klasy związku Właściwości elektronowe: rozkład gęstości elektronowej moment dipolowy potencjał elektrostatyczny Właściwości spektroskopowe (widma IR, UV, NMR, ESR) Struktury związków przejściowych i kompleksów aktywnych

44 Mechanika kwantowa: wynikiOrbitale cząsteczkowe: Orbitale graniczne – określają gdzie są najbardziej reaktywne elektrony (HOMO) i gdzie podążą (LUMO) niezbędne do wyjaśnienia chemicznej reaktywności

45

46 Teoria funkcjonału gęstościGęstość elektronowa i funkcja falowa mogą być używane alternatywnie to opisu stanu podstawowego układu Istnieje funkcjonał gęstości dla określonej liczby elektronów wyznaczający energię układu Czas CPU jest proporcjonalny do n3, gdzie n odpowiada liczbie funkcji bazowych

47

48 Metody półempiryczne Pominięcie bezpośredniego rozpatrywania elektronów niewalencyjnych całki rezonansowe Hcore są empirycznymi funkcjami parametrycznymi Parametry są dobierane aby najlepiej odtwarzać wartości eksperymentalne (np. ciepło tworzenia) Czas CPU jest proporcjonalny do n2 , gdzie n odpowiada liczbie funkcji bazowych

49 Porównanie metod półempirycznych

50

51 Mechanika molekularnaMetoda przewidywania stabilnej konfiguracji poprzez minimalizację elektronowej energii układu cząsteczek Pole siłowe – to wyrażenie opisujące energię układu jako położenia jej jąder (1) (2) (3)

52 Pole siłowe Amber Prosta funkcja harmoniczna Stała siłowaDługość wiązania słabnie

53 Pole siłowe Amber Rozważamy tylko stany podstawoweWybór „typu atomu” (parametryzacja) jest podstawowy dla wyników obliczeń: np. AMBER ma 5 typów tlenu: karbonylowy, alkoholowy, kwasowy, estrowy/eterowy i występujący w wodzie Netropsyna

54 oddziaływania niewiążącePole siłowe Amber: oddziaływania niewiążące Oddziaływania Van der Waalsa Oddziaływania elektrostatyczne Wartości qi i qj wyznaczamy metodami ab initio i półempirycznymi

55 Porównanie metod obliczeniowychmetody ab initio i DFT Obliczenia dla propanu Poziom teorii metody półempiryczne mechanika molekularna Wielkość układu

56 Porównanie metod: czas obliczeń i jakość wynikówWłaściwości propanu

57 Zakład Syntezy i Technologii Środków LeczniczychUniwersytet Medyczny w Białymstoku Mickiewicza 2, Białystok tel. (085) fax (48-85)