Corneal Permeability and Uptake of Twenty-Five Drugs: Species Comparison and Quantitative Structure–Permeability Relationships

The purpose of this study was to determine corneal permeability and uptake in rabbit, porcine, and bovine corneas for twenty-five drugs using an N-in-1 (cassette) approach and relate these parameters to drug physicochemical properties and tissue thickness through quantitative structure permeability relationships (QSPRs). A twenty-five-drug cassette containing β-blockers, NSAIDs, and corticosteroids in solution at a micro-dose was exposed to the epithelial side of rabbit, porcine, or bovine corneas mounted in a diffusion chamber, and the corneal drug permeability and tissue uptake were monitored using an LC-MS/MS method. Data obtained were used to construct and evaluate over 46,000 quantitative structure–permeability (QSPR) models using multiple linear regression, and the best-fit models were cross-validated by Y-randomization. Drug permeability was generally higher in rabbit cornea and comparable between bovine and porcine corneas. Permeability differences between species could be explained in part by differences in corneal thickness. Corneal uptake between species correlated with a slope close to 1, indicating generally similar drug uptake per unit weight of tissue. A high correlation was observed between bovine, porcine, and rabbit corneas for permeability and between bovine and porcine corneas for uptake (R(2) ≥ 0.94). MLR models indicated that drug characteristics such as lipophilicity (LogD), heteroatom ratio (HR), nitrogen ratio (NR), hydrogen bond acceptors (HBA), rotatable bonds (RB), index of refraction (IR), and tissue thickness (TT) are of great influence on drug permeability and uptake. When data for all species along with thickness as a parameter was used in MLR, the best fit equation for permeability was Log (% transport/cm(2)·s) = 0.441 LogD − 8.29 IR + 8.357 NR − 0.279 HBA − 3.833 TT + 10.432 (R(2) = 0.826), and the best-fit equation for uptake was Log (%/g) = 0.387 LogD + 4.442 HR + 0.105 RB − 0.303 HBA − 2.235 TT + 1.422 (R(2) = 0.750). Thus, it is feasible to explain corneal drug delivery in three species using a single equation.