Искусственный интеллект в разработке лекарств

Материал из MachineLearning.

(Различия между версиями)

Версия 21:32, 16 июля 2026

Содержание

1 Математическая репрезентация химических данных
2 Графовые нейронные сети (GNN) и виртуальный скрининг
3 Генеративный дизайн молекул (De Novo Design)
- 3.1 Вариационные автоэнкодеры (VAE) и латентная оптимизация
- 3.2 Эквивариантные диффузионные модели (E(n)-EDM)
4 Предсказание 3D-структуры белков: Архитектура AlphaFold 2
5 От in silico к клинике: связь с доказательной медициной
6 Ограничения и вызовы
7 Литература
8 См. также

Статья написана с использованием LLM и проверена участником Dmitrii Vishovan 01:29, 17 июля 2026 (MSD)

Искусственный интеллект в разработке лекарств (Шаблон:Lang-en, дизайн in silico) — это междисциплинарная область на стыке машинного обучения, хемоинформатики и структурной биологии. Целью направления является автоматизация и ускорение процесса поиска, генерации и тестирования новых фармакологически активных соединений с использованием алгоритмов ИИ.

Традиционный цикл разработки нового лекарственного препарата (от поиска мишени до выхода на рынок) занимает в среднем 10–12 лет и стоит более 2,5 млрд долларов США, при этом более 90 % препаратов-кандидатов проваливаются на стадии клинических испытаний. Применение искусственного интеллекта позволяет радикально сократить пространство поиска, состоящее из $10^{60}$ возможных малых молекул, предсказывать свойства соединений до их физического синтеза и генерировать молекулы de novo, не имеющие аналогов в природе.

Математическая репрезентация химических данных

Фундаментальной проблемой применения методов машинного обучения в хемоинформатике является задача векторизации — преобразования дискретной химической структуры в непрерывное математическое представление (эмбеддинг), пригодное для дифференцируемой оптимизации. В современной вычислительной химии выделяют три базовые парадигмы репрезентации молекул:

Строковые линейные нотации (1D-представления): Наиболее распространенным стандартом является SMILES (Simplified Molecular-Input Line-Entry System) и его структурно-строгая альтернатива SELFIES (Self-Referencing Embedded Strings). Молекула описывается как одномерная последовательность символов, полученная путем обхода молекулярного графа (обычно в глубину). Данная парадигма позволяет напрямую адаптировать мощный аппарат NLP (в частности, архитектуру Трансформер) для задач языкового моделирования химии. Главный математический недостаток SMILES — проблема нестрогой инъективности: одна и та же молекула может иметь множество валидных строковых представлений, что требует применения канонизации графа или алгоритмов аугментации данных при обучении.
Молекулярные графы (2D-топология): Представляют собой топологически строгую абстракцию, где молекула описывается как неориентированный атрибутированный граф $G = (\mathcal{V}, \mathcal{E})$ . Множество вершин $\mathcal{V}$ соответствует атомам, а множество ребер $\mathcal{E}$ — ковалентным связям. Каждому узлу $v_i$ сопоставляется начальный вектор признаков $h_i \in \mathbb{R}^d$ (включающий атомарный номер, тип гибридизации, частичный заряд по Гастейгеру), а каждому ребру $e_{ij}$ — свойства химической связи. В отличие от строк, графовое представление инвариантно к перестановкам индексов атомов (permutation invariance), что делает его естественным входом для графовых нейронных сетей (GNN).
Пространственные репрезентации (3D-конформации): Описывают молекулу через набор непрерывных декартовых координат атомов $\mathbf{R} \in \mathbb{R}^{N \times 3}$ в трехмерном пространстве (например, в виде облаков точек или воксельных сеток). Поскольку физические взаимодействия (связывание лиганда с карманом белка-мишени) определяются пространственной комплементарностью и силами Ван-дер-Ваальса, 3D-представления критически важны для молекулярного докинга и оценки свободной энергии связывания. Обучение на таких данных требует использования архитектур, обладающих свойством $E(3)$ -эквивариантности (устойчивости предсказаний к трансляциям, ротациям и отражениям координатного пространства), таких как эквивариантные графовые сети (EGNN).

Графовые нейронные сети (GNN) и виртуальный скрининг

Виртуальный скрининг — это вычислительная процедура фильтрации сверхбольших химических библиотек (свыше $10^9$ соединений) для идентификации молекул-кандидатов (хитов), обладающих высокой вероятностью связывания с целевой биологической мишенью.

Современным стандартом, пришедшим на смену классическим методам QSAR (количественные соотношения «структура-свойство») на базе фиксированных молекулярных отпечатков (например, Morgan fingerprints), являются Графовые нейронные сети (GNN). В отличие от жестких векторов, GNN обучают репрезентации графа $G = (\mathcal{V}, \mathcal{E})$ сквозным образом (end-to-end), используя парадигму передачи сообщений (Message Passing Neural Networks, MPNN).

Архитектура Message Passing

В молекулярных графах критически важны не только свойства атомов, но и свойства связей (одинарная, двойная, ароматическая). На каждом слое $l$ скрытое представление (эмбеддинг) атома $h_v$ обновляется путем агрегации «сообщений» от его соседей $\mathcal{N}(v)$ с обязательным учетом признаков ребер $e_{vu}$ .

Математически шаг передачи сообщения формализуется следующим образом:

$m_v^{(l+1)} = \sum_{u \in \mathcal{N}(v)} M^{(l)} \left( h_v^{(l)}, h_u^{(l)}, e_{vu} \right)$

$h_v^{(l+1)} = U^{(l)} \left( h_v^{(l)}, m_v^{(l+1)} \right)$

Где:

$M^{(l)}$ — обучаемая функция вычисления сообщений (обычно многослойный персептрон, MLP).
$m_v^{(l+1)}$ — агрегированное сообщение (в качестве функции агрегации чаще всего используется сумма, так как она позволяет нейросети различать изоморфные графы лучше, чем функции mean или max).
$U^{(l)}$ — функция обновления узла (часто реализуется через рекуррентные блоки, такие как GRU, или механизмы остаточных связей ResNet).

Операция Readout

После $L$ итераций обмена сообщениями (где $L$ определяет радиус рецептивного поля графа) локальные атомные эмбеддинги необходимо свернуть в единый вектор молекулы $h_G$ . Это выполняется с помощью функции Readout, которая обязана быть инвариантной к перестановкам индексов атомов (Permutation Invariant):

$h_G = R \left( \{ h_v^{(L)} \mid v \in G \} \right)$

Для этой цели применяются как простые пулинги (Sum/Mean Pooling), так и параметризованные механизмы глобального внимания (Global Attention Pooling или Set2Set), которые позволяют сети «фокусироваться» на конкретных функциональных группах молекулы, ответственных за её биологическую активность.

Предсказание ADMET

Вектор $h_G$ передается на полносвязные слои классификатора или регрессора. В индустрии фармакологии профилирование ADMET (Absorption, Distribution, Metabolism, Excretion, Toxicity) традиционно формулируется как задача многозадачного обучения (Multi-Task Learning, MTL). Одна GNN обучается предсказывать сразу десятки свойств: растворимость в воде (logS), липофильность (logP), проницаемость гематоэнцефалического барьера и ингибирование ферментов печени (цитохромы P450). Разделение параметров на нижних слоях (Shared Representation) позволяет сети извлекать более устойчивые химические закономерности.

Исторический прецедент: Ярким примером мощи MPNN стало открытие антибиотика галицина (Halicin) исследователями из MIT (Stokes et al., 2020). Направленная графовая сеть (Directed MPNN) была обучена на наборе из 2335 молекул с бинарными метками (ингибирует ли молекула рост E. coli). Затем модель была применена для виртуального скрининга библиотеки Repurposing Hub (107 млн соединений), выявив структурно уникальную молекулу с мощнейшим бактерицидным эффектом, которая не была бы найдена классическими методами эвристического поиска.

Генеративный дизайн молекул (De Novo Design)

В отличие от виртуального скрининга, ограниченного конечными размерами предсинтезированных библиотек, генеративные алгоритмы решают обратную задачу (inverse design) — направленное конструирование молекулярных структур de novo (с нуля) путем стохастического сэмплирования из выученного распределения химического пространства.

Вариационные автоэнкодеры (VAE) и латентная оптимизация

В пионерской работе Gómez-Bombarelli и соавт. (2018) была предложена архитектура вариационного автоэнкодера для одномерных представлений молекул (SMILES). Математически VAE выполняет вероятностное отображение дискретной топологии молекулы $x$ в непрерывное низкоразмерное латентное многообразие (latent manifold).

Энкодер моделирует приближенное апостериорное распределение $q_\phi(z|x)$ , параметризуя среднее $\mu$ и дисперсию $\sigma^2$ многомерного гауссовского распределения, в то время как декодер $p_\theta(x|z)$ вероятностно восстанавливает исходную структуру. Обучение сети сводится к максимизации вариационной нижней оценки (Evidence Lower Bound, ELBO), где дивергенция Кульбака-Лейблера (KL-divergence) выступает в роли регуляризатора, принуждая латентное пространство быть непрерывным и гладким.

Оптимизация (Latent Space Optimization): Фундаментальное преимущество такого подхода заключается в переходе от сложной комбинаторной оптимизации на дискретных графах к дифференцируемой оптимизации в $\mathbb{R}^d$ . Параллельно с VAE обучается суррогатная нейросеть-предиктор $f(z)$ , оценивающая физико-химические свойства. Применяя методы байесовской оптимизации или градиентного подъема (gradient ascent) по поверхности $f(z)$ , алгоритм находит новый латентный вектор $z^*$ , максимизирующий целевую функцию (например, аффинность или растворимость), после чего декодирует его в совершенно новую, структурно валидную молекулу.

Эквивариантные диффузионные модели (E(n)-EDM)

Доминирующей парадигмой генерации молекулярных конформаций в последние годы стали вероятностные диффузионные модели (Denoising Diffusion Probabilistic Models, DDPM). В отличие от 1D/2D-генерации, современные архитектуры работают непосредственно с 3D-координатами атомов $\mathbf{R} \in \mathbb{R}^{N \times 3}$ и их категориальными признаками $\mathbf{h}$ (типами элементов).

Процесс генерации формализуется как марковская цепь. Прямой диффузионный процесс постепенно разрушает структуру молекулы, добавляя гауссовский шум к координатам атомов согласно заданному расписанию дисперсии. Обучаемая нейронная сеть (как правило, графовая) аппроксимирует функцию оценки (score function) распределения данных — $\nabla_{\mathbf{R}} \log p_t(\mathbf{R})$ . Решение обратного стохастического дифференциального уравнения (SDE) позволяет итеративно «очищать» шум (denoising), формируя молекулу из случайного облака точек.

При генерации молекулы непосредственно внутри рецепторного «кармана» белка-мишени (Structure-Based Drug Design), архитектура модели обязана обладать свойством эквивариантности по отношению к группе евклидовых преобразований $E(3)$ . Это гарантирует, что при любом вращении или параллельном переносе белка в пространстве, сгенерированная молекула-лиганд повернется синхронно, сохраняя физическую стабильность конформации и термодинамически выгодную свободную энергию связывания (Binding Free Energy, $\Delta G$ ).

Предсказание 3D-структуры белков: Архитектура AlphaFold 2

Белки — это макромолекулы, трехмерная (3D) структура которых строго детерминирует их биологическую функцию и форму активного центра (рецепторного «кармана»). Задача алгоритмического предсказания 3D-конформации белка по его 1D-аминокислотной последовательности (проблема фолдинга) оставалась нерешенной более 50 лет из-за колоссального комбинаторного пространства возможных состояний (парадокс Левинталя).

В 2020 году система AlphaFold 2 от Google DeepMind совершила фундаментальный прорыв, решив эту задачу с атомарной точностью (среднее отклонение предсказанных координат от реальных составило менее 1 Ангстрема). Архитектура AlphaFold 2 представляет собой сложный сквозной (end-to-end) пайплайн, который перевел задачу фолдинга из плоскости симуляции молекулярной динамики в плоскость глубокого обучения на графах и многообразиях.

Двухколейная архитектура: Evoformer

Ключевая инновация модели — архитектурный блок Evoformer (Evolutionary Transformer), который заменяет классические сверточные сети двухколейным (dual-track) механизмом непрерывного обмена информацией. На вход модуль получает два потока:

Множественное выравнивание последовательностей (Multiple Sequence Alignment, MSA): 2D-матрица, где собраны эволюционно родственные белки из генетических баз данных. Она позволяет сети улавливать коэволюцию: если два аминокислотных остатка часто мутируют синхронно в процессе эволюции организмов, они, с высокой вероятностью, находятся в физическом контакте в свернутом белке.
Матрица парных взаимодействий (Pair Representation): Тензор размерности $N \times N \times d$ (где $N$ — длина белка), кодирующий гипотетические расстояния между всеми парами остатков.

Внутри Evoformer эти представления обновляют друг друга. Для обработки матрицы расстояний применяется механизм Triangle Attention (треугольное мультипликативное обновление). Он математически принуждает сеть соблюдать правило неравенства треугольника (расстояние $d(i, j) \le d(i, k) + d(k, j)$ ), отсекая геометрически невозможные конформации графа на самых ранних этапах логического вывода.

Структурный модуль (Structure Module) и IPA

После 48 слоев Evoformer абстрактные тензоры передаются в Структурный модуль, который выполняет прямой рендеринг 3D-координат.

Вместо предсказания абсолютных декартовых координат, модуль оперирует локальными системами отсчета. Каждый аминокислотный остаток моделируется как твердое тело (rigid body), описываемое вектором трансляции и матрицей поворота в группе $SE(3)$ . Обновление координат происходит с помощью механизма Invariant Point Attention (IPA). Это специализированный механизм внимания, который по своей физической сути инвариантен к глобальным вращениям и сдвигам всего белка в пространстве.

Для градиентной оптимизации используется строгая метрика — функция потерь FAPE (Frame Aligned Point Error). Она вычисляет расстояния между атомами не в глобальной системе координат, а в локальных системах отсчета каждого остатка, штрафуя модель исключительно за искажение относительной геометрии химических связей.

Влияние на дизайн лекарств

За создание алгоритма разработчикам AlphaFold была присуждена Нобелевская премия по химии 2024 года. В контексте Drug Discovery алгоритм совершил революцию: он открыл доступ к сверхточному виртуальному скринингу и генеративному дизайну лигандов для сотен тысяч белков-мишеней (включая рецепторы орфанных заболеваний), для которых ранее не существовало экспериментальных рентгеноструктурных или криоэлектронных моделей. Это окончательно замкнуло цикл разработки лекарств in silico.

От in silico к клинике: связь с доказательной медициной

Разработка молекулы с помощью алгоритмов генеративного ИИ — это лишь первый этап. После успешного прохождения доклинических испытаний (на клетках и животных) препарат выходит на стадию клинических исследований на людях (R&D pipeline).

На этом этапе возникает потребность в разработке систем-компаньонов — прогностических моделей, которые будут предсказывать, ответит ли конкретный пациент на новое ИИ-лекарство (персонализированная медицина). Для оценки того, стоит ли внедрять такой прогностический алгоритм в реальную практику, традиционно используется Анализ кривых решений (Decision Curve Analysis, DCA). DCA позволяет математически доказать, что тестирование популяции перед назначением нового препарата приносит пациентам большую «чистую пользу» (Net Benefit), чем стратегии слепого назначения дорогостоящего лекарства всем подряд или полного отказа от терапии. Таким образом, цепочка ИИ замыкается: от генерации молекул с помощью графовых сетей до обоснования клинической полезности с помощью биостатистики.

Ограничения и вызовы

Синтезируемость: Генеративные модели часто предлагают молекулы со сложной геометрией, которые невозможно или экономически невыгодно синтезировать в реальной химической лаборатории. В современные loss-функции внедряются штрафы за сложность синтеза (Synthetic Accessibility Score).
Обобщение вне распределения (OOD): Модели, обученные на известных химических базах (например, ZINC или ChEMBL), плохо предсказывают свойства совершенно новых классов соединений (out-of-distribution generalization).
Нехватка качественных данных: Данные in vivo (результаты на людях и животных) крайне скудны, зашумлены и часто скрыты коммерческой тайной фармацевтических компаний.

Литература

Stokes J.M. et al. A Deep Learning Approach to Antibiotic Discovery // Cell. — 2020. — Т. 180. — № 4. — С. 688-702.
Gómez-Bombarelli R. et al. Automatic Chemical Design Using a Data-Driven Continuous Representation of Molecules // ACS Central Science. — 2018. — Т. 4. — № 2. — С. 268-276.
Jumper J. et al. Highly accurate protein structure prediction with AlphaFold // Nature. — 2021. — Т. 596. — С. 583-589.
Meyers J. et al. De novo design of targeted protein degraders with generative diffusion models // Nature Computational Science. — 2023.
Walters W.P., Barzilay R. Critical assessment of AI in drug discovery // Expert Opinion on Drug Discovery. — 2020. — Т. 16. — № 9. — С. 937-947.

См. также

Источник — «http://www.machinelearning.ru/wiki/index.php?title=%D0%98%D1%81%D0%BA%D1%83%D1%81%D1%81%D1%82%D0%B2%D0%B5%D0%BD%D0%BD%D1%8B%D0%B9_%D0%B8%D0%BD%D1%82%D0%B5%D0%BB%D0%BB%D0%B5%D0%BA%D1%82_%D0%B2_%D1%80%D0%B0%D0%B7%D1%80%D0%B0%D0%B1%D0%BE%D1%82%D0%BA%D0%B5_%D0%BB%D0%B5%D0%BA%D0%B0%D1%80%D1%81%D1%82%D0%B2»

Категории: Машинное обучение в биологии | Глубокое обучение | Прикладные задачи машинного обучения

@@ Строка 59: / Строка 59: @@
 При генерации молекулы непосредственно внутри рецепторного «кармана» белка-мишени (Structure-Based Drug Design), архитектура модели обязана обладать свойством '''эквивариантности''' по отношению к группе евклидовых преобразований <tex>E(3)</tex>. Это гарантирует, что при любом вращении или параллельном переносе белка в пространстве, сгенерированная молекула-лиганд повернется синхронно, сохраняя физическую стабильность конформации и термодинамически выгодную свободную энергию связывания (Binding Free Energy, <tex>\Delta G</tex>).
-== Предсказание структуры белков (AlphaFold) ==
+== Предсказание 3D-структуры белков: Архитектура AlphaFold 2 ==
-Белки — это молекулярные машины, трехмерная структура которых определяет их функцию (и то, какое лекарство к ним подойдет). Задача предсказания 3D-структуры белка по его аминокислотной последовательности (фолдинг) оставалась нерешенной более 50 лет.
-В 2020 году система '''AlphaFold 2''' от компании Google DeepMind совершила революцию, решив эту задачу с атомарной точностью. Модель использует архитектуру [[Трансформер (модель)|Evoformer]], которая одновременно обрабатывает эволюционные паттерны (множественные выравнивания последовательностей) и пространственные графы остатков. За эту разработку создателям была присуждена Нобелевская премия по химии в 2024 году. Знание точной структуры белков открыло путь для применения генеративных моделей ИИ к мишеням, для которых ранее не существовало рентгеноструктурных данных.
+Белки — это макромолекулы, трехмерная (3D) структура которых строго детерминирует их биологическую функцию и форму активного центра (рецепторного «кармана»). Задача алгоритмического предсказания 3D-конформации белка по его 1D-аминокислотной последовательности (проблема фолдинга) оставалась нерешенной более 50 лет из-за колоссального комбинаторного пространства возможных состояний (парадокс Левинталя).
+В 2020 году система '''AlphaFold 2''' от Google DeepMind совершила фундаментальный прорыв, решив эту задачу с атомарной точностью (среднее отклонение предсказанных координат от реальных составило менее 1 Ангстрема). Архитектура AlphaFold 2 представляет собой сложный сквозной (end-to-end) пайплайн, который перевел задачу фолдинга из плоскости симуляции молекулярной динамики в плоскость глубокого обучения на графах и многообразиях.
+=== Двухколейная архитектура: Evoformer ===
+Ключевая инновация модели — архитектурный блок '''Evoformer''' (Evolutionary Transformer), который заменяет классические сверточные сети двухколейным (dual-track) механизмом непрерывного обмена информацией. На вход модуль получает два потока:
+# '''Множественное выравнивание последовательностей (Multiple Sequence Alignment, MSA):''' 2D-матрица, где собраны эволюционно родственные белки из генетических баз данных. Она позволяет сети улавливать коэволюцию: если два аминокислотных остатка часто мутируют синхронно в процессе эволюции организмов, они, с высокой вероятностью, находятся в физическом контакте в свернутом белке.
+# '''Матрица парных взаимодействий (Pair Representation):''' Тензор размерности <tex>N \times N \times d</tex> (где <tex>N</tex> — длина белка), кодирующий гипотетические расстояния между всеми парами остатков.
+Внутри Evoformer эти представления обновляют друг друга. Для обработки матрицы расстояний применяется механизм '''Triangle Attention''' (треугольное мультипликативное обновление). Он математически принуждает сеть соблюдать правило неравенства треугольника (расстояние <tex>d(i, j) \le d(i, k) + d(k, j)</tex>), отсекая геометрически невозможные конформации графа на самых ранних этапах логического вывода.
+=== Структурный модуль (Structure Module) и IPA ===
+После 48 слоев Evoformer абстрактные тензоры передаются в Структурный модуль, который выполняет прямой рендеринг 3D-координат.
+Вместо предсказания абсолютных декартовых координат, модуль оперирует локальными системами отсчета. Каждый аминокислотный остаток моделируется как твердое тело (rigid body), описываемое вектором трансляции и матрицей поворота в группе <tex>SE(3)</tex>. Обновление координат происходит с помощью механизма '''Invariant Point Attention (IPA)'''. Это специализированный механизм внимания, который по своей физической сути инвариантен к глобальным вращениям и сдвигам всего белка в пространстве.
+Для градиентной оптимизации используется строгая метрика — функция потерь '''FAPE''' (Frame Aligned Point Error). Она вычисляет расстояния между атомами не в глобальной системе координат, а в локальных системах отсчета каждого остатка, штрафуя модель исключительно за искажение относительной геометрии химических связей.
+=== Влияние на дизайн лекарств ===
+За создание алгоритма разработчикам AlphaFold была присуждена Нобелевская премия по химии 2024 года. В контексте Drug Discovery алгоритм совершил революцию: он открыл доступ к сверхточному виртуальному скринингу и генеративному дизайну лигандов для сотен тысяч белков-мишеней (включая рецепторы орфанных заболеваний), для которых ранее не существовало экспериментальных рентгеноструктурных или криоэлектронных моделей. Это окончательно замкнуло цикл разработки лекарств ''in silico''.
 == От in silico к клинике: связь с доказательной медициной ==