Prediction of blood-brain partitioning: a model based on ab initio calculated quantum chemical descriptors

Computer Assisted Models for Blood Brain Barrier Permeation of 1, 5-Benzodiazepines

AIM

To generate and validate predictive models for blood-brain permeation (BBB) of CNS molecules using the QSPR approach.

BACKGROUND

Prediction of molecules crossing BBB remains a challenge in drug delivery. Predictive models are designed for the evaluation of a set of preclinical drugs which may serve as alternatives for determining BBB permeation by experimentation.

OBJECTIVE

The objective of the present study was to generate QSPR models for the permeation of CNS molecules across BBB and its validation using existing in-house leads.

METHODS

The present study envisaged the determination of the set of molecular descriptors which are considered significant correlative factors for BBB permeation property. Quantitative Structure- Property Relationship (QSPR) approach was followed to describe the correlation between identified descriptors for 45 molecules and highest, moderate and least BBB permeation data. The molecular descriptors were selected based on drug-likeness, hydrophilicity, hydrophobicity, polar surface area, etc. of molecules that served the highest correlation with BBB permeation. The experimental data in terms of log BB were collected from available literature, subjected to 2D-QSPR model generation using a regression analysis method like Multiple Linear Regression (MLR).

RESULTS

The best QSPR model was Model 3, which exhibited regression coefficient as R²= 0.89, F = 36; Q²= 0.7805 and properties such as polar surface area, hydrophobic hydrophilic distance, electronegativity, etc., which were considered key parameters in the determination of the BBB permeability. The developed QSPR models were validated with in-house 1,5-benzodiazepines molecules and correlation studies were conducted between experimental and predicted BBB permeability.

CONCLUSION

The QSPR model 3 showed predictive results that were in good agreements with experimental results for blood-brain permeation. Thus, this model was found to be satisfactory in achieving a good correlation between selected descriptors and BBB permeation for benzodiazepines and tricyclic compounds.

Analyzing the substitution effect on the CoMFA results within the framework of density functional theory (DFT)

Though QSAR was originally developed in the context of physical organic chemistry, it has been applied very extensively to chemicals (drugs) which act on biological systems, in this idea one of the most important QSAR methods is the 3D QSAR model. However, due to the complexity of understanding the results it is necessary to postulate new methodologies to highlight their physical-chemical meaning. In this sense, this work postulates new insights to understand the CoMFA results using molecular quantum similarity and chemical reactivity descriptors within the framework of density functional theory. To obtain these insights a simple theoretical scheme involving quantum similarity (overlap, coulomb operators, their euclidean distances) and chemical reactivity descriptors such as chemical potential (μ), hardness (ɳ), softness (S), electrophilicity (ω), and the Fukui functions, was used to understand the substitution effect. In this sense, this methodology can be applied to analyze the biological activity and the stabilization process in the non-covalent interactions on a particular molecular set taking a reference compound.

A Generally Applicable Computer Algorithm Based on the Group Additivity Method for the Calculation of Seven Molecular Descriptors: Heat of Combustion, LogPO/W, LogS, Refractivity, Polarizability, Toxicity and LogBB of Organic Compounds; Scope and Limits of Applicability

A generally applicable computer algorithm for the calculation of the seven molecular descriptors heat of combustion, logP_{octanol/water}, logS (water solubility), molar refractivity, molecular polarizability, aqueous toxicity (protozoan growth inhibition) and logBB (log (c_blood/c_brain)) is presented. The method, an extendable form of the group-additivity method, is based on the complete break-down of the molecules into their constituting atoms and their immediate neighbourhood. The contribution of the resulting atom groups to the descriptor values is calculated using the Gauss-Seidel fitting method, based on experimental data gathered from literature. The plausibility of the method was tested for each descriptor by means of a k-fold cross-validation procedure demonstrating good to excellent predictive power for the former six descriptors and low reliability of logBB predictions. The goodness of fit (Q²) and the standard deviation of the 10-fold cross-validation calculation was >0.9999 and 25.2 kJ/mol, respectively, (based on N = 1965 test compounds) for the heat of combustion, 0.9451 and 0.51 (N = 2640) for logP, 0.8838 and 0.74 (N = 1419) for logS, 0.9987 and 0.74 (N = 4045) for the molar refractivity, 0.9897 and 0.77 (N = 308) for the molecular polarizability, 0.8404 and 0.42 (N = 810) for the toxicity and 0.4709 and 0.53 (N = 383) for logBB. The latter descriptor revealing a very low Q² for the test molecules (R² was 0.7068 and standard deviation 0.38 for N = 413 training molecules) is included as an example to show the limits of the group-additivity method. An eighth molecular descriptor, the heat of formation, was indirectly calculated from the heat of combustion data and correlated with published experimental heat of formation data with a correlation coefficient R² of 0.9974 (N = 2031).

A Genetic Algorithm Based Support Vector Machine Model for Blood-Brain Barrier Penetration Prediction

Blood-brain barrier (BBB) is a highly complex physical barrier determining what substances are allowed to enter the brain. Support vector machine (SVM) is a kernel-based machine learning method that is widely used in QSAR study. For a successful SVM model, the kernel parameters for SVM and feature subset selection are the most important factors affecting prediction accuracy. In most studies, they are treated as two independent problems, but it has been proven that they could affect each other. We designed and implemented genetic algorithm (GA) to optimize kernel parameters and feature subset selection for SVM regression and applied it to the BBB penetration prediction. The results show that our GA/SVM model is more accurate than other currently available log BB models. Therefore, to optimize both SVM parameters and feature subset simultaneously with genetic algorithm is a better approach than other methods that treat the two problems separately. Analysis of our log BB model suggests that carboxylic acid group, polar surface area (PSA)/hydrogen-bonding ability, lipophilicity, and molecular charge play important role in BBB penetration. Among those properties relevant to BBB penetration, lipophilicity could enhance the BBB penetration while all the others are negatively correlated with BBB penetration.

Computer modeling of blood brain barrier permeability of physiologically active compounds

At present work discusses the current level of computer modeling the relationship structure of organic compounds and drugs and their ability to penetrate the BBB. All descriptors that influence to this permeability within classification and regression QSAR models are generalized and analyzed. The crucial role of H-bond in processes both passive, and active transport across BBB is observed. It is concluded that further research should be focused on interpretation the spatial structure of a full-size P-glycoprotein molecule with high resolution and the creation of QSAR models describing the quantitative relationship between structure and active transport of substances across BBB.

A 2D-QSPR approach to predict blood-brain barrier penetration of drugs acting on the central nervous system

Drugs acting on the central nervous system (CNS) have to cross the blood-brain barrier (BBB) in order to perform their pharmacological actions. Passive BBB diffusion can be partially expressed by the blood/brain partition coefficient (logBB). As the experimental evaluation of logBB is time and cost consuming, theoretical methods such as quantitative structure-property relationships (QSPR) can be useful to predict logBB values. In this study, a 2D-QSPR approach was applied to a set of 28 drugs acting on the CNS, using the logBB property as biological data. The best QSPR model [n = 21, r = 0.94 (r² = 0.88), s = 0.28, and Q² = 0.82] presented three molecular descriptors: calculated n-octanol/water partition coefficient (ClogP), polar surface area (PSA), and polarizability (α). Six out of the seven compounds from the test set were well predicted, which corresponds to good external predictability (85.7%). These findings can be helpful to guide future approaches regarding those molecular descriptors which must be considered for estimating the logBB property, and also for predicting the BBB crossing ability for molecules structurally related to the investigated set.

Prediction of blood-brain partitioning: a model based on molecular electronegativity distance vector descriptors

The objective of this paper is to build a reliable model based on the molecular electronegativity distance vector (MEDV) descriptors for predicting the blood-brain barrier (BBB) permeability and to reveal the effects of the molecular structural segments on the BBB permeability. Using 70 structurally diverse compounds, the partial least squares regression (PLSR) models between the BBB permeability and the MEDV descriptors were developed and validated by the variable selection and modeling based on prediction (VSMP) technique. The estimation ability, stability, and predictive power of a model are evaluated by the estimated correlation coefficient (r), leave-one-out (LOO) cross-validation correlation coefficient (q), and predictive correlation coefficient (R(p)). It has been found that PLSR model has good quality, r=0.9202, q=0.7956, and R(p)=0.6649 for M1 model based on the training set of 57 samples. To search the most important structural factors affecting the BBB permeability of compounds, we performed the values of the variable importance in projection (VIP) analysis for MEDV descriptors. It was found that some structural fragments in compounds, such as -CH(3), -CH(2)-, =CH-, =C, triple bond C-, -CH<, =C<, =N-, -NH-, =O, and -OH, are the most important factors affecting the BBB permeability.

New predictive models for blood-brain barrier permeability of drug-like molecules

PURPOSE

The goals of the present study were to apply a generalized regression model and support vector machine (SVM) models with Shape Signatures descriptors, to the domain of blood-brain barrier (BBB) modeling.

MATERIALS AND METHODS

The Shape Signatures method is a novel computational tool that was used to generate molecular descriptors utilized with the SVM classification technique with various BBB datasets. For comparison purposes we have created a generalized linear regression model with eight MOE descriptors and these same descriptors were also used to create SVM models.

RESULTS

The generalized regression model was tested on 100 molecules not in the model and resulted in a correlation r2 = 0.65. SVM models with MOE descriptors were superior to regression models, while Shape Signatures SVM models were comparable or better than those with MOE descriptors. The best 2D shape signature models had 10-fold cross validation prediction accuracy between 80-83% and leave-20%-out testing prediction accuracy between 80-82% as well as correctly predicting 84% of BBB+ compounds (n = 95) in an external database of drugs.

CONCLUSIONS

Our data indicate that Shape Signatures descriptors can be used with SVM and these models may have utility for predicting blood-brain barrier permeation in drug discovery.