Using Descriptor Cutoffs and Multiparameter Optimization for Preparation of CNS Targeted Libraries

Expert-driven In Silico Drug Discovery Solutions
1 September 2020
Svitlana Kondovych
Senior Researcher

Disorders of the central nervous system (CNS) have been constantly puzzling the global health research community due to their plurality, diversity, and, above all, intractability. The human brain is an extremely well-protected organ, which is undoubtedly beneficial to us; however, this complex security system turns up to be a double-edged sword. It is hard to damage the brain. It is not any easier to cure it.

The most significant factor that constrains the discovery of new CNS drugs is the limited blood-brain barrier (BBB) penetrability. More than 98 % of small-molecule drugs do not cross BBB, whereas larger molecules essentially never succeed [1]. Thus, increasing BBB-targeting efficiency is one of the key issues in the design of novel CNS-active medications. Recent advances in CNS drug research rest upon medicinal chemistry [2-4] and nanotechnology-based approaches [3] to go through the BBB.

To optimize the study of potentially CNS-active compounds, researchers attempt to define the physicochemical characteristics of the successful CNS drugs and derive the favorable selection rules. Numerous reviews and statistical studies (see, e.g., [4] and references therein) aim to systematize available BBB permeation data, extract correlations, and suggest strategies on choosing the working drug.

The main parameters that determine a drug permeability include the molecular weight (MW) of a compound, its topological polar surface area (TPSA), charge, and lipophilicity. The latter is usually characterized by the calculated octanol-water partition coefficient (ClogP), as well as by its value at physiological pH = 7.4 (ClogD). Overall, CNS-targeted drugs appear to have high lipophilicity and low flexibility, they have fewer hydrogen-bond donors (HbD < 3) and acceptors (HbAc < 7), reduced molecular weight (150 – 400 Da), fewer formal charges (particularly, negative charges) and lower polar surface area (TPSA ≤ 65 Å2). The exact cut-off values for each of these parameters make up the subject of continuous search, elaboration, and discussion in the scientific community [1,2,4-18].

In most cases, to explain and predict CNS-related pharmacokinetic properties, one needs to solve a complex multi-parameter puzzle, taking into account not only values for individual physicochemical parameters but also considering their interplay. A review of the multivariate screening methods and their application is given in [4,19]. Among the most widely employed techniques are the CNS multiparameter optimization (MPO) [17] and the “Golden Triangle” rule [18,20].

Fig. 1. Multiparameter optimization approach comprises six physicochemical parameters

Fig. 2. The “Golden Triangle” rule utilizes the correlation between MW and lipophilicity

The former method implies the simultaneous estimation of the set of six parameters, namely: MW, TPSA, ClogP, ClogD, HbD, and the acid-base dissociation constant, pKa (Fig. 1). The algorithm returns the MPO desirability score, a number from 0 to 6, which may serve as a reference point for CNS drug design ideas. As a rule, the threshold value of MPO score ≥ 4 identifies the potential CNS-active compounds [17].

The latter approach – the “Golden Triangle” rule – involves joint estimation of molecular weight and lipophilicity ClogD of potential drugs. In fact, the limit of ClogD depends on MW [18]: greater lipophilicity corresponds to lower MWs, resulting in a triangle shape of the graph (Fig. 2) for acceptable drug-like compounds [20]. Considering this correlation may provide better predictions than treating the parameters separately.

Life Chemicals follows the most up-to-date criteria to offer its collection of compounds with desired CNS properties. The Life Chemicals CNS Screening Library comprises over 9,900 small organic molecules, carefully singled out from the database using the optimal threshold values for the main physicochemical parameters. Fig. 3 illustrates the resulting distributions and the calculated MPO desirability score.

To prevent the failure in BBB transport we additionally optimize the number of S, Cl and N atoms in the CNS-targeted compounds, ensure the absence of quaternary ammonium groups and allow not more than one carboxylic acid group [1]. All the compounds are subject to PAINS and toxicophore filters.

Please visit our Website or contact us at for more information.

 Distribution of the main physicochemical parameters and calculated MPO score  in the Life Chemicals CNS Library

Fig. 3. Distribution of the main physicochemical parameters and calculated MPO score 
in the Life Chemicals CNS Library



[1] Pardridge W. M. Blood–brain barrier delivery. Drug Discovery Today2007, 12 (1/2); 54–61.

[2] Pajouhesh H., Lenz G. R. Medicinal chemical properties of successful central nervous system drugs. NeuroRx., 2005, 2; 541–553.

[3] Vlieghe P., Khrestchatisky M. Medicinal chemistry based approaches and nanotechnology-based systems to improve CNS drug targeting and delivery. Medicinal Research Reviews2013, 33(3):457-516.

[4] Rankovic, Z. CNS Drug Design: Balancing Physicochemical Properties for Optimal Brain Exposure. Journal of Medicinal Chemistry2015, 58(6), 2584–2608.

[5] Banks, W. A., & Greig, N. H. Small molecules as central nervous system therapeutics: old challenges, new directions, and a philosophic divide. Future Medicinal Chemistry2019, 11(6):489-493.

[6] Didziapetris, R.; Japertas, P.; Avdeef, A.; Petrauskas, A. Classification analysis of P-glycoprotein substrate specificity. J. Drug Targeting2003, 11, 391−406.

[7] Hansch C, Leo A.J. Substituent constant for correlation analysis in chemistry and biology. New York: Wiley, 1979.

[8] Desai, P. V.; Sawada, G. A.; Watson, I. A.; Raub, T. J. Integration of in silico and in vitro tools for scaffold optimization during drug discovery: predicting P-glycoprotein efflux. Mol. Pharmaceutics2013, 10, 1249−1261.

[9] Ghose, A. K.; Herbertz, T.; Hudkins, R. L.; Dorsey, B. D.; Mallamo, J. P. Knowledge-based, central nervous system (CNS) lead selection and lead optimization for CNS drug discovery. ACS Chem. Neurosci.2012, 3, 50−68.

[10] Van de Waterbeemd, H.; Camenisch, G.; Folkers, G.; Chretien, J. R.; Raevsky, O. A. Estimation of blood−brain barrier crossing of drugs using molecular size and shape, and H-bonding descriptors. J. Drug Targeting1998, 6, 151−165.

[11] Hitchcock, S. A.; Pennington, L. D. Structure−brain exposure relationships. J. Med. Chem.2006, 49, 7559−7583

[12] Hitchcock, S. A. Structural modifications that alter the P-glycoprotein efflux properties of compounds. J. Med. Chem.,2012, 55, 4877−4895.

[13] Mao, F., Ni, W., Xu, X., Wang, H., Wang, J., Ji, M., & Li, J. Chemical Structure-Related Drug-Like Criteria of Global Approved Drugs. Molecules2016, 21(1), 75.

[14] Ilardi, E. A., Vitaku, E., & Njardarson, J. T. Data-Mining for Sulfur and Fluorine: An Evaluation of Pharmaceuticals To Reveal Opportunities for Drug Design and Discovery. Journal of Medicinal Chemistry2014, 57(7), 2832–2842.

[15] Smith, B. R., Eastman, C. M., & Njardarson, J. T. Beyond C, H, O, and N! Analysis of the Elemental Composition of U.S. FDA Approved Drug Architectures. Journal of Medicinal Chemistry, 2014, 57(23), 9764–9773.

[16] Wager, T. T., Chandrasekaran, R. Y., Hou, X., Troutman, M. D., Verhoest, P. R., Villalobos, A., & Will, Y. Defining Desirable Central Nervous System Drug Space through the Alignment of Molecular Properties, in Vitro ADME, and Safety Attributes. ACS Chemical Neuroscience2010, 1(6), 420–434.

[17] Wager, T. T., Hou, X., Verhoest, P. R., & Villalobos, A. Moving beyond Rules: The Development of a Central Nervous System Multiparameter Optimization (CNS MPO) Approach To Enable Alignment of Druglike Properties. ACS Chemical Neuroscience2010, 1(6), 435–449.

[18] Waring, M.J. Defining optimum lipophilicity and molecular weight ranges for drug candidates molecular weight dependent lower logD limits based on permeability. Bioorg. Med. Chem. Lett.2009, 19, 2844–2851.

[19] Zerroug, A., Belaidi, S., BenBrahim, I., Sinha, L., & Chtita, S. Virtual screening in drug-likeness and structure/activity relationship of pyridazine derivatives as Anti-Alzheimer drugs. Journal of King Saud University – Science2019, 31, 595–601.

[20] Johnson T. W., Dress K. R., & Edwards, M. Using the Golden Triangle to optimize clearance and oral absorption. Bioorg. Med. Chem. Lett., 2009, 19, 5560-5564.

1 September 2020, 12:11 Svitlana Kondovych Computational Chemistry

Comments ()