QSAR Modeling of Histamine H3R Antagonists/inverse Agonists as Future Drugs for Neurodegenerative Diseases.

BACKGROUND: Histamine H3 receptor (H3R) is associated with several neuropsychological diseases, and thus it is an important target involved in several CNS disorders, such as narcolepsy, attention deficit hyperactivity disorder and schizophrenia. Since QSAR modeling is a feasible approach to explain the role of the molecular substituents in the biological activity, it can help in improving the design of better H3R ligands for these conditions.
METHODS: This article reviews papers previously published in literature to show the current status of the contribution from QSAR modeling to reach H3R antagonists/inverse agonists.
RESULTS: Classical and 3D-QSAR models were retrieved, showing that the steric and hydrophobic properties of the H3R ligands are most important to reach good affinity.
CONCLUSION: Although QSAR methods are valuable to design better H3R antagonists/inverse agonists, pharmacokinetics should also be considered in future models to ensure good CNS penetration. Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.org.

INTRODUCTION

The histaminergic system in the CNS is mainly modulated by histamine, a biogenic amine involved in several pathophysiological effects. The effects of histamine are produced through interaction with histamine G protein coupled receptors (GPCRs). There are four known subtypes of histamine receptors, H1 to H4 (H1R-H4R), differing in localization and mechanism of cellular signaling [1].

The H3R was discovered in 1983 by Arrang and colleagues, by an experimental observation that H1R antagonists did not show any response in this target and the H2R antagonists exhibited variable affinities not correlated to the H2R affinities, suggesting considerable differences between the classic receptors and H3R [2, 3]. Peripherally, H3R can be found in nerve endings of the gastrointestinal tract and the heart [4]. However, the H3R is found predominantly in the central nervous system (CNS) and is highly expressed in the cerebral cortex, basal ganglia [5], hippocampus [6], nucleus accumbens and substantia nigra [7]. The H3R is located in the presynaptic neurons, where it acts both as autoreceptor (modulating the synthesis and release of histamine in histaminergic neurons), and heteroreceptor in non-histaminergic neurons, regulating release of other neurotransmitters such as acetylcholine, norepinephrine, dopamine and serotonin [8]. H3R has different signaling pathways such as, inhibition of Na+/H+ exchanger activity, modulation of the MAPK pathway, activation of the AKT/GSK-3β axis and activation of phospholipase A2 [9]. However, signal transduction by H3R is primarily mediated by Gi/o protein, leading to reduction in intracellular cAMP concentration and influx of calcium into neurons, so it acts as inhibitory controller of neurotransmitters release [10].

Nowadays, there are at least six human isoforms for H3R, although its exact physiological role is still unclear. However, the CNS distribution of each isoform is considerably different, possibly leading to different pharmacological profile [11]. In addition, there are important differences in the binding affinities of H3R antagonists among species that are attributable to differences in some amino acids. While first-generation H3R antagonists were generally more potent at rodent receptors (including imidazole and non-imidazole compounds), more recent non-imidazole compounds are much more potent at human receptors [12]. These differences must be considered during the design process, since it can lead to human activity profile far from the predicted by the animal receptors.

Considering the large expression of H3R in the CNS and its ability to control the release of several distinct neurotransmitters, H3R has become an interesting target for bioactive molecules to treat neurological and psychiatric disorders such as Alzheimer's (AD) and Parkinson’s (PD) disease, epilepsy, hyperactivity, attention deficit disorder and schizophrenia [1, 13].

Since H3R acts through a negative feed-back mechanism, H3R antagonists/inverse agonists can be useful in increasing the release of neurotransmitters, thereby helping in the treatment of conditions involved in reducing neurotransmitter activity [14]. Recently, the H3R inverse agonist pitolisant or tripolisant (Wakix®, Fig. ) has been approved in the United States and the European Union for the treatment of narcolepsy with or without cataplexy. In addition, pitolisant is also being evaluated in clinical trials with indications for treatment of other disorders such as schizophrenia and some types of dementia [15]. With regards to neurodegenerative diseases, pitolisant has shown promising results in clinical trials for PD, especially in improving the excessive daytime sleepiness in PD patients [16].

Several antagonists/inverse agonists have also demonstrated pro-cognitive activity in cognitive deficit models. Administration of H3R antagonists in hypothalamic tuberomammillary nucleus increased the release of acetylcholine, dopamine and norepinephrine into the prefrontal cortex [17] and may cause both improvement and increased cognitive functions [18-20]. H3R antagonists have also shown to elevate acetylcholine levels in cortex and hippocampus, enhancing memory. However, preclinical studies have shown that H3R antagonists activate signaling pathways that may improve cognitive efficacy and disease-modifying effects beyond symptomatic alleviation in AD. For instance, administration of ABT-239 (a H3R antagonist, Fig. ) in normal mice showed increased cortical cAMP response element-binding (CREB) and glycogen synthase kinase 3β (GSK3β) phosphorylation, producing cognitive efficacy independently of increasing acetylcholine release [12]. In addition, some H3R antagonists such as ABT-239 and thioperamide (Fig. ), demonstrated neuroprotective action in vitro and in vivo in neurotoxicity models [20].

Clinical studies have shown that H3R antagonists may present therapeutic potential for neurodegenerative diseases. The H3R antagonist GSK-239512 (Fig. ) was developed for the treatment of various types of dementia with cognitive impairment. In phase I trials, it was evaluated with indication for AD with mild to moderate symptoms [21] and in phase II studies with indication for schizophrenia [22]. PF-03654746 (Fig. ) was evaluated for efficacy and safety in volunteers with excessive daytime sleepiness (EDS) associated with narcolepsy and presented significant improvement in symptoms, as pitolisant [23], and also completed phase I clinical trials in patients with schizophrenia, showing significant improvement in cognitive symptoms [24]. Other H3R antagonist, SAR-110894, was evaluated in phase II clinical trials for treatment of AD, in association with donepezil [25]. Finally, the compound MK-0249 was evaluated for the treatment of AD associated with cognitive deficits [26], and also in phase II studies for patients with paranoid schizophrenia (Fig. ) [27].

In summary, numerous advantages have been observed for H3R antagonists/inverse agonists that may translate into promising AD and PD enhancing agents. Considering all these therapeutic potentials, H3R ligands are constantly a focus in search of new chemical entities, and thus the structural requirements necessary for appropriate interaction with the receptor and the design of more active compounds have been constantly reported. Considering this, the quantitative structure-activity relationships (QSAR) strategy can be highlighted. In addition to the appropriate pharmacodynamic requirements, the QSAR approach may also assist in predicting the pharmacokinetic profile of the designed molecules [28]. Pharmacokinetics is among the most important bottlenecks in the drug discovery, since several compounds that fail in the preclinical and clinical phases do not fulfill the ADME (absorption, distribution, metabolism and excretion) requirements or present inappropriate significant toxicity [20, 29].

In this paper, we performed a review of the main structural requirements to improve H3R antagonist/inverse agonist activity as suggested by quantitative structure-activity relationships (QSAR) models reported in literature. This data is quite important to guide a more rational design of new H3R ligands with better activity profile. The reports were raised in PubMed and Web of Science databases, using the terms “QSAR” and “H3 ligands” as keywords.

QSAR

Several H3R ligands were identified and evaluated as potent H3R antagonists and different approaches trying to explore structural determinants for H3R antagonistic activity. The QSAR approach is a feasible technique for identifying the important structural determinants to a defined biological activity. QSAR is a statistical model that has been successfully used in agrochemical, toxicological and environmental studies beyond drug discovery [30]. Accordingly, QSAR approach was also explored to determine H3R affinity, activity and selectivity.

The QSAR approach is based on the assumption that similar molecules have similar activities, and considering that similarity refers to chemical characteristics, the activity of a molecule is dependent of certain physicochemical properties, called as descriptors. In medicinal chemistry, this assumption is known as structure-activity relationships (SAR).

Through statistically validated mathematical models, the QSAR modeling allows to predict and quantify the SAR variables, present in a series of analogues with known quantitative activity, responsible for the affinity and activity in a specific target. In these mathematical models, the different structural, physicochemical and conformational properties of the molecules are expressed by structural descriptors (independent variables). The different types of descriptors are obtained experimentally or calculated in silico, and evaluated for correlation with previously known activity of the compounds (dependent variable) [31, 32].

The QSAR studies started with the work of Corwin Hansch group through the investigation of the role of hydrophobicity in the biological activity of compounds. In later works, they also added the influence of electronic and steric properties on the model. The earlier QSAR models are known as Hansch's analysis [33].

To present reliable results, it is necessary that these models should be statistically validated. Basically, the validations of the QSAR models are subjected to statistical model validation parameters such as degree of adjustment, degree of significance and degree of predictability [34]. The analysis of adjustment of the model can be performed by calculating the correlation coefficient (R), standard deviation (s) [31] and other parameters such as standard error (SE) of prediction (standard error, SEP) or estimated (standard error of estimate, SEE) and its deviation (standard deviation error in prediction, SDEP). The statistical significance of a model can be evaluated using, per example, the multiple determination coefficient (R) and Fisher's test (F) [34]. The predictability can be determined by cross-validation processes, through evaluation of cross-correlation coefficient (Q), leave-N-out, scrambling of dependent variable and also through prediction of the activity of a test set [32].

In addition to the classical QSAR modeling (2D-QSAR), QSAR studies can enhance the steric factors through a three-dimensional approach (3D-QSAR), such as comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA). These approaches allow to obtain descriptors (such as Coulomb and Lennard-Jones potentials, hydrogen bonding and others) in a three-dimensional environment, to identify possible points of interaction of a ligand in the space, and possibly to the pharmacological target responsible for the biological activity [35].

As described, the QSAR modeling allows to obtain simple and fast models that can be useful in explaining the dependence of the activity on certain descriptors, and more importantly, in predicting the activity of new and non-tested compounds. However, as any statistical model, mispredictions can happen and it is important to understand the applicability domain of each model. Therefore, the use of dataset obtained by different methodologies or animal species, choice of redundant descriptors, employment of non-sense descriptors, inadequate number of variables in each model and others mistakes should be observed [33].

QSAR STUDIES ON H3R ANTAGONISTS

Table presents the descriptors, with the meanings and nature, addressed throughout the review.

There are few reports on QSAR studies with H3R antagonists in literature. Among them, one of the sets evaluated was comprised of acylated histamine derivatives (pKi 5.7-7.3) previously reported by Stark et al. [36]. Agrawal and colleagues presented a classical QSAR model using topological descriptors and obtained through multiple regression analysis (MRA). The MRI descriptor was considered essential to attribute activity and to prediction power of the model. The addition of a qualitative indicator parameter (Ip) in the presence of a benzene ring in R4 in the acylated histamine structure of the compounds (Fig. ) also resulted in improvement in the obtained model. Thus, it was proposed that the presence of an aromatic ring in these kinds of compounds is essential to good activity [37]. However, the number of molecules used to build the models was relatively low (n = 13), leading to quite simple models. Since QSAR modeling is a ligand-based approach, representative number of molecules is indispensable to statistical significance. To avoid correlations by coincidence using QSAR approach, it is suggested to use at least 5 compounds per descriptor in the final model obtained through classical MRA [34]. The obtained statistical parameters for the best model were: SD = 0.2199, R = 0.9281, F = 12.430, R = 0.7921 and Q = 0.4221.

3D-QSAR models (CoMFA e CoMSIA) were built using a set of 144 H3R antagonists (118 in the training set, 26 in the test set) with general structure consisting on an heterocycle linked to a phenyl group and a heteroaromatic group (Fig. ). The compounds were aligned based on the center of the heterocycle to obtain de CoMFA models. The 3D descriptors used were of electrostatic, hydrophobic and donor of hydrogen bonding nature, in addition to added descriptors such as HOMO and LUMO. The use of additional descriptors was important for the improvement of 3D-QSAR models (CoMFA: q = 0.721, r = 0.931, SEE = 0.236 e CoMSIA: q = 0.700, r = 0.921, SEE = 0.252) and better predictive capacity on the contribution of the descriptors in H3R affinity. The CoMFA results showed that the presence of bulky groups near the heterocycle and especially, near the phenyl group, is favorable for the activity. Moreover, electronegative and positively charged groups in the heterocycle and in the phenyl group, respectively, may increase the affinity for H3R. The CoMSIA models showed that hydrophilicity in the heteroaromatic region is favorable to activity, whereas in the heterocyclic region, the H3R affinity may be increased through hydrogen bonding acceptors and hydrophobic groups (Fig. ) [38]. In fact, hydrophobic groups may constitute an additional point of interaction due to the presence of a lipophilic pocket between the TM3 and TM6 (transmembrane domains) in the H3R [39], and compound presented in Fig. ( is an example of high activity compound.

Fig. (4)

Example of H3R ligand used by Chen [38] for CoMFA and CoMSIA studies.

A set of 58 (44 in the training set) arylbenzofuran derivatives (pKi 7.2-10.1) previously reported by Gfesser et al. [40] including ABT-239 was also used to built QSAR models which could help to improve the affinity for both human and rat H3R [41]. The final model was fully (internal and external) validated. The model is dependent on four descriptors (EHOMO, logDpH7.4, Mor19V e Mor30M) to explore the structural requirements of this derivatives to achieve adequate human H3R affinity and presented relatively satisfactory statistics (r = 0.754, F = 40.7, SE = 0.317, q2LOO = 0.71). The negative correlation between the descriptor EHOMO and pKi indicates that lower EHOMO values lead to higher H3R affinity. The researchers observed that EHOMO was highly correlated to the electronic density in the aromatic system, suggesting that it can perform specific interactions in the receptor. Moreover, more hydrophilic substituents may lead to higher affinity molecules, as indicated by logDpH7.4. Molecular docking studies were used to explore and validate the obtained models. Human H3R homology model from the bovine rhodopsin structure was built and the compounds were docked using known interaction points between H3R and histamine. The selected interactions were between Asp114 (the most important in biogenic amines receptors) and the aliphatic NH2 group, and between Glu206 and the Nt imidazole ring of histamine. The results suggested that the compounds may interact through charge transfer between the phenyl or benzofuran rings of the ligands and the hydrophobic π system comprised by Ala122, Phe211, Val214, Trp371, Phe367, Phe398 and Leu401 amino acids (Fig. ). The QSAR model also showed unfavorable correlation of lipophilicity with molecular mass (Mor30M) and positive correlation with molecular volume (Mor19V), suggesting that bulky hydrophilic substituents may improve the H3R affinity of the ligands [41]. It must be stressed that although hydrophilicity is favorable to H3R affinity, it may be detrimental to pharmacokinetics, since it is expected that the compounds reach the CNS and they must be capable of crossing the blood-brain barrier (BBB) indeed. In addition, more hydrophilic compounds can have poor oral absorption, leading to low bioavailability in vivo. Fortunately, all the compounds from the set were quite lipophilic, but a threshold for hydrophilicity of the substituents should be observed in the design of new molecules. The QSAR model to rat H3R affinity was even better (r = 0.840, F = 69.8, SE = 0.288, q2LOO = 0.81), however the descriptors present in both human and rat models were different, showing the importance of specie-related pKi value. In summary, when performing QSAR studies, species variability regarding pKi values must be considered as key information for prediction of affinity for H3R.

Da Costa and Trsic [14] reported QSAR models to a set of 28 4-phenoxypiperidine derivatives (pKi 5.27-9.20) previously published by Dvorak et al. (Fig. ) [42]. Quantum-chemical descriptors were used to investigate the contribution of electronic characteristics in the binding affinity to H3R. The calculated descriptors have presented high correlation with the antagonist activity, classifying the compounds into two distinct groups with higher and lower activity by using hierarchical cluster analysis (HCA) and principal component analysis (PCA). The presence of εFERMO (frontier effective-for-reaction molecular orbital energies) in the model showed important statistic contribution (R = 0.927, F = 80.11, SEP = 0.141) with high prediction power (88%) of the pKi values. The descriptors HOMO and LUMO were also evaluated in some built models [43], however these descriptors showed poor correlation to the activity than εFERMO descriptor, possibly due to limitations of the HOMO and LUMO model of electron transfer. The model showed that molecules with high nucleophilic characteristic may have higher H3R affinity, since εFERMO value showed favorable contribution to pKi. On the other hand, the electrophilicity (ω descriptor) showed unfavorable contribution to the affinity. In addition, the authors demonstrated that higher electronic density in the piperidine nitrogen atom (N1) showed higher activity (pKi), showing importance of such atom for interaction with the receptor.

Several H3R ligands show the imidazole ring present in histamine. Imidazole nitrogen interacts with H3R through the proton transfer [38], being crucial for the activation of H3R. Since the blockade of the receptor is more relevant than activation to improve neurodegenerative diseases, this feature is undesirable to the development of H3R ligands with good pharmacological profile. Piperidinylurea derivatives (Fig. ) from a set of 15 compounds were employed in QSAR studies. The activity was evaluated as pIC50, and the biparametric QSAR models were validated by different techniques, presenting good statistical quality (R = 0.8820, R = 0.7780, AdjR = 0.7409, F(2,12,0.05) = 21.0210, SEE = 0.2988, t(12,0.0005) = 13.7890, X(0.05) = 0.1438 and p = 0.0000). The negative correlation of the vsurf_DD13 and vsurf_Wp4 values indicated that small distance between polar/hydrophobic groups and the atom of the target molecule, as well as superficial hydrophobicity of the compounds, are favorable to H3R antagonistic activity [44]. Moreover, considering these compounds also have considerable human ether-a-go-go-related gene (hERG) channel blocking activity and that the properties involved in hERG activity are quite different from those involved in H3R activity, the obtained model can be considered quite useful in improving H3R affinity with lower hERG blockade activity. The hERG channel blockade is involved in dangerous cardiac arrhythmias which means high toxicity to humans, and thus must be avoided. In summary, lower distance between polar and aromatic/hydrophobic groups should be small to avoid hERG activity, as well as more hydrophobic molecules can have better H3R affinity than hERG activity. Polar groups on the surface area also increase hERG activity.

The role of lipophilicity (logP) and basicity (pKa) were evaluated with a set of 11 2-aminobenzimidazole (Fig. ) derivatives [45]. It is important to stress that lipophilicity also increases the binding to plasma proteins, affecting the distribution of the compounds in vivo. The pKa can also be correlated to ionization state, that can influence in both pharmacokinetics and receptor affinity. The pKi values of compounds (pKi 6.91-9.37) were measured through binding assays in rat receptors from cerebral cortex synaptosomes. The obtained QSAR models were not considered statistically satisfactory by using classical MRA. When using the partial least square (PLS) analysis, better but not sufficient model was obtained. Although the QSAR data was not conclusive, some SAR information regarding these compounds was obtained. Protonation of 2-aminobenzimidazole nucleus is unfavorable to basicity, suggesting these compounds may interact in the neutral form with the receptor. In addition, the lipophilicity given by the substituent in the benzimidazole nucleus is favorable to the H3R affinity. The work also suggests that optimal lipophilicity is around logP 2.4, as indicated by the second-order descriptor logP2 [45].

In another work, a set of 38 H3R antagonists (pKi 4.05-9.89) was evaluated through 3D-QSAR (CoMFA and CoMSIA methods) which employs steric, electrostatic, lipophilic, hydrogen bonding acceptor and hydrogen bonding donor descriptors. The compounds had as common structure a polar heterocycle group (containing an imidazole or thiazole group) attached to an imidazole ring by an alkyl spacer. Initially the models were built using the imidazole, the terminal apolar group and a hydrogen-bond donor group of thioperamide in the alignment. However, this led to poor models, suggesting the role of hydrogen-bond is different. Using the nitrogen atom with a lone electron pair present in all the compounds instead of the donor group, the models improved significantly and then this alignment was used. The best model was obtained (CoMSIA: n = 38, LVs = 3, Q = 0.85, SDEP = 0.45, R = 0.91, s = 0.38) using only the steric field as descriptor, proving to be determinant for H3R affinity. It is possible to verify that increasing volume in the polar region of the heterocycle can improve the affinity, but the substitution in the benzene ring with bulky groups should be avoided (Fig. ). In addition, the presence of donor groups and hydrogen bonding in the polar region of the molecule did not demonstrate statistical significance in the H3R affinity [46]. It is interesting that ABT-239 (Fig. ) is quite similar to the compounds used in this set, but it presents a bulky 4-cyanobenzene group in this prohibited region and also present high affinity. Maybe benzofurans can be considered outliers to this model, but it limits its application domain.

A set of monosubstituted benzyl analogues of thioperamide (Fig. ) was synthesized and also evaluated through classical QSAR modeling to obtain better tools for positron emission tomography (PET) and single photon emission computed tomography (SPECT) applications, and halogens are suitable atoms to reach this objective [47]. In this work, the intrinsic activity (pA2 value) was used as dependent variable for H3R activity. The synthesized thioperamide analogues have shown clear influence of the substituent in the benzyl group in the activity, correlated to steric and electronic factors. To evaluate this quantitatively, steric (the dihedral angle φ) and electronic (the electron density δ) descriptors were calculated. The best model (n = 13, r = 0.93, s = 0.28, F = 31.57) exhibited influence of both descriptors. The angle between the phenyl group and the isothiourea group (φ) showed minor contribution to explain the activity, whereas the electron density of the substituted carbon of the phenyl ring (δ) was highly correlated to the pA2 values. The results suggest that substituents which lead to more negative charge on the carbon linked to them (i.e. when linked to iodine) may have better affinity for H3R, and ortho-substitution can also improve the activity due to its steric influence on φ. Although the model is limited due to the small size of the dataset, the residual values from the predicted activities were quite low, and thus the model brought important information allowing some preliminary evaluation to this set of compounds.

H3R has considerable homology with H4R (~ 37% total, 68% within TM) [48], and so obtaining molecules with selectivity towards each is a hard task. Due to the similarity between them, it is likely that several compounds that bind at one receptor may have considerable affinity for the other. The proposed pharmacophore model for H3R ligands may also be applicable for H4R ligands [1], making it difficult to achieve selective compounds. In fact, several ligands designed to bind to H3R present considerable affinity for H4R, such as thioperamide and clobenpropit (Fig. ). Lim et al. [48] used a set of 22 clobenpropit analogues (pKi 6.5–8.6) to determine the main characteristics that play the role in the affinity for both H3R and H4R. Electronic, hydrophobic and steric descriptors were calculated and used to search valid QSAR models. Two final models were obtained, each for H3R affinity (r = 0.982, R = 0.964, S = 0.099, F = 115.091, F = 2.965, q = 0.946) and for H4R affinity (r = 0.946, R = 0.894, S = 0.166, F = 35.994, F = 2.965, q = 0.801). It was verified that the descriptors which correlate with the affinity for H3R also show high correlation with H4R affinity and vice-versa, proving that it is really difficult to design clobenpropit analogues with considerable selectivity towards either receptors. The only descriptor with positive correlation to H4R affinity that had shown inverse correlation with H3R affinity was the energy of bond stretch-bend (E_stb), although its contribution alone is not sufficient to achieve such objective. Lower E_stb values may favor affinity for H4R without increasing the affinity for H3R [48].

Finally, another study aiming determinants of selectivity was also performed [49]. In this work, QSAR models were constructed for a series of 14 compounds derived from clobenpropit with an additional lipophilic moiety linked to the thiourea group. The authors mixed classical and 3D-QSAR models, as well as data obtained from docking of clobenpropit and VUF5228 in both H3R and H4R. The obtained models showed that the higher topological diameter increases the affinity for H3R, corroborating with data from other QSAR models, and also BCUT_SMR_0 descriptor (r = 0.900, R = 0.810, S = 0.312, F = 23.454, F = 3.982, q = 0.699). The mixed classical and 3D model (r = 0.955, R = 0.912, S = 0.277, F = 57.249, F = 3.982, q = 0.862) furnished some information regarding selectivity towards H3R, by using the difference of activity between both receptors (ΔpKi = pKi H3R - pKi H4R). The selectivity can be explained by 3D interaction probes which may represent different interaction points in the receptors (represented by compound VUF5228, Fig. ).

These probes (DRY) represent hydrophobic interactions with 4-chlorophenyl group and the additional cyclohexyl group, and these groups led to higher H3R selectivity. By combining the results from 3D QSAR with data from homology modeling of H4R, the result suggests that Asn4.57 and Glu5.46 residues are determinants of H4R selectivity. However, the 3D models are not useful in explaining which substitution pattern may lead to higher H3R selectivity, although the additional cyclohexyl led to higher affinity for H3R.

CONCLUSION

Many researchers have been directing efforts to obtain models that explain how molecules can be modified to achieve higher H3R affinity, antagonistic activity and selectivity over other targets. Among the several reports in literature presenting compounds which were designed and evaluated through SAR analysis, mainly piperidine, benzofuran, piperidinylurea and imidazole series were deeply studied through QSAR approaches. The data from these reports suggest that the steric and hydrophobic factors play the most important role in the H3R activity. Important features regarding the interaction with H3R have been described, but it must be adequately explored to improve both the affinity and the selectivity. However, considering the wide range of compounds reported, it is noted that the number of QSAR studies in the literature is still scarce. In addition, with some exceptions, such models focus only on the pharmacodynamic aspect of the biological activity. Considering that the H3R is mainly found on CNS, pharmacokinetic and toxicological aspects such as BBB penetration, metabolic stability, distribution coefficients and promiscuity should be used as dependent variable in QSAR studies. Pharmacokinetics is among the main bottlenecks to reach clinical trials, and it must be predicted in the earlier stages of drug discovery process to maximize the chances to reach the patient bedside.

Consent for Publication

Not applicable.

Table 1

List of descriptors from reported QSAR models for H3R activities and their meaning.

Descriptor	Meaning	Nature
MRI	Molecular redundancy index: indicates the capacity and symmetry of a molecule	Topological
Ip1	Indicator parameter for the presence or absence of one benzene ring	-
EHOMO	Highest occupied molecular orbital energy	Thermodynamic
logDpH7.4	Distribution coefficient in pH 7.4	Hydrophobic
Mor19V	3D atomic coordinates by the transform used in electron diffraction studies weighted by volume	3D-Steric
Mor30M	3D atomic coordinates by the transform used in electron diffraction studies weighted by mass	3D-Steric
εFERMO	Frontier effective-for-reaction molecular orbital energy	Thermodynamic
ω	Electrophilicity index	Electronic
vsurf_DD13	Distance between polar (H2O probe) and hydrophobic (DRY probe) groups to the receptor	3D-Mixed
vsurf_Wp4	Polarizability on the van der Waals surface	3D-Mixed
φ	Dihedral angle	Steric
δ	Electron density	Electronic
E_stb	The bond stretch-bend cross-term potential energy descriptor calculated from stored 3D conformations	Thermodynamic
DRY	Interaction energy value of hydrophobic probe	3D-Mixed
BCUT_SMR_0	BCUT descriptors using atomic contribution to molar refractivity	Steric-electronic