Q S AR = QUANTITATIVE STRUCTURE - ACTIVITY RELATIONSHIPS KVANTITATIVNÍ VZTAHY MEZI (CHEMICKOU) STRUKTUROU A (BIOLOGICKOU) AKTIVITOU © Oldřich Farsa 2018 Hledáme vztah, v němž je biologická aktivita funkcí struktury, resp. parametrů, které se strukturou souvisejí A= f (struktura) 2 základní přístupy QSAR •regresní analýza - hledá matematický popis funkce, většinou s použitím lineární apod. regrese •empirické metody - hledají pouze extrémy (maxima a minima) dané funkce Regresní analýza hledá rovnici ve tvaru A = a+ ax+ax+...ax, 0 11 2 2 n n kde A je biologická aktivita, x jsou parametry vyplývající ze struktury sloučeniny, a , b jsou regresní koeficienty (aQ ...absolutní člen), získané výpočtem. V případě tzv. Hanschovy metody jsou x fyzikálně- chemické parametry vyplývající ze struktury, v případě tzv. Free-Wilsonova přístupu parametry x vyjadřují přítomnost nebo nepřítomnost určitého substituentu nebo strukturního fragmentu v molekule. Hanschova metoda regresní analýzy A = a + a x + ax +...ax 0 11 2 2 n n A ... biologická aktivita, často kvůli linearizaci vztahu v převrácené hodnotě nebo logaritmu (1/MIC ... převrácená hodnota minimální inhibiční koncentrace u antimikrobních látek, log LD50, logaritmy konstant vyjadřujících afinitu látky k receptoru, také ale některé fa rm ako kin etické parametry, např. log BB, vyjadřující schopnost látky pronikat hematoencefalickou bariérou aj. al ... an... regresní koeficienty, tj. koeficienty získané výpočtem s použitím např. lineární regrese "Klasické" parametry xi ... xn •hydrofobní •elektronové •sterické a) Hydrofobní parametry - v rovnici často v kvadratickém tvaru - vyjadřují poměr rozpustnosti látky ve vodě a v lipidech; mají pro účinek často zásadní význam, zejména pro průnik bariérovými systémy organismu; např. log P(oktanol/voda), log P(cyklohexan/voda) aj., parametr Rm z rozdělovači TLC na tzv. reverzní fázi (stacionární fáze lipofilní, mobilní fáze hydrofilní): ^ = 108(^-1), dále logaritmus tzv. kapacitního faktoru log k' z rozdělovači chromatografie plynové (GC) nebo vysokoúčinné kapalinové (HPLC) log£' = log( tr to ), kde tr je retenční čas dané látky a t0 je tzv. mrtvý retenční čas, tj. retenční čas látky, která se na koloně nezadržuje (např. u HPLC na reverzní fázi s oktadecylovaným silikagelem se pro stanovení t0 používá dusitan sodný). (Hanschův) lipofilní parametr n - pro řady (série) látek, obsahující na stejném strukturním fragmentu (nejčastěji benzenovém jádře) různé substituenty Px T^=l0g —= l0gPx-l0gP// * H kde Px je rozdělovači koeficient substituované látky a PH rozděl, koeficient látky nesubstituované. Vypočtené hydrofobní parametry Kromě experimentálně stanovených hydrofobních parametrů se dnes stále častěji používají odhady těchto parametrů získaných výpočtem podle různých algoritmů, z nichž nejjednodušší jsou patrně výpočetní postupy pro odhad log P (oktan o l/vod a) součtem log P inkrementů podle Rekkera a Nysse i loggia,-/,, kde f. je log P příslušného fragmentu a a. je počet, kolikrát se tento fragment vyskytuje, nebo podle Hansche (a Lea) vzorcem i j logP=X a,/,+X bjfj, kde f. je fragmentová konstanta, f. je korekční faktor a a.ab. je odpovídající počet výskytů daného parametru. Takto vypočtený parametr se běžně označuje ClogP. Často se používají i složitější postupy vyžadující počítače a vhodný software, který často umožňuje i optimalizaci struktury metodami molekulové mechaniky a výpočty dalších parametrů pro výpočty QSAR (pro PC např. Molgen, HyperChem). Shoda vypočteného odhadu log P, např. Clog P, s experimentálně stanovenou hodnotou, bývá pro různé výpočetní postupy velmi odlišná, často však pro výpočty postačí, existuje-li lineární závislost mezi log P a Clog P. Mlog P (Moriguchi et al.) •odvozen statistickou analýzou log P and strukturních dat 1230 sloučenin •definován vzorcem log P= 1.244(CXf6- l.OniNOft+OAMPKX (/=60.5) (í=53.5) (í=33.3) - 0.145(Í/B)° "s + 0.51L HB -+ 0,268POL -2,2\UMP + 0,912 A L K - 0.192RNG - 3,6840^ 4- 0.474A 02 (f-9.5) (f=!3.l) (f=22.l) (f^lOS) + l.582NCS+OJ7$BLM-1041 (/=16.4) (/=5.0) Jt=l230t r= 0.952, f-0.411, ^(13,1216)^900.4 Mlog P (pokracovani) kde Parameter Typcot Description CX MO Ufi HB POL AMP ALK RNG S02 SCS BLM N \ N N D N N N N \T n Summation of numbers of carbon and halogen atoms weighted by C: 1.0, F: 0.5,0:1.0, Br: 1.5, and 1: 2.0 Total number of N and O atoms Proximity effect of N,0; X-Y: 2.0, X-A-Y: 1.0 (Xt Y: N/O, A: Q S, or P) with a correction (-1) for carboxamide/sulfonamide Total number of unsaturated bondsexcept those in NOj Dummy variable for (lie presence of intramolecular hydrogen bond as ortho -OH and -CO- R, -OH and -NH 3t -N H 2 and -COOH, or 8-OH/NHj in qumolines, 5 or 8-OM/NH3 in quinoxah'nes, etc. Number of aromatic polarsubstituents (aromatic substituents excluding Ar-CX,- and Ar-CX = C<. X: CorH) Amphoteric property; a-aminoacid: 1.0, aminobenzoic acid: 0,5, pyridinecarboxylic acid: O.S D ummy varia ble for alka ne, al k e nc, cycloal kanc, or cycloaJkene (hydrocarbons with 0 or I double bond) Dummy variable for the presence of ring structures except benzene and itscondenccd rings (aromatic, heteroaromatic, and hydrocarbon rings) + Quaternary nitrogen: >N<, 1.0; N oxide, 0,5 Number of nitro groups Isothiocyanato (-N=C=S), L0; thiocyanato (-S-CsN), 0.5 Dummy variable for the presence of /Madam □) N, numerical variable: D, dummy variable. Variace na Hansche a Lea VG Viswanadhan, V. N.; Ghose, A. K.; Revankar, G. R.; Robins, R. K., J. Chem. Inf. Comput. Sci., 1989, 29, 163-172 •separace log P na „atomové hydrofobicity", lišící se podle toho, v jakém seskupení je atom umístěn •celkem 120 „druhů atomů" •odvozeny analýzou souboru 893 látek o znaných lipofilitách KLOP Klopman, G.; Li, Ju-Yun.; Wang, S.; Dimayuga, M.: J.Chem.Inf.Comput.Sci., 1994, 34, 752 •multivariační regresní analýza souboru log P 1663 organických látek vzniklo 98 seskupení s parciálními log P a 98 korekčních faktorů Příklad QSAR závislosti pouze s hydrofobním parametrem Účinnost fenolů jako induktorů apoptosy nádorových buněk Hansch, C. etal.: Bioorg. Med. Chem. 11, 617(2003) TU k 1. I>iii far QtA H 3 Ĺoanrjd Loj lfC ťredJ&j L/C Cb|F ■ Km rj'j.-.'. 2.7* 2JÍ -fcJQ* J.Ti lir 1.41 -4.14 |J7 4-<:nH.,.:0 A 2JŠ4 Jíl •: •:• 1*7 2*Í : ■- i: 1JŮ 4-i": ľ4-jií!=>:■] L 44 I.4Í □Udí ■ ■■ 1 Jwi ft> ta Jhůi jaT1 2*9 -4.77 4.M 4- l-jĹha>j] 2JW 22* -ĎĽÍI 2*0 ■' Hedo!1 ■44 I.7Í : 4* ll 4-MeS-jJŕLO> ■in 1.72 ÚlSS 2Jl 2.12 2.41 OH estradiol (1) log 1.C ,, 4.0 3.0 2.0 1.0 0.0-1.0 log 1/C = 0,67(±0,21)ClogP + 0.37(±0.63) n = 8, ŕ = 0,910, s = 0,201, q2 = 0,863 CH, HO- \\ // // \ // OH H3C diethylstilbestrol (8) Regresní krivka - Inhibitory apoptosy 1 C log P 2.0 3.0 4.0 b) Elektronové parametry -souvisejí s elektronovým obalem molekuly • Hammetova konstanta a - pro m- a p-substituované deriváty benzenu; vyjadřuje elektrondonorové (+M, +1) nebo elektronakceptorové (-M, -I) vlastnosti substituentu, resp. jeho vliv na elektronovou hustotu substituentu v pol. 1; definována na základě Hammetovy rovnice, odvozené pomoci rychlostních konstant hydrolýzy ethylesterů substituovaných benzoových kyselin uváděné obvykle ve tvaru log k= log kH+ p.a, kde k je rychlostní konstanta reakce, p tzv. reakční konstanta, charakterizující podmínky reakce (zejména rozpouštědlo a teplota) a kH rychlostní konstanta hydrolýzy nesubstituované látky, tj. ethylbenzoátu. (Pro o-deriváty tato lineární závislost neplatí, uplatňují se zde ve zvýšené míře sterické vlivy). Pro Hammetovu konstantu a tedy platí logk-logkH P a vyjadřuje míru elektrondonorového nebo elektronakceptorového působení substituentu, charakterizuje tedy souhrnně jeho mezomerní a indukcí efekt. Pro elekrondonorové substituenty nabývá záporných hodnot (a < 0), pro elektronakceptorové kladných (a > 0). Běžně se rozlišují následující typy Hammetovy konstanty: am konstanta pro substituent v m-poloze gp konstanta pro substituent v p-poloze Tyto konstanty (a lineární tvar Hammetovy rovnice) platí obecně pro reakce na druhém nebo dalším atomu od benzenového jádra a pro reakce na první atomu, který není silně konjugován, pro reakční centra na silně konjugovaném atomu (dusík anilinu, kyslík fenolu) byly pro substituenty v poloze para zavedeny tzv. duální konstanty, pro akceptorové substituenty ap-, použitelné např. pro disociace na dusíku anilinu nebo kyslíku fenolu a pro donorové ap+, používané u reakcí, kde na a-atomu vůči jádru vzniká aspoň částečný kladný náboj. •odvozené konstanty: a , a, a*, podobné Swain-Luptonovy konstanty ^ 9? •parametry ze spekter a jiných fyzikálních měření - chem. posuny ô z NMR, vlnová délka maxima absorbance Xmax z UV-VIS spekter, vlnočet v výrazného absorpčního pásu v IČ spektrech, půlvlnný potenciál Ei/2z polarografie •vypočtené elektronové parametry: polarita, polarizovatelnost, parciální náboj na určitém atomu c) Sterické parametry - vyjadřují "celkovou objemnost" molekuly •van der Waalsovy poloměry vp Taftová sterická konstanta Esodvozená pomocí rychlostních konstant hydrolýzy esterů alkanových kyselin E,= logf, Kh kde k je rychlostní konstanta hydrolýzy esteru příslušné alkanové kyseliny RCOOFT a k obdobná konstanta pro ester kyseliny octové CH COOR - standard. E není čistě sterickým parametrem, zahrnuje o S částečně i elektronové vlivy (+1). E (CH ) = 0, objemnější substituenty E < 0, méně objemné E > 0 So S S Sterická konstanta v odvozená z van der Waalsových poloměrů [pm] v[pm] = r[pm] - 120pm 120 pm je poloměr vodíku •pro symetrické substituenty získány z modelů, pro nesymetrické rovněž z kyselé hydrolýzy esterů Hodnoty sterických konstant Es a v pro některé substituenty Substituent -H -CH3 -C2H5 -C3H7 -CH(CH3)2 -C4H9 -CH2CH(CH3)2 -CH(CH3)C2H5 -C5HII cyklo-C5H9 cyklo-C6H Es +1,24 0 -0,07 -0,36 -0,47 -0,39 -0,93 -1,13 -0,40 -0.51 0,52 0,56 0,68 0,76 0,68 0,98 1,02 0,68 1.25 Další parametry užívané v QSAR ^většinou vypočtené >zahrnují často 2 i 3 typy vlivů (hydrof.+elektron.+ster.) "Klasické" •parachor = M_ (1/4) r d y , kde y je povrchové napětí, M molární hmotnost a d hustota. •molární refrakce (= molekulární refraktivita) MR (též CMR); def. vztah označován Lorentzova-Lorenzova rovnice 0 arD 2 1 , 2M (n-\)M MR-n--- +-=-L---—, n d (n +2) d kde n je index lomu. "Neklasické" •solvatační energie - je-li pro vodu, pak hydratační energie AGQW •povrchy molekuly různého typu - polární van der Waalsův, nepolární, přístupný vodě, dynamický polární (DPSA), topologický polární (TPSA) aj. •objemy molekuly - polární, přístupný vodě aj. Free- Wilsonova metoda regresní analýzy •hledá závislost biologické aktivity na přítomnosti nebo nepřítomnosti určitých substituentů nebo strukturních fragmentů v molekule. Jde o vlastně o statistickou separaci aktivity na příspěvky určitých částí molekul, tj. předpokládá se aditivita vlivu substituentů nebo jiných částí molekuly. Metoda vede k řešení soustav rovnic o větším počtu neznámých, které se v jednoduchých případech dají řešit maticovým počtem, jinak pomocí statistických programů umožňujících multilineární regresi (MLR). •obě metody lze též kombinovat. Část nezávislých proměnných pak vyjadřuje fyzikálně-chemické vlastnosti látek, další, nazývané „indikátorové proměnné" (symbol I), udávají přítomnost nebo nepřítomnost určitých molekulových fragmentů. Indikátorových proměnných bývá obvykle malý počet, často pouze jediná. Free- Wilsonova metoda •aby poskytovala prognózy, jsou nutná nejméně 2 obměňovaná místa •každá látka popsána rovnicí •indexy značíme jednak substituenty, jednak obměňovaná místa molekuly => každné nezávislé proměnné a každému regresnímu koeficientu přiřazeny 2 indexy, z nichž první označuje příslušnost k obměňované poloze (např. indexy 1,2 ... až j) a druhý určitý substituent či jiný obměňovaný fragment (např. indexy 1,2 ... až m u obměn na prvním obměň, místě; 1,2 ... až n na druhém obměň, místě; 1,2 ... až p na třetím obměň, místě atd. až např. po z) •biologické aktivity y jednotlivých látek rovněž označíme indexy V1 = 30 + ai,1X1,1 + ai,2X1,2 + + Ví,™ + a2,1X2,1 + a2,2X2,2 + - + Wn + ^Al + + - + 33Ap + + a x J,z J,Z v=a+a x +a x +......+ a x +a x +a x +...+a x +a x +a x +...+a x +.......+ y2 0 1,1 1,1 1,2 1,2 1,m 1,m 2,1 2,1 2,2 2,2 2,n 2,n 3,1 3,1 3,2 3,2 3,p 3,p a x J,z J,Z v=a+a x +a x +......+ a x +a x +a x +...+a x +a x +a x +...+a x +.......+ y3 0 1,1 1,1 1,2 1,2 1,m 1,m 2,1 2,1 2,2 2,2 2,n 2,n 3,1 3,1 3,2 3,2 3,p 3,p a x J,z J,Z v = a+ ax +ax +......+ a x +ax +ax + ... + a x +ax +ax + ... + a x +.......+ y\ 0 1,1 1,1 1,2 1,2 1,m 1,m 2,1 2,1 2,2 2,2 2,n 2,n 3,1 3,1 3,2 3,2 3,p 3,p a x J,z J,Z •skutečná soustava rovnic jednodušší, řada členů rovna 0 •přesto v daném stavu neřešitelná; proměnné x nejsou nezávislé => nutnost zjednodušujícícho předpokladu Řešení dle Freeho a Wilsona •předpoklad: součet příspěvků vlivů jednotlivých substituentů na biol. aktivitu ve sledované skupině pro každé obměňované místo je nulový => symetrické rovnice •b ... počet výskytů daného fragmentu; indexy mají stejný význam jako v předchozí soustavě rovnic ba + b a +.......+ ba =0 1,1 1,1 1,2 1,2 1,m 1,m ba + b a +.......+ ba =0 2,1 2,1 2,2 2,2 2,m 2,m b. a + b. a +.......+ b.a.. = 0 1,1 1,1 i,2 i,2 i,j i,j •pomocí symetrických rovnic se na každém obměňovaném místě jeden substituent nahradí (vyjádří) kombinací ostatních a soustava se stává řešitelnou Příklad - biologicky aktivní thiobenzanilidy R R1 1. thiobenzanilid H 2. 4'-chlor-4-methylthiobenzanilid CH3 3. 4-chlor-4'-methylthiobenzanilid Cl 4. 4,4'-dichlorthiobenzanilid 5. 4'-methylthiobenzanilid •v souboru sloučenin 2 obměňovaná místa Cl H R2 H Br CH, Cl * CH. Matice struktury Název thiobenzanilid 4'-brom-4-methylthiobenzanilid 4-chlor-4'-methylthiobenzanilid 4,4'-dichlorthiobenzanilid 4'-methylthiobenzanilid 1,1 (H) 1 0 0 0 1 2 bi,i 1,2 (CH,) 0 1 0 o 0 1 b1,2 1,3 (Cl) o 0 1 1 o 2 bi,s 2,1 (H) 1 0 0 0 0 1 b2,1 2,2 (CH,) 0 0 1 0 1 2 b2,2 2,3 (Cl) 0 1 0 1 o 2 b2,3 •první obměňovanou polohou pol. 4; v matici první 3 sloupce •druhou obměňovanou polohohou pol. 4'; v matici další 3 sloupce •pod tabulkou uvedeno, kolikrát se daný substituent vyskytuje (b) •pro každé substituované místo lze napsat symetrickou rovnici: 2a +a +2a =0 1,1 1,2 1,3 aoi + 2 a + 2 a = 0 2.1 2,2 2,3 •ze symetrických rovnic lze libovolný regresní koeficient v každé obměňované poloze vyjádřit pomocí ostatních, např. a = - 2 a - 2 a 1.2 1,1 1,3 a = - 2 a - 2 a 2,1 2,2 2,3 •vztahů se využije k úpravě (zjednodušení) matice proměnných x, které již lze považovat za nezávislé Název \i \s X2,2 X2,3 (H) (Cl) (CH,) (Cl) thiobenzanilid 1 0 -2 0 4'-chlor-4-methylthiobenzanilid -2 -2 0 1 4-chlor-4'-methylthobenzanilid 0 1 1 0 4,4'-dichlorthiobenzanilid 0 1 0 1 4'-methylthiobenzanilid 1 0 1 0 v matici na každém místě o jeden substituent méně •přepíšeme-li do soustavy rovnic, dostáváme y = a + a -2a -2a yi 0 1,1 2,2 2,3 v = a -2a -2a + a '2 0 1,1 1,3 2,3 v = a + a + a y3 0 1,3 2,2 v = a + a + a Jř4 0 1,3 2,3 v = a + a + a y5 0 1,1 2,1 •soustava 5 rovnic o 5 neznámých je řešitelná (metoda dosazovací, sčítací, maticový počet, multilineární regrese na počítači) •za y dosadíme příslušné hodnoty biologické aktivity jednotlivých látek •získáme hodnoty regresních koeficientů a„, a_ a_ „, a . a„„ 3 3 0' 1,1' 1,3' 2,2' 2,3 •hodnoty ai 2 a a22 se získají dosazením do předchozích rovnic, ve kterých byly vyjádřeny pomocí ostatních regresních koeficientů •aQ zde vyjadřuje průměrnou hodnotu studované aktivity souboru, ostatní hodnoty a s příslušnými indexy udávají, jakým způsobem obměňované substituenty v jednotlivých polohách tuto aktivitu ovlivňují (tj. jsou to kladné nebo záporné odchylky od průměrné aktivity) • pomocí získaných hodnot regresních koeficientů lze výpočtem odhadnout aktivity látek s kombinacemi použitých substituentů, které ve výchozí matici nebyly zahrnuty (tj. 4-chlorthiobenzanilid, 4'-chlorthiobenzanilid a 4,4'-dimethylthiobenzanilid) •Aktivity: antituberkulotické, vyjádřené jako log MIC proti Mycobacterium kansasii Řešení dle Fujity a Bana •jednodušší: místo symetrických rovnic se vynechá jeden ze sloupců pro každé ze substitucí obměňovaných míst •vhodné vynechat sloupce odpovídající substituci dané polohy vodíkem => zavedení nulové hodnoty pro každé substituční místo; hodnota odpovídá substituentu ve vynechaném sloupci •příspěvky odpovídající vlivům substituentu na biol. aktivitu jsou pak relativní vůči této hodnotě => výhodné mít v souboru nesubst. látku Zjednodušená matice struktury Název X X X X 1,2 1,3 2,2 2,3 (CH,) (Cl) (CH,) (Cl) thiobenzanilid 0 0 0 0 4'-brom-4-methylthiobenzanilid 1 0 0 1 4-chlor-4'-methylthiobenzanilid 0 1 1 0 4,4'-dichlorthiobenzanilid 0 1 0 1 4'-methylthiobenzanilid 0 0 1 0 •soustava rovnic k řešení pak vypadá následovně: v = a '1 o •na výstupu hodnota aQ odpovídá aktivitě nesubst. látky (zde thiobenzanilid) •ostatní hodnoty a odpovídají příspěvkům jednotlivých substituentů k aktivitě Empirické metody QSAR •používány s výhodou tam, kde nelze jednoduše zjistit matem, popis funkce A = f(struktura) •hledají pouze extrémy (maxima a/nebo minima) dané funkce; matem, popis přitom zůstává "černou skříňkou" •při jejich aplikaci syntetik usměrňuje volbu látky k syntéze na základě biol. hodnocení látky předchozí Fibonacci - Leonardo Pisánský (cca 1175 - 1250) •přinesl do Evropy arabský číselný systém •Fibonacciho posloupnost: nekonečná posloupnost přirozených čísel, začínající 0, 1, 1, 2, 3, 5, 8, 13, 21, ... (čísla nacházející se ve Fibonacciho posloupnosti jsou někdy nazývána Fibonacciho čísla), kde každé číslo je součtem dvou předchozích, •rekurzivní zápis: 0, pro n-0; F(n) -\\pron-1; F(n -1) + F(n - 2) jinak. •explicitní vyjádření - vztah pro n-tý člen: Optimalizace podle jednoho strukturního parametru - Fibonacciho optimalizace Látky se seřadí podle rostoucí hodnoty strukturního parametru, u něhož se přepokládá, že aktivitu významně ovlivňuje. Počet látek musí odpovídat počtu bodů v některém z Fibonacciho intervalů (viz tab. 1), není-li tomu tak, některá z krajních látek, u nichž je malá pravděpodobnost, že by byly nejúčinnější, se vypustí nebo se naopak přidá fiktivní krajní látka. K syntéze se vyberou látky, které mají v intervalu pořadí uvedené ve druhém sloupci tabulky. Provede se biologické hodnocení těchto látek a podle jeho výsledků se část intervalu od jednoho z krajních bodů směrem k méně aktivní sloučenině vyloučí. Výsledný soubor je dalším Fibonacciho intervalem. Výběr se opakuje tak dlouho, dokud se nedospěje k nejúčinnější sloučenině. Metoda umožňuje podstatně snížit počet syntetizovaných a testovaných látek, např. namísto 589 látek, které by bylo nutno připravit a otestovat, abychom mezi nimi našli nejúčinnější, stačí připravit a stanovit biol. aktivitu u pouhých 13 sloučenin (viz sloupec C tab. 1). Tab. 1. Fibonacciho optimalizace Legenda: A ... počet látek příslušného Fibonacciho intervalu B ... pořadí látek určených pro syntézu a testování v rámci intervalu C ... počet látek potřebný k dosažení optimalizace A B C A B C A B C 2 la2 2 20 8 a 13 6 143 55 a 89 10 4 2a3 3 33 13 a21 7 222 89 a 144 11 7 3a5 4 54 21 a 34 8 366 144 a 233 12 12 5a8 5 88 34 a 55 9 589 233 a 377 13 Optimalizace podle více strukturních parametrů Simplexová metoda Každou látku lze charakterizovat jako bod v n-rozměrném prostoru, ve kterém biologická aktivita je jednou souřadnicí a zbývajících n - 1 souřadnic představují fyzikální a fyzikálně-chemické vlastnosti sloučeniny, o nichž se předpokládá, že aktivitu ovlivňují. Pokud pracujeme v (klasickém) trojrozměrném prostoru, tj. optimalizujeme-li pouze dva parametry, můžeme optimalizaci provést i graficky na milimetrovém papíře. Pracujeme vlastně v průmětu do roviny vlastností. K (syntéze a) hodnocení se vyberou 3 látky, které v rovině souřadnic vlastností nejsou od sebe příliš vzdálené s jejichž souřadnice v této rovině tvoří trojúhelník, nejlépe rovnostranný. Porovnáme aktivity těchto 3 sloučenin a z bodu, patřícího látce o nejnižší aktivitě, vedeme polopřímku středem spojnice dvou bodů s vyššími aktivitami (alternativně bodem, který vznikne rozdělením této spojnice v převráceném poměru aktivit) a na ní nalezneme bod stejně vzdálený, ale opačně orientovaný, než bod s nejnižší aktivitou. Pokud tento bod neodpovídá žádné sloučenině, použijeme k hodnocení látku nejbližší. Z tohoto bodu a dvou předchozích vytvoříme další trojúhelník, který podrobíme stejnému optimalizačnímu pochodu. Postup opakujeme, dokud nacházíme látky o rostoucí aktivitě. Jakmile aktivita začne klesat, můžeme látku s dosaženou nejvyšší aktivitou považovat za optimální. Simplexová metoda Optimalizace podle více strukturních parametrů Optimalizační schémata •postup racionálních myšlenkových pochodů farm. chemika •zohledňují hydrofobní, elektronové a sterické parametry •nejsou univerzální; pro určitý typ obměn na určité struktuře bývá nutné sestavit nové Schéma obměn substituentů na fenylu (Topliss 1972) 3-CI mene 4-OCH, <=" 4-CI 3-CF3) 4-CI; 3-CF3, 4-N02 3,4 CL I 3-CH3, 4-N(CH3)2 stejné 4-N(CH3)2 4-CH, Ny^nebo 4-CF3, (4-Br, 4-I), 2,4-CI2 »4-C(CH3)3; 3,4-(CH3)2 4-NH2; 4-OH; 3-CH3, 4-OCH3 S-NÍCH^; 3-NH2 <= 3-C! 3-CF3; 3-Br; 3-I; 3,5-CI, 3-CH3 - 2-CI; 2-CH3; 2-OCH3 Komentář ke schématu obměn substituentu na fenylu Při optimalizaci podle uvedeného schématu postupujeme tak, že nejprve připravíme nesubstituovanou sloučeninu a její 4-chlorderivát. Substituce chlorem sníží elektronovou hustotu v poloze 1 a zároveň zvýší lipofilitu (4-CI: a = 0,23; 71 = 0,71); je-li 4-chlorderivát aktivnější, lze lipofilitu a elektronakceptorové vlastnosti zvýšit další substitucí chlorem. Je-li 4-chlorderivát méně aktivní, můžeme předpokládat, že pokles elektronové hustoty ovlivnil aktivitu negativně, a připravíme 4-methoxyderivát, který má lipofilitu prakticky stejnou jako nesubst. látka, el. hustota v pol. 1 je však vyšší (4-OCH3: a = -0,27,71 = 0,02). Není-li zřetelný rozdíl mezi aktivitou nesubst. látky a 4-chlorderivátu; můžeme předpokládat, že vliv el. hustoty a lipofility působí proti sobě a zkusíme připravit 4-methylderivát (4-CH3: a = -0,17,71 = 0,56). Je-li u všech látek substituovaných v pol. 4 aktivita nižší něž u nesubstituované látky, je zřejmé, že substituce v poloze 4 je stericky nevýhodná, a budeme připravovat látky substituované v polohách 2 a 3. V jednotlivých větvích tohoto schématu je možné postupovat dále, dokud nedojdeme ke sloučenině s maximální aktivitou.