Rhythmomimetic drug delivery: modeling, analysis and numerical simulation

TL;DR

This paper models and analyzes a glucose-driven rhythmic drug delivery device based on chemomechanical interactions in a gel membrane, demonstrating oscillatory behavior suitable for hormone therapy applications.

Contribution

It introduces a mathematical model capturing oscillations in a drug delivery device, linking bifurcation analysis to practical hormone release timing.

Findings

Existence of oscillatory solutions in the model.

Frequency of oscillations can be tuned to match hormone release.

Limit cycle analysis reveals the system's rhythmic behavior.

Abstract

We develop, analyse and numerically simulate a model of a prototype, glucose-driven, rhythmic drug delivery device, aimed at hormone therapies, and based on chemomechanical interaction in a polyelectrolyte gel membrane. The pH-driven interactions trigger volume phase transitions between the swollen and collapsed states of the gel.For a robust set of material parameters, we find a class of solutions of the governing system that oscillate between such states, and cause the membrane to rhythmically swell, allowing for transport of the drug, fuel and reaction products across it, and collapse, hampering all transport across it. The frequency of the oscillations can be adjusted so that it matches the natural frequency of the hormone to be released. The work is linked to extensive laboratory experimental studies of the device built by Siegel's team. The thinness of the membrane and its clamped boundary, together with the homogeneously held conditions in the experimental apparatus, justify neglecting spatial dependence on the fields of the problem. Upon identifying the forces and energy relevant to the system, and taking into account its dissipative properties, we apply Rayleigh's variational principle to obtain the governing equations. The material assumptions guarantee the monotonicity of the system and lead to the existence of a three dimensional limit cycle. By scaling and asymptotic analysis, this limit cycle is found to be related to a two-dimensional one that encodes the volume phase transitions of the model. The identification of the relevant parameter set of the model is aided by a Hopf bifurcation study of steady state solutions.