The oxidoreductase PYROXD1 uses NAD(P)+ as an antioxidant to sustain tRNA ligase activity in pre-tRNA splicing and unfolded protein response.

The tRNA ligase complex (tRNA-LC) splices precursor tRNAs (pre-tRNA), and Xbp1-mRNA during the unfolded protein response (UPR). In aerobic conditions, a cysteine residue bound to two metal ions in its ancient, catalytic subunit RTCB could make the tRNA-LC susceptible to oxidative inactivation. Here, we confirm this hypothesis and reveal a co-evolutionary association between the tRNA-LC and PYROXD1, a conserved and essential oxidoreductase. We reveal that PYROXD1 preserves the activity of the mammalian tRNA-LC in pre-tRNA splicing and UPR. PYROXD1 binds the tRNA-LC in the presence of NAD(P)H and converts RTCB-bound NAD(P)H into NAD(P)+, a typical oxidative co-enzyme. However, NAD(P)+ here acts as an antioxidant and protects the tRNA-LC from oxidative inactivation, which is dependent on copper ions. Genetic variants of PYROXD1 that cause human myopathies only partially support tRNA-LC activity. Thus, we establish the tRNA-LC as an oxidation-sensitive metalloenzyme, safeguarded by the flavoprotein PYROXD1 through an unexpected redox mechanism.

Reversible Control of RNA Splicing by Photoswitchable Small Molecules.

Dynamics are intrinsic to both RNA function and structure. Yet, the available means to precisely provide RNA-based processes with spatiotemporal resolution are scarce. Here, our work pioneers a reversible approach to regulate RNA splicing within primary patient-derived cells by synthetic photoswitches. Our small molecule enables conditional real-time control at mRNA and protein levels. NMR experiments, together with theoretical calculations, photochemical characterization, fluorescence polarization measurements, and living cell-based assays, confirmed light-dependent exon inclusion as well as an increase in the target functional protein. Therefore, we first demonstrated the potential of photopharmacology modulation in splicing, tweaking the current optochemical toolkit. The timeliness on the consolidation of RNA research as the driving force toward therapeutical innovation holds the promise that our approach will contribute to redrawing the vision of RNA.

Expanding the Accessible Chemical Space of SIRT2 Inhibitors through Exploration of Binding Pocket Dynamics.

Considerations of binding pocket dynamics are one of the crucial aspects of the rational design of binders. Identification of alternative conformational states or cryptic subpockets could lead to the discovery of completely novel groups of the ligands. However, experimental characterization of pocket dynamics, besides being expensive, may not be able to elucidate all of the conformational states relevant for drug discovery projects. In this study, we propose the protocol for computational simulations of sirtuin 2 (SIRT2) binding pocket dynamics and its integration into the structure-based virtual screening (SBVS) pipeline. Initially, unbiased molecular dynamics simulations of SIRT2:inhibitor complexes were performed using optimized force field parameters of SIRT2 inhibitors. Time-lagged independent component analysis (tICA) was used to design pocket-related collective variables (CVs) for enhanced sampling of SIRT2 pocket dynamics. Metadynamics simulations in the tICA eigenvector space revealed alternative conformational states of the SIRT2 binding pocket and the existence of a cryptic subpocket. Newly identified SIRT2 conformational states outperformed experimentally resolved states in retrospective SBVS validation. After performing prospective SBVS, compounds from the under-represented portions of the SIRT2 inhibitor chemical space were selected for in vitro evaluation. Two compounds, NDJ18 and NDJ85, were identified as potent and selective SIRT2 inhibitors, which validated the in silico protocol and opened up the possibility for generalization and broadening of its application. The anticancer effects of the most potent compound NDJ18 were examined on the triple-negative breast cancer cell line. Results indicated that NDJ18 represents a promising structure suitable for further evaluation.

Synthesis, in silico, and in vitro studies of novel dopamine D2 and D3 receptor ligands.

Dopamine is an important neurotransmitter in the human brain and its altered concentrations can lead to various neurological diseases. We studied the binding of novel compounds at the dopamine D2 (D2 R) and D3 (D3 R) receptor subtypes, which belong to the D2 -like receptor family. The synthesis, in silico, and in vitro characterization of 10 dopamine receptor ligands were performed. Novel ligands were docked into the D2 R and D3 R crystal structures to examine the precise binding mode. A quantum mechanics/molecular mechanics study was performed to gain insights into the nature of the intermolecular interactions between the newly introduced pentafluorosulfanyl (SF5 ) moiety and D2 R and D3 R. A radioligand displacement assay determined that all of the ligands showed moderate-to-low nanomolar affinities at D2 R and D3 R, with a slight preference for D3 R, which was confirmed in the in silico studies. N-{4-[4-(2-Methoxyphenyl)piperazin-1-yl]butyl}-4-(pentafluoro-λ6-sulfanyl)benzamide (7i) showed the highest D3 R affinity and selectivity (pKi values of 7.14 [D2 R] and 8.42 [D3 R]).

Synthesis and Biological Activity of a Cytostatic Inhibitor of MLLr Leukemia Targeting the DOT1L Protein.

Histone methyltransferase DOT1L catalyzes mono-, di- and trimethylation of histone 3 at lysine residue 79 (H3K79) and hypermethylation of H3K79 has been linked to the development of acute leukemias characterized by the MLL (mixed-lineage leukemia) rearrangements (MLLr cells). The inhibition of H3K79 methylation inhibits MLLr cells proliferation, and an inhibitor specific for DOT1L, pinometostat, was in clinical trials (Phase Ib/II). However, the compound showed poor pharmacological properties. Thus, there is a need to find new potent inhibitors of DOT1L for the treatment of rearranged leukemias. Here we present the design, synthesis, and biological evaluation of a small molecule that inhibits in the nM level the enzymatic activity of hDOT1L, H3K79 methylation in MLLr cells with comparable potency to pinometostat, associated with improved metabolic stability and a characteristic cytostatic effect.

Modulating Protein-Protein Interactions with Visible-Light-Responsive Peptide Backbone Photoswitches.

Life relies on a myriad of carefully orchestrated processes, in which proteins and their direct interplay ultimately determine cellular function and disease. Modulation of this complex crosstalk has recently attracted attention, even as a novel therapeutic strategy. Herein, we describe the synthesis and characterization of two visible-light-responsive peptide backbone photoswitches based on azobenzene derivatives, to exert optical control over protein-protein interactions (PPI). The novel peptidomimetics undergo fast and reversible isomerization with low photochemical fatigue under alternatively blue-/green-light irradiation cycles. Both bind in the nanomolar range to the protein of interest. Importantly, the best peptidomimetic displays a clear difference between isomers in its protein-binding capacity and, in turn, in its potential to inhibit enzymatic activity through PPI disruption. In addition, crystal structure determination, docking and molecular dynamics calculations allow a molecular interpretation and open up new avenues in the design and synthesis of future photoswitchable PPI modulators.

Ultra-performance liquid chromatography tandem mass spectrometry for the rapid, simultaneous analysis of ziprasidone and its impurities.

The separation and characterization of the unknown degradation product of second-generation antipsychotic drug ziprasidone are essential for defining the genotoxic potential of the compound. The aim of this study was to develop a simple UHPLC method coupled with tandem mass spectrometry (MS/MS) for chemical characterization of an unknown degradant, and the separation and quantification of ziprasidone and its five main impurities (I-V) in the raw material and pharmaceuticals. Chromatographic conditions were optimized by experimental design. The MS/MS fragmentation conditions were optimized individually for each compound in order to obtain both specific fragments and high signal intensity. A rapid and sensitive UHPLC-MS/MS method was developed. All seven analytes were eluted within the 7 min run time. The best separation was obtained on the Acquity UPLC BEH C18 (50 × 2.1 mm × 1.7 µm) column in gradient mode with ammonium-formate buffer (10 mm; pH 4.7) and acetonitrile as mobile phase, with the flow rate of 0.3 mL min-1 and at the column temperature of 30°C. The new UHPLC-MS/MS method was fully validated and all validation parameters were confirmed. The fragmentation pathways and chemical characterization of an unknown degradant were proposed and it was confirmed that there are no structural alerts concerning genotoxicity.

Severe hepatotoxicity induced by efavirenz in a treatment-naïve, low body mass index HIV-infected, female patient with no hepatitis and other virus co-infections.

Combined antiretroviral therapy (cART) consisting of two nucleoside reverse transcriptase inhibitors (NRTIs) and one non-nucleoside reverse transcriptase inhibitor (NNRTI), such as efavirenz, is still the first-line treatment in resource-limited settings. However, efavirenz has shown strong prominence of disadvantages with variance in plasma concentration and central nervous side effects. Our study presents HIV infected, drug naïve, female patient with relatively low BMI, CYP2B6 516G>T (rs3745274) genotype with high efavirenz plasma concentration. In this case report, the patient was admitted at the hospital 6 months after cART initiation with drug-induced severe hepatotoxicity. Furthermore, pathophysiological findings proved confluent parenchymal necrosis after aspiration liver biopsy, with mild to moderate inflammation in portal tracts with focal interface hepatitis. All other possible causes were excluded. Thus, we conclude that efavirenz has a potential harmful effect in patients with low BMI, specific genotyping and interindividual pharmacokinetics affecting high plasma concentration.

Synthesis, In Silico, and In Vitro Evaluation of Anti-Leishmanial Activity of Oxadiazoles and Indolizine Containing Compounds Flagged against Anti-Targets.

Due to the lack of approved vaccines against human leishmaniasis and the limitations of the current chemotherapy inducing side effects and drug resistance, development of new, effective chemotherapeutic agents is essential. This study describes the synthesis of a series of novel oxadiazoles and indolizine-containing compounds. The compounds were screened in silico using an EIIP/AQVN filter followed by ligand-based virtual screening and molecular docking to parasite arginase. Top hits were further screened versus human arginase and finally against an anti-target battery to tag their possible interactions with proteins essential for the metabolism and clearance of many substances. Eight candidate compounds were selected for further experimental testing. The results show measurable in vitro anti-leishmanial activity for three compounds. One compound with an IC50 value of 2.18 µM on Leishmania donovani intramacrophage amastigotes is clearly better positioned than the others as an interesting molecular template for further development of new anti-leishmanial agents.

Comparison of local versus general anesthesia in patients undergoing transcatheter aortic valve replacement: A meta-analysis.

BACKGROUND: Transcatheter aortic valve replacement (TAVR) is typically performed under general anesthesia (GA). However, there is increasing data supporting the safety of performing TAVR under local anesthesia/conscious sedation (LA). We performed a meta-analysis to gain better understanding of the safety and efficacy of LA versus GA in patients with severe aortic stenosis undergoing TAVR. METHODS AND RESULTS: We comprehensively searched EMBASE, PubMed, and Web of Science. Effect sizes were summarized using risk ratios (RRs) difference of the mean (DM), and 95% CIs (confidence intervals) for dichotomous and continuous variables respectively. Twenty-six studies and 10,572 patients were included in the meta-analysis. The use of LA for TAVR was associated with lower overall 30-day mortality (RR, 0.73; 95% CI, 0.57-0.93; P = 0.01), use of inotropic/vasopressor drugs (RR, 0.45; 95% CI, 0.28-0.72; P < 0.001), hospital length of stay (LOS) (DM, -2.09; 95% CI, -3.02 to -1.16; P < 0.001), intensive care unit LOS (DM, -0.18; 95% CI, -0.31 to -0.04; P = 0.01), procedure time (DM, -25.02; 95% CI, -32.70 to -17.35; P < 0.001); and fluoroscopy time (DM, -1.63; 95% CI, -3.02 to -0.24; P = 0.02). No differences were observed between LA and GA for stroke, cardiovascular mortality, myocardial infarction, permanent pacemaker implantation, acute kidney injury, paravalvular leak, vascular complications, major bleeding, procedural success, conduction abnormalities, and annular rupture. CONCLUSION: Our meta-analysis suggests that use of LA for TAVR is associated with a lower 30-day mortality, shorter procedure time, fluoroscopy time, ICU LOS, hospital length of stay, and reduced need for inotropic support.

Evaluation of Biological Activity and Computer-Aided Design of New Soft Glucocorticoids.

Soft glucocorticoids are compounds that are biotransformed to inactive and non-toxic metabolites and have fewer side effects than traditional glucocorticoids. A new class of 17ß-carboxamide steroids has been recently introduced by our group. In this study, local anti-inflammatory activity of these derivatives was evaluated by use of the croton oil-induced ear edema test. Glucocorticoids with the highest maximal edema inhibition (MEI) were pointed out, and the systemic side effects of those with the lowest EC50 values were significantly lower in comparison to dexamethasone. A 3D-QSAR model was created and employed for the design of 27 compounds. By use of the sequential combination of ligand-based and structure-based virtual screening, three compounds were selected from the ChEMBL library and used as a starting point for the design of 15 derivatives. Molecular docking analysis of the designed derivatives with the highest predicted MEI and relative glucocorticoid receptor binding affinity (20, 22, 24-1, 25-1, 27, VS7, VS13, and VS14) confirmed the presence of interactions with the glucocorticoid receptor that are important for the activity.

The Forty-Sixth Euro Congress on Drug Synthesis and Analysis: Snapshot †.

The 46th EuroCongress on Drug Synthesis and Analysis (ECDSA-2017) was arranged within the celebration of the 65th Anniversary of the Faculty of Pharmacy at Comenius University in Bratislava, Slovakia from 5-8 September 2017 to get together specialists in medicinal chemistry, organic synthesis, pharmaceutical analysis, screening of bioactive compounds, pharmacology and drug formulations; promote the exchange of scientific results, methods and ideas; and encourage cooperation between researchers from all over the world. The topic of the conference, "Drug Synthesis and Analysis," meant that the symposium welcomed all pharmacists and/or researchers (chemists, analysts, biologists) and students interested in scientific work dealing with investigations of biologically active compounds as potential drugs. The authors of this manuscript were plenary speakers and other participants of the symposium and members of their research teams. The following summary highlights the major points/topics of the meeting.

A combined ligand- and structure-based approach for the identification of rilmenidine-derived compounds which synergize the antitumor effects of doxorubicin.

The clonidine-like central antihypertensive agent rilmenidine, which has high affinity for I1-type imidazoline receptors (I1-IR) was recently found to have cytotoxic effects on cultured cancer cell lines. However, due to its pharmacological effects resulting also from α2-adrenoceptor activation, rilmenidine cannot be considered a suitable anticancer drug candidate. Here, we report the identification of novel rilmenidine-derived compounds with anticancer potential and devoid of α2-adrenoceptor effects by means of ligand- and structure-based drug design approaches. Starting from a large virtual library, eleven compounds were selected, synthesized and submitted to biological evaluation. The most active compound 5 exhibited a cytotoxic profile similar to that of rilmenidine, but without appreciable affinity to α2-adrenoceptors. In addition, compound 5 significantly enhanced the apoptotic response to doxorubicin, and may thus represent an important tool for the development of better adjuvant chemotherapeutic strategies for doxorubicin-insensitive cancers.

The influence of bile salts on the distribution of simvastatin in the octanol/buffer system.

INTRODUCTION: Distribution coefficient (D) is useful parameter for evaluating drugs permeability properties across biological membranes, which are of importance for drugs bioavailability. Given that bile acids are intensively studied as drug permeation-modifying and -solubilizing agents, the aim of this study was to estimate the influence of sodium salts of cholic (CA), deoxycholic (DCA) and 12-monoketocholic acids (MKC) on distribution coefficient of simvastatin (SV) (lactone [SVL] and acid form [SVA]) which is a highly lipophilic compound with extremely low water solubility and bioavailability. METHODS: LogD values of SVA and SVL with or without bile salts were measured by liquid-liquid extraction in n-octanol/buffer systems at pH 5 and 7.4. SV concentrations in aqueous phase were determined by HPLC-DAD. Chem3D Ultra program was applied for computation of physico-chemical properties of analyzed compounds and their complexes. RESULTS: Statistically significant decrease in both SVA and SVL logD was observed for all three studied bile salts at both selected pH. MKC exerted the most pronounced effect in the case of SVA while there were no statistically significant differences between observed bile salts for SVL. The calculated physico-chemical properties of analyzed compounds and their complexes supported experimental results. CONCLUSIONS: Our data indicate that the addition of bile salts into the n-octanol/buffer system decreases the values of SV distribution coefficient at both studied pH values. This may be the result of the formation of hydrophilic complexes increasing the solubility of SV that could consequently impact the pharmacokinetic parameters of SV and the final drug response in patients.

Predicting targets of compounds against neurological diseases using cheminformatic methodology.

Recently developed multi-targeted ligands are novel drug candidates able to interact with monoamine oxidase A and B; acetylcholinesterase and butyrylcholinesterase; or with histamine N-methyltransferase and histamine H3-receptor (H3R). These proteins are drug targets in the treatment of depression, Alzheimer's disease, obsessive disorders, and Parkinson's disease. A probabilistic method, the Parzen-Rosenblatt window approach, was used to build a "predictor" model using data collected from the ChEMBL database. The model can be used to predict both the primary pharmaceutical target and off-targets of a compound based on its structure. Molecular structures were represented based on the circular fingerprint methodology. The same approach was used to build a "predictor" model from the DrugBank dataset to determine the main pharmacological groups of the compound. The study of off-target interactions is now recognised as crucial to the understanding of both drug action and toxicology. Primary pharmaceutical targets and off-targets for the novel multi-target ligands were examined by use of the developed cheminformatic method. Several multi-target ligands were selected for further study, as compounds with possible additional beneficial pharmacological activities. The cheminformatic targets identifications were in agreement with four 3D-QSAR (H3R/D1R/D2R/5-HT2aR) models and by in vitro assays for serotonin 5-HT1a and 5-HT2a receptor binding of the most promising ligand (71/MBA-VEG8).

The stress stability of olanzapine: studies of interactions with excipients in solid state pharmaceutical formulations.

Stress stability testing represents an important part of the drug development process. It is used as an important tool for the identification of degradation products and degradation pathways, as well as for the assessment of changes in physical form of drug molecules. The impact of excipients on the stability of olanzapine confirms that levels of impurities and degradants are limiting parameters and are therefore used for stability evaluation. The major degradation product of olanzapine was identified as 2-methyl-5,10-dihydro-4H-thieno[2,3-b][1,5]benzodiazepine-4-one (III). The structure of III was determined by using LC-MS, IR and NMR. Compatibility and stress stability results demonstrated that tablet formulations of olanzapine are sensitive to temperature and moisture. In samples protected from moisture, the increase in concentration of III was shown to be highly temperature dependent and the degradation followed zero-order kinetics. In addition, studies of olanzapine with excipients and in formulated tablets revealed polymorphic phase changes in some samples, influenced by a combination of stress temperature and humidity conditions. Polymorphic transitions were monitored using x-ray powder diffraction (XRPD) analysis and exhibited no correlation between the phase change (appearance of a new polymorph) and the degradation process.

Application of parallel artificial membrane permeability assay technique and chemometric modeling for blood-brain barrier permeability prediction of protein kinase inhibitors.

Aim: This study aims to investigate the passive diffusion of protein kinase inhibitors through the blood-brain barrier (BBB) and to develop a model for their permeability prediction. Materials & methods: We used the parallel artificial membrane permeability assay to obtain logPe values of each of 34 compounds and calculated descriptors for these structures to perform quantitative structure-property relationship modeling, creating different regression models. Results: The logPe values have been calculated for all 34 compounds. Support vector machine regression was considered the most reliable, and CATS2D_09_DA, CATS2D_04_AA, B04[N-S] and F07[C-N] descriptors were identified as the most influential to passive BBB permeability. Conclusion: The quantitative structure-property relationship-support vector machine regression model that has been generated can serve as an efficient method for preliminary screening of BBB permeability of new analogs.

[Box: see text].

Machine Learning-Based Prediction of Glioma IDH Gene Mutation Status Using Physio-Metabolic MRI of Oxygen Metabolism and Neovascularization (A Bicenter Study).

The mutational status of the isocitrate dehydrogenase (IDH) gene plays a key role in the treatment of glioma patients because it is known to affect energy metabolism pathways relevant to glioma. Physio-metabolic magnetic resonance imaging (MRI) enables the non-invasive analysis of oxygen metabolism and tissue hypoxia as well as associated neovascularization and microvascular architecture. However, evaluating such complex neuroimaging data requires computational support. Traditional machine learning algorithms and simple deep learning models were trained with radiomic features from clinical MRI (cMRI) or physio-metabolic MRI data. A total of 215 patients (first center: 166 participants + 16 participants for independent internal testing of the algorithms versus second site: 33 participants for independent external testing) were enrolled using two different physio-metabolic MRI protocols. The algorithms trained with physio-metabolic data demonstrated the best classification performance in independent internal testing: precision, 91.7%; accuracy, 87.5%; area under the receiver operating curve (AUROC), 0.979. In external testing, traditional machine learning models trained with cMRI data exhibited the best IDH classification results: precision, 84.9%; accuracy, 81.8%; and AUROC, 0.879. The poor performance for the physio-metabolic MRI approach appears to be explainable by site-dependent differences in data acquisition methodologies. The physio-metabolic MRI approach potentially supports reliable classification of IDH gene status in the presurgical stage of glioma patients. However, non-standardized protocols limit the level of evidence and underlie the need for a reproducible framework of data acquisition techniques.

Quantitative structure-mobility relationship analysis of imidazoline receptor ligands in CDs-mediated CE.

The performed quantitative structure-mobility relationship (QSMR) study has investigated relative migration times of 11 guanidine/imidazoline derivatives, imidazoline receptor ligands, in CE system containing one of CDs, α-, ß-, or γ-CD, using linear and nonlinear modeling methods. The analyzed ligands and their inclusion complexes with CDs were fully examined and optimized at semiempirical parametrized model 3 level. The density functional theory, such as B3LYP/6-31G+(d,p)/3-21G(d)/STO-3G(d,p)/STO-3G(d), and ab initio theory, such as HF/3-21G(d)/STO-3G(d), were applied for molecular descriptors computation of the optimized ligands and their complexes. Predictive performances of the developed QSMR models were tested by use of the cross-validation and external test set prediction. Obtained results for Q(2) values (0.869, 0.911, and 0.966 for CE system containing α-, ß-, and γ-CD, respectively) and root mean squared error of prediction (0.239, 0.242, and 0.288 for α-, ß-, and γ-CD, respectively) were proved high predictive power of the proposed models. Finally, multitarget QSMR model, using the ligands descriptors (X) and the relative migration time in presence of α-CD (Y1), ß-CD (Y2), and γ-CD (Y3), has been created. The multitarget QSMR model can be used as initial screening predictive tool for CE migration behavior of other related guanidine/imidazoline derivatives in presence of α-, ß-, and γ-CD.

Partial least square and hierarchical clustering in ADMET modeling: prediction of blood-brain barrier permeation of α-adrenergic and imidazoline receptor ligands.

PURPOSE. Rate of brain penetration (logPS), brain/plasma equilibration rate (logPS-brain), and extent of blood-brain barrier permeation (logBB) of 29 α-adrenergic and imidazoline-receptors ligands were examined in Quantitative-Structure-Property Relationship (QSPR) study. METHODS. Experimentally determined chromatographic retention data (logKw at pH 4.4, slope (S) at pH 4.4, logKw at pH 7.4, slope (S) at pH 7.4, logKw at pH 9.1, and slope (S) at pH 9.1) and capillary electrophoresis migration parameters (µeff at pH 4.4, µeff at pH 7.4, and µeff at pH 9.1), together with calculated molecular descriptors, were used as independent variables in the QSPR study by use of partial least square (PLS) methodology. RESULTS. Predictive potential of the formed QSPR models, QSPR(logPS), QSPR(logPS-brain), QSPR(logBB), was confirmed by cross- and external validation. Hydrophilicity (Hy) and H-indices (H7m) were selected as significant parameters negatively correlated with both logPS and logPS-brain, while topological polar surface area (TPSA(NO)) was chosen as molecular descriptor negatively correlated with both logPS and logBB. The principal component analysis (PCA) and hierarchical clustering analysis (HCA) were applied to cluster examined drugs based on their chromatographic, electrophoretic and molecular properties. Significant positive correlations were obtained between the slope (S) at pH 7.4 and logBB in A/B cluster and between the logKw at pH 9.1 and logPS in C/D cluster. CONCLUSIONS. Results of the QSPR, clustering and correlation studies could be used as novel tool for evaluation of blood-brain barrier permeation of related α-adrenergic/imidazoline receptor ligands.This article is open to POST-PUBLICATION REVIEW. Registered readers (see "For Readers") may comment by clicking on ABSTRACT on the issue's contents page.PURPOSE. Rate of brain penetration (logPS), brain/plasma equilibration rate (logPS-brain), and extent of blood-brain barrier permeation (logBB) of 29 α-adrenergic and imidazoline-receptors ligands were examined in Quantitative-Structure-Property Relationship (QSPR) study. METHODS. Experimentally determined chromatographic retention data (logKw at pH 4.4, slope (S) at pH 4.4, logKw at pH 7.4, slope (S) at pH 7.4, logKw at pH 9.1, and slope (S) at pH 9.1) and capillary electrophoresis migration parameters (µeff at pH 4.4, µeff at pH 7.4, and µeff at pH 9.1), together with calculated molecular descriptors, were used as independent variables in the QSPR study by use of partial least square (PLS) methodology. RESULTS. Predictive potential of the formed QSPR models, QSPR(logPS), QSPR(logPS-brain), QSPR(logBB), was confirmed by cross- and external validation. Hydrophilicity (Hy) and H-indices (H7m) were selected as significant parameters negatively correlated with both logPS and logPS-brain, while topological polar surface area (TPSA(NO)) was chosen as molecular descriptor negatively correlated with both logPS and logBB. The principal component analysis (PCA) and hierarchical clustering analysis (HCA) were applied to cluster examined drugs based on their chromatographic, electrophoretic and molecular properties. Significant positive correlations were obtained between the slope (S) at pH 7.4 and logBB in A/B cluster and between the logKw at pH 9.1 and logPS in C/D cluster. CONCLUSIONS. Results of the QSPR, clustering and correlation studies could be used as novel tool for evaluation of blood-brain barrier permeation of related α-adrenergic/imidazoline receptor ligands. This article is open to POST-PUBLICATION REVIEW. Registered readers (see "For Readers") may comment by clicking on ABSTRACT on the issue's contents page.