Сетевой фармакологический и биоинформатический анализы основывались на общедоступных базах данных и не предусматривали участия людей или животных. Эксперименты in vitro проводились с использованием установленной клеточной линии и, согласно институциональным требованиям, не требовали этического одобрения.
Скрининг потенциальных компонентов и поиск мишеней
Потенциальные компоненты и соответствующие им мишени были получены из базы данных и аналитической платформы системной фармакологии традиционной китайской медицины (TCMSP) по поисковому запросу «dongchongxiacao» 13 июля 2025 г. Были отобраны соединения с прогнозируемой пероральной биодоступностью не менее 30% и лекарственной пригодностью (drug-likeness) не менее 0,18. Оба критерия скрининга соответствовали семь потенциальных компонентов. Для описательного сравнения из TCMSP были записаны молекулярная масса, пероральная биодоступность, лекарственная пригодность и другие доступные фармакокинетические дескрипторы.
Информация о связях «соединение-мишень» для отобранных компонентов была загружена непосредственно из базы данных TCMSP. При формировании недублированного набора мишеней повторяющиеся записи удалялись вручную, в то время как исходные связи «соединение-мишень» сохранялись для последующего построения сети. Поскольку в TCMSP названия мишеней были представлены в полном виде, а не в форме стандартизированных символов генов, поиск каждой мишени проводился по базе данных UniProt (дата доступа: 14 июля 2025 г.). Выбиралась соответствующая запись для человека, и название мишени преобразовывалось в официальный символ гена. Стандартизированный набор данных «соединение-мишень» использовался для анализа пересечения «заболевание-мишень» и анализа сети «соединение-мишень».
Сбор генов, ассоциированных с гломерулонефритом, и поиск пересечений с мишенями
Гены, связанные с заболеванием, были получены из баз данных GeneCards, DisGeNET, PharmGKB и Therapeutic Target Database по поисковому запросу «glomerulonephritis» 14 июля 2025 года. Для GeneCards были отобраны записи с показателем релевантности выше медианного значения, равного 10. Подходящие записи из других баз данных включались в соответствии с их курируемыми аннотациями «болезнь-мишень». Названия генов были стандартизированы согласно официальным символам генов человека, а дубликаты удалены. В результате этой процедуры было выявлено 1 134 уникальных генов, связанных с гломерулонефритом. Поскольку при поиске использовался более широкий термин «glomerulonephritis», полученный набор данных был рассмотрен как релевантный молекулярным процессам, связанным с хроническим гломерулонефритом, а не специфичный для одного гистопатологического подтипа.
Мишени, связанные с соединениями, полученные из TCMSP, были пересечены с набором генов, связанных с заболеванием, с помощью программы R. Было выявлено двадцать три общих мишени. Диаграммы Венна были построены с использованием пакета R venn, а сеть перекрывающихся мишеней соединений была сконструирована и визуализирована с помощью Cytoscape. Взаимосвязи между потенциальными компонентами и общими мишенями обобщены в Таблице 1.
Анализ функционального обогащения и белок-белковых взаимодействий
Анализ обогащения по Gene Ontology (GO) и Kyoto Encyclopedia of Genes and Genomes (KEGG) проводили в среде R с использованием пакетов clusterProfiler и org.Hs.eg.db. Термины Gene Ontology анализировали по категориям биологических процессов, клеточных компонентов и молекулярных функций. Значения P корректировали для множественного тестирования по методу Бенджамини-Хохберга; были отобраны термины или пути со скорректированными значениями P ниже 0.05. Обогащенные термины и пути представляли в виде независимо построенных столбчатых и точечных диаграмм.
23 перекрывающихся мишени были внесены в базу данных STRING с выбором организма Homo sapiens и минимальным показателем взаимодействия более 0,4. Полученная сеть белок-белковых взаимодействий была импортирована в Cytoscape. Топология сети оценивалась с помощью CytoNCA на основе показателей степени, центральности по посредничеству, центральности по близости, центральности по собственному вектору, локальной средней связности и центральности сети. Узловые мишени определялись как узлы, для которых значения степени, центральности по посредничеству, центральности по близости, центральности по собственному вектору, локальной средней связности и центральности сети превышали соответствующие медианные значения сети. В результате данного анализа было выделено девять узловых мишеней для последующей интерпретации.
Молекулярный докинг и валидация путем редокинга
Молекулярный докинг проводили с использованием структуры циклооксигеназы-2 в комплексе с целекоксибом (PDB ID 3LN1). В качестве рецептора были оставлены цепь A и связанная с ней группа гема, при этом в процессе подготовки рецептора кристаллографические молекулы воды и сокристаллизованный лиганд целекоксиб были удалены. Нативная конформация целекоксиба была извлечена из цепи A отдельно для валидации путем повторного докинга. С помощью AutoDockTools были добавлены полярные атомы водорода и заряды Гастейгера, после чего рецептор был сохранен в формате PDBQT. PDB ID 3LN1 представляет собой структуру мышиной COX-2, связанной с целекоксибом, с разрешением 2.40 Å. Структуры β-ситостерола и линолеил ацетата были получены из PubChem и подготовлены с помощью Open Babel. Были добавлены атомы водорода и парциальные заряды Гастейгера, определены вращающиеся связи, и каждый лиганд был преобразован в формат PDBQT. Докинг целекоксиба, β-ситостерола и линолеил ацетата выполняли с помощью AutoDock Vina с использованием одного и того же рецептора и параметров поиска.
Бокс докинга был центрирован по нативному региону связывания целекоксиба с использованием следующих координат: center_x = 30.9886 center_y = -22.2836 center_z = -16.5072 size_x = 26 Å size_y = 26 Å size_z = 26 Å. Значение полноты поиска (exhaustiveness) было установлено на уровне 16, количество выходных мод — 10, а энергетический диапазон — 3 kcal/mol. Для визуализации была выбрана наиболее высокоранжированная поза с наименьшей расчетной энергией связывания. Для валидации протокола был проведен редокинг нативного целекоксиба при тех же условиях, а среднеквадратичное отклонение (RMSD) тяжелых атомов между кристаллографической и редокированной позами было рассчитано с помощью PyMOL. Значение RMSD ниже 2 Å считалось приемлемым воспроизведением референсной позы. Остатки, расположенные в пределах 4 Å от каждого задокированного лиганда, были определены в PyMOL. Показатели докинга и позы интерпретировались как оценки структурной совместимости, а не как прямое доказательство связывания с мишенью или ингибирования фермента.
Культивирование клеток и приготовление соединений
Гломерулярные мезангиальные клетки крысы HBZY-1 культивировали в модифицированной среде Игла в модификации Дульбекко с высоким содержанием глюкозы, дополненной 10% фетальной бычьей сыворотки и 1% пенициллина-стрептомицина, при 37 °C в увлажненных условиях с 5% CO₂. Для всех экспериментов использовали клетки в пассажах 2–7; перед использованием подтверждалось отсутствие микоплазмы. В качестве химически определенных соединений использовали β-ситостерин и целекоксиб. β-ситостерин растворяли в абсолютном этаноле для приготовления стокового раствора с концентрацией 10 mM (приблизительно 4,15 mg/mL). При необходимости для получения визуально прозрачного раствора проводили кратковременную соникацию и нагрев до 50–60 °C. Целекоксиб растворяли в диметилсульфоксиде для приготовления стокового раствора с концентрацией 10 mM (приблизительно 3,81 mg/mL). Стоковые растворы разделяли на аликвоты и хранили при −20 °C. Перед обработкой стоковые растворы разбавляли непосредственно в предварительно нагретой культуральной среде и немедленно перемешивали.
Концентрации растворителя были согласованы внутри каждого эксперимента. В анализе жизнеспособности клеток конечная концентрация этанола была доведена до 0.2% v/v во всех группах, так как максимальная протестированная концентрация β-sitosterol составила 20 µM. В эксперименте по индуцированию воспаления с помощью LPS конечные концентрации этанола и диметилсульфоксида во всех группах были доведены до 0.1% и 0.01% v/v соответственно. В эксперименте по индуцированию профибротического эффекта с помощью TGF-β1 конечная концентрация этанола была доведена до 0.1% v/v во всех группах. В группы Control, LPS и TGF-β1 добавляли эквивалентные объемы соответствующих растворителей.
Анализ жизнеспособности клеток
Клетки высевали в 96-луночные планшеты с плотностью 5 × 103 клеток на лунку в 100 µL полной среды и оставляли для адгезии на 24 ч. Для оценки диапазона концентраций β-sitosterol клетки инкубировали с 0,5, 1, 5, 10 или 20 µM β-sitosterol в течение 24 ч. Для анализа жизнеспособности клеток при воздействии LPS клетки предварительно обрабатывали β-sitosterol в концентрациях 1, 5 или 10 µM в течение 2 ч, после чего подвергали воздействию 1 µg/mL липополисахарида в течение дополнительных 24 ч.
Группы обработки включали: контрольную группу, β-sitosterol в концентрации 0.5 µM, β-sitosterol в концентрации 1 µM, β-sitosterol в концентрации 5 µM, β-sitosterol в концентрации 10 µM, β-sitosterol в концентрации 20 µM, LPS, LPS с β-sitosterol в концентрации 1 µM, LPS с β-sitosterol в концентрации 5 µM и LPS с β-sitosterol в концентрации 10 µM. В каждую лунку добавляли по 10 µL реагента Cell Counting Kit-8, после чего планшеты инкубировали в течение 2 h. Затем измеряли оптическую плотность при 450 nm. Жизнеспособность клеток определяли путем вычитания значения холостой пробы без клеток и выражали относительно контрольной группы. Для каждого условия использовали шесть технических лунок в каждом из шести независимых биологических экспериментов. Значения технических лунок из каждого эксперимента усредняли для получения одного значения биологического повтора для статистического анализа.
Воспалительная активация, индуцированная LPS
Для экспериментов по воспалительной активации клетки HBZY-1 высевали в 6-луночные планшеты с плотностью 2 × 105 клеток на лунку и подвергали воздействию после достижения примерно 70%–80% конфлюэнтности. Клетки содержали в среде с добавлением 1% фетальной бычьей сыворотки в течение 12 ч до начала обработки. β-Ситостерол или целекоксиб вводили за 2 ч до стимуляции LPS. Затем клетки подвергали воздействию 1 µg/mL LPS в течение 24 ч.
Экспериментальные группы включали контрольную группу, β$\beta$-ситостерол при 10 µM, ЛПС, ЛПС плюс β$\beta$-ситостерол при 1 µM, ЛПС плюс β$\beta$-ситостерол при 5 µM, ЛПС плюс β$\beta$-ситостерол при 10 µM, и LPS с целекоксибом при 1 µMЦелекоксиб был включен в качестве функционального положительного контроля для ингибирования COX-2-зависимого синтеза простагландинов. По завершении обработки собирали супернатанты культур для измерения уровня PGE2. Клетки собирали отдельно для анализа методом количественной ПЦР в реальном времени на Tnf, Il6и Ptgs2 и для вестерн-блот-анализа COX-2.
Профибротическая активация, индуцированная TGF-β1
Для оценки профибротической активации использовали отдельную модель, индуцированную TGF-β1. Клетки HBZY-1 высевали в 6-луночные планшеты с плотностью 2 × 105 клеток на лунку и культивировали в среде с добавлением 1% эмбриональной бычьей сыворотки в течение 12 ч перед обработкой. Клетки предварительно обрабатывали β-ситостеролом в течение 2 ч, после чего стимулировали рекомбинантным TGF-β1 в концентрации 10 ng/mL в течение 48 ч.
Группы лечения включали контрольную группу, β-ситостерол (10 µM), TGF-β1, TGF-β1 + β-ситостерол (5 µM) и TGF-β1 + β-ситостерол (10 µM). По окончании периода обработки клетки собирали для проведения количественной ПЦР в реальном времени для анализа экспрессии Col1a1 и Acta2, а также для вестерн-блот-анализа альфа-гладкомышечного актина (α-SMA).
Количественная ПЦР в реальном времени
Общую РНК выделяли из обработанных клеток с использованием реагента для экстракции РНК на основе фенола и гуанидиния в соответствии с протоколом производителя. Концентрацию и чистоту РНК определяли с помощью спектрофотометра. Для обратной транскрипции использовали только образцы с соотношением A260/A280 от 1.8–2.0. Для каждого образца 1 µg общей РНК подвергали обратной транскрипции в комплементарную ДНК в конечном объеме реакции 20 µL. Полученную комплементарную ДНК перед количественной ПЦР разводили в пять раз водой, свободной от нуклеаз. Количественную ПЦР проводили в общем объеме 20 µL, содержащем 10 µL 2× SYBR Green master mix, по 0.4 µL прямых и обратных праймеров (10 µM), 2 µL разведенной комплементарной ДНК и 7.2 µL воды, свободной от нуклеаз. Условия амплификации включали этап начальной денатурации при 95 °C в течение 30 s, после чего следовали 40 циклов денатурации при 95 °C в течение 5 s и отжига-элонгации при 60 °C в течение 30 s. Затем проводили анализ кривых плавления в диапазоне от 65 °C до 95 °C для проверки специфичности амплификации.
Измерялись уровни экспрессии генов Tnf, Il6, Ptgs2, Col1a1, Acta2 и Gapdh. Ген Gapdh использовали в качестве внутреннего референсного гена. Относительную экспрессию генов определяли методом 2−ΔΔCt, используя группу Control в качестве калибратора. Последовательности праймеров приведены в Дополнительной таблице 1. Каждый образец анализировали в трех технических повторностях; было проведено шесть независимых биологических экспериментов. Значения, полученные в технических повторностях, усредняли перед проведением статистического анализа.
Измерение уровня PGE2 методом иммуноферментного анализа
Супернатанты культур собирали через 24 ч после стимуляции LPS и центрифугировали при 1,000 × g в течение 10 мин при 4 °C для удаления клеточного дебриса. Концентрации PGE2 определяли с помощью имеющегося в продаже набора для иммуноферментного анализа в соответствии с инструкциями производителя.
Стандарты и образцы наносили на планшет для анализа в дубликатах. После завершения этапов связывания антител и проявления цвета регистрировали оптическую плотность при 450 nm с помощью микропланшетного ридера. Концентрации PGE2 определяли по четырехпараметрической логистической стандартной кривой и выражали в pg/mL. Анализировали супернатанты из шести независимых экспериментов с клетками. Перед статистическим анализом значения дубликатов для каждого биологического образца усредняли.
Вестерн-блоттинг
Обработанные клетки лизировали в буфере для радиоиммунопреципитационного анализа с добавлением смеси ингибиторов протеаз. Лизаты инкубировали на льду в течение 30 min и центрифугировали при 12,000 × g в течение 15 min при 4 °C. Концентрацию белка в супернатантах определяли с помощью бицинхонинового метода анализа белков. Равные количества белка (30 µg на лунку) смешивали с буфером для нанесения и денатурировали при 95 °C в течение 5 min. Белки в экспериментах с LPS разделяли на 10% полиакриламидных гелях с додецилсульфатом натрия, тогда как в экспериментах с TGF-β1 использовали 12% гели. Разделенные белки затем переносили на поливинилиденфторидные мембраны.
Мембраны блокировали 5% раствором обезжиренного молока в трис-солевом буфере, содержащем 0.1% Tween 20, в течение 1 ч при комнатной температуре, после чего инкубировали в течение ночи при 4 °C с первичными антителами к COX-2 (1:1000), α-SMA (1:1000) и GAPDH (1:5000). Анализ COX-2 проводили в образцах из эксперимента с LPS, в то время как α-SMA анализировали в образцах из эксперимента с TGF-β1. После промывки мембраны инкубировали с вторичными антителами, конъюгированными с пероксидазой хрена (1:5000), в течение 1 ч при комнатной температуре. Полосы белков визуализировали с помощью усиленного хемилюминесцентного субстрата и фиксировали с использованием системы хемилюминесцентной визуализации. Время экспозиции составило примерно 10 s для COX-2, 10 s для α-SMA и 2 s для GAPDH.
Интенсивность полос количественно определяли с помощью ImageJ. Сигналы COX-2 и α-SMA были нормализованы по соответствующему сигналу GAPDH. Вестерн-блоттинг проводили с использованием образцов белка, полученных в шести независимых экспериментах по культивированию и обработке клеток.
Статистический анализ
Статистический анализ проводили с помощью программы GraphPad Prism. Данные представлены в виде среднего значения ± стандартного отклонения. Для каждого количественного анализа было проведено шесть независимых экспериментов по культивированию и обработке клеток, при этом каждый независимый эксперимент представлял собой одну биологическую повторность. Измерения технических повторов, полученные в рамках одного и того же биологического эксперимента, усреднялись перед статистическим анализом и не рассматривались как независимые наблюдения.
Для анализа жизнеспособности клеток в каждом биологическом эксперименте использовали по шесть технических лунок для каждого условия. Реакции количественной ПЦР проводили в трех технических повторностях для каждого биологического образца, а измерения уровня PGE2 — в двух повторностях. Количественный анализ вестерн-блоттинга основывался на шести независимо подготовленных образцах белка. Таким образом, каждый количественный столбец содержит значения шести индивидуальных биологических повторов.
Различия между несколькими группами анализировали с помощью однофакторного дисперсионного анализа с последующим использованием критерия множественных сравнений Тьюки. Двустороннее значение P < 0.05 считалось статистически значимым. Значения отдельных биологических повторностей представлены вместе со средним значением и стандартным отклонением.