Предсказание и скрининг потенциальных мишеней PLB
Информация о химической структуре PLB была получена из базы данных PubChem (https://pubchem.ncbi.nlm.nih.gov/). Молекулярные мишени PLB были предсказаны на основе лиганд-ориентированного двухмерного и трехмерного сходства с использованием баз данных SwissTarget (https://swisstargetprediction.ch/index.php), SEA (https://sea.bkslab.org/), TargetNet (http://targetnet.scbdd.com/calcnet/index/), PharmMapper (https://www.lilab-ecust.cn/pharmmapper/index.html), Comparative Toxicogenomics Database (CTD, https://ctdbase.org/), Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP, https://www.tcmsp-e.com/index.php) и HERB (http://herb.ac.cn/). Для обеспечения надежности и воспроизводимости данных применялись следующие специфические для каждой базы данных пороги фильтрации: в TCMSP были сохранены мишени с биодоступностью при пероральном приеме (OB) ≥ 30% и показателем лекарственно-подобности (DL) ≥ 0.18; в SwissTargetPrediction были оставлены записи с вероятностью предсказания ≥ 0.5; в PharmMapper были отобраны мишени с нормализованным показателем соответствия (fit score) ≥ 4.0; в SEA в качестве статистически значимых были включены только мишени с E-value < 0.01; в TargetNet были сохранены мишени с вероятностью предсказания > 0.5; в CTD были включены только мишени с курируемыми доказательствами (выведенными из взаимодействий химическое вещество–ген) и баллом взаимодействия > 0.3; в HERB были сохранены мишени с баллом подтверждения в литературе ≥ 0.4. Все идентификаторы белков, полученные из вышеуказанных баз данных, были приведены к официальным символам генов человека HGNC с использованием базы данных UniProt (https://www.uniprot.org/) с ограничением по виду Homo sapiens. После удаления дубликатов из всех источников для последующего анализа было получено в общей сложности 50 уникальных мишеней, связанных с PLB.
Поиск генов, ассоциированных с раком предстательной железы
Гены, связанные с раком предстательной железы, были получены из баз данных GeneCards (https://www.GeneCardss.org/), DrugBank (https://go.drugbank.com/), CTD и HERB. Для обеспечения достоверности данных применялись специфические для каждой базы критерии включения. В GeneCards были отобраны гены с показателем релевантности ≥ 0.5, так как этот порог позволяет выделить гены с умеренными или сильными доказательствами их связи с искомым заболеванием. В DrugBank были включены только записи, имеющие экспериментальное подтверждение (например, одобренные FDA или исследуемые препараты для лечения рака предстательной железы), и исключены вычислительно предсказанные или теоретические взаимодействия. В CTD были сохранены только записи с уровнями доказательности, классифицируемыми как «маркер» или «механизм», на основе курируемых взаимодействий «химическое вещество — ген — болезнь». В HERB были включены мишени с подтвержденным литературой показателем достоверности ≥ 0.4 для обеспечения достаточных экспериментальных данных или результатов текстового анализа. Все символы генов были приведены к номенклатуре HGNC с помощью базы данных UniProt (https://www.uniprot.org/), вид был ограничен Homo sapiens, а дубликаты были удалены, что позволило получить 1 19 уникальных мишеней, связанных с раком предстательной железы, для последующего анализа.
Идентификация перекрывающихся мишеней и коррекция аннотаций конфликтов
Пересечение 50 мишеней PLB и 1 19 мишеней рака предстательной железы было выполнено с использованием стандартизированных символов HGNC в программе Venny 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/index.html), в результате чего было выявлено 151 перекрывающихся потенциальных мишень. Для разрешения конфликтов аннотаций между базами данных был внедрен трехэтапный конвейер коррекции: (i) несоответствия псевдонимов генов были унифицированы с помощью инструмента сопоставления ID UniProt; (ii) избыточные паралогичные гены были исключены с использованием CD‑HIT с порогом сходства последовательностей > 0.4; и (iii) противоречивые функциональные аннотации сохранялись только в том случае, если они были подтверждены как минимум двумя независимыми базами данных, в то время как уникальные конфликтующие описания из одного источника отбраковывались.
Построение сети и идентификация хабов
Сеть белок-белковых взаимодействий (PPI) была построена путем ввода перекрывающихся мишеней в базу данных STRING (версия 12.0, https://cn.string-db.org/) с ограничением поиска по виду Homo sapiens. Для обеспечения высокой надежности были сохранены только те взаимодействия, совокупный показатель достоверности которых составлял ≥ 0.70; все узлы, не имеющие связей, были исключены из сети. Построенная PPI-сеть была визуализирована и проанализирована с точки зрения топологии в программе Cytoscape (версия 3.10.2), после чего с помощью плагина cytoHubba были определены хабовые гены, при этом в качестве основного алгоритма ранжирования использовался метод MCC (Maximal Clique Centrality), а в качестве вспомогательной метрики перекрестной проверки — Degree.
Анализ функционального обогащения
Функциональная аннотация согласно Gene Ontology (GO) и анализ обогащения путей Kyoto Encyclopedia of Genes and Genomes (KEGG) для дифференциально экспрессируемых генов были выполнены с использованием Bioinformatics Online Analysis Platform (https://www.bioinformatics.com.cn, дата последнего доступа: 4 мая 2026 г.). Платформа использует точный тест Фишера для расчета необработанных P-значений, а для коррекции множественного тестирования с целью контроля частоты ложноположительных результатов (FDR) применялся метод Бенджамини — Хохберга (BH). Термины GO (биологический процесс, клеточный компонент и молекулярная функция) и пути KEGG со значением FDR < 0.05 считались статистически значимо обогащенными.
Молекулярный докинг
Кристаллические структуры выбранных белков-мишеней, представляющие активные конформации, были получены из RCSB Protein Data Bank (https://www.rcsb.org/). Подготовка белков проводилась в программе Discovery Studio и включала удаление молекул воды и гетероатомов, удаление избыточных цепей, восстановление недостающих остатков, добавление атомов водорода при pH 7.4, присвоение зарядов Гастейгера и минимизацию энергии. Связывающий карман был определен на основе координат совместно кристаллизованного лиганда с использованием сетки, охватывающей все критические остатки, участвующие в связывании субстрата. Молекулярный докинг PLB выполнялся с помощью веб-сервера SwissDock (http://www.swissdock.ch/) с использованием модуля Attracting Cavities. Для оценки аффинности рассчитывались свободные энергии связывания, и конформация с наименьшей энергией для каждой мишени была выбрана в качестве окончательной позы докинга. Результаты моделирования молекулярного докинга были визуализированы для их валидации.
Моделирование молекулярной динамики
Моделирование молекулярной динамики (МД) проводили с помощью GROMACS 202.2. Для описания белка использовали силовое поле Amber14SB, сольватация системы осуществлялась в воде TIP3P. Параметры плюмбагина, включая частичные заряды AM1-BCC и типы атомов GAFF2, были сгенерированы с помощью Antechamber с последующим преобразованием топологии с помощью ACPYPE и назначением параметров ионов по Юнгу-Читему. Каждый комплекс белок-лиганд помещали в усеченный додекаэдрический расчетный бокс с минимальным расстоянием от белка до границы бокса 1.2 nm, сольватировали молекулами воды TIP3P и нейтрализовали добавлением 0.15 M Na⁺/Cl⁻. После минимизации энергии методом наискорейшего спуска (Fmax < 1,00 kJ·mol⁻1·nm⁻1) систему подвергали последовательному NVT- и NPT-уравновешиванию по 20 ps каждое при 298 K. Затем проводили продуктивное моделирование в течение 20 ns в условиях NPT с шагом интегрирования 2 fs, используя схему обрезания Verlet. Дальнодействующие электростатические взаимодействия рассчитывали методом PME (particle mesh Ewald), при этом для электростатических взаимодействий и взаимодействий Ван-дер-Ваальса использовали радиус обрезания 1.2 nm. Для ограничения связей, содержащих водород, применяли алгоритм LINCS. Температуру поддерживали на уровне 298 K с помощью термостата Нозе-Гувера, а давление — на уровне 1 bar с помощью баростата Парринелло-Рамана. Координаты сохраняли каждые 10 ps для последующего анализа. Анализ траекторий и визуализацию выполняли с помощью утилит GROMACS, VMD и PyMOL, а расчеты свободной энергии связывания методом MM-PBSA проводили с помощью gmx_MMPBSA в соответствующих случаях.
Оценка уравновешивания траекторий молекулярной динамики и протокол расчета свободной энергии методом MM-PBSA
Для каждого комплекса белок–плумбагин с помощью GROMACS было проведено полноатомное производство молекулярно-динамического моделирования продолжительностью 20 ns. Временная точка 10 ns была установлена в качестве границы уравновешивания: сегмент 0–10 ns был определен как фаза конформационной релаксации, в течение которой остов белка и карман связывания лиганда подвергались непрерывным конформационным перестройкам; сегмент 10–20 ns был идентифицирован как термодинамически стабильное плато, о чем свидетельствовало отсутствие однонаправленного дрифта RMSD, радиуса гирации (Rg), порезидульного RMSF, скрытой SASA лиганда, межмолекулярных водородных связей и энергий взаимодействия лиганда с рецептором, которые демонстрировали лишь незначительные стационарные флуктуации. Все количественные кинетические параметры и свободные энергии связывания MM/PBSA рассчитывались исключительно на основе равновесной фазы 100–200 ns, из которой кадры равномерно извлекались каждые 10 ps, что дало 1 0 равновесных снимков для каждой системы в качестве входных данных для gmx_MMPBSA. Первые 10 ns траекторий релаксации были отсечены, чтобы исключить влияние конформационного дрифта на расчеты свободной энергии.
Прогностический анализ целевых генов при раке предстательной железы
Прогностический анализ проводился с использованием наборов данных по раку предстательной железы из The Cancer Genome Atlas (TCGA). Данные секвенирования РНК (STAR-counts) и соответствующая клиническая информация были получены с портала Genomic Data Commons (https://portal.gdc.cancer.gov). Количество экспрессированных генов было переведено в единицы транскриптов на миллион (TPM) и нормализовано с использованием log2(TPM + 1). После исключения образцов с неполной клинической информацией в анализ было включено 498 случаев. Медианный уровень экспрессии каждого гена использовался для разделения пациентов на группы с высокой и низкой экспрессией. Для оценки общей выживаемости (OS) для SRC и выживаемости без прогрессирования (PFS) для MAPK3 был проведен анализ выживаемости Каплана-Мейера с использованием лог-рангового критерия и одномерной регрессии пропорциональных рисков Кокса; были рассчитаны отношения рисков (HRs) и 95% доверительные интервалы (CIs). Прогностическая эффективность дополнительно оценивалась путем построения временно-зависимых рабочих характеристик приемника (ROC-кривых) для периодов в 1, 3 и 5 лет. Статистический анализ выполнялся в программе R версии 4.0.3, при этом значение P < 0.05 считалось статистически значимым.