ISSN 2070-4909 (print) ISSN 2070-4933 (online)
CoBpeivieHHafl $apMa^KOHOMUKa / $apMa^nMfleMMonomfl
^
0 :r
s
^
o 0 T Ci 0
o m
*
1
f
m o o _o
o
o
■v
J
Ii
0
1
T TO I
CO TO
Üf
w
0 X
/
rn
ib" co
m
„-rosvuttw anaW G\J^co
OVA
www.pharmacoeconomics.ru
<n g
'iE o c o o
o o
TO OT
iE co ra ^
.c m
> ^
> Is-§ -lh
S^ <D
TO I—
TO O 0
TO O
TO I TO T
ra
^
o co
_Q
lO
s .0 I— TO
m .0
S
T
>
c o c
0
1
*
o
TO I— I
—
FARMAKOEKONOMIKA
Modern Pharmacoeconomics and Pharmacoepidemiology
m
0
2022 Vol. 15 No. 2
ftw 15
0
eis
2022
https://doi.org/10.17749/2070-4909/farmakoekonomika.2022.142
ISSN 2070-4909 (print) ISSN 2070-4933 (online)
Биофизическое моделирование указывает на высокое сродство этиловых эфиров омега-3 полиненасыщенных жирных кислот к ферментам провоспалительного каскада арахидоновой кислоты
Торшин И.Ю., Громова О.А.
Федеральный исследовательский центр «Информатика и управление» Российской академии наук (ул. Вавилова, д. 44, корп. 2, Москва 119333, Россия)
Для контактов: Громова Ольга Алексеевна, e-mail: unesco.gromova@gmail.com
РЕЗЮМЕ
Актуальность. Омега-3 полиненасыщенные жирные кислоты (ю-3 ПНЖК) широко используются в современной кардиологии. Противовоспалительный эффект эйкозапентаеновой кислоты (ЭПК) и докозагексаеновой кислоты (ДГК) осуществляется, в частности, за счет модуляции активности каскада арахидоновой кислоты (АРК).
Цель: провести биофизический анализ взаимодействий различных форм ПНЖК (АРК, ЭПК, ДГК, этиловые эфиры ЭПК и ДГК) с тремя таргетными белками - ферментами каскада АРК (циклооксигеназа-1 (ЦОГ-1), ЦОГ-2, 5-липоксигеназа (5-ЛПОГ)).
Материал и методы. Минимизацию энергии комплексов «лиганд-рецептор» при различных начальных расположениях молекулы относительно кармана рецептора и молекулярные энергетические профили белка (энергии переноса боковой цепи из растворителя в белок) вычисляли посредством пакета ЕСММ^.
Результаты. Этиловые эфиры ЭПК и ДГК характеризовались значительно большими абсолютными значениями энергий специфического связывания таргетных белков, чем АРК или просто ЭПК/ДГК. Например, в случае ЦОГ-2 ЭПК, ДГК и этиловый эфир ЭПК имели близкие значения ААб (-3,0...-3,1 ккал/моль), в то время как комплекс ЦОГ-2 с этиловым эфиром ДГК был несколько более стабилен (ААв = -3,4 ккал/моль). В случае фермента 5-ЛПОГ наиболее стабильным был комплекс с этиловым эфиром ЭПК (ААб = -1,62 ккал/моль). Более высокие абсолютные значения энергий связывания указывают на более высокое сродство этиловых эфиров ЭПК и ДГК к ферментам каскада АРК.
Заключение. Полученные результаты позволяют предположить, что высокое сродство этиловых эфиров ЭПК и ДГК к таргетным белкам каскада АРК является одной из молекулярных основ терапевтической эффективности препаратов ю-3 ПНЖК с высокой степенью стандартизации.
КЛЮЧЕВЫЕ СЛОВА
Аффинность к рецепторам, омега-3 полиненасыщенные жирные кислоты, ПНЖК, молекулярная фармакология. Статья поступила: 11.05.2022 г.; в доработанном виде: 09.06.2022 г.; принята к печати: 30.06.2022 г. Конфликт интересов
Авторы заявляют об отсутствии необходимости раскрытия конфликта интересов в отношении данной публикации. Финансирование
Работа выполнялась по государственному заданию № 0063-2019-0003 «Математические методы анализа данных и прогнозирования» с использованием инфраструктуры Центра коллективного пользования «Высокопроизводительные вычисления и большие данные» ФИЦ ИУ РАН (г. Москва).
Вклад авторов
Авторы сделали эквивалентный вклад в подготовку публикации. Для цитирования
Торшин И.Ю., Громова О.А. Биофизическое моделирование указывает на высокое сродство этиловых эфиров омега-3 полиненасыщенных жирных кислот к ферментам провоспалительного каскада арахидоновой кислоты. ФАРМАКОЭКОНОМИКА. Современная фармакоэкономика и фармакоэпидемиология. 2022; 15 (2): 221-229. https://doi.org/10.17749/2070-4909/farmakoekonomika.2022.142.
к ц
ф :г
s
^
о ф
т
с^ ф
о
ш к
S
т. го ш
о
СО _о
о
о
к ц
ч
о
т.
ф
т га т. п го
it
(Л
ф X
ш ю
О)
ю
(Л
о
Е о с о
о ф
о о
ГО ОТ с ^ Е со го ^ л ю о. от
5 ^
> Is-§ +
о. Ф
го н
го о
го
I
го т
го
^
о
со ц
_Q
ю
S _0 н го
к
0
с^ ф
ш
н
ф
1
-О S
т
>
ц
о с
0
1
о
Ü о
£
™ о
I ■& ГО I CIS
Review articles
FWÉÉMÉi
:
Biophysical modeling indicates a high affinity of ethyl esters of omega-3 polyunsaturated fatty acids to the enzymes of the pro-inflammatory arachidonic acid cascade
Torshin I.Yu., Gromova O.A.
Federal Research Center"Informatics and Management', Russian Academy of Sciences (44 corp. 2 Vavilov Str, Moscow 2119333, Russia)
Corresponding author: Olga A. Gromova, e-mail: unesco.gromova@gmail.com
SUMMARY
Background. Omega-3 polyunsaturated fatty acids (ro-3 PUFAs) are widely used in modern cardiology. The anti-inflammatory effect of eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) is realized, in particular, by modulating the activity of the arachidonic acid (ARA) cascade.
Objective: to conduct biophysical analysis of the interactions of various forms of PUFAs (ARA, EPA, DHA, EPA and DHA ethyl esters) with three target proteins-enzymes of the ARA cascade (cyclooxygenase-1 (COX-1), COX-2, 5-lipoxygenase (5-LPOG)).
Material and methods. The minimization of the energy of the ligand-receptor complexes at various initial positions of the molecule relative to the receptor pocket and the molecular energy profiles of the protein (the energy of side chain transfer from solvent to protein) were calculated using the ECMMS package.
Results. EPA and DHA ethyl esters were characterized by significantly higher absolute values of specific binding energies of target proteins than ARA or simply EPA/DHA. For example, in the case of COX-2, EPA, DHA and EPA ethyl ester had similar AAG values (-3.0...-3.1 kcal/mol), while the COX-2 complex with DHA ethyl ester was somewhat more stable (AAG = -3.4 kcal/mol). In the case of the 5-LPOG enzyme, the complex with EPA ethyl ester was the most stable (AAG = -1.62 kcal/mol). Higher absolute binding energies indicate a higher affinity of EPA and DHA ethyl esters to enzymes of the ARA cascade.
Conclusion. The results suggest that the high affinity of EPA and DHA ethyl esters to the enzymes of the ARA cascade is one of the molecular bases of therapeutic efficacy of ro-3 PUFAs with high degree of standardization.
KEYWORDS
Receptor affinity, omega-3 polyunsaturated fatty acids, PUFAs, molecular pharmacology. Received: 11.05.2022; in the revised form: 09.06.2022; accepted: 30.06.2022 Conflict of interests
The authors declare they have nothing to disclose regarding the conflict of interests with respect to this manuscript. Funding
The research was conducted under the government contract No. 0063-2019-0003 "Mathematical methods of data analysis and forecasting" using the infrastructure of the Shared Research Facilities "High Performance Computing and Big Data" of Federal Research Center "Informatics and Management", Russian Academy of Sciences (Moscow).
Authors' contribution
The authors contributed equally to this article. For citation
Torshin I.Yu., Gromova O.A. Biophysical modeling indicates a high affinity of ethyl esters of omega-3 polyunsaturated fatty acids to the enzymes of the pro-inflammatory arachidonic acid cascade. FARMAKOEKONOMIKA. Sovremennaya farmakoekonomika i farmakoepidemio-logiya / FARMAKOEKONOMIKA. Modern Pharmacoeconomics and Pharmacoepidemiology. 2022; 15 (2): 221-229 (in Russ.). https://doi. org/10.17749/2070-4909/farmakoekonomika.2022.142.
ВВЕДЕНИЕ/ INTRODUCTION
Увеличение потребления омега-3 полиненасыщенных жирных кислот (ю-3 ПНЖК) достоверно ассоциировано с дозозависи-мым снижением смертности от всех причин [1]. Детальный аудит доказательных публикаций по ю-3 ПНЖК, проведенный посредством современных математических методов анализа больших данных [2], подтвердил существование обширной клинической базы, указывающей на обоснованность использования стандартизированных препаратов на основе ю-3 ПНЖК для снижения сердечно-сосудистой смертности.
Аудит доказательных публикаций [2] указал на необходимость изучения количественного состава препаратов на основе ю-3 ПНЖК. Сравнительное исследование состава 10 образцов ю-3 ПНЖК с хроматографическим определением более 40 метаболи-
тов жирных кислот и других соединений позволило разработать количественные маркеры для идентификации более стандартизированных образцов ю-3 ПНЖК - «эйкозапентаеновая кислота (ЭПК) + докозагексаеновая кислота (ДГК)», «ю-6 + ю-11», коэффициент стандартизации и др. [3].
Анализ количественного и качественного состава фармацевтических образцов ю-3 ПНЖК показал, что в них входят разные формы ю-3 ПНЖК (триглицериды, карбоновые кислоты, этиловые эфиры) (рис. 1, а). Каждая из этих форм имеет химические особенности, отражающиеся на их фармакологических свойствах. Молекулы триглицеридов (сложные эфиры с глицерином) достаточно громоздки: остаток глицерина - многоатомный спирт, он включает три гидроксильные (-ОН) группы (рис. 1, Ь). Для взаимодействия с таргетными белками триглицериды должны
га g
_û
ю
s _0 I—
га
т
>
о с
0
1
*
о
§ £
° i ? ф
H Œ Ï О Œ о
ê!
£
? о
I -&
га i eis
Основные моменты
Что уже известно об этой теме?
► Стандартизированные препараты на основе омега-3 полиненасыщенных жирных кислот (ю-3 ПНЖК), воздействующие на синтез эйкозаноидов и докозаноидов в каскаде арахидоновой кислоты (АРК), успешно используются в терапии кардиоваскулярных заболеваний
► Препараты ю-3 ПНЖК существенно различаются по количественному и качественному составу
► Анализы состава препаратов указывают на перспективы применения сложных эфиров ю-3 ПНЖК
Что нового дает статья?
► Сложные эфиры ю-3 ПНЖК характеризуются значительно большими энергиями специфического связывания ферментов каскада АРК, чем свободные формы ПНЖК
► Абсолютные значения энергии связывания сложных эфиров ю-3 ПНЖК с таргетными белками (циклооксигеназы, липоксигеназа) возрастают с увеличением длины цепи сложного эфира
► Более высокие значения энергий связывания указывают на более высокое сродство сложных эфиров эйкозапентаеновой кислоты (ЭПК) и докозагексаеновой кислоты (ДГК) к таргетным белкам
Как это может повлиять на клиническую практику в обозримом будущем?
► Более стабильные комплексы с таргетными белками соответствуют более интенсивному образованию противовоспалительных эйкозаноидов из ю-3 ПНЖК
► Препараты, содержащие сложные эфиры ю-3 ПНЖК, могут более эффективно взаимодействовать с различными ферментами, чем препараты, содержащие только ю-3 ПНЖК
► Более высокое сродство сложных эфиров ЭПК и ДГК к таргетным белкам каскада АРК может повышать терапевтическую эффективность соответствующих препаратов
Отрицательный заряд
скомпенсирован/-The negative charge is compensated
Арахидоновая кислота / Arachidonic acid
5 8 ^ 11 v 14 .........Отрицательный заряд (в щелочной среде) / The negative charge (in alkaline medium)
Эйкозапентаеновая кислота (ЭПК) / Eicosapentaenoic acid (EPA) и
Докозагексаеновая кислота (ДГК) / Docosahexaenoic acid (DHA)
^^ _
4 7 10 13 16 19
Highlights
What is already known about the subject?
► Standardized drugs based on omega-3 polyunsaturated fatty acids (ra-3 PUFAs), which affect the synthesis of eicosanoids and docosanoids in the arachidonic acid (ARA) cascade, are successfully used in the treatment of cardiovascular diseases
► Medicines of ra-3 PUFAs differ significantly in quantitative and qualitative composition
► Analyzes of the composition of medicines indicate the prospects for the use of ra-3 PUFA esters
What are the new findings?
► ra-3 PUFA esters are characterized by significantly higher specific binding energies for the enzymes of the ARA cascade than free forms of PUFAs
► The absolute values of the binding energy of ra-3 PUFA esters with target proteins (cyclooxygenases, lipoxygenase) increase with increasing the ester chain length
► Higher binding energies indicate a higher affinity of eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) esters to target proteins
How might it impact the clinical practice in the foreseeable future?
► More stable complexes with target proteins correspond to more intense formation of anti-inflammatory eicosanoids from ra-3 PUFAs
► Medicines containing ra-3 PUFA esters may interact more effectively with various enzymes than medicines containing only ra-3 PUFAs
► The higher affinity of EPA and DHA esters to target proteins of the ARA cascade may increase the therapeutic efficacy of the respective medicines
Остаток глицерина - многоатомного спирта с тремя -ОН группами / The residue of glycerin - polyatomic alcohol with three -OH groups
H,C-0
Могут быть представлены различные типы жирных кислот / Different types of fatty acids can be represented
НС-О
HoC^O
Объем «хвоста» очень мал - всего два атома углерода / The volume of the "tail" is very small - only two carbon atoms
Этиловый эфир ЭПК / EPA ethyl ester 5 8 ^^14
Этиловый эфир ДГК / DHA ethyl ester
a
Рисунок 1. Полиненасыщенные жирные кислоты и их производные, изученные в настоящем исследовании: а - химические структуры; b - сравнение структур триглицеридов и этиловых эфиров Figure 1. The studied polyunsaturated fatty acids and their derivatives: a - chemical structures; b - comparison of the structures of triglycerides and ethyl esters
предварительно гидролизоваться до свободных ПНЖК, что неизбежно сказывается на фармакокинетике образцов ю-3 ПНЖК, содержащих триглицериды.
Свободные ПНЖК, как правило, представлены в организме в диссоциированной форме и характеризуются отрицательным
зарядом в карбоксильной группе (-СОО-). Заряд карбоксила меняется в зависимости от рН и может воспрепятствовать взаимодействию с таргетными белками. Этиловые эфиры, в отличие от триглицеридов, могут взаимодействовать с таргетными белками без предварительного гидролиза эфира до ПНЖК и спирта
к (D
S
О (D Т Œ (D
О ^
m
о: s
X СО m о
СО
.û ^
о с о s
о: ^
О X (D Т СО X
со СО X
41 F^ О
С ч—
-1- СО
а т
w g
О L-
15
8 ¡? (D TJ-
8
СО Oi
F ^ Е со
со ^
SZ ю Œ d
Ét
w 5
с^ (D
Т
О с
о
X
*
о
СО
СО X СО т
СО ^
о со s -Û ю
со
I— —
0 X ГС СО
1 Ê
° i &|
? ё
I— Œ
ï °
£ а ? о
х -а
СО х
b
(см. рис. 1, b). В отличие от свободных жирных кислот в этиловых эфирах ю-3 ПНЖК некомпенсированный отрицательный заряд карбоксильной группы отсутствует.
Приведенные выше аргументы носят качественный характер, так что их следует дополнить количественными оценками взаимодействия рассматриваемых молекул с таргетными белками. Как известно, ю-3 ПНЖК характеризуются значительными противовоспалительными и вазодилаторными свойствами в связи с их влиянием на метаболизм простагландинов через каскад арахи-доновой кислоты (АРК) [4]. В ходе биотрансформаций в каскаде АРК ПНЖК трансформируются циклооксигеназами (ЦОГ) в про-станоиды (простагландины, простациклины, тромбоксаны), а посредством липоксигеназ (ЛПОГ) - в лейкотриены. В отличие от АРК ю-3 ПНЖК и их производные проявляют гораздо меньшее провоспалительное действие - наоборот, предотвращают синтез провоспалительных простаноидов [5, 6] и способствуют синтезу противовоспалительных [7], нейропротекторных [8] и вазодила-торных простаноидов [9]. Более того, ю-3 ПНЖК необходимы для биосинтеза медиаторов разрешения воспаления - резолвинов, нейропротектинов, маресинов и др. [10].
Как видно из приведенного выше краткого экскурса в биохимию ю-3 ПНЖК, наиболее важными белками - рецепторами ю-3 ПНЖК являются ЦОГ и ЛПОГ. Взаимодействия различных форм ю-3 ПНЖК перспективно исследовать именно с этими ферментами. Расчет значений энергий взаимодействия исследуемых таргет-ных молекул с вышеупомянутыми таргетными белками ранее не проводился. Оценка сродства различных форм ПНЖК к этим таргетным белкам позволяет выявить те формы ю-3, от которых следует ожидать наибольшего терапевтическего эффекта.
Цель - провести биофизический анализ взаимодействий различных форм ПНЖК (АРК, ЭПК, ДГК, этиловые эфиры ЭПК и ДГК) с тремя таргетными белками - ферментами каскада АРК (ЦОГ-1, ЦОГ-2, 5-ЛПОГ).
МАТЕРИАЛ И МЕТОДЫ / MATERIAL AND METHODS
Исследование проводилось современными методами биофизического моделирования и биоинформатики, включающими расчет энергий взаимодействий в комплексах «молекула-рецептор».
Химические структуры исследованных ПНЖК приведены на рисунке 1. Выбор молекул для исследования очевиден: АРК -основной субстрат каскада АРК, ЭПК и ДГК - формы «морских» ю-3 ПНЖК с наибольшей базой клинико-эпидемиологических исследований. Этиловые эфиры ДГК и ЭПК - действующие начала препарата ю-3 ПНЖК, который характеризуется наибольшей доказательной базой среди препаратов ю-3 ПНЖК.
В работе исследовались взаимодействия этих ПНЖК с ЦОГ-1, ЦОГ-2 и 5-ЛПОГ. Модели пространственных структур указанных ферментов были приготовлены на основе структур, имеющихся в базе данных Protein Data Bank (PDB)1, содержащей все известные к настоящему времени структуры белков, определенные посредством диффракционных экспериментов или ядерного магнитного резонанса (табл. 1). Структуры из PDB были отобраны с учетом: 1) наличия в соответствующих файлах информации о сайтах связывания лигандов (жирных кислот и коферментов); 2) качества эксперимента по установлению пространственной структуры. В случае ЦОГ коферментами являлись протопорфирин с ионом кобальта, а в случае 5-ЛПОГ - ион трехвалентного железа.
На основе отобранных из базы данных PDB-файлов были приготовлены модели структур ферментов человека (методика опи-
Таблица 1. Использованные в исследовании файлы базы данных Protein Data Bank (PDB) по пространственным структурам ферментов каскада арахидоновой кислоты
Table 1. The Protein Data Bank (PDB) database files on the spatial structures of enzymes of the arachidonic acid cascade used in №е study
Фермент / Enzyme PDB-файлы / PDB-files
Циклооксигеназа-1 I Cyclooxygenase-1 1diy 1igx
Циклооксигеназа-2 I Cyclooxygenase-2 3hs5 3hs6 3hs7
5-липоксигеназа I 5-lipoxygenase 1ik3
сана в [11]). Затем были исследованы возможности размещения молекул ю-3 ПНЖК в связывающих карманах ферментов. Для этого проводилась минимизация энергии комплексов «лиганд-рецептор» при различных начальных расположениях молекулы относительно кармана рецептора.
Метод молекулярных энергетических профилей / The method of molecular energy profiles
В данной работе для анализа структур комплексов ферментов с ПНЖК использовался метод молекулярных энергетических профилей. Молекулярные энергетические профили являются новой концепцией в биоинформатике и биофизике белка. Они используются для определения аминокислотных остатков в структуре белка, которые обладают экстремальными (наиболее высокими или наиболее низкими) значениями энергии переноса боковой цепи из растворителя в белок (AG0CT). Остатки с экстремальными значениями AG0CT часто являются наиболее важными для реализации биологических функций белка [11].
Молекулярные энергетические профили представляет собой список значений AG0CT для каждого остатка в белке. Список упорядочен в соответствии с порядковыми номерами остатков в аминокислотной последовательности белка-рецептора. Для заданной трехмерной структуры белка значения AG0CT рассчитываются с использованием методов классической молекулярной механики. Процедура расчета состоит из двух этапов: минимизации энергии молекулярной системы и собственно расчета AG0CT. Ниже приводится детальное описание этого вычислительного метода.
Минимизация энергии
Минимизация полной энергии осуществлялась посредством алгоритмической комбинации, которая включает перемежающиеся стадии минимизации энергии посредством конъюгирован-ных градиентов и коротких периодов молекулярной динамики. Подобного рода комбинация (названная «интерактивной оптимизацией») позволяет избежать «застревания» системы в локальных минимумах. Интерактивная оптимизация является итеративной процедурой. Один цикл минимизации состоит из 50...200 шагов минимизации конъюгированных градиентов, за которой следует 100...1000 фемтосекунд молекулярной динамики. Применение метода конъюгированных градиентов в молекулярной механике заключается в смещении каждого атома в направлении действующих на этот атом сил. Для молекулы из N атомов сила, действующая на каждый атом, определяется как:
F (Xj) = V XU полн,
где ип0лн - полная внутренняя энергия молекулярной системы,
j = 1...N.
к ц
ф :г
s
о ф т
œ ф
0
ш
к s
1
га ш
о
СО _о
о
о
к ц
ч
0
1
ф
т га i
со га
if
(Л
ф X
о ю
О)
ю
(Л
о
Е о с о
о ф
о о
ГО ОТ
с ^
Е со
го ^
л ю
о. от % ^
> Is-§ +
о. Ф
го
I—
го о
го
I
го т го
о
со g
_û
ю
s _0 I—
го 15
ГС
0 ф
ш
I—
ф
1
_0
s
т
>
ц
о с
0
1
X о
1 www.rcsb.org.
Œ о ê!
S
? о
X -& ГО I CIS
Для г'-го шага минимизации по конъюгированным градиентам новые координаты каждо го атом а определяются как:
х\ = +'-1 +kxh'. j j 1 -
h ) = F(x '.) + ■
'r k ''h h-
Fkx?) I2 7
Сполн молекуляриой с;истеглы, соде ржащей N атомов, В кова-лентных связей (связывающих термов), /2 угловых термов и T торсионных термов, определяется как:
B ( T N N
j = а _паа ej +L Ej +а а в* ,
оолн / f связ / 1 угл / t торс несвяз '
j=) е=) j=) j=) k=j+)
< = ^ (0j - 0Ц с2,
1 n
j =~аа vm (i+c°s(m—j --о сс,
2 m=1
где £увяз - терм энергииу-й ковалентной связи; Еуугл - терм энергии /-го углового изгиба двух ковалентных связей; £{орс - терм энергии вращения /-го диэдрал ьного угла.
Терм несвязывающих взаимодействий определяется как сумма электростатическогг и Вандервнальсового термов:
т71,к _ pj,k + ej,к -^несвяз ^ эл т В ДВ-
Электростатич еские взаи модействия о писываются законом Кулона:
цЦ2 = qjqк/4яее0гм ,
а Вандерваальсовы - посредством потенциала Леннарда-Джонса:
М A
Fj,k _-мк + AA_kL
Ц]Е!ДВ Г6 гп кМк rj0
Силовые константы ККЬ, К» Ут, равновесные значения связей и углов г0, е0, ф0, а также константь\Аин,В]к определены в соответствии с потенциалом 3. Значения частичных атомных зарядов qJ, дк определены на основе потенциала /\IMMBER при составлении которого использовалпсь квантово-механические расчеты. Вычисления проводили посеедством пакета Е(СМ1\5Е! [10].
Молекулярная динамика
Процедур;! молекулярной динамики основана на втором законе Ньютона:
Р(хи . = »ту =»2.^(0, К и ^ ) = Х, (1 ).
Значения ускор(5иий= ( индиви дуальных г^т^мов находятся на основе иполн (см. выше), а затем численно интегрируются с исполь-
зованием алгоритма Верле. Скорости и ководинаты в следующий момент времени г + Зз.
Х (1 + .1) = ХН1) + 7(1 )А1 + 0,У 5(1 )(А1 )2, .(г м Аг) = м 0,УАП (а) м 5(1 м А))).
Начальны е скорости атомов были рандомизи рованы в соответствии с распределен)ем МЛакс^елл^:
D(v) = j /2 KkTf/2 4nv2 e-
Температура системы фзикси(эовалась в соответствии с равенством:
1 N 3
- YmVf = —NkBT.
v (t) J 5 T B
н ХкК
Для расчетов использовались Т = 300К и е = 1.
Расчет эне|]гий переноса
Расчеты энергий! переноса боковой цепи из растворителя в белок (ДОои) проводили на о с нове нес в язываюнциго потенциала, определяемого наосновании вышеприведенных уравнений молекулярной механики и термв э нергии со льватации. Таким образом, ДОост г'-го аминлкислотного остатиа, состоящнго изк атомов:
Ав* = АЕн + АА О1 .
ост несвяз нольв
В данном случае Д£НесВн представляет собой .азность между несвязываюАлй энер=ией взаимодействия атомов аминокислотного остатка со всеми атомами белка и энерглей взаимодействия атомов внутри аминокислотного остатка:
M i
= Е' — Е' х£"
несв,белок несв,вакк несв,белок
и определяется как:
к ОТ
е' = "V V о е/3 + лЛз)
-^оенв,бел.ж йзП эл 13ДВ/ '
у=л з=л,зг^у
к к
Е' = \ \3Е]'3 + Еу'3Л
несв,вакк эл т -^ВДВ/ '
/=.1 з=у+й
где у = 1...Н соответствует л;уммированию по атомам г'-го аминокислотного остатка [10].
ДДОСольв вычислялась кcвк рзао^тость между термами сольватации внутри и вне белка:
АА^Сольв =ʧуASAi,
/=л
где 5у - атомный сольватационный параметр, полученный на основании экспериментов по определению LogP в системе «октанол-вода», а - изменение доступности поверхности растворителю при переносе остатка из вакуума в белок.
РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ / RESULTS AND DISCUSSION
Получены полноатомные пространственные структуры комплексов исследованных молекул с таргетными белками. Затем проведен расчет энергий связывания и измерены молекулярные энергетические профили. В последующих разделах полученные результаты расчетов предваряются краткими введениями по биологии соответствующих ферментов.
к ц
ф
:г
s
^
о ф т
œ ф
0
m
гс s
1 га ш о
СО _о
о
о
ГС
ц
0
1
ф
т га i
со га
ы
ф X
m ю
О)
ю
(Л
о
'1= о с о
о ф
о о ГО От с ^ E C!D ГО ^ л ю а. от
5 ^^
> г-§ нн
о. Ф
го
I—
го о
го
I
го т го
«J
со g
-Û
ю
.û н го
ГС
0 ф
ш
I—
ф
1
-Û
s
т
>
ц
о с
0
1
*
о
Œ о
? о i ■& ГО I CIS
Основным параметром, позволяющим оценить сродство той или иной молекулы к определенному белку-рецептору, является AAG -изменение свободной энергии (AG) образования комплекса «молекула-рецептор» минус изменение свободной энергии молекулы рецептора при отсутствии связанного лиганда. Для стабильных комплексов значения AAG отрицательны. Чем меньше значение AAG (т.е. чем больше абсолютное значение AAG), тем больше стабильность комплекса и тем выше сродство соответствующего лиганда к таргетному белку.
Результаты расчетов AAG комплексов «лиганд-белок» / Results of calculations of ligand-protein complexes AAG
В таблице 2 и на рисунке 2 суммированы результаты расчетов AAG по всем вычислительным экспериментам, проведенным в ходе исследования. Как видно из расчета энергий связывания, ю-3 ПНЖК и их производные характеризовались большими по сравнению с АРК абсолютными значениями энергий связывания AAG с белками-рецепторами.
Циклооксигеназы-1 и -2: основные ферменты синтеза простагландинов / Cyclooxygenases-1 and -2: the main enzymes of prostaglandin synthesis
Существует две разновидности циклооксигеназ (также известных как простагландин Н2 синтетазы) - ЦОГ-1 и ЦОГ-2. Обе разновидности на 65% идентичны по аминокислотной последователь-
ности и практически идентичны по пространственной структуре. Ткани организма отличаются друг от друга по уровням экспрессии ЦОГ-1 и ЦОГ-2. ЦОГ-1 экспрессируется в большинстве клеток и является необходимым ферментом для биосинтеза простагландинов. ЦОГ-2 - индуцибельный фермент, который экспрессируется в активированных макрофагах и других клетках при воспалении.
Результаты расчета энергий связывания (см. рис. 2) показывают, что ЭПК, ДГК и их этиловые эфиры характеризуются более высокими абсолютными значениями энергий связывания с ЦОГ-1 и ЦОГ-2 по сравнению с ю-6 АРК. За исключением комплекса ЭПК-ЦОГ-1, энергии связывания ЭПК и ДГК на 0,4...1,8 ккал/моль превышали абсолютные значения энергии стабилизации комплексов АРК-ЦОГ-1 и АРК-ЦОГ-2. В случае ЦОГ-2 ЭПК, ДГК и этиловый эфир ЭПК характеризовались близкими значениями ААG (-3,0...-3,1 ккал/моль), в то время как комплекс ЦОГ-2 с этиловым эфиром ДГК был несколько более стабилен (ААG = -3,4 ккал/моль).
Анализ молекулярно-энергетических профилей (рис. 3) показал, что этот выигрыш в энергии стабилизации комплексов обусловлен не изменениями в характере взаимодействий каких-либо отдельных аминокислотных остатков, а более эффективной (по сравнению с АРК) аккомодацией молекул ю-3 ПНЖК и их производных в структуре таргетных белков ЦОГ-1 и ЦОГ-2.
Анализ молекулярно-энергетических профилей, представленных на рисунке 3, позволил сделать вывод, что более высокая
Таблица 2. Энергии связывания исследованных молекул полиненасыщенных жирных кислот с ферментами каскада арахидоновой кислоты, ккал/моль Table 2. Binding energies of the studied polyunsaturated fatty acids molecules with enzymes of the arachidonic acid cascade, kcal/mol
Фермент / молекула ПНЖК // Enzyme / PUFA-molecule Циклооксигеназа-1 / Cyclooxygenase-1 Циклооксигеназа-2 / Cyclooxygenase-2 5-липоксигеназа / 5-lipoxygenase
АРК/ARA -1,65 -2,44 -0,69
ЭПК / EPA -1,01 -3,02 -0,95
ДГК / DHA -2,67 -3,00 -1,21
Этиловый эфир ЭПК / EPA ethyl ester -3,42 -3,10 -1,62
Этиловый эфир ДГК / DHA ethyl ester -3,34 -3,42 -1,23
Примечание. ПНЖТ - полиненасыщенная жирная кислота; АРК - арахидоновая кислота; ЭПК - эйкозапентаеновая кислота; ДГК - докозагексаеновая кислота. Расчет был проведен на основе приготовленных моделей комплексов «ПНЖК-белок». Приведены значения относительных энергий образования комплекса «молекула-рецептор». Жирным шрифтом выделены молекулы с наименьшими значениями энергий связывания (т.е. с наибольшим сродством к таргетным белкам).
Note. PUFA - polyunsaturated fatty acid; ARA - arachidonic acid; EPA - eicosapentaenoic acid; DHA - docosahexaenoic acid. The calculation was carried out on the basis of the prepared models of PUFA-protein complexes. Values of relative energies of formation of the molecule-receptor complex are given. Molecules with the lowest binding energies (i.e., with the highest affinity to target proteins) are highlighted in bold.
Циклооксигеназа-1 / Cyclooxygenase-1 Циклооксигеназа-2 / Cyclooxygenase-2
5-липоксигеназа / 5-lipoxygenase
0,00
ol
I -0,50 al c
w" -1,00
= -1,50 ол
-E -2,00
л а
"2,50
<
-3,00 -3,50
АРК / ARA
ЭПК / EPA
ДГК / DHA
Этиловый эфир ЭПК / EPA ethyl ester ■ Этиловый эфир ДГК / DHA ethyl ester
Рисунок 2. Энергии связывания исследованных молекул полиненасыщенных жирных кислот с ферментами каскада арахидоновой кислоты.
АРК - арахидоновая кислота; ЭПК - эйкозапентаеновая кислота; ДГК - докозагексаеновая кислота
Figure 2. Binding energies of the studied polyunsaturated fatty acids molecules with the enzymes of the arachidonic acid cascade.
ARA - arachidonic acid; EPA - eicosapentaenoic acid; DHA - docosahexaenoic acid
TC ■Ц
Ф
:r
s
^
о ф
T
œ
Ф
о
m
тс s т. го m о
СО _о
о
о
тс ц
ч
о
т.
ф
т га т.
со га
if
(Л
ф X
о ю
ОТ
ю
(Л
о
Е о с о
о ф
о о
ГО ОТ с ^ Е со го ^ л ю о. от
5 ^
> Is-§ +
о. Ф
го
I—
го о
го
I
го т
го
^
о
со g
_û
ю
s _0 I—
го 15
тс
0 ф
ш
I—
ф
1
_0
s
т
>
ц
о с
0
1
*
о
Œ о ê!
£
? о
I -& ГО I CIS
0.2 0.1 H
1 0 га
.3 -0.1 -I
es
J -0.2£ -0.3 -
о
il -0.4-
га ж
-0-5"
^ -0.6 --0.7
0.2
_ 0.1
о
Е
5 0
о
и" -0.1
=1
| -0.3 g -0.4
| -о, -0.6 -0.7
АРК - ЦОГ-1 / ARA - COX-1
Аминокислотные остатки ЦОГ-1 / Amino acid residues of COX-1
Этиловый эфир ДГК - ЦОГ-1 / DHA ethyl ester - COX-1
0.2
l 0.1 ol
l/mo 0
a c
S- -01" -0.2 -
ï -0.3-м
ал/ -0.4
-0.5-- -0.6 -0.7
Рисунок 3. Молекулярно-энергетические профили комплексов полиненасыщенных жирных кислот с циклооксигеназой-1 (ЦОГ-1): а - комплекс арахидоновой кислоты (АРК) с ЦОГ-1; b - комплекс докозагексаеновой кислоты (ДГК) с ЦОГ-1; c - комплекс этилового эфира ДГК с ЦОГ-1
Figure 3. Molecular energy profiles of polyunsaturated fatty acids complexes with cyclooxygenase-1 (COX-1): a - arachidonic acid (ARA) with COX-1 ; b - docosahexaenoic acid (DHA) with COX-1; c - DHA ethyl ester with COX-1
Аминокислотные остатки ЦОГ-1 / Amino acid residues of COX-1
стабильность структур этиловых эфиров ЭПК и ДГК обусловлена взаимодействиями этилового «хвоста» с гидрофобным карманом в структуре ЦОГ (рис. 4). Например, в структуре ЦОГ-1 гидрофобный карман сформирован аминокислотами валин-116, арги-нин-120, аланин-527, лейцин-531. Этиловый «хвост» молекулы эфира увеличивает площадь ее поверхности по сравнению со свободной ПНЖК, тем самым увеличивая площадь поверхности взаимодействия этиловых эфиров ЭПК и ДГК с рецептором. Поскольку взаимодействия гидрофобные, увеличение площади взаимодействия молекул с ЦОГ приводит к увеличению стабильности комплексов «лиганд-ЦОГ».
Арахидонат 5-липоксигеназа: основной фермент биосинтеза лейкотриенов / Arachidonate 5-lipoxygenase: the main enzyme of leukotriene biosynthesis
Арахидонат 5-ЛПОГ (ген Alox5) преобразует ПНЖК в лейкотрие-ны. Липоксигеназный путь синтеза лейкотриенов существует в различных лейкоцитах, включая мастоциты (тучные клетки), эозинофилы, нейтрофилы, моноциты и базофилы. Трансформации АРК посредством 5-ЛПОГ приводят к образованию провоспа-лительных лейкотриенов 4-й группы (LTB4, LTC4, LTD4, LTE4). При действии 5-ЛПОГ на ЭПК синтезируются противовоспалительные лейкотриены 5-й группы (LTB5, LTC5, LTD5, LTE5).
тс ц
ф :г
s
^
о ф т
œ ф
о
ш
тс s т. го ш о
СО _о
о
о
тс ц
ч
о
т.
ф
т га т.
со га
If
(Л
ф X
о ю
ОТ
ю
(Л
о
'1= о с о
о ф
о о
ГО ОТ
с -t
Е СО
ГО ^
л ю
о. от > ^
> Г-§ +
О. Ф
го
I—
го о
го
I
го т
го
^
о
со g
_û
ю
s -û I— ГО
тс
0 ф
ш
I—
ф
1
_о
s
т
>
ц
о г:
о
ï
*
о
Œ о
£
? О I ■в" ГО I CIS
Рисунок 4. Пространственная структура комплекса этилового эфира эйкозапентаеновой кислоты (ЭПК) с циклооксигеназой-2 (ЦОГ-2). Показан гидрофобный карман, в котором расположен этиловый «хвост» молекулы. Структура комплекса ЭПК с ЦОГ-1 практически идентична Figure 4. Spatial structure of the complex of eicosapentaenoic acid (EPA) ethyl ester with cyclooxygenase-2 (COX-2). A hydrophobic pocket is shown, in which the ethyl "tail" of the molecule is located. The structure of the EPA with COX-1 complex is almost identical
Расчеты, представленные в настоящей работе, показывают, что комплексы ю-3 ПНЖК с 5-ЛПОГ намного более стабильны (AAG = -0,95...-1,62 ккал/моль), чем комплекс АРК с 5-ЛПОГ (AAG = -0,69 ккал/моль). Наиболее стабильным был комплекс 5-ЛПОГ с этиловым эфиром ЭПК (AAG = -1,62 ккал/моль). Анализ молекулярно-энергетических профилей комплексов 5-ЛПОГ с ДГК и этиловым эфиром ДГК продемонстрировал, что обе молекулы имеют близкие значения энергии связывания (AAG = -1,21... -1,23 ккал/моль) вследствие сравнительно малого объема связывающего сайта ПНЖК в структуре фермента 5-ЛПОГ.
Меньший объем сайта связывания ДГК приводит к тому, что добавление этилового «хвоста» не только обеспечивает увеличение поверхности взаимодействия с белком-рецептором, но и усиливает стерические препятствия (пространственные напряжения) в структуре фермента. Поэтому увеличение поверхности молекулы этилового эфира ДГК компенсируется возрастанием стерических препятствий, так что не происходит значительного увеличения энергии связывания с 5-ЛПОГ. В то же время объем молекулы ЭПК меньше, чем объем ДГК на 5...10 А3, поэтому этиловый эфир ЭПК образует более стабильный комплекс с 5-ЛПОГ (рис. 5).
ЗАКЛЮЧЕНИЕ / CONCLUSION
В работе представлены результаты биофизического анализа взаимодействий АРК, ЭПК, ДГК и этиловых эфиров ЭПК и ДГК с ЦОГ-1/2 и 5-ЛПОГ. Исследование привело к нескольким выводам, весьма интересным с точки зрения молекулярной фармакологии ЭПК и ДГК.
ЛИТЕРАТУРА:
1. Li W.L., Zhang N.H., Ge S.W., Xu G. Dietary omega-3 fatty acid intake and mortality in CKD population: a 1999-2014 NHANES analysis. Am J Nephrol. 2021; 52 (12): 909-18. https://doi. org/10.1159/000520027.
2. Торшин И.Ю., Громова О.А., Кобалава Ж.Д. О репрессиях ю-3 полиненасыщенных жирных кислот адептами доказательной медицины. ФАРМАКОЭКОНОМИКА. Современная фармако-экономика и фармакоэпидемиология. 2019; 12 (2): 91-114. https:// doi.org/10.17749/2070-4909.2019.12.2.91-114.
Рисунок 5. Пространственная структура комплекса этилового эфира эйкозапентаеновой кислоты с 5-липоксигеназой. Показан гидрофобный карман связывания молекулы омега-3 полиненасыщенной жирной кислоты Figure 5. Spatial structure of the eicosapentaenoic acid ethyl ester complex with 5-lipoxygenase. Hydrophobic binding pocket of omega-3 polyunsaturated fatty acid molecule is shown
Во-первых, ЭПК и ДГК могут образовывать намного более стабильные комплексы с ферментами ЦОГ-1/2 и 5-ЛПОГ, чем АРК. С одной стороны, этот результат соответствует вытеснению АРК из активных центров таргетных ферментов, что приводит к замедлению образования провоспалительных эйкозаноидов. С другой стороны, более стабильные комплексы с таргетными белками соответствуют более интенсивному образованию противовоспалительных эйкозаноидов из ЭПК и ДГК.
Во-вторых, этиловые эфиры ЭПК и ДГК образуют более стабильные комплексы с ферментами каскада АРК, чем свободные ЭПК и ДГК. Результаты биофизического моделирования показали, что это обусловлено увеличением площади поверхности молекулы при присоединении этилового «хвоста» в составе сложных эфиров ЭПК и ДГК.
В-третьих, разница в молекулярном объеме между ЭПК и ДГК (ЭПК - меньший объем, ДГК - больший объем) и их этиловыми эфирами обеспечивает разницу в связывании белками-рецепторами с различными объемами сайта связывания лигандов. Например, более объемный этиловый эфир ДГК лучше связывается с большим сайтом в ферментах ЦОГ-1/2, а этиловый эфир ЭПК с меньшим объемом лучше связывается с меньшим по объему сайтом связывания 5-ЛПОГ. Поэтому препараты, содержащие этиловые эфиры и ЭПК, и ДГК, могут более эффективно взаимодействовать с различными ферментами, чем препараты, содержащие этиловые эфиры только ЭПК или только ДГК.
Большее сродство этиловых эфиров ЭПК и ДГК к молекулам ферментов каскада АРК является одной из фундаментальных основ высокой терапевтической эффективности препаратов ю-3 ПНЖК с высокой степенью стандартизации.
3. Торшин И.Ю., Громова О.А., Зайчик Б.Ц., Ружицкий А.О. Комплексное исследование состава экстрактов жира рыб и количественные критерии для различения стандартизированных экстрактов омега-3 полиненасыщенных жирных кислот. Кардиология. 2020; 60 (5): 47-56. https://doi.org/10.18087/ cardio.2020.5.n1053.
4. Simopoulos A.P. Omega-3 fatty acids in inflammation and autoimmune diseases. J Am Coll Nutr. 2002; 21 (6): 495-505. https://doi. org/10.1080/07315724.2002.10719248.
rc c; ф zï
о ф
T Œ Ф
0
m
к s
1 го m о to .о с; о
о s
к с;
0
1
ф т го
I
со го
I
& о
С ч—
ГО
ÜH
О
с; о
ю
ОТ
ю
Е о с о
о ф
о
0 го от
F ^ Е со го ^ .£= Ю Œ ОТ
1 S « S
Ci Ф
го
го о
го
I
го т го
о
го
Б _û
ю
го
т
>
с; о
EZ
0
1
S о
о X
к
го
° s
£ I
? ф
I— Œ
£ О
55 Ü ? о
ü -а
ГО х
qs
5. Gonzâlez-Périz A., Planagumà A., Gronert K., et al. Docosahexaenoic acid (DHA) blunts liver injury by conversion to protective lipid mediators: protectin D1 and 17S-hydroxy-DHA. FASEB J. 2006; 20 (14): 2537-9. https://doi.org/10.1096/fj.06-6250fje.
6. Khalfoun B., Thibault F., Watier H., et al. Docosahexaenoic and eicosapentaenoic acids inhibit in vitro human endothelial cell production of interleukin-6. Adv Exp Med Biol. 1997; 400B: 589-97.
7. Wada M., DeLong C.J., Hong Y.H., et al. Enzymes and receptors of prostaglandin pathways with arachidonic acid-derived versus eicosapentaenoic acid-derived substrates and products. J Biol Chem. 2007; 282 (31): 22254-66. https://doi.org/10.1074/jbc.M703169200.
8. Serhan C.N., Chiang N., Van Dyke T.E. Resolving inflammation: dual
REFERENCES:
1. Li W.L., Zhang N.H., Ge S.W., Xu G. Dietary omega-3 fatty acid intake and mortality in CKD population: a 1999-2014 NHANES analysis. Am J Nephrol. 2021; 52 (12): 909-18. https://doi.org/10.1159/000520027.
2. Torshin I.Yu., Gromova O.A., Kobalava Zh.D. Concerning the "repression" of ro -3 polyunsaturated fatty acids by adepts of evidence-based medicine. FARMAKOEKONOMIKA. Sovremennaya farmakoekonomika i farma-koepidemiologiya / FARMAKOEKONOMIKA. Modern Pharmacoecono-mics and Pharmacoepidemiology. 2019; 12 (2): 91-114. (in Russ.). https://doi.org/10.17749/2070-4909.2019.12.2.91-114.
3. Torshin I.Yu., Gromova O.A., Zaychik B.Ts., Ruzhitsky A.O. Comprehensive study of the composition of fish fat extracts and quantitative criteria for distinguishing standardized omega-3 polyunsaturated fatty acids extracts. Kardiologiia. 2020; 60 (5): 47-56 (in Russ.). https://doi.org/10.18087/cardio.2020.5.n1053.
4. Simopoulos A.P. Omega-3 fatty acids in inflammation and autoimmune diseases. J Am Coll Nutr. 2002; 21 (6): 495-505. https://doi. org/10.1080/07315724.2002.10719248.
5. Gonzalez-Periz A., Planaguma A., Gronert K., et al. Docosahexaenoic acid (DHA) blunts liver injury by conversion to protective lipid mediators: protectin D1 and 17S-hydroxy-DHA. FASEB J. 2006; 20 (14): 2537-9. https://doi.org/10.1096/fj.06-6250fje.
anti-inflammatory and pro-resolution lipid mediators. Nat Rev Immunol. 2008; 8 (5): 349-61. https://doi.org/10.1038/nri2294.
9. Mori T.A. Omega-3 fatty acids and blood pressure. Cell Mol Biol (Noisy-le-grand). 2010; 56 (1): 83-92.
10. Serhan C.N., Libreros S., Nshimiyimana R. E-series resolvin metabolome, biosynthesis and critical role of stereochemistry of specialized pro-resolving mediators (SPMs) in inflammation-resolution: preparing SPMs for long COVID-19, human clinical trials, and targeted precision nutrition. Semin Immunol. 2022: 101597. https://doi. org/10.1016/j.smim.2022.101597.
11. Torshin I.Yu. Bioinformatics in the post-genomic era: the role of biophysics. NY: Nova Science Publishers; 2013: 255 pp.
6. Khalfoun B., Thibault F., Watier H., et al. Docosahexaenoic and eicosapentaenoic acids inhibit in vitro human endothelial cell production of interleukin-6. Adv Exp Med Biol. 1997; 400B: 589-97.
7. Wada M., DeLong C.J., Hong Y.H., et al. Enzymes and receptors of prostaglandin pathways with arachidonic acid-derived versus eicosapentaenoic acid-derived substrates and products. J Biol Chem. 2007; 282 (31): 22254-66. https://doi.org/10.1074/jbc. M703169200.
8. Serhan C.N., Chiang N., Van Dyke T.E. Resolving inflammation: dual anti-inflammatory and pro-resolution lipid mediators. Nat Rev Immunol. 2008; 8 (5): 349-61. https://doi.org/10.1038/nri2294.
9. Mori T.A. Omega-3 fatty acids and blood pressure. Cell Mol Biol (Noisy-le-grand). 2010; 56 (1): 83-92.
10. Serhan C.N., Libreros S., Nshimiyimana R. E-series resolvin metabolome, biosynthesis and critical role of stereochemistry of specialized pro-resolving mediators (SPMs) in inflammation-resolution: preparing SPMs for long COVID-19, human clinical trials, and targeted precision nutrition. Semin Immunol. 2022: 101597. https://doi. org/10.1016/j.smim.2022.101597.
11. Torshin I.Yu. Bioinformatics in the post-genomic era: the role of biophysics. NY: Nova Science Publishers; 2013: 255 pp.
тс
Ф
:r
s
^
о ф
T
œ
Ф
о
m
тс s x со m о
СО _о
о
о
тс ц
ч
о
х
ф
т га х
со га
it
(Л
ф X
о ю
ОТ
ю
(Л
о
Е о с о
о ф
о о
ГО ОТ с ^ Е со го ^ л ю а. от
5 ^
> Is-% +
Сведения об авторах
Торшин Иван Юрьевич - к.ф-м.н., к.х.н., старший научный сотрудник Института фармакоинформатики ФИЦ «Информатика и управление» РАН. ORCID ID: https://orcid.org/0000-0002-2659-7998; WoS ResearcherlD: C-7683-2018; Scopus Author ID: 7003300274; РИНЦ SPIN-код: 1375-1114.
Громова Ольга Алексеевна - д.м.н., профессор, научный руководитель Института фармакоинформатики ФИЦ «Информатика и управление» РАН. ORCID ID: https://orcid.org/0000-0002-7663-710X; WoS ResearcherlD: J-4946-2017; Scopus Author ID: 7003589812; РИНЦ SPIN-код: 6317-9833. E-mail: unesco.gromova@gmail.com.
About the authors
Ivan Yu. Torshin - PhD (Phys. Math.), PhD (Chem.), Senior Researcher, Institute of Pharmacoinformatics, Federal Research Center "Informatics and Management", RAS (Moscow, Russia). ORCID ID: https://orcid.org/0000-0002-2659-7998; WoS ResearcherID: C-7683-2018; Scopus Author ID: 7003300274; RSCI SPIN-code: 1375-1114.
Olga A. Gromova - Dr. Med. Sc., Professor, Research Supervisor, Institute of Pharmacoinformatics, Federal Research Center "Informatics and Management", RAS (Moscow, Russia). ORCID ID: https://orcid.org/0000-0002-7663-710X; WoS ResearcherID: J-4946-2017; Scopus Author ID: 7003589812; RSCI SPIN-code: 6317-9833. E-mail: unesco.gromova@gmail.com.
« 5
о. <D
ГО I—
ro о
ro
X
ro
T
ro
^
о
CO
g
_û
ю
s _0 I—
со
тс о
Œ ф
m
I—
ф
X
_0 s
T
>
ц
о с
о
X
*
о
Œ о ê!
£
? О X -& ГО х CIS