Аналіз зв’язку структура-властивість на основі молекулярних ґраткових моделей

Размещено на http://allbest.ru

Національна Академія Наук України

Фізико-хімічний Інститут ім. О.В. Богатського

Артеменко Анатолій Григорович

УДК 547.1:541.6

Автореферат дисертації

на здобуття наукового ступеня кандидата хімічних наук

Аналіз зв'язку структура-властивість на основі Молекулярних ґраткових моделей

02.00.03 - органічна хімія

Одеса-2000

Дисертацією є рукопис.

Робота виконана на кафедрi органiчної хiмiї Одеського національного унiверситету iм. I.I.Мечникова Міністерства освіти і науки України та у Фізико-хімічному інституті ім. О.В.Богатського НАН України.

Науковий керівник - кандидат хімічних наук, доцент Кузьмін Віктор Євгенович, Фізико-хімічний інститут ім. О.В. Богатського НАН України, старший науковий співробітник

Офіційні опоненти:доктор хімічних наук, професор Ільченко Андрій Якович, Інститут органічної хімії НАН України, провідний науковий співробітник - кандидат хімічних наук Сидельнікова Тетяна Андріївна, Одеський державний медичний університет Міністерства охорони здоров'я України, старший викладач кафедри медичної хімії

Провідна установа Інститут біоорганічної хімії та нафтохімії НАН України, м. Київ.

Захист відбудеться “ 19 ” січня 2001 р. о 14 годині на засіданні спеціалізованої вченої ради Д 41.219.02 в Фізико-хімічному Інституті ім. О.В. Богатського НАН України, 65080, Одеса, вул. Люстдорфська дорога, 86.

З дисертацією можна ознайомитись в бібліотеці Фізико-хімічного інституту ім. О.В. Богатського НАН України, 65080, Одеса, вул. Люстдорфська дорога, 86.

Автореферат розісланий “ 14 ” грудня 2000 р.

Вчений секретар

спеціалізованої вченої ради,

кандидат хімічних наук,

старший науковий співробітник Л.О. Литвинова

ЗАГАЛЬНА ХАРАКТЕРИСТИКА РОБОТИ

Актуальність теми. У наш час досить широкого розвитку набули дослідження в області кількісних співвідношень між структурою та властивостями органічних сполук (QSAR/QSPR QSAR(QSPR) - quantitative structure-activity(property) relationship (кількісний зв'язок “структура - активність (властивість)”).).^.Зростаючий інтерес до цих проблем пояснюється насамперед тим, що виявлення подібних зв'язків сприяє кращому розумінню механізмів функціонування нових матеріалів, реагентів і препаратів. Використання комп'ютерних технологій дозволяє на базі виявлених закономірностей провести не тільки попереднє відсівання (скринінг), але й молекулярний дизайн органічних речовин, що виявляють ту або іншу властивість. Все це дозволяє істотно знизити необхідні для розробки нового матеріалу, реагенту, препарату колосальні затрати часу та ресурсів.

Існує досить велика кількість методів і моделей для вирішення задач “структура-властивості” від самих простих до вельми складних. Однак майже всі вони використовують однопланову структурну інформацію. Найбільш прості моделі засновані на використанні констант замісників, топологічних індексів молекул (програми ЭММА, CODESSA), інформації про кількість структурних фрагментів того або іншого типу (метод HQSAR Hologram QSAR - метод QSAR, заснований на знятті “відбитків” (молекулярных голограмм) різних структурних фрагментів з молекулярної топологічної структури.). У більш складних 3D (тривимірних) QSAR методах (наприклад, CoMFA Comparative Molecular Field Analysis (порівняльний аналіз молекулярних полів) - один з найбільш популярних методів 3D QSAR.) основну роль відіграють характеристики, що відображають особливості міжмолекулярної взаємодії досліджуваних сполук. На жаль, більшість відомих моделей QSAR/QSPR ігнорують стереохімічні особливості молекул. Зрозуміло, що використання таких методів не завжди дозволяє досягнути надійного прогнозування цільової властивості.

У зв'язку з цим, актуальною є розробка методів, що використовують всі типи структурних параметрів (інтегральних, локальних, польових) в комплексі, що дає можливість отримати практично вичерпну інформацію про структуру молекули. Звісно, при цьому необхідно використовувати спеціальні методи для встановлення зв'язку структура-властивість, здатні обробляти великі масиви структурної інформації без істотної втрати надійності моделі.

Зв'язок роботи з науковими програмами, планами, темами. Робота виконувалась у вiдповiдностi з темою державного бюджету “Розвиток концепції структурно-хімічної подібності для вирішення задач “структура-властивість” з метою спрямованого синтезу нових речовин та матеріалів - мезоморфогенів, альтернативних хладонів, комплексонів, біологічно активних сполук” (№ державної реєстрації 0195 V 015695 1995-1998рр.), а також в рамках міжнародних проектів INTAS “Розвиток та застосування алгоритму просторового навчання штучних нейронних мереж в дослідженнях 3D QSAR” (грант INTAS-UA 95-0060) та “Антиракові та антивірусні властивості макроциклічних піридинофанів: зв'язок структура - активність та дизайн лікарських препаратів” (грант INTAS-UA 97-31528).

Мета і задачі дослідження. Розробка і реалізація підходу, що забезпечує генерацію однорідної системи структурних параметрів, які відображають локальні, інтегральні і польові характеристики молекулярної структури. Вивчити зв'язок між структурою і:

біологічною активністю стероїдів і 1,4-бенздіазепінів;

здатністю утворення комплексів з переносом заряду монозаміщених бензолів;

люмінесцентними властивостями -дикетонатів рідкоземельних елементів (РЗЕ);

здатністю (-)-ментонів створювати спіральне упорядкування в нематичній мезофазі.

Виявити структурні характеристики молекул, сприяючі і перешкоджаючі вияву досліджуваних властивостей.

Наукова новизна одержаних результатів. Запропонована інваріантна система кодування просторової структури молекули на основі побудови ламаної спіральної кривої в рамках ґраткової моделі. Показано, що молекулярний код може бути використаний для оцінки хіральності та структурної подібності молекул.

Розроблена система опису локальних характеристик структури, що враховує в тій або іншій області простору молекули не тільки наявність/ відсутність атому або групи атомів, але й їх властивості. Для функціональних груп, що є потенційними фармакофорами, запропоновані характеристики, що відображають доступність захаращуваність, а також взаємне положення та орієнтацію цих груп.

Як польові параметри молекулярної структури запропоновані імовірності заселення різних областей навколомолекулярного простору атомами (групами), що моделюють функціональну (біологічну) мішень (псевдорецептор).

На основі гармонічного аналізу локальних і польових характеристик молекул запропоновані нові структурні інваріанти, що відображають характер розподілу відповідних властивостей атомів, фрагментів або молекулярних полів відносно просторової структури молекул.

Практичне значення одержаних результатів. Розроблено комплекс комп'ютерних програм “MolLatt”, що реалізує послідовність процедур для вирішення різноманітних задач QSAR/QSPR:

побудова ґраткових моделей молекул і визначення їх молекулярного коду;

обчислення параметрів структурної подібності і хіральності різних молекул;

розрахунок більше ніж 20 тис. локальних, інтегральних та польових структурних характеристик молекулярної структури;

аналіз зв'язку структура-властивості методами множинного регресійного аналізу, регресії часткових найменших квадратів (PLS), тренд-вектору;

візуалізація всіх етапів обробки молекулярної структури і вирішення задачі “структура-властивості”.

Вказаний комплекс програм дозволяє вирішувати широкий спектр задач “структура-властивості”.

Особистий внесок здобувача. Ідеологія ґраткового представлення молекул, розрахунку структурних параметрів та моделювання псевдорецептору розроблена автором роботи спільно з науковим керівником. Особистий внесок: реалізація методу за допомогою комплексу комп'ютерних програм і вирішення конкретних задач QSAR.

Апробація результатів дисертації. Основні результати роботи докладалися та обговорювалися на I, II, III Міських семінарах “Застосування обчислювальної техніки і математичного моделювання в прикладних наукових дослідженнях" (Одеса, 1994-1996), 7-му Міжнародному симпозіумі по хіральній дискримінації (Єрусалим, 1995), 11-й Міжнародній конференції по зв'язку між математичними, хімічними і комп'ютерними науками “MATH/CHEM/COMP` 96" (Дубровник, Хорватія, 1996), 7-й Міжнародній конференції по математичній хімії (Жирона, Іспанія, 1997), I-й Всеросійській конференції “Молекулярне моделювання” (Москва, 1998), на двох симпозіумах по міжмолекулярній взаємодії і конформації молекул (Казань, 1997 р. і Твер, 1999 р.), I-й конференції молодих вчених-хіміків Південного регіону України (Одеса, 1998), Міжнародній конференції “Успіхи сучасної органічної та біоорганічної хімії”(Одеса, 1999), 7-му Міжнародному симпозіумі по молекулярним аспектам хіміотерапії (Гданськ, Польща, 1999), II з'їзді онкологів країн СНД (Київ, 2000).

Публікації. За матеріалами дисертації опубліковано 5 статей і 11 тез доповідей.

Структура та обсяг роботи. Дисертаційна робота складається із вступу, чотирьох розділів, висновків, списку використаних джерел і трьох додатків. Робота викладена на 141 сторінці, містить 34 малюнки, 25 таблиць, список використаних джерел із 112 найменувань.

ОСНОВНИЙ ЗМІСТ РОБОТИ

1. Розвиток моделей 3D QSAR

Проведений аналіз існуючих методів вирішення задач 3D QSAR/QSPR. Детально розглянуті типи структурних параметрів, що використовуються для ефективного моделювання біологічної активності та інших властивостей молекул, зумовлених їх структурою.

З особливою увагою розглянуті методи моделювання рецепторів, що дозволяють здійснювати не тільки прогноз біологічної активності лікарських препаратів, але й формувати гіпотези про механізм їх біологічної дії. Сформульовані основні проблеми, які виникають при моделюванні “біологічної мішені”, обговорюються можливі шляхи їх вирішення. Проаналізовані переваги і нестачі ряду найбільш відомих методів 3D QSAR.

У кінці розділу зазначено, що серйозним недоліком більшості методів є використання однопланової структурної інформації та ігнорування стереохімічних особливостей молекул. Зроблений висновок про можливість універсалізації всіх типів структурних параметрів (інтегральних, локальних, польових) в рамках ґраткового моделювання молекул.

2. ґраткові моделі молекулярної структури

Розділ присвячений викладенню ґраткового підходу до моделювання молекулярної структури, опису інтегральних, локальних і польових структурних параметрів, що розраховуються в його рамках. Суттю методу є розміщення молекули з фіксованою Просторову структуру молекул одержують через експеримент (наприклад, РСА) або шляхом розрахунку методами молекулярної механіки, квантової хімії і т.п. просторовою структурою в ґратках, складених з кубічних комірок (розміром 1.8-2.0 A Розмір комірки є настроювальним параметром моделі і може варіюватися.) (рис. 1).



Размещено на http://allbest.ru

Рис. 1. Приклад молекули бутану, розміщеної в кубічних ґратках

Рис. 2. Схема побудови правої кривої спіральної лінії в просторових кубічних ґратках

У рамках підходу в залежності від задачі, що вирішувалась, можливе суміщення молекул досліджуваного ряду по центру маси та головним осям інерції або по групі атомів, яка лежить в основі спільного жорсткого структурного фрагменту. Відомості про кожну комірку отриманих тривимірних ґраток, а саме, наявність в ній того чи іншого атому (фрагменту), значення тієї чи іншої фізико-хімічної властивості атому (фрагменту) або польових характеристик навколомолекулярного простору, являють собою джерело структурної інформації і використовуються як структурні дескриптори для вирішення задач “структура-властивості”.

Далі в просторі ґраток будується ламана лінія (рис. 2), що бере початок в центрі маси молекули і проходить через центри всіх комірок ґраток, так щоб “спіральність” цієї лінії була постійною (всі повороти тільки “ліворуч” або тільки “праворуч”). Фактично ця ламана являє собою складну криву, що складається з вкладених один в одний коаксіальних фрагментів спіралі. Ламана спіральна лінія використовується для побудови молекулярного коду (МК) кожної із досліджуваних сполук.

ПМК: 2.35.2.9.1.6.0.1.2.17

Рис. 3. Приклад побудови молекулярного коду для молекули СHFClBr

Структура коду досить проста: він складається з послідовності чисел b_i, що є характеристиками атомів Звичайно використовуються характеристики атомів: порядковий номер, маса, рефракція, поляризовність, ван-дер-ваальсовий радіус, ліпофільність, заряд. Можливе також використання характеристик фрагментів типу CH₂, CH₃, NH₂ і т.п. та чисел a_i, що визначають відстань між атомами “i” та “i+1” вдовж вказаної вище ламаної, виражену цілими числами (кількість пустих комірок) (рис. 3). Такий код містить в стислому вигляді всю інформацію про форму та стереохімію досліджуваної сполуки. Для кожної молекули можна будувати як лівий (ЛМК), так і правий код (ПМК), в залежності від типу спіральної лінії. На основі МК Для розрахунків використовується ПМК. з'являється можливість оцінки структурної подібності / відмінності (СП/СВ) сполук. Для цієї мети запишемо ПМК для двох довільних молекул:

; ,

де n₁n₂. Тут більш коротка послідовність комірок доповнюється нулями, так щоб сумарна кількість елементів в обох кодах була однакова. СВ молекул М₁ і М₂ може бути розрахована, як Евклідова відстань R в просторі параметрів (а_i, а'_i) і (b_i, b'_i), відповідно. Якщо не враховувати відмінності між природою величин a_i та b_i, і записати ПМК двох структур у вигляді послідовностей ПМК(M₁)=X₁X_2...X_i...X_N, ПМК(M₂) =Y₁Y_2...Y_i...Y_N, де N=2n₂, то з'являється можливість оцінити СП/СВ одним параметром, якщо використати коефіцієнт парної лінійної кореляції (R) або відстань Камберра:

Крім параметрів структурної подібності, що розраховуються на основі МК, в рамках моделі розраховуються параметри, що відображають СП/СВ фізико-хімічних властивостей (ліпофільності, заряду, електростатичного потенціалу та інших) атомів у вузлах ґраток, що накладаються.

Параметри хіральності можна розраховувати на основі оцінки структурної подібності енантиомерів.

Було досліджено поведінку розроблених характеристик структурної подібності та хіральності для різних серій модельних сполук зі структурою, що систематично змінюється (довжиною зв'язку, масою замісника, торсійним кутом, положенням замісника в ланцюгу). Було показано, що запропоновані характеристики досить коректно кількісно оцінюють рівень хіральності та структурної подібності молекул і можуть бути використані для вирішення задач QSAR/QSPR.

Приведено докладний опис структурних характеристик, що розраховуються в рамках ґраткового підходу.

Інтегральні: Головні моменти інерції, дипольний момент, ліпофільність, молярна рефракція, параметри структурної подібності і хіральності, що розраховуються на основі МК, гармонічні параметри молекулярної структури (див. нижче).

Локальні. Наявність / відсутність в комірці атому певного типу, значення фізико-хімічної властивості атому (фрагменту); характеристики потенційних фармакофорів (в межах ґраткової моделі розраховується три типи характеристик фармакофорів):

параметр доступності / захаращуваності кожного із фармакофорів вздовж шести напрямків координатних осей, відповідно до головних осей інерції;

відстані між центрами пар фармакофорних груп;

характеристики трійок фармакофорних груп:

площа трикутника, утвореного центрами фармакофорних груп;

змішаний добуток векторів, направлених від центра маси до центрів фармакофорних груп.

Польові (характеристики навколомолекулярного простору):

електростатичний потенціал (всі зарядові характеристики розраховуються по методу вирівнювання електронеґативностей [Кузьмін В.Є., Берестецька О.Л., 1983];

потенціал поля ліпофільності [Croizet F., 1990];

імовірність заселення пустої комірки різними атомами (“щупами”), а також імовірність того, що комірка залишається вакантною (розраховуються на основі енергії взаємодії Використане універсальне силове поле UFF [Rappe A.K., 1992]. між молекулою і відповідним атомом (“щупом") в комірці ґраток);

наявність / відсутність передбачуваного донору (або акцептору) водневого зв'язку в комірці. Мається на увазі, що водневий зв'язок може бути утворений між молекулою, що аналізується, і функціональною мішенню.

На основі ґраткової моделі з'являється можливість ефективного використання всього комплексу параметрів одночасно. Гармонічний аналіз функцій розподілу локальних і польових характеристик (Р) вздовж спіральної кривої дозволяє виділити інтегральні структурні параметри молекул: від окремих фрагментів (високочастотні гармоніки) до всієї молекули (низькочастотні гармоніки):

,

де k-номер гармоніки, q_k -амплітуда, (k -фазовий кут. Як параметри для вирішення задач “стуктура-властивості” використовуються значення амплітуд (q_k).

Структурні параметри, отримані внаслідок гармонічного аналізу локальних і польових характеристик молекул, хоч і важко інтерпретувати, але, проте, очевидно, що вони відображають широкий спектр розподілів відповідних властивостей атомів, фрагментів або молекулярних полів відносно просторової структури молекул. Саме тому можлива побудова рівнянь зв'язку структура-властивість тільки на їх основі, що продемонстровано у третьому розділі роботи.

Слід підкреслити, що запропонована система параметрів дає майже вичерпну (а часто надмірну) інформацію про структуру молекули. Саме через це для встановлення зв'язку структура-властивість необхідно використовувати спеціальні методи, здатні обробляти великі масиви структурної інформації без істотної втрати надійності моделі. Все це дозволяє використовувати структурну інформацію для дослідження практично будь-яких властивостей сполук. Крім того, аналіз локальних і польових характеристик молекул дозволяє конструювати функціональну (біологічну) мішень (“псевдорецептор”) як “зліпок” з конгломерату активних і неактивних сполук навчальної вибірки.

3. Вирішення задач QSAR/QSPR за допомогою ґраткової моделі

Описано застосування розробленого підходу для вирішення різних задач “структура-властивості”, зокрема для моделювання біологічної активності стероїдів і 1,4-бенздіазепінів, здатності утворення комплексів з перенесенням заряду монозаміщених бензолів, люмінесцентних властивостей -дикетонатів РЗЕ, здатності (-)-ментонів створювати спіральне упорядкування в нематичній мезофазі.

У всіх задачах досліджувані молекули оптимізувалися методами молекулярної механіки Просторова структура молекул будувалась за допомогою програми HyperChem, безкоштовна “trial version” якої отримана з http://www.hypercube.com.. Пошук кількісних співвідношень структура-властивість проводився за допомогою методів множинного лінійного регресійного аналізу (МЛР), модернізованого і класичного тренд-вектору Процедура тренд-вектору дає можливість не конкретизувати вигляд відповідної залежності і використовувати велику кількість структурних параметрів. Платою за це є обмеження можливості вирішення зворотної задачі (прогнозування величини властивості), яка може бути вирішена тільки в ранговій шкалі. (МТВ/КТВ) [Витюк Н.В., Кузьмін В.Є., 1994], а також методу регресії часткових найменших квадратів (PLS), що широко використовується для побудови рівнянь QSAR/QSPR при великій кількості дескрипторів. У всіх задачах розглядалося декілька моделей для однієї властивості, що дало можливість більш повно і різносторонньо інтерпретувати вплив структури на досліджувані властивості.

3.1 Вивчення спорідненості стероїдів до кортікостероїд- та тестостерон- зв'язуючого глобуліну

Аналізувався вплив структури на спорідненість молекул стероїдів (31 сполука) до кортикостероїд- і тестостерон - зв'язуючого глобуліну (КЗГ/ТЗГ), згідно з даними роботи [Cramer R.D., 1988]. Об'єкти дослідження були вибрані не випадково. Дана серія сполук є добре вивченою і часто використовується для тестування працездатності розроблених методів QSAR:

R₁=-OH, =O; R₂=-H, -OH, =O, -COMe; R₃=-H, -OH, =O, -COCH₂OH, -COCH₃; R₄=H, -OH; R₅=-OH, =O; R₆=-H.

Внаслідок аналізу було отримано 4 приблизно рівноцінних моделі (табл. 1).

Таблиця 1

Характеристики моделей, отриманих в результаті вирішення задачі “структура-спорідненість до КЗГ” для молекул стероїдів.

Модель	Структурні параметри	Кількість параметрів	Метод аналізу	Характеристики моделі
I	гармонічні	5	МТВ	R=0.972; CVR CVR - коефіцієнт кореляції, обчислений при ковзному контролі. Ковзний контроль дозволяє оцінити досліджувану властивість молекули “в умовах екзамену”, тобто коли розглядувана молекула видаляється з навчальної вибірки і регресійна модель будується без неї. Потім аналізована властивість прогнозується для цієї молекули. Вказана процедура виконується для всіх досліджуваних сполук.=0.936
II	локальні й інтегральні	7	КТВ	R=0.974; CVR=0.950
III	польові	5	КТВ	R=0.973; CVR=0.924
IV	всі види параметрів	17	PLS	R=0.980; CVR=0.972;А Кількість латентних (прихованих) змінних.=1

Було виявлено, що найбільш важливим фактором для вияву спорідненості стероїдів до КЗГ є наявність sp²-гібридного атому оксигену, пов'язаного з атомом карбону С₂. Наявність ОН- групи в цьому ж положенні дає протилежний ефект. Потрібно зазначити також неґативний вплив на зв'язування з КЗГ груп =О, -ОН (в псевдоекваторіальному положенні), пов'язаних з карбоном С₁₅ і позитивний вплив гідроксильної групи в екваторіальному положенні біля атому С₁₁. Вплив інших параметрів не такий значний.

Отримана в цій задачі краща модель (IV) за статистичними показниками (R=0.980, CVR=0.972) значно переважає модель, отриману методом CoMFA (R=0.947, CVR=0.814) в роботі [Cramer R.D., 1988]

Внаслідок аналізу впливу структури на зв'язування стероїдів з ТЗГ також були одержані 4 моделі (табл. 2).

Таблиця 2

Характеристики моделей, отриманих в результаті вирішення задачі “структура-спорідненість до ТЗГ” для молекул стероїдів.

Модель	Структурні параметри	Кількість параметрів	Метод аналізу	Характеристики моделі
V	гармонічні	8	КТВ	R=0.993; CVR=0.908
VI	локальні й інтегральні	5	КТВ	R=0.970; CVR=0.949
VII	польові	6	КТВ	R=0.989; CVR=0.940
VIII	всі види параметрів	19	PLS	R=0.986; CVR=0.916; А=2

Аналіз структурних параметрів показав важливість для вияву активності характеристик структурної подібності за зарядами. Як і чекалось, чим більш подібна молекула за зарядами до однієї з активних молекул, тим більша активність їй присутня. Із локальних характеристик сприяє вияву активності наявність гіроксильної групи в псевдоекваторіальному положенні атома С₁₅. Цікаво зазначити, що в попередньому випадку (зв'язування КЗГ) вплив цього замісника прямо протилежний (див. вище). Наявність ОН- групи в аксіальному положенні атому карбону С₂ трохи зменшує активність.

Аналіз гармонічних параметрів показав, що в обох випадках (КЗГ і ТЗГ) характер взаємного просторового розташування атомів оксигену стероїдних молекул істотно впливає на їх спорідненість до відповідних рецепторів; можна зробити припущення про участь цих атомів в утворенні водневих зв'язків з білковими молекулами.

При моделюванні афінітету стероїдів до ТЗГ отримана нами краща модель (VIII) за статистичними показниками (R=0.986, CVR=0.916) значно переважає модель, отриману методом CoMFA (R= 0.934, CVR= 0.745) в роботі [Cramer R.D., 1988].

3.2 Дослідження константи утворення комплексів з переносом заряду для ряду похідних бензолу

Враховуючи, що в ряді випадків біологічна активність зумовлена утворенням міжмолекулярних комплексів з переносом заряду (див. [Livingstone D.J., 1992]), був досліджений вплив структури монозаміщених бензолів (38 молекул) на константу утворення комплексів з переносом заряду з 1,3,5-тринітробензолом згідно з даними роботи [Livingstone D.J., 1992]:

	R= -H, -alc, -N(CH3)3, -Ph, -CHO, -CO-alc, -O-alc, -SCH3, -CF3, -CN, -Hal, -NO2,-NH2, -NH-alc, -CH2CN, -N(alc)2, -CH2OH, -CO-N(alc)2, -SO2-N(alc)2, -O-Ph, -O-CH2-Ph ; alc=CH3, C2H5, C3H7

Внаслідок статистичного аналізу було отримано 5 моделей (табл.3).

Таблиця 3

Характеристики моделей, отриманих в результаті вирішення задачі структура-властивість для молекул монозаміщених бензолів.

Модель	Структурні параметри	Кількість параметрів	Метод аналізу	Характеристики моделі
IX	гармонічні	8	КТВ	R=0.964; CVR=0.959
X	локальні та інтегральні	9	КТВ	R=0.961; CVR=0.843
XI	польові	7	КТВ	R=0.979; CVR=0.867
XII	інтегральні	3	МЛР	R=0.967;F F - критерій Фішера.=151;SE SE - стандартна помилка передбачення.=0.12
XIII	всі види параметрів	23	PLS	R=0.986; CVR=0.943; А=2

Аналіз відібраних параметрів показав, що процесу комплексоутворення, як і чекалося, перешкоджають електронеґативні замісники, що зменшують електронодонорні властивості ароматичного кільця. Замісники, що легко поляризуються, сприяють реакційній здатності досліджуваних молекул. Інтегральні параметри також показують (модель XII), що, чим менше ліпофільність і більше поляризовність всієї молекули, тим вище її здатність до утворення комплексів з переносом заряду. Виявлено, що, чим ближче досліджувані сполуки до активних молекул з точки зору полярності (електронеґативності атомів) структури, тим вище їх реакційна здатність.

Результати аналізу цілком логічно відображають фізико-хімічні особливості реакції утворення комплексів з переносом заряду і узгоджуються з виводами [Livingstone D.J., 1992], однак, отримана в теперішній роботі краща модель (XIII) за статистичними показниками (R=0.986) трохи переважає модель, отриману цими авторами (R= 0.975) на основі квантово-механічних розрахунків.

3.3 Дослідження психотропної активності 1,4-бенздіазепінів

На наступному етапі був досліджений ряд 1,4-бенздіазепінів, добре вивчених психотропних препаратів:

	X=-H, -Me Y=-H, -Me,-Hal, -NO2, -CF3, -S-Me, -SO-Me, -SO2Me, -CN, -N(Me)2, -S-Et Z=-H, -Me, -Hal, -CF3, -OMe, -NO2

Незважаючи на велику кількість інформації про біологічну активність похідних 1,4-бенздіазепінів, в роботі цілеспрямовано була використана обмежена вибірка (38 молекул), для якої автори статті [Miyashita Y., 1981] не змогли побудувати задовільну модель QSAR.

Внаслідок нашого вирішення задачі “структура - психотропна активність” були отримані 3 приблизно рівноцінних моделі (табл.4).

Таблиця 4

Характеристики моделей, отриманих в результаті вирішення задачі структура-психотропна активність для ряду похідних 1,4-бенздіазепінів

Модель	Структурні параметри	Кількість параметрів	Метод аналізу	Характеристики моделі
XIV	гармонічні	11	КТВ	R=0.966; CVR=0.931
XV	локальні та інтегральні	10	КТВ	R=0.933; CVR=0.731
XVI	польові та інтегральні	14	КТВ	R=0.945; CVR=0.766

Аналіз структурних параметрів показав, що вияву психотропної активності сприяє наявність електроноакцепторних угрупувань (наприклад, -NO₂, CF₃) в положенні замісника Y. Наявність вуглеводневих радикалів в цьому положенні чинить протилежну дію.

Замісник Z також впливає на психотропну активність: її збільшенню сприяють невеликі за розмірами групи, які несуть частковий неґативний заряд, наприклад, такі як, F. Цілком ймовірно, що групи >С=O та >N=C- беруть участь у взаємодії з бенздіазепіновим рецептором, оскільки їх захаращуваність впливає на психотропну активність.

За результатами моделювання можна передбачити наявність високої психотропної активності для сполуки із замісниками СН₃, СF₃ і F в положеннях X, Y і Z, відповідно, якої немає у вибірці і дані про його активність невідомі.

Таким чином, краща модель (XIV) за статистичними показниками (R=0.966) набагато переважає модель, отриману авторами оригінальної статті (R<0.6).

3.4 Дослідження інгібування бенздіазепінового рецептору

Як випливає з опису запропонованого методу, він цілком може бути застосовний для дослідження компліментарності взаємодії ліганд - рецептор.

У зв'язку з цим, була зроблена спроба в рамках ґраткової моделі дослідити вплив структури 1,4-бенздіазепін-2-онів (57 молекул, включаючи ряд комерційних препаратів) на константу інгібування бензодіазепінового сайту (BzR) сімейства GABA_A рецепторів (LogIC50 IC₅₀ в нмоль/л), за даними роботи [Winkler D.A., 1998]:

Размещено на http://allbest.ru

R1 =-H, -CH3, -CH2CF3, -CH2OCH3, -C(CH3)3, -(CH2)2OH, -C(Me)3; R2'=-H, -Hal, -NO2, -CF3 ; R3 =-H, -CH3, -OH, -Cl ; R6'=-H, -F, -Cl; R7 =-H, -Me, -Hal, -CF3, -NO2, -N3, -NH2, NHOH, -CN, -Et, CH=CH2, CHO, -COMe ; R8 =-H, -Me, -Cl .

Внаслідок моделювання були отримані наступні результати (табл. 5):

Таблиця 5

Характеристики моделей, отриманих в результаті вирішення задачі структура- LogIC50 для ряду похідних 1,4-бенздіазепінів.

Модель	Структурні параметри	Кількість параметрів	Метод аналізу	Характеристики моделі
XVII	гармонічні	6	МЛР	R = 0.935;CVR = 0.913; F = 52.0; SE = 0.28
XVIII	локальні та польові	6	МЛР	R = 0.955; CVR = 0.917; F = 77.7; SE =0.24
XIX	всі види параметрів	29	PLS	R=0.974;CVR=0.941;A=2



Рис. 4. Польові параметри, які визначають величину LogIC₅₀ молекул 1,4-бенздиазепінів

Локальні та польові характеристики зосереджені в області замісників і відображають вплив кожного з них на величину LogIC₅₀. Так, імовірності присутності азотів (можливі фрагменти рецептору) в комірках 868, 145 і 213 (рис. 4) чутливі до впливу замісника в положенні R₁: наявність об'ємних замісників (наприклад, трет-бутілу) збільшує величину LogIC₅₀. Виявлено також, що сприяти зв'язуванню молекул з ГАМК-рецептором будуть великі значення ліпофільності замісника R₇, а також наявність електроноакцепторних угрупувань (наприклад, -NO₂, -CF₃, -F) в положеннях R₇ і R₂'. Наявність атома фтору в положенні замісника R₆' трохи зменшує величину LogIC₅₀.

Аналіз параметрів структурної відмінності показує, що чим менше відмінності структури (зокрема, по потенціалу ліпофільності та по наявності атомів галогену) по відношенню до молекул з низкими значеннями LogIC₅₀, тим менше значення параметру інгібування.

Отримані моделі, що відображають інгібування бенздіазепінового рецептору, істотно перевищують за статистичними характеристиками результати авторів [Winkler D.A., 1998] (R<0.75).

3.5 Аналіз зв'язку інтенсивності люмінесценції -дикетонатів лантаноїдів (III) з їх структурою

Велика кількість і різноманітність структурних параметрів дозволяє з оптимізмом дивитися на можливості ґраткової моделі. Тому, в рамках дисертаційної роботи, була проаналізована ефективність розробленого підходу, крім задач “структура-біологічна активність”, для дослідження інших властивостей. Так, зокрема, проведено аналіз зв'язку інтенсивності люмінесценції - дикетонатів лантаноїдів (III) з їх структурою. Як об'єкти дослідження розглядалися трис-хелатні -дикетонатні комплекси Автори виражають щиру вдячність ст.н.с. Мешковій С.Б. (ФХІ ім. О.В. Богатського НАН України) за надані дані по люмінесценції -дикетонатів лантаноїдів. Sm³⁺, Eu³⁺, Nd³⁺, Yb³⁺ (15 молекул):

	М= Sm, Eu, Nd, Yb R1= -Me, -CF3, -C3F7, -C6F13, -C8F17, CF3-O-CF2-, CF3-O-CF2-CF2-, C3F7-O-CF(CF3)-, CF3-O-(CF2)3-O-CF(CF3)- R2= -Me, -Ph, -C4H3S (-тієніл)

Виявилося, що найбільш тісно з інтенсивністю люмінесценції пов'язані характеристики поля ліпофільності. Зокрема, вирішальний вплив на люмінесценцію всіх -дикетонатів РЗЕ має потенціал поля ліпофільності в комірці 517 (рис. 5, табл. 6).

Таблиця 6

Регресійні моделі зв'язку структура-інтенсивність люмінесценції для комплексів лантаноїдів з -дикетонами.

Y=a + bF_Lip(517)	a	b	R	CVR	SE	F
Iлюм (Sm)	-49.5	67.4	0.971	0.960	9.2	216.6
Iлюм (Eu)	-95.2	115.2	0.959	0.944	19.0	150.6
Iлюм (Nd)	-35.2	54.2	0.968	0.957	7.9	190.5
Iлюм (Yb)	-23.4	50.4	0.956	0.941	8.5	138.2



Рис. 5. Комірка 517 навколомолеку-лярного простору -дикетонатів РЗЕ

ґратковий просторовий структура молекула

Важливо зазначити, що ця комірка якраз відповідає області простору в районі фторалкільних груп. Чим вище вказаний потенціал, тим інтенсивніше люмінесценція. Таким чином, отримані результати повністю відповідають гіпотезі про створення фторалкільними групами гідрофобного бар'єру, перешкоджаючого проникненню молекул води у внутрішню координаційну сферу комплексу, що протидіє гасінню люмінесценції.

3.6 Аналіз закручувальної здатності хіральних добавок, індукуючих спіральне упорядкування в нематичній фазі

Показано застосування розроблених характеристик хіральності молекул для ряду похідних (-)-ментону [(-)-транс-2S, 5R-2-ізопропил-5-метилциклогексанону], що індукують спіральне упорядкування в нематичній мезофазі:

		R=-H, - F, -Cl, -Br, -NO2, -CN, -OMe, -OEt, -Ph.

Внаслідок статистичного аналізу PLS-методом з використанням інтегральних параметрів вдалося отримати рівняння, що забезпечує задовільну якість зв'язку структура - закручувальна здатність (R=0.930; CVR=0.839; A=4).

Виявилося, що здатність закручувальних добавок визначається не тільки хіральністю молекул, але й їх формою, сумарною електронеґативністю атомів і молярною рефракцією. Ймовірно останні два параметри відображають здатність молекул до міжмолекулярних взаємодій, що, безсумнівно, є важливим фактором, який впливає на закручувальну здатність. Можна передбачити, що сумарна електронеґативність в певній мірі характеризує здатність молекул до електростатичних взаємодій, а молярна рефракція відображає здатність молекул до дисперсійних взаємодій.

Зазначено, що при вирішенні всіх розглянутих задач QSAR/QSPR головною метою була демонстрація ефективності розробленого нами підходу, а не вичерпне дослідження тієї або іншої властивості. Зрозуміло, що в кожному конкретному випадку можливе проведення більш глибокого дослідження, із збільшенням об'ємів і різноманітності навчальних і екзаменаційних вибірок і з більш зосередженою увагою до тих або інших типів структурних характеристик. Однак все це може бути предметом окремих самостійних досліджень.

4. Комплекс програм “MolLatt”

Розділ присвячений опису функціонування розробленого програмного комплексу “MolLatt”^©, за допомогою якого проводилася реалізація вищевикладених задач. Графічна візуалізація молекул, їх ґраткових моделей, полів і т.п. була виконана з використанням графічного стандарту фірми Silicon Graphics Open Graphic Language (OpenGL).

Зазначено, що програмний комплекс “MolLatt” складається з двох взаємопов'язаних програм “Lattice Model” і “DesAn” (Descriptors Analysis), що реалізовують розрахунок і обробку параметрів, відповідно. Обмін даними між програмами реалізовано за допомогою DDE DDE (Dynamic Data Exchange) - протокол обміну даними між різними програмами, які працюють у середовищі MS Windows..

Також описані використані в програмі “DesAn” прискорені методи побудови моделей зв'язку структура-властивість (алгоритми “генетичного типу”).

ВИСНОВКИ

Розроблена схема побудови ламаної спіральної кривої, що проходить через систему концентричних кубів молекулярної ґраткової моделі, адекватно описує конформаційні і конфігураційні особливості просторової структури молекули і дозволяє виразити їх у вигляді молекулярного коду.

Процедура тренд-вектору є ефективним інструментом для попереднього відбору структурних параметрів, що використовуються надалі для вирішення задач “структура-властивість” PLS-методом.

Наявність sp²-гібридного атому оксигену, пов'язаного з атомом карбону С₂ стероїдного скелету, є важливим фактором для вияву спорідненості стероїдів до кортикостероїд - зв'язуючого глобуліну. Наявність гідроксильної групи в псевдоекваторіальному положенні атому С₁₅ сприяє зв'язуванню стероїдів з тестостерон - зв'язуючим глобуліном. В обох випадках характер взаємного просторового розташування атомів оксигену стероїдних молекул істотно впливає на їх спорідненість до відповідних рецепторів.

Показано, що в рамках ґраткового підходу без проведення квантовохімічних розрахунків вдається побудувати адекватні статистичні моделі, що відображають вплив структури монопохідних бензолу на їх здатність до утворення комплексів з переносом заряду.

Для ряду заміщених 1,4-бенздіазепін-2-онів вдалося побудувати адекватні і досить надійні статистичні моделі, які пов'язують характеристики структури досліджуваних сполук з психотропною активністю (R=0.966) і здатністю інгібувати ГАМК-рецептор (R=0.974). Крім того, вдалося передбачити сполуку з потенційно високою психотропною активністю - 1-метил-5-(о-фтор)феніл-7-трифторметил-1,2-дигідро-3Н-бенздіазепін-2-он.

Збільшенню інтенсивності люмінесценції - дикетонатів лантаноїдів (III) сприяють фторалкільні замісники, що створюють високий потенціал поля ліпофільності в області хелатних циклів. Ймовірно, створений цими групами гідрофобний бар'єр, перешкоджає проникненню молекул води до іонів РЗЕ і, таким чином, протидіє гасінню люмінесценції.

Для похідних (-)-ментону виявлено, що три групи параметрів, які визначають хіральність, форму молекул, а також їх здатність до міжмолекулярних взаємодій дозволяють кількісно оцінити ефективність цих сполук як легуючих добавок, що забезпечують спіральне упорядкування в нематичній мезофазі.

СПИСОК ОПУБЛІКОВАНИХ ПРАЦЬ ЗА ТЕМОЮ ДИСЕРТАЦІЇ

Кузьмин В.E., Артеменко A.Г. Анализ хиральности молекул в рамках решеточной модели пространственной структуры // Ж. структ. химии.- 1998.- Т. 39, № 3.-С. 537-542.

Мешкова С.Б., Шапиро Ю.Е., Кузьмин В.Е., Артеменко А.Г., Русакова Н.В., Пыхтеева Е.Г., Большой Д.В. Конформационный анализ комплексов ланатана(Ш) и неодима (Ш) с перфтор--дикетонами Коорд. химия.-1998.-Т.24, №9.- С.714-718.

Кузьмин В.Е., Артеменко А.Г., Ковдиенко Н.А., Желтвай А.И. Решеточные модели молекул для решения задач QSAR // Хим.-фарм. журн.-1999.- № 9.- С.14-20.

Kuz'min V.E., Artemenko A.G, Kovdienko N.A., Tetko I.V., Livingstone D.J. Lattice Model for QSAR Studies J. Mol. Modeling.-2000.-№ 6.-P.517-526.

Кузьмин В.Е., Мешкова С.Б., Артеменко А.Г. Анализ связи интенсивности люминесценции -дикетонатов лантаноидов (III) с их структурой в рамках решеточной модели // Коорд. химия.-2000.- Т.26, №9.- С.705-710.

Хромов А.И., Кузьмин В.Е., Артеменко А.Г. Решеточно-топологическое моделирование молекул Материалы Ш семинара “Моделирование в прикладных научных исследованиях”.-Одесса.-1996.- С.15-17.

Кузьмин В.Е., Хромов А.И., Артеменко А.Г., Ковдиенко Н.А. Решеточно-топологическая модель молекул Тезисы докладов городского семинара “Применение вычислительной техники и математического моделирования в прикладных научных исследованиях”.-Одесса.-1994.- С. 8.

Kuz'min V.E., Artemenko A.G. A Lattice Model of the Molecule. The Chirality Analysis and Drug Design Abstracts of 7^th International Symposium on Chiral Discrimination (ISCD 7).- Jerusalem (Israel).-1995.-P.68.

Кузьмин В.E., Артеменко A.Г. Решеточная модель для исследования связей “структура-свойство” Тезисы докладов II городского семинара “Применение вычислительной техники и математического моделирования в прикладных научных исследованиях”.-Одесса.-1995.- С. 14.

Artemenko A.G., Kruglyak Y.A., Kuz'min V.E., Slone J.E. Complex 3D Representation of a System of Structural Parameters for Solving “Structure-Property” Problems // Abstracts of 7^th International Conference on Mathematical Chemistry.- Girona (Spain).-1997.-P.95.

Artemenko A.G., Kuz'min V.E., Ognichenko L.A., Onischenko S.I. New Informational-Field Characteristics of Molecular Recognition for Problem Solving “Structure-Property” // Abstracts of 7^th International Conference on Mathematical Chemistry.- Girona (Spain).-1997.-P.99.

Кузьмин В.E., Артеменко A.Г., Онищенко С.В., Огниченко Л.Н. Моделирование информационного поля молекул Тезисы докладов I Всероссийской конференции “Молекулярное моделирование”.-Москва.-1998.

Артеменко A.Г., Кузьмин В.E. Визуализация многомерного признакового пространства молекул для интерактивного решения задач 3D-QSAR Тезисы докладов I Всероссийской конференции “Молекулярное моделирование”.-Москва.-1998.

Артеменко A.Г. Решение задач QSAR в рамках решеточной модели Тезисы докладов I конференции молодых ученых и студентов-химиков южного региона Украины.-Одесса.-1998.- С. 7.

Artemenko A.G., Kuz'min V.E. Lattice Models for QSAR Tasks // Abstracts of 7^th International Symposium on Molecular Aspects of Chemotherapy.- Gdaтsk (Poland).-1999.-P.170.

Кузьмин В.E., Артеменко A.Г., Ковдиенко Н.А., Алиханиди С.Э., Огниченко Л.Н., Хромов А.И. Решеточные молекулярные модели для исследования связи “структура-противораковая активность” Тезисы докладов II съезда онкологов стран СНГ.-Киев.-2000. - Т.22, №261.

АНОТАЦІЯ

Артеменко А.Г. Аналіз зв'язку структура-властивість на основі молекулярних ґраткових моделей. - Рукопис.

Дисертація на здобуття наукового ступеня кандидата хімічних наук за спеціальністю 02.00.03 - органічна химия. - Фізико-хімічний інститут ім. О. В. Богатського НАН України, Одеса, 2000.

Розроблено і реалізовано ґратковий підхід до моделювання просторової структури молекул. Такий підхід забезпечує генерацію однорідної системи структурних параметрів, що відображають локальні, інтегральні та польові характеристики молекулярної структури для ефективного вирішення різних задач “структура-властивості”. Запропонована інваріантна система кодування просторової структури молекул. Отриманий молекулярний код використаний для оцінки хіральності і структурної подібності молекул. Модернізована система опису локальних і польових характеристик молекулярної структури.

Розроблений підхід застосовано для вирішення різноманітних задач “структура-властивості”. Виявлено, що спорідненість стероїдних молекул до кортикостероїд- і тестостерон- зв'язуючого глобуліну в значній мірі зумовлена просторовим розташуванням атомів оксигену. Показано, що адекватний опис впливу структури монопохідних бензолу на їх здатність до утворення комплексів з переносом заряду можливий в рамках ґраткової моделі без проведення квантовохімічних розрахунків.

Для ряду заміщених 1,4-бенздіазепін-2-онів вдалося передбачити сполуку з потенційно високою психотропною активністю - 1-метил-5-(о-фтор)феніл-7-трифторметил-1,2-дигідро-3Н-бенздіазепін-2-он. Показано, що фторалкільні замісники, що створюють високий потенціал поля ліпофільності в області хелатних циклів, сприяють збільшенню інтенсивності люмінесценції - дикетонатів лантаноїдів (III). Виявлено, що використовуючи три групи параметрів, які визначають хіральність, форму молекул, а також їх здатність до міжмолекулярних взаємодій, можливо кількісно оцінити закручувальну здатність похідних (-)-ментону в нематичній мезофазі.

На основі ґраткового підходу розроблено комплекс комп'ютерних програм, в межах якого можливе вирішення різноманітних задач QSAR/QSPR.

Ключові слова: задачі “структура-властивості”, ґратковий підхід, структурні параметри, просторова структура молекули, молекулярний код, спорідненість, активність, інтенсивність люмінесценції, закручувальна здатність.

Аннотация

Артеменко А.Г. Анализ связи структура-свойство на основе молекулярных решеточных моделей. - Рукопись.

Диссертация на соискание ученой степени кандидата химических наук по специальности 02.00.03 - органическая химия. - Физико-химический институт им. А. В. Богатского НАН Украины, Одесса, 2000.

Разработан и реализован решеточный подход к моделированию пространственной структуры молекул. Такой подход обеспечивает генерацию однородной системы структурных параметров, отражающих локальные, интегральные и полевые характеристики молекулярной структуры для эффективного решения различных задач “структура-свойства”. Предложена инвариантная система кодирования пространственной структуры молекул. Полученный молекулярный код использован для оценки хиральности и структурного подобия молекул. Модернизирована система описания локальных и полевых характеристик молекулярной структуры.

Разработанный подход применен для решения различных задач “структура-свойства”. Обнаружено, что сродство стероидных молекул к кортикостероид- и тестостерон- связывающему глобулину в значительной мере обусловлено пространственным расположением атомов кислорода. Показано, что адекватное описание влияния структуры монопроизводных бензола на их способность к образованию комплексов с переносом заряда возможно в рамках решеточной модели без проведения квантово-химических расчетов.

Для ряда замещенных 1,4-бенздиазепин-2-онов удалось предсказать соединение с потенциально высокой психотропной активностью - 1-метил-5-(о-фтор)фенил-7-трифторметил-1,2-дигидро-3Н-бенздиазепин-2-он. Показано, что фторалкильные заместители, создающие высокий потенциал поля липофильности в области хелатных циклов, способствуют увеличению интенсивности люминесценции -дикетонатов лантаноидов (III). Обнаружено, что, используя три группы параметров, определяющих хиральность, форму молекул, а также их способность к межмолекулярным взаимодействиям, возможно количественно оценить закручивающую способность производных (-)-ментона в нематической мезофазе.

На основе решеточного подхода разработан комплекс компьютерных программ, в рамках которого возможно решение самых разных задач QSAR/QSPR.

Ключевые слова: задачи “структура-свойства”, решеточный подход, структурные параметры, пространственная структура молекулы, молекулярный код, сродство, активность, интенсивность люминесценции, закручивающая способность.

ANNOTAtion

Artemenko A.G. Analysis of structure-property relationship on the base of molecular lattice models. -Manuscript.

Thesis for a candidate's degree in Chemical Sciences on speciality 02.00.03. - organic chemistry.- A.V.Bogatsky Phisico-Chemical Institute, National Academy of Sciences of Ukraine, Odessa, Ukraine, 2000.

Lattice approach to modeling of a spatial molecular structure is developed and realized. This approach provides a generation of homogeneous system of structural parameters, which reflect local, integral and field characteristics of molecular structure for effective solution of diverse “structure-property” tasks.

Placing a molecule with a fixed spatial structure into a lattice of cubic cells (the size of a cell 1.8-2.0Е) is essence of method. The invariant disposition of the molecule in the lattice is achieved by superposition of the centre of mass of the molecule with the origin of the coordinates. In addition the principal axes of inertia of the molecule are also superimposed with the coordinate axes of the lattice. If the analyzed structures contain a large common structural fragment, their alignment is carried out mainly according to this fragment. An information about each cell of the spatial lattice, namely, presence or absence of atom (fragment) in it, a value physical-chemical property of atom (fragment) or field characteristics of environment space, are the source of structural information and used as structural descriptors for “structure-property” studies.

An invariant system of encoding of spatial molecular structure is offered on the base of construction of broken spiral curve within the frameworks of lattice model. The molecular code used for estimation chirality and structural similarity of molecules.

A system of description of local structure characteristics, taking into accounts not only presence or absence of atom (group) in that or another region of molecular space, but their diverse physical-chemical properties is offered. Besides, а new characteristics are developed for potential pharmacophores (functional groups), reflecting their accessibility /congestion, and also mutual position and orientation.

A probabilities of an occupancy of a different environment space regions by atoms (groups) which modeling functional (biological) target (pseudoreceptor) are offered.

A system of local and field parameters of molecular structure which is calculated on the base of lattice model allowed to create a scheme of constructing of fragments of biological target (pseudoreceptor) using as “matrix” the molecules of work set of compounds under investigation.

The new structural invariants that reflect the distribution of atoms, fragments or molecular fields relatively to spatial structure of molecules are offered on the base of harmonic analysis of local and field characteristics of molecules.

Lattice approach is used for solution of diverse “structure-property” tasks. The structural descriptors of molecules that promote and prevent from displaying of researched properties are detected. The multiple linear regression analysis, trend-vector procedure and partial least squares (PLS) method (widely used for QSAR/QSPR studies with great number of descriptors) were applied to detect “structure-property” relationships.

...