Predicting Transdermal Fentanyl Delivery Using Mechanistic Simulations for Tailored Therapy

Transdermal drug delivery is a key technology for administering drugs. However, most devices are “one-size-fits-all”, even though drug diffusion through the skin varies significantly from person-to-person. For next-generation devices, personalization for optimal drug release would benefit from an augmented insight into the drug release and percutaneous uptake kinetics. Our objective was to quantify the changes in transdermal fentanyl uptake with regards to the patient’s age and the anatomical location where the patch was placed. We also explored to which extent the drug flux from the patch could be altered by miniaturizing the contact surface area of the patch reservoir with the skin. To this end, we used validated mechanistic modeling of fentanyl diffusion, storage, and partitioning in the epidermis to quantify drug release from the patch and the uptake within the skin. A superior spatiotemporal resolution compared to experimental methods enabled in-silico identification of peak concentrations and fluxes, and the amount of stored drug and bioavailability. The patients’ drug uptake showed a 36% difference between different anatomical locations after 72 h, but there was a strong interpatient variability. With aging, the drug uptake from the transdermal patch became slower and less potent. A 70-year-old patient received 26% less drug over the 72-h application period, compared to an 18-year-old patient. Additionally, a novel concept of using micron-sized drug reservoirs was explored in silico. These reservoirs induced a much higher local flux (µg cm^-2 h^-1) than conventional patches. Up to a 200-fold increase in the drug flux was obtained from these small reservoirs. This effect was mainly caused by transverse diffusion in the stratum corneum, which is not relevant for much larger conventional patches. These micron-sized drug reservoirs open new ways to individualize reservoir design and thus transdermal therapy. Such computer-aided engineering tools also have great potential for in-silico design and precise control of drug delivery systems. Here, the validated mechanistic models can serve as a key building block for developing digital twins for transdermal drug delivery systems.

Amino Acid and Oligopeptide Effects on Calcium Carbonate Solutions

Biological organisms display sophisticated control of nucleation and crystallization of minerals. In order to mimic living systems, deciphering the mechanisms by which organic molecules control the formation of mineral phases from solution is a key step. We have used computer simulations to investigate the effects of the amino acids arginine, aspartic acid, and glycine on species that form in solutions of calcium carbonate (CaCO₃) at lower and higher levels of supersaturation. This provides net positive, negative, and neutral additives. In addition, we have prepared simulations containing hexapeptides of the amino acids to consider the effect of additive size on the solution species. We find that additives have limited impact on the formation of extended, liquid-like CaCO₃ networks in supersaturated solutions. Additives control the amount of (bi)carbonate in solution, but more importantly, they are able to stabilize these networks on the time scales of the simulations. This is achieved by coordinating the networks and assembled additive clusters in solutions. The association leads to subtle changes in the coordination of CaCO₃ and reduced mobility of the cations. We find that the number of solute association sites and the size and topology of the additives are more important than their net charge. Our results help to understand why polymer additives are so effective at stabilizing dense liquid CaCO₃ phases.

Simulation of Calcium Phosphate Prenucleation Clusters in Aqueous Solution: Association beyond Ion Pairing

Classical molecular dynamics simulations and free energy methods have been used to obtain a better understanding of the molecular processes occurring prior to the first nucleation event for calcium phosphate biominerals. The association constants for the formation of negatively charged complexes containing calcium and phosphate ions in aqueous solution have been computed, and these results suggest that the previously proposed calcium phosphate building unit, [Ca(HPO4)3]4-, should only be present in small amounts under normal experimental conditions. However, the presence of an activation barrier for the removal of an HPO4 2- ion from this complex indicates that this species could be kinetically trapped. Aggregation pathways involving CaHPO4, [Ca(HPO4)2]2-, and [Ca(HPO4)3]4- complexes have been explored with the finding that dimerization is favorable up to a Ca/HPO4 ratio of 1:2.

Prediction of Steam Burns Severity using Raman Spectroscopy on ex vivo Porcine Skin

Skin burns due to accidental exposure to hot steam have often been reported to be more severe than the ones occurring from dry heat. While skin burns due to flames or radiant heat have been thoroughly characterized, the mechanisms leading to steam burns are not well understood and a conundrum still exists: can second degree burns occur without destruction of the epidermis, i.e. even before first degree burns are detected? Skin permeability is dependent both on temperature and on the kinetic energy of incoming water molecules. To investigate the mechanism underlying the injuries related to steam exposure, we used porcine skin as an ex vivo model. This model was exposed to either steam or dry heat before measuring the skin hydration via confocal Raman microspectroscopy. The results show that during the first minute of exposure to steam, the water content in both the epidermis and dermis increases. By analyzing different mechanisms of steam diffusion through the multiple skin layers, as well as the moisture-assisted bio-heat transfer, we provide a novel model explaining why steam burns can be more severe, and why steam can penetrate deeper and much faster than an equivalent dry heat.

Using simulation to understand the structure and properties of hydrated amorphous calcium carbonate

The protocols used give a range of behaviours comparable to previous experiments and an insight into the structure of ACC.