ADME Evaluation in Drug Discovery. 3. Modeling Blood-Brain Barrier Partitioning Using Simple Molecular Descriptors

In this paper, QSPR models were developed for in vivo blood-brain partitioning data (logBB) of a large data set consisting of 115 diverse organic compounds. The best model is based on three descriptors: n-octanol/water partition coefficient calculated using the SLOGP approach, logP; high-charged polar surface areas based on the Gasteiger partial charges, HCPSA, and the excessive molecular weight larger than 360, MW(360). The model bears good statistical significance, n = 78, r = 0.88, q = 0.86, s = 0.36, F = 81.5. The actual prediction potential of the model was validated through two external validation sets of 37 diverse compounds. The predicted results demonstrate that the model bears better prediction potential than many other models and can be used for logBB estimations for drug and drug-like molecules. Comparison of several logP calculation approaches suggests that logP calculated by SLOGP can be used as a significant descriptor for the prediction of molecular transport properties because SLOGP gives the most similar results with CLOGP. The QSPR model indicates that larger polar surface areas have a more negative contribution to logBB, but the absolute partial charges on the atoms surrounded by the polar surfaces should be larger than 0.10|e|. Meanwhile, tight junction membranes limit the size of hydrophilic molecules that can cross the membrane with a molecular weight of approximately 360, because when a molecule's weight is larger than 360 it shows a negative contribution to logBB. The computations of molecular surface, partial charges, logP, and logBB have been accomplished using a program called Drug-BB. Moreover, to improve the efficiency of the computations of logP, we made an extensive reparametrization of SLOGP, and the newly developed SLOGP model is only based on simple atomic addition. Further, we developed a set of parameters to calculate the topological polar surface area (TPSA), thus the high-charged topological polar surface area (HCTPSA) could be estimated from the 2D connection information of a molecule. Adopting the new strategies, the estimations of logP, HCTPSA, and logBB are only based on the topological structure of a molecule and therefore, can be used for fast screening of virtual libraries having millions of molecules.

[1]  F. Lombardo,et al.  Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. , 2001, Advanced drug delivery reviews.

[2]  D. E. Clark,et al.  Rapid calculation of polar molecular surface area and its application to the prediction of transport phenomena. 2. Prediction of blood-brain barrier penetration. , 1999, Journal of pharmaceutical sciences.

[3]  M. Yazdanian,et al.  In vitro blood-brain barrier permeability of nevirapine compared to other HIV antiretroviral agents. , 1998, Journal of pharmaceutical sciences.

[4]  J. Taskinen,et al.  Relationship between immobilised artificial membrane chromatographic retention and the brain penetration of structurally diverse drugs. , 1997, Journal of pharmaceutical and biomedical analysis.

[5]  Bernard Testa,et al.  Lipophilicity and hydrogen-bonding capacity of H1-antihistaminic agents in relation to their central sedative side-effects , 1994 .

[6]  M. Feher,et al.  A simple model for the prediction of blood-brain partitioning. , 2000, International journal of pharmaceutics.

[7]  Tingjun Hou,et al.  ADME Evaluation in Drug Discovery. 2. Prediction of Partition Coefficient by Atom-Additive Approach Based on Atom-Weighted Solvent Accessible Surface Areas , 2003, J. Chem. Inf. Comput. Sci..

[8]  F. Lombardo,et al.  Computation of brain-blood partitioning of organic solutes via free energy calculations. , 1996, Journal of medicinal chemistry.

[9]  Yiannis N. Kaznessis,et al.  Prediction of blood-brain partitioning using Monte Carlo simulations of molecules in water , 2001, J. Comput. Aided Mol. Des..

[10]  G. Caldwell.  Compound optimization in early- and late-phase drug discovery: acceptable pharmacokinetic properties utilizing combined physicochemical, in vitro and in vivo screens. , 2000, Current opinion in drug discovery & development.

[11]  Tingjun Hou,et al.  ADME evaluation in drug discovery , 2002, Journal of molecular modeling.

[12]  J H Lin,et al.  Blood-brain barrier permeability and in vivo activity of partial agonists of benzodiazepine receptor: a study of L-663,581 and its metabolites in rats. , 1994, The Journal of pharmacology and experimental therapeutics.

[13]  M. Sanner,et al.  Reduced surface: an efficient way to compute molecular surfaces. , 1996, Biopolymers.

[14]  Philip L. Smith,et al.  In vitro models to predict blood-brain barrier permeability , 1997 .

[15]  A. Ghose,et al.  Atomic Physicochemical Parameters for Three‐Dimensional Structure‐Directed Quantitative Structure‐Activity Relationships I. Partition Coefficients as a Measure of Hydrophobicity , 1986 .

[16]  J. Michael.  Textbook of Medical Physiology , 2005 .

[17]  Glen Eugene Kellogg,et al.  HINT: A new method of empirical hydrophobic field calculation for CoMFA , 1991, J. Comput. Aided Mol. Des..

[18]  Harpreet S. Chadha,et al.  Hydrogen-bonding. Part 36. Determination of blood brain distribution using octanol-water partition coefficients. , 1995, Drug design and discovery.

[19]  Thomas A. Halgren.  Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[20]  David J. Begley,et al.  Potential of Immobilized Artificial Membranes for Predicting Drug Penetration Across the Blood−Brain Barrier , 1998, Pharmaceutical Research.

[21]  J. Gasteiger,et al.  ITERATIVE PARTIAL EQUALIZATION OF ORBITAL ELECTRONEGATIVITY – A RAPID ACCESS TO ATOMIC CHARGES , 1980 .

[22]  U Norinder,et al.  Theoretical calculation and prediction of brain-blood partitioning of organic solutes using MolSurf parametrization and PLS statistics. , 1998, Journal of pharmaceutical sciences.

[23]  Harpreet S. Chadha,et al.  Hydrogen bonding. 33. Factors that influence the distribution of solutes between blood and brain. , 1994, Journal of pharmaceutical sciences.

[24]  A. Ghose,et al.  Prediction of Hydrophobic (Lipophilic) Properties of Small Organic Molecules Using Fragmental Methods: An Analysis of ALOGP and CLOGP Methods , 1998 .

[25]  C R Ganellin,et al.  Predicting the brain-penetrating capability of histaminergic compounds. , 1994, Drug design and discovery.

[26]  T. Halgren.  Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[27]  A. Leo.  CALCULATING LOG POCT FROM STRUCTURES , 1993 .

[28]  Robert P. Sheridan,et al.  PATTY: A Programmable Atom Typer and Language for Automatic Classification of Atoms in Molecular Databases. , 1994 .

[29]  Juan M. Luco,et al.  Prediction of the Brain-Blood Distribution of a Large Set of Drugs from Structurally Derived Descriptors Using Partial Least-Squares (PLS) Modeling , 1999, J. Chem. Inf. Comput. Sci..

[30]  Bruce L. Bush,et al.  PATTY: A programmable atom type and language for automatic classification of atoms in molecular databases , 1993, J. Chem. Inf. Comput. Sci..

[31]  P. Selzer,et al.  Fast calculation of molecular polar surface area as a sum of fragment-based contributions and its application to the prediction of drug transport properties. , 2000, Journal of medicinal chemistry.

[32]  R Griffiths,et al.  Development of a new physicochemical model for brain penetration and its application to the design of centrally acting H2 receptor histamine antagonists. , 1988, Journal of medicinal chemistry.

[33]  Hugh Davson,et al.  Physiology of the CSF and Blood-Brain Barriers , 1996 .

[34]  Gordon M. Crippen,et al.  Prediction of Physicochemical Parameters by Atomic Contributions , 1999, J. Chem. Inf. Comput. Sci..

[35]  P. Carrupt,et al.  Molecular fields in quantitative structure–permeation relationships: the VolSurf approach , 2000 .

[36]  S Sarre,et al.  Brain, liver and blood distribution kinetics of carbamazepine and its metabolic interaction with clomipramine in rats: a quantitative microdialysis study. , 1995, The Journal of pharmacology and experimental therapeutics.

[37]  Michał J. Markuszewski,et al.  Brain/blood distribution described by a combination of partition coefficient and molecular mass , 1996 .

[38]  C. Decker,et al.  Prediction of pharmacokinetic properties using experimental approaches during early drug discovery. , 2001, Current opinion in chemical biology.

[39]  Denis M. Bayada,et al.  Polar Molecular Surface as a Dominating Determinant for Oral Absorption and Brain Penetration of Drugs , 1999, Pharmaceutical Research.