Using QSARs for predictions in drug delivery

Drug delivery research is an inherently empirical process, however high-throughput approaches could take advantage of understanding drug/material interactions such as from electrostatic, hydrophobic, or other non-covalent interactions between therapeutic molecules and a drug delivery polymer. Cyclodextrin polymers have been investigated for drug delivery specifically due to their capacity to exploit this affinity interaction to change the rate of drug release. Testing drug candidates; however, for affinity is time-consuming, making computational predictions more effective. One option, molecular "docking" programs, provide predictions of affinity, but lack reliability, as their accuracy with cyclodextrin remains unverified experimentally. Alternatively, quantitative structure-activity relationship models (QSARs), which analyze statistical relationships between molecular properties, appear more promising. Previously constructed QSARs for cyclodextrin are not publicly available, necessitating an openly accessible model. Around 600 experimental affinities between cyclodextrin and guest molecules were cleaned and imported from published research. The software PaDEL-Descriptor calculated over 1000 chemical descriptors for each molecule, which were then analyzed in R to create several QSARs with different statistical methods. These QSARs proved highly time efficient, calculating in minutes what docking programs would take hours to accomplish. Additionally, on test sets, QSARs reached R2 values of around 0.7-0.8. The speed, accuracy, and accessibility of these QSARs improve evaluation of individual drugs and facilitate screening of large datasets for potential candidates in cyclodextrin affinity-based delivery systems. An app was built to rapidly access model predictions for end users using the "shiny" library in R. To demonstrate the usability for drug release planning, the QSAR predictions were coupled with a mechanistic model of diffusion within the app. Integrating new modules should provide an accessible approach to use other cheminformatic tools in the field of drug delivery.

[1]  K. A. Connors,et al.  Population characteristics of cyclodextrin complex stabilities in aqueous solution. , 1995, Journal of pharmaceutical sciences.

[2]  R. Tibshirani,et al.  Lasso and Elastic-Net Regularized Generalized Linear Models [R package glmnet version 4.0-2] , 2020 .

[3]  Gerta Rücker,et al.  y-Randomization and Its Variants in QSPR/QSAR , 2007, J. Chem. Inf. Model..

[4]  M. Cordeiro,et al.  Convenient QSAR model for predicting the complexation of structurally diverse compounds with beta-cyclodextrins. , 2009, Bioorganic & medicinal chemistry.

[5]  Danail Bonchev,et al.  Statistical modelling of molecular descriptors in QSAR/QSPR , 2012 .

[6]  H. V. von Recum,et al.  Providing sustained transgene induction through affinity-based drug delivery. , 2016, Journal of biomedical materials research. Part A.

[7]  Kurt Hornik,et al.  Misc Functions of the Department of Statistics, ProbabilityTheory Group (Formerly: E1071), TU Wien , 2015 .

[8]  Y. Inoue,et al.  Complexation Thermodynamics of Cyclodextrins. , 1998, Chemical reviews.

[9]  T. Mano,et al.  Development of machine learning models of β-cyclodextrin and sulfobutylether-β-cyclodextrin complexation free energies. , 2011, International journal of pharmaceutics.

[10]  C. O’Sullivan,et al.  Amperometric detection of antibodies in serum: performance of self-assembled cyclodextrin/cellulose polymer interfaces as antigen carriers. , 2011, Organic & biomolecular chemistry.

[11]  Jahan B. Ghasemi,et al.  QSPR approaches to elucidate the stability constants between β-cyclodextrin and some organic compounds: Docking based 3D conformer , 2016 .

[12]  Jahan B. Ghasemi,et al.  3D-QSAR and docking studies of the stability constants of different guest molecules with beta-cyclodextrin , 2014, Journal of Inclusion Phenomena and Macrocyclic Chemistry.

[13]  Jie Xu,et al.  Quantitative structure–property relationship study of β-cyclodextrin complexation free energies of organic compounds , 2015 .

[14]  G. Nikolic,et al.  In silico prediction of the β-cyclodextrin complexation based on Monte Carlo method. , 2015, International journal of pharmaceutics.

[15]  Takahiro Suzuki,et al.  A Nonlinear Group Contribution Method for Predicting the Free Energies of Inclusion Complexation of Organic Molecules with - and -Cyclodextrins , 2001, J. Chem. Inf. Comput. Sci..

[16]  A. Tropsha,et al.  Beware of q2! , 2002, Journal of molecular graphics & modelling.

[17]  Peter Wolschann,et al.  Predicting Complexation Thermodynamic Parameters of β-Cyclodextrin with Chiral Guests by Using Swarm Intelligence and Support Vector Machines , 2009, International journal of molecular sciences.

[18]  Jahan B. Ghasemi,et al.  An alignment independent 3D-QSAR study for predicting the stability constants of structurally diverse compounds with β-cyclodextrin , 2011 .

[19]  CHUN WEI YAP,et al.  PaDEL‐descriptor: An open source software to calculate molecular descriptors and fingerprints , 2011, J. Comput. Chem..

[20]  M. Shoichet,et al.  Mathematical model accurately predicts protein release from an affinity-based delivery system. , 2015, Journal of controlled release : official journal of the Controlled Release Society.

[21]  Dan C. Fara,et al.  Quantitative Structure-Property Relationship Modeling of beta-Cyclodextrin Complexation Free Energies , 2004, J. Chem. Inf. Model..

[22]  Arthur J. Olson,et al.  AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading , 2009, J. Comput. Chem..

[23]  Dan C. Fara,et al.  Quantitative Structure—Property Relationship Modeling of β-Cyclodextrin Complexation Free Energies. , 2004 .

[24]  H. V. von Recum,et al.  Affinity-based drug delivery. , 2011, Macromolecular bioscience.

[25]  Supratik Kar,et al.  On a simple approach for determining applicability domain of QSAR models , 2015 .

[26]  Alexander Tropsha,et al.  Best Practices for QSAR Model Development, Validation, and Exploitation , 2010, Molecular informatics.

[27]  R. Setzer,et al.  Solvers for Initial Value Problems of Differential Equations(ODE, DAE, DDE) , 2015 .

[28]  Reed B. Jacob,et al.  High-Throughput Virtual Screening Molecular Docking Software for Students and Educators , 2018 .

[29]  Arthur J Olson,et al.  Small-molecule library screening by docking with PyRx. , 2015, Methods in molecular biology.

[30]  G. Saidel,et al.  Experimental Studies and Modeling of Drug Release from a Tunable Affinity-Based Drug Delivery Platform , 2011, Annals of Biomedical Engineering.