Finite Element Analysis for Predicting Skin Pharmacokinetics of Nano Transdermal Drug Delivery System Based on the Multilayer Geometry Model
Authors Gu Y, Gu Q, Yang Q, Yang M, Wang S, Liu J
Received 8 May 2020
Accepted for publication 30 July 2020
Published 12 August 2020 Volume 2020:15 Pages 6007—6018
Checked for plagiarism Yes
Review by Single anonymous peer review
Peer reviewer comments 2
Editor who approved publication: Dr Linlin Sun
Yongwei Gu,1– 3,* Qing Gu,4,* Qing Yang,1,2 Meng Yang,3 Shengzhang Wang,5 Jiyong Liu1– 3
1Department of Pharmacy, Fudan University Shanghai Cancer Center, Shanghai 200032, People’s Republic of China; 2Department of Oncology, Shanghai Medical College, Fudan University, Shanghai 200032, People’s Republic of China; 3Department of Pharmacy, Changhai Hospital, Second Military Medical University, Shanghai 200433, People’s Republic of China; 4Department of Pharmacy, Jingan District Zhabei Central Hospital, Shanghai 200070, People’s Republic of China; 5Institute of Biomechanics, Department of Aeronautics and Astronautics, Fudan University, Shanghai 200433, People’s Republic of China
*These authors contributed equally to this work
Correspondence: Shengzhang Wang
Department of Aeronautics and Astronautics, Institute of Biomechanics, Fudan University, Shanghai 200433, People’s Republic of China
Department of Pharmacy, Fudan University Shanghai Cancer Center, Shanghai 200032, People’s Republic of China; Department of Oncology, Shanghai Medical College, Fudan University, Shanghai 200032, People’s Republic of China
Background: Skin pharmacokinetics is an indispensable indication for studying the drug fate after administration of transdermal drug delivery systems (TDDS). However, the heterogeneity and complex skin structured with stratum corneum, viable epidermis, dermis, and subcutaneous tissue inevitably leads the drug diffusion coefficient (Kp) to vary depending on the skin depth, which seriously limits the development of TDDS pharmacokinetics in full thickness skin.
Methods: A multilayer geometry skin model was established and the Kp of drug in SC, viable epidermis, and dermis was obtained using the technologies of molecular dynamics simulation, in vitro permeation experiments, and in vivo microdialysis, respectively. Besides, finite element analysis (FEA) based on drug Kps in different skin layers was applied to simulate the paeonol nanoemulsion (PAE-NEs) percutaneous dynamic penetration process in two and three dimensions. In addition, PAE-NEs skin pharmacokinetics profile obtained by the simulation was verified by in vivo experiment.
Results: Coarse-grained modeling of molecular dynamic simulation was successfully established and the Kp of PAE in SC was 2.00× 10− 6 cm2/h. The Kp of PAE-NE in viable epidermis and in dermis detected using penetration test and microdialysis probe technology, was 1.58× 10− 5 cm2/h and 3.20× 10− 5 cm2/h, respectively. In addition, the results of verification indicated that PAE-NEs skin pharmacokinetics profile obtained by the simulation was consistent with that by in vivo experiment.
Discussion: This study demonstrated that the FEA combined with the established multilayer geometry skin model could accurately predict the skin pharmacokinetics of TDDS.
Keywords: nano transdermal drug delivery system, skin pharmacokinetics, finite element analysis, multilayer geometry model, diffusion coefficient, paeonol nanoemulsion