기술 요약: PILL-CoDe – 자동 미분을 통한 폴리필(Polypill)의 역설계

문제 정의

폴리필(Polypills)—여러 활성 제약 성분(API)과 부형제를 결합한 단일 경구 투여 형태—은 조화로운 다상 방출(multi-phase release)을 통해 개인 맞춤형 고정 용량 조합 요법을 실현할 수 있는 잠재력을 제공한다. 적층 제조(Additive Manufacturing, AM)는 복잡한 다중 재료 기하학적 구조의 물리적 구현을 가능하게 하지만, 기존의 제형 워크플로우는 임의적인 파라미터 스윕(parameter sweeps)에 의존하고 있다. 이러한 방식은 형상, 조성 및 방출 거동의 고차원 설계 공간을 체계적으로 탐색하기에는 부적합하다.

기존의 역설계 접근 방식은 세 가지 주요 한계점을 가진다:

1. 결합되지 않은 최적화(Decoupled Optimization): 이들은 일반적으로 기하학적 형상과 재료 조성을 분리하여 최적화하며, 이로 인해 필(pill)의 외부 형상과 여러 활성 성분 및 부형제의 공간적 분포를 공동 최적화하는 데 실패한다.
2. 단순화된 물리 모델: 많은 모델이 용해(dissolution)를 단순한 기하학적 전면 전파(front-propagation) 문제(Eikonal 방정식에 의해 지배됨)로 취급하며, 이는 국소적 농도 축적, 포화 효과 및 확산 제한 수송과 같은 핵심 현상을 간과한다. 이러한 현상은 고전적 용해 이론(예: Noyes-Whitney 및 Nernst-Brunner 모델)의 중심이다.
3. 계산적 장벽: 복잡하게 결합된 물리를 최적화 루프에 포함시키는 과정은 종종 번거로운 수동 아드조인트 감도(adjoint sensitivities) 유도를 요구하며, 물리적 충실도가 높아질수록 이 과정은 점점 더 까폐해진다.

방법론

저자들은 규정된 약물 방출 역학을 일치시키기 위해 폴리필의 기하학적 형상과 부형제 분포를 동시에 공동 설계(co-design)하는 프레임워크인 PILL-CoDe를 제시한다. 이 프레임워크는 JAX로 구현되어 전체 파이프라인이 엔드 투 엔드(end-to-end)로 미분 가능하다.

1. 하이브리드 설계 표현 (Hybrid Design Representation)

제조 가능성(위상적 연결성)과 설계 자유도(재료 불균질성) 사이의 균형을 맞추기 위해, 본 프레임워크는 하이브리드 파라미터화를 채택한다:

• 기하학적 형상 (Supershapes): 폴리필의 전역적 기하 구조는 슈퍼셰이프(Gielis curves)를 사용하여 파라미터화된다. 이는 결과물인 투여 형태가 하나의 연결된 단일 도메인이 되도록 보장한다. 기하 구조는 이동($c_x, c_y$), 회전($\theta$), 반축 스케일링($a, b$), 사각형도($n$), 그리고 회전 대칭성($m$)의 7개 파라미터로 정의된다.
• 재료 분포 (Neural Implicit): $S$개의 서로 다른 부형제의 내부 분포는 좌표 기반 신경망으로 표현된다. 이 네트워크는 공간 좌표를 입력으로 받아 각 부형제의 부피 분율을 출력한다. Softmax 활성화 함수는 단위 분할 제약 조건($\sum \gamma_s = 1$)을 강제한다. 이 방식은 설계 해상도를 시뮬레이션 메쉬로부터 분리하여 고해상도의 연속적인 재료 필드를 생성할 수 있게 한다.

2. 물리 정보 기반 순방향 모델 (Physics-Informed Forward Model)

용해 과정은 결합된 편미분 방정식(PDE) 시스템에 의해 지배된다:

• 계면 진화: 수정된 Allen-Cahn 위상장(phase-field) 방정식이 움직이는 고체-용매 계면을 모델링한다. 여기에는 Nernst-Brunner 역학에서 유도된 용해 강제 항이 포함되어, 위상의 진화와 국소 용질 농도를 결합한다.
• 질량 수송: **피크 확산(Fickian diffusion)**은 용매 내 용해된 종의 수송을 지배한다. 확산을 위한 소스 항은 위상장 계면과 결합되어 질량 보존을 보장한다.
• 경계 조건: 시스템은 이상적인 in vitro 싱크 조건(경계에서의 농도 0)과 위상장에 대한 제로 플럭스(zero-flux) 조건을 가정한다.

3. 최적화 프레임워크

역설계 문제는 시뮬레이션된 순시 질량 방출률과 규정된 목표 프로파일 사이의 평균 제곱 오차(MSE)를 최소화하는 것으로 정식화된다.

• 제약 조건: 최적화는 다음을 강제한다:
  • 그레이성 제약(Grayness Constraint): 연속적인 재료 분율을 이진 상태(0 또는 1)로 유도하는 지속 전략을 사용하여 물리적 실현 가능성을 확보한다.
  • 부피 제약: 스무스 미니멈(smooth minimum) 공식을 사용하여 특정 부형제의 최소 부피 분율을 강제한다.
  • 박스 제약(Box Constraints): 물리적 유의성을 보장하기 위해 기하학적 파라미터의 범위를 제한한다.
• 감도 분석: 프레임워크는 JAX를 통한 **자동 미분(AD)**을 활용하여 정확한 감도를 계산한다. 과도(transient) 비선형 솔버의 계산 비용을 관리하기 위해 저자들은 다음을 사용한다:
  • Newton-Raphson 반복 과정을 우회하기 위한 내재 함수 정리(Implicit Function Theorem, IFT).
  • 긴 과도 시뮬레이션과 관련된 메모리 병목 현상을 완화하기 위한 체크포인팅(Checkpointing).

주요 기여

본 논문은 네 가지 구체적인 기여를 식별한다:

1. 하이브리드 파라미터화: 전역 기하 구조를 위한 슈퍼셰이프와 내부 다중 재료 분포를 위한 좌표 기반 신경망을 결합한 새로운 설계 공간.
2. 미분 가능한 용해-확산 시뮬레이터: 수정된 Allen-Cahn 위상장 방정식과 Fickian 확산을 결합하여, 전면 전파 모델에는 없는 국소 농도 축적 및 포화 효과를 포착하는 물리 정보 기반 모델.
3. 엔드 투 엔드 역설계 프레임워크: 지배 방정식이 최적화 루프 내에 직접 삽립되며, 감도가 이산화된 솔버를 통해 자동으로 계산되어 수동 아드조인트 유도를 피하는 정식화.
4. 확장된 도달 가능 집합(Reachable Set): 기하학적 형상만을 최적화할 때보다 기하학적 형상과 부형제 분포를 공동 설계하는 것이 더 넓은 범위의 방출 프로파일을 달성할 수 있음을 입증하였으며, 이를 통해 복잡하고 비단조적인 역학을 매칭할 수 있음을 보여줌.

결과

프레임워크는 75×75 메쉬와 2개 층의 신경망을 사용하는 2D 도메인에서의 수치 실험을 통해 검증되었다.

• 가설 검증:
  • 단조 타겟(Monotonic Targets): 단일 재료 최적화(기하학만 고려)는 단순한 단조 방출 프로파일을 성공적으로 매칭했다.
  • 비단조 타겟(Non-Monotonic Targets): 단일 재료 설계는 초기 느린 방출, 급격한 가속, 급격한 감소가 필요한 복잡한 프로파일을 매칭하는 데 실패했다(정규화된 MSE: $3.01 \times 10^{-2}$). 반면, 다중 재료 공동 설계는 이러한 프로파일을 성공적으로 매칭했다(정규화된 MSE: $9.52 \times 10^{-5}$).
  • 메커니즘: 최적화기는 초기 방출을 억제하기 위해 느리게 용해되는 부형제를 주변부에 배치하고, 껍질 침식 시 폭발적 방출을 일으키기 위해 빠르게 용해되는 부형제를 내부부에 배치하였다.
• 초기화 강건성: 프레임워크는 다양한 초기 기하 구조(예: "spike", "circle", "sunflower")로부터 동등한 기능적 성능으로 수렴하였으며, 이는 로컬 미니마에 대한 강건성을 입증하며 설계자가 제조 가능성이나 섭취 용이성에 따라 형상을 선택할 수 있게 함을 보여준다.
• 열화 사례 연구: 프레임키는 우주 비행 조건 등을 시뮬레이션하기 위해 열화된 원료를 포함하는 시나리오에 적용되었다. 신선한, 중간 정도 노화된, 그리고 심하게 노화된 부형제에 대한 최소 부피 제약을 변화시킴으로써, 최적화기는 타겟 방출 프로파일을 맞추기 위해 재료 분포를 성공적으로 조정하였으며, 이는 자원이 제한된 환경에서의 유용성을 입증한다.

의의 및 주장

본 논문은 PILL-CoDe가 제어 방출 폴리필의 체계적 설계를 위한 핵심적인 격차를 해결한다고 주장한다. 단순한 전면 전파 모델을 넘어 결합된 용해-확산 프레임워크로 이동함으로써, 이 방법은 다중 재료 시스템에서 방출을 지배하는 물리적 메커니즘(농도 구배, 포화)을 포착한다.

저자들은 자신들의 접근 방식이 제조 제약 조건 하에서 형상과 조성을 체계적으로 공동 설계할 수 있게 한다고 주장하며, 이는 복잡한 물리 법칙에 대한 감도 유도의 어려움으로 인해 이전에는 불가능했던 능력이다. 결과는 기하학적 형상만으로는 달성할 수 없는 복잡하고 환자 맞춤적인 타겟에 맞춰 방출 역학을 조절하기 위해 다중 재료 역설계가 필수적임을 시사한다.

저자들은 다음과 같은 점을 언급하며 주지의 범위를 겸허하게 제한한다:

• 현재 모델은 이상적인 싱크 조건 하의 용해 지배 영역으로 제한된다.
• 도달 가능한 방출 프로파일 집합은 기하학, 재료 분포, 그리고 물리적 제약(예: 비음수 방출률) 사이의 상호작용에 의해 경계가 정해진다.
• 실제 구현을 위해서는 실험적으로 검증된 재료 데이터셋, 3D 확장, 그리고 적층 제조 제약(예: 최소 특징 크기) 및 in vivo 생리학적 요인의 통합을 포함하는 향후 연구가 필요하다.