Abstract:Yongwei Gu, 1– 3, * Qing Gu, 4, * Qing Yang, 1, 2 Meng Yang, 3 Shengzhang Wang, 5 Jiyong Liu 1– 3 1 Department of Pharmacy, Fudan University Shanghai Cancer Center, Shanghai 200032, People's Republic of China; 2 Department of Oncology, Shanghai Medical College, Fudan University, Shanghai 200032, People's Republic of China; 3 Department of Pharmacy, Changhai Hospital, Second Military Medical University, Shanghai 200433, People's Republic of China; 4 Department of Pharmacy, Jingan District Zhabei Central Hospital, Shanghai 200070, People's Republic of China; 5 Institute 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 Tel/Fax +86-21-65647825 Email szwang@fudan.edu.cn Jiyong Liu 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 Tel/Fax +86-21-64175590 Email liujiyong@fudan.edu.cn 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 ( K p) 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 K p 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 K ps 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 K p of PAE in SC was 2.00× 10 − 6 cm 2 /h. The K p of PAE-NE in viable epidermis and in dermis detected using penetration test and microdialysis probe technology, was 1.58× 10 − 5 cm 2 /h and 3.20× 10 − 5 cm 2 /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

Penetration pathways, influencing factors and predictive models for dermal absorption of exobiotic molecules: A critical review

New trends in mechanistic transdermal drug delivery modelling: Towards an accurate description of skin microstructure

Methods for evaluating penetration of drug into the skin: A review

Penetration of nanoparticles into human skin.

Characterisation of skin penetration pathways using stimulated Raman scattering microscopy

Epidermal Penetration of Gold Nanoparticles and Its Underlying Mechanism Based on Human Reconstructed 3D Episkin Model

Progress in Topical and Transdermal Drug Delivery Research-Focus on Nanoformulations

An Investigation into the Transdermal Behavior of Active Ingredients by Combination of Experiments and Multiscale Simulations

Penetration profile and human cadaver skin distribution of finasteride from vesicular nanocarriers.

Impact of solvent dry down, vehicle pH and slowly reversible keratin binding on skin penetration of cosmetic relevant compounds: I. Liquids

Finite Element Analysis for Predicting Skin Pharmacokinetics of Nano Transdermal Drug Delivery System Based on the Multilayer Geometry Model

A Data-Mining Approach for the Quantitative Assessment of Physicochemical Properties of Molecular Compounds in the Skin Flux

Skin Structure, Physiology, and Pathology in Topical and Transdermal Drug Delivery

Xenobiotic-metabolizing enzymes in the skin of rat, mouse, pig, guinea pig, man, and in human skin models

Impact of solvent dry down, phase change, vehicle pH and slowly reversible keratin binding on skin penetration of cosmetic relevant compounds: II. Solids

Revisiting techniques to evaluate drug permeation through skin

Mechanism investigation of ethosomes transdermal permeation.

Novel Pharmaceutical Strategies for Enhancing Skin Penetration of Biomacromolecules

Effect of vehicle on the in vitro penetration and metabolism of genistein and daidzein in ex vivo skin explants and the Phenion full‐thickness skin model

Methods for the characterisation of dermal uptake: Progress and perspectives for organophosphate esters

Enhanced Penetration Strategies for Transdermal Delivery