Analiza QSPR leków stosowanych w leczeniu niewydolności nerek i jej powikłań z użyciem indeksów stopnia i zmodyfikowanego odwrotnego stopnia

Dlaczego struktura leków nerkowych ma znaczenie

Niewydolność nerek to narastający globalny problem zdrowotny, a wiele osób z uszkodzonymi nerkami polega na złożonych kombinacjach leków do kontroli ciśnienia krwi, hormonów, minerałów i równowagi płynów. Jednak projektowanie i ulepszanie takich leków jest powolne i kosztowne, jeśli każdą obiecującą cząsteczkę trzeba wytworzyć i przetestować w laboratorium. W badaniu pokazano, jak kształty i połączenia w cząsteczkach leków na przewlekłą chorobę nerek można przekształcić w liczby, które wiarygodnie prognozują kluczowe właściwości fizyczne, pomagając badaczom przesiewać i dopracowywać związki komputerowo, zanim trafią do probówki.

Od narysowania cząsteczki do przewidywania jej zachowania

Chemicy często szkicują cząsteczki jako sieci atomów połączonych wiązaniami. W tej pracy te szkice traktowane są jak grafy matematyczne: atomy stają się węzłami, a wiązania — krawędziami. Autorzy koncentrują się na dziewiętnastu rzeczywistych lekach stosowanych w zarządzaniu przewlekłą chorobą nerek i jej powikłaniami, w tym lekach obniżających poziom cukru we krwi, cholesterolu, ciśnienia krwi lub parathormonu oraz lekach leczących anemię lub wysoki poziom fosforanów we krwi. Dla każdego leku przekształcają rysunek chemiczny w graf, a następnie obliczają kilka opisowych liczb, które uchwytują stopień połączeń i rozgałęzienia cząsteczki. Te deskryptory, znane jako indeksy topologiczne, służą jako strukturalne odciski palców, które można porównywać między wieloma związkami.

Proste reguły z dalekosiężnym wpływem

Kluczową ideą w badaniu jest „stopień” atomu, czyli po prostu liczba wiązań z nim związanych. W oparciu o to badacze definiują rodziny indeksów opartych na stopniu, które podsumowują wzorce, takie jak częstotliwość łączenia atomów o określonej łączności. Wprowadzają też „zmodyfikowane indeksy odwrotnego stopnia”, które przemapowują te łączności za pomocą prostej reguły kontrolowanej przez niewielki parametr. Poprzez regulację tego parametru generują z tego samego grafu molekularnego kilka pokrewnych odcisków, akcentujących różne aspekty układu cząsteczki. Korzystając z metody partycjonowania krawędzi i oprogramowania do algebry komputerowej, systematycznie obliczają te indeksy dla wszystkich dziewiętnastu leków nerkowych.

Łączenie struktury z właściwościami rzeczywistymi

Aby sprawdzić, czy te odciski rzeczywiście są użyteczne, autorzy porównują je z mierzonymi właściwościami fizykochemicznymi pobranymi z publicznych baz danych. Należą do nich masa molowa, powierzchnia zdolna do interakcji z wodą, liczba ciężkich atomów, ogólna złożoność strukturalna, temperatura wrzenia, jak silnie cząsteczka załamuje światło (molowa refraktywność), jak łatwo ulega polaryzacji elektronów (polaryzowalność) oraz zajmowana przez nią objętość molowa. Następnie dopasowują trzy typy modeli statystycznych — liniowy, sześcienny krzywoliniowy i logarytmiczny — między każdym indeksem a każdą właściwością. Siła dopasowania jest mierzona współczynnikami korelacji, które wskazują, jak dobrze formuła oparta na indeksie potrafi odtworzyć dane eksperymentalne.

Wyszukiwanie najlepszych numerycznych odcisków

Analiza pokazuje, że niektóre indeksy wyraźnie przewyższają inne. Tradycyjny indeks oparty na stopniu zwany przeformułowanym Zagreb-1 przewiduje liczbę ciężkich atomów niemal idealnie za pomocą prostego równania liniowego. Dla większości pozostałych właściwości lepiej sprawdzają się jednak krzywoliniowe relacje sześcienne niż linie proste czy formy logarytmiczne. W tych przypadkach szczególnie wyróżniają się zmodyfikowane indeksy odwrotnego stopnia. Przy określonym ustawieniu parametru modyfikacji (l = 2) indeks łączności atom–wiązanie ściśle śledzi masę molową, a indeks harmoniczno–geometryczny oddaje, jaka część powierzchni cząsteczki może wchodzić w interakcje z wodą. Podobnie inne zmodyfikowane indeksy odwrotnego stopnia przy różnych wartościach parametru najlepiej opisują złożoność, temperaturę wrzenia i objętość molową, podczas gdy indeks arytmetyczno–geometryczny oparty na stopniu skutecznie łączy strukturę z molową refraktywnością i polaryzowalnością.

Co to oznacza dla przyszłych terapii nerek

Dla czytelnika niezwiązanego ze specjalistyczną dziedziną wniosek jest taki, że autorzy zbudowali zestaw matematycznych skrótów: patrząc tylko na to, jak atomy są połączone w proponowanym leku nerkowym, mogą dokładnie oszacować wiele istotnych cech fizycznych bez pomiarów laboratoryjnych. Wśród testowanych narzędzi zmodyfikowane indeksy odwrotnego stopnia — szczególnie przy jednej preferowanej wartości parametru — okazały się najbardziej wszechstronne, dając najlepsze prognozy dla większości właściwości. Modele te nie odnoszą się bezpośrednio do choroby nerek, ale mogą znacznie przyspieszyć wczesne etapy odkrywania i optymalizacji leków, koncentrując wysiłki eksperymentalne na najbardziej obiecujących kandydatach i ostatecznie pomagając nowym terapiom dla niewydolności nerek szybciej trafić do pacjentów.

Cytowanie: Godlin, J.J.J., Radha, S. QSPR analysis of the drugs used to treat renal failure and its complications using degree and modified reverse degree indices. Sci Rep 16, 8889 (2026). https://doi.org/10.1038/s41598-026-37586-9

Słowa kluczowe: leki w przewlekłej chorobie nerek, modelowanie struktury molekularnej, indeksy topologiczne, predykcja QSPR, projektowanie leków na niewydolność nerek