Gemini 4 Carbon: odtwórz referencyjny zestaw pKa do odkrywania leków z surowych danych laboratoryjnych
Eigent
Od surowych danych laboratoryjnych do referencyjnego zestawu do odkrywania leków
pKa to wartość pH, przy której cząsteczka przyłącza lub oddaje proton. Wpływa na to, jak lek się rozpuszcza, przenika przez błony i wiąże się ze swoim celem molekularnym, dlatego jest jedną z pierwszych właściwości, jakie mierzą badacze leków.
Sam pomiar to dla aparatu rutyna. Przekształcenie surowych zapisów w wiarygodne wartości referencyjne to już praca eksperta. Wyszkolony analityk koryguje odczyty elektrody, dopasowuje widma UV w całym zakresie pH, przelicza pomiary w mieszaninach metanol–woda z powrotem na czystą wodę i decyduje, którym przebiegom można ufać. Przebieg po przebiegu wymaga to mnóstwa czasu i staranności.
Dlatego oddaliśmy analizę Gemini 4 Carbon w Eigent. Jeden brief i 97 surowych plików z aparatu – każdy obliczony wynik usunięty, bez internetu, bez podpowiedzi – zamieniły się w kompletny referencyjny zestaw pKa dla 35 związków lekopodobnych. 99.7% pojedynczych pomiarów było poprawnych: to najlepszy wynik spośród wszystkich modeli, jakie dotąd przetestowaliśmy na tym zadaniu.
Przepływ pracy w skrócie
Jeden brief poprowadził całość od surowych plików do gotowego zestawu referencyjnego. Każdy etap zostawił coś, co mogliśmy sprawdzić.
| Etap | Co się wydarzyło | Co powstało |
|---|---|---|
| Dekodowanie | Rozszyfrowanie nieudokumentowanego formatu plików aparatu i kalibracji jego elektrody | pH i skorygowane pH dla każdego odczytu |
| Dopasowanie | Dopasowanie widm UV w całym zakresie pH, przy wielu długościach fali jednocześnie | pKa dla każdego z 291 miareczkowań |
| Ekstrapolacja | Przeniesienie wyników z mieszanin metanol–woda z powrotem na czystą wodę | pKa w wodzie dla każdego z 97 przebiegów |
| Ocena | Akceptacja lub odrzucenie przebiegów na podstawie dowodów – wytrącenie, błędne widma, niefizyczne nachylenia | Średnie dla związków wyłącznie z wiarygodnych przebiegów |
| Dostarczenie | Zapisanie odpowiedzi i skryptu, który ją wygenerował | Zestaw referencyjny i jeden skrypt analityczny w Pythonie |
| Sprawdzenie | Ocena odpowiedzi i potwierdzenie, że przebieg pozostał offline | Wyniki na każdym poziomie i ślad audytowy |
Oto jak Gemini 4 Carbon to zrobił – łącznie z tym, gdzie nie dał rady, i dlaczego.
Co dostarczył
- Kompletny referencyjny zestaw pKa – wartość pKa i zawartość metanolu dla każdego miareczkowania, pKa w wodzie dla każdego przebiegu oraz średnia, odchylenie standardowe i lista zaakceptowanych przebiegów dla każdego związku.
- Analiza w postaci kodu – jeden skrypt w Pythonie, który odtwarza całą odpowiedź z surowych plików. Po ponownym uruchomieniu daje tę samą odpowiedź, identyczną co do bajtu.
- Pisemne podsumowanie metody – jak skalibrowano odczyty, jak dopasowano widma, jak ekstrapolowano przebiegi oraz które przebiegi pominięto i dlaczego.
Konfiguracja
- Eigent 1.0.7, w jednej sesji w trybie Single Agent, z Gemini 4 Carbon jako modelem.
- Lokalny folder zawierający wyłącznie surowe pliki z automatycznego titratora, listę związków wraz z ich strukturami, stałe fizyczne dla mieszanin metanol–woda oraz środowisko Pythona z typowymi bibliotekami naukowymi i RDKit.
- Offline i samodzielnie. Brief wykluczał korzystanie z internetu, a w trakcie przebiegu nikt nie odpowiadał na pytania. Każda decyzja – któremu przebiegowi zaufać, którego ustawienia się trzymać – należała do agenta.
- Wyłącznie surowe sygnały. Każdą wartość obliczoną przez oprogramowanie aparatu usunięto, więc nic w plikach nie zdradzało odpowiedzi.
- Żadnych wskazówek co do metody. Brief opisywał pomiary i to, co należy podać, ale nie sposób ich analizy. Wybór odpowiedniej chemii fizycznej był częścią zadania.
Jak pracował Gemini 4 Carbon
Gemini 4 Carbon realizował zadanie w terminalu jako pojedynczy agent, bez udziału człowieka. Przebieg podzielił się na pięć faz.
-
Rozpoznanie. Model przeczytał listę związków i stałe kalibracyjne, sparsował XML aparatu dla wszystkich typów oznaczeń i opublikował pierwszą aktualizację statusu: "now analyzing instrument recordings, medium composition, and electrode/UV equations."
-
Inżynieria wsteczna aparatu. Ustalił, jak titrator przelicza miliwolty z elektrody na pH i jak czteroparametrowa kalibracja elektrody przekłada tę wartość na rzeczywiste stężenie jonów wodorowych. Z dozowanych objętości wyprowadził algebrę stojącą za jedną ze stałych składu –
a = 1 + b·(0.2 + 0.025/V)– zgodną z wartością w pliku z dokładnością do 3 × 10⁻⁵. Następnie przejrzał ustawienia analityka we wszystkich 97 plikach: okna pH i długości fali, limity zmętnienia oraz liczbę przejść. -
Budowa i testowanie analizy. Zbudował dopasowania widm dla wielu długości fali jednocześnie, wypróbował je na jednym związku, a następnie uruchomił na wszystkich 291 miareczkowaniach. Poprawkę na rozcieńczenie sprawdził bezpośrednio na plateau pojedynczych form, a za pomocą RDKit sklasyfikował każdy związek jako kwas, zasadę lub amfolit. Znalazł też każdy skan, w którym analityk ręcznie nadpisał ustawienie – i zastosował się do tych decyzji.
-
Kontrola jakości jak w laboratorium. Druga aktualizacja statusu: "now examining edge cases, turbid sweeps, multiprotic zwitterions, and outlier runs." Co wychwycił:
- C21_4 – związek wytrącił się z wody. Widma tego przebiegu zgadzały się z pozostałymi przebiegami tylko na poziomie 0.20, więc wyłączył go ze średniej dla związku.
- C27_1 – fiolka o niskiej absorbancji z błędnym widmem, wykryta dzięki osobnemu dopasowaniu dwóch pasm UV.
- C29_2 – niefizyczne nachylenie ekstrapolacji (+200 wobec +98 i +114 w przebiegach siostrzanych), więc pominął ten przebieg.
- Ujemna absorbancja przy 500 nm – jej przyczyną okazał się pęcherzyk powietrza w widmie referencyjnym.
Przeanalizował chronologię eksperymentów w laboratorium i powtórne przebiegi ze stycznia–lutego 2022, aby zrozumieć, dlaczego niektóre związki mają dodatkowe przebiegi. Przetestował też własną metodę na niezależnych danych: porównał ekstrapolację ważoną i nieważoną na 15 związkach, które zmierzono również bezpośrednio w wodzie. Wariant ważony miał mniejszy błąd – 0.0604 wobec 0.0623 – więc przy nim pozostał.
-
Dostarczenie. Napisał jeden skrypt analityczny liczący 451 wierszy, uruchomił go i zweryfikował schemat, typy i zakresy wyników.
Łącznie: 76 wywołań modelu i 85 poleceń w terminalu, bez ani jednego nieudanego wywołania i bez udziału człowieka.
Wynik
W folderze pojawiły się dwie rzeczy: zestaw referencyjny – 291 miareczkowań, 97 przebiegów i 35 związków – oraz skrypt, który go zbudował. Końcowa wiadomość Gemini 4 Carbon podsumowała metodę – skład medium i kalibrację elektrody, dopasowanie widm UV przy wielu długościach fali, ekstrapolację do wody i selekcję przebiegów – oraz wymieniła każdy pominięty przebieg wraz z uzasadnieniem.
Jak wypadł w ocenie
Odpowiedź oceniliśmy względem zestawu referencyjnego, poza folderem agenta. Ocena jest rygorystyczna: każde miareczkowanie, przebieg i związek to osobny test typu zaliczony/niezaliczony, a miareczkowanie zalicza test tylko wtedy, gdy jego pKa mieści się w granicach niepewności statystycznej dopasowania referencyjnego, a zawartość metanolu jest poprawna. Przebieg zalicza test tylko wtedy, gdy zaliczą go wszystkie trzy jego miareczkowania, a związek – tylko wtedy, gdy zalicza go każdy użyty przebieg, więc jeden błąd może pociągnąć za sobą kolejne.
| Poziom | Poprawne |
|---|---|
| Miareczkowania | 290 / 291 |
| Przebiegi | 91 / 97 |
| Związki | 32 / 35 |
To 99.7% poprawnych pojedynczych pomiarów, każdy w granicach tolerancji statystycznej dopasowania referencyjnego.
Sprawdziliśmy też sam przebieg, przeglądając wywołania narzędzi w Dev mode. W Eigent 1.0.7 sesje Single Agent zawsze mają narzędzia do wyszukiwania w sieci, przeglądarki i pobierania stron, więc ta kontrola ma znaczenie: przebieg, który cokolwiek wyszukał, się nie liczy.
| Kontrola | Co ustaliliśmy |
|---|---|
| Nie użył żadnych danych zewnętrznych | Żadnych wywołań wyszukiwania w sieci, przeglądarki ani pobierania stron. Żadne z 85 poleceń w terminalu nie sięgnęło do sieci, instalatora pakietów ani żadnej ścieżki poza folderem zadania. |
| Nie naruszył danych wejściowych | Dane wejściowe są identyczne co do bajtu z oryginalnymi plikami. |
| Odpowiedź to jego własna analiza | Ponowne uruchomienie skryptu odtwarza odpowiedź identyczną co do bajtu. |
Gdzie zawiódł – i dlaczego
Gemini 4 Carbon nie zaliczył wszystkich testów. Wszystkie błędy wynikają z dwóch wyborów metodycznych:
- Jak ważyć ekstrapolację. Model ważył ekstrapolację z mieszanin metanol–woda niepewnością pomiaru – ten wybór odbiega od konwencji przyjętej w referencji. Kosztowało to 5 przebiegów i 3 związki. Gemini 4 Carbon dokonał tego wyboru na podstawie dowodów – dopasowanie ważone miało mniejszy błąd względem związków zmierzonych w wodzie – a brief nie określał metody.
- Której skali pH użyć. Zastosował zapisane przez analityka okno pH na innej skali pH niż referencja. Kosztowało to 1 miareczkowanie i 1 przebieg.
Żaden z tych wyborów nie jest luką w samej analizie. Oba to konwencje, które rozpozna ekspert w tej dziedzinie – taki rodzaj wewnętrznej reguły, którą laboratorium zwykle wpisuje do swoich instrukcji.
Był to też pojedynczy przebieg. Wyniki w tak trudnym zadaniu mogą się różnić między przebiegami.
Czego się nauczyliśmy
- Surowe dane wymuszają uczciwość agenta. Skoro wartości obliczone przez oprogramowanie producenta usunięto, każda wartość musiała wynikać z sygnałów – a agent musiał pokazać swoją pracę.
- Opisz pomiar, nie metodę. Gemini 4 Carbon sam wybrał chemię fizyczną i uzasadnił każdy wybór dowodami z danych.
- Kod jest produktem końcowym. Skrypt, który odtwarza odpowiedź, naukowiec może ponownie uruchomić, przeczytać i przejrzeć.
- Błędy to konwencje, a nie pomyłki. Tam, gdzie zawiódł, dokonał uzasadnionego wyboru, który rozstrzygnęłyby wewnętrzne reguły laboratorium.
Surowe pliki z aparatury nie muszą czekać w kolejce na analizę. Przekaż je Gemini 4 Carbon w Eigent – a potem oceń rozumowanie, nie tylko wynik.
Podziękowania: nagranie powstało z użyciem Gemini 4 Carbon i Eigent. Muzyka: "Dance With You", Everet Almond.
