Next Article in Journal
Updates on Aptamer Research
Next Article in Special Issue
Quantitative Structure–Activity Relationships for Structurally Diverse Chemotypes Having Anti-Trypanosoma cruzi Activity
Previous Article in Journal
The Effects of Proteasome Inhibitors on Telomerase Activity and Regulation in Multiple Myeloma Cells
Previous Article in Special Issue
Structure-Based Approach for the Prediction of Mu-opioid Binding Affinity of Unclassified Designer Fentanyl-Like Molecules
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Three-Dimensional Quantitative Structure-Activity Relationships (3D-QSAR) on a Series of Piperazine-Carboxamides Fatty Acid Amide Hydrolase (FAAH) Inhibitors as a Useful Tool for the Design of New Cannabinoid Ligands

1
Escuela de Quimica y Farmacia, Facultad de Medicina, Universidad Andres Bello, Quillota 980, Viña del Mar 2531015, Chile
2
Pharmacy Department, Faculty of Chemistry, Pontificia Universidad Católica de Chile, Vicuña Mackenna 4860, Santiago 7820436, Chile
3
Departamento de Ciencias Farmacéuticas, Facultad de Ciencias, Universidad Católica del Norte, Avenida Angamos 0610, Antofagasta 1270709, Chile
4
Instituto de Química y Bioquímica, Facultad de Ciencias, Universidad de Valparaíso, Av. Gran Bretaña 1111, Valparaíso 2360102, Chile
5
Centro de Investigación Farmacopea Chilena (CIFAR), Universidad de Valparaíso, Santa Marta 183, Valparaíso 2360134, Chile
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(10), 2510; https://doi.org/10.3390/ijms20102510
Submission received: 21 February 2019 / Revised: 27 April 2019 / Accepted: 7 May 2019 / Published: 21 May 2019
(This article belongs to the Special Issue QSAR and Chemoinformatics Tools for Modeling)

Abstract

:
Fatty Acid Amide Hydrolase (FAAH) is one of the main enzymes responsible for endocannabinoid metabolism. Inhibition of FAAH increases endogenous levels of fatty acid ethanolamides such as anandamide (AEA) and thus consitutes an indirect strategy that can be used to modulate endocannabinoid tone. In the present work, we present a three-dimensional quantitative structure-activity relationships/comparative molecular similarity indices analysis (3D-QSAR/CoMSIA) study on a series of 90 reported irreversible inhibitors of FAAH sharing a piperazine-carboxamide scaffold. The model obtained was extensively validated (q2 = 0.734; r2 = 0.966; r2m = 0.723). Finally, based on the information derived from the contour maps we designed a series of 10 new compounds with high predicted FAAH inhibition (predicted pIC50 of the best-proposed compounds = 12.196; 12.416).

Graphical Abstract

1. Introduction

The endocannabinoid system (ECS) remains a highly relevant topic in the scientific community as it is involved in several regulatory actions and pathophysiological conditions [1]. Current available knowledge suggests that the ECS is constituted by the cannabinoid receptors, type 1 and 2, the main endogenous ligands anandamide (AEA) and 2-arachidonyl glycerol (2-AG) as well as the enzymes that participate in their biosynthesis (N-acyl phosphatidylethanolamine phospholipase D, or NAPE-PLD) and degradation (Fatty Acid Amide Hydrolase or FAAH and Monoacylglycerol Lipase, or MAGL) [2].
FAAH is an integral membrane protein of ~60 kDa (579 amino acids) that belongs to the amidase family of enzymes [3]. It exists as a dimer in its membrane-associated form [3] and is highly expressed in the brain, liver, kidney, and testis. While most other mammalian serine hydrolases possess a Ser-His-Asp catalytic triad, FAAH has a distinctive Ser-Ser-Lys triad [4].
In vitro studies show that FAAH can inactivate various lipid amides and allows modifications in both the amide headgroup and lipid acyl chain. In addition to AEA, other bioactive lipid amide substrates include related N-acylethanolamine (NAE), oleamide (cis-9,10-octadecanoamide), and N-acyl taurines (NATs) [4].
Endogenous FAAH substrates such as AEA play key regulatory functions in the body and have been implicated in a variety of pathological disorders, including pain, inflammation, anxiety [5], depression and vascular hypertension [6]. Therefore, inhibition of FAAH represents a rational therapeutic approach to treat conditions where higher endocannabinoid activity can be beneficial. Furthermore, as opposed to direct cannabinoid activation, enzyme inhibition offers spatio-temporal control, increasing endocannabinoid activity only at the sites where lipid signaling molecules are being produced. Accordingly, animal studies showed that FAAH inhibitor URB597 elevated endocannabinoid tone without producing motor side effects [7]. For this reason, diverse FAAH inhibitors have been developed [6].
The first-generation of FAAH inhibitors were designed in order to covalently bind to the catalytic residue Ser241 [3,4]. Despite their ability to block FAAH in in vitro pharmacological assays, they remained poor candidates for preclinical studies due to their lack of selectivity [4]. Subsequently, FAAH inhibitors with significantly improved selectivity were developed, including carbamates (ORG-231295), α-ketoheterocycles (OL-135) carbamoyl tetrazoles (LY-2183240), benzothiazole derivatives and piperidine/piperazine ureas [4] (PF-3845, PF-04457845) (Figure 1).
Some quantitative structure-activity relationships (QSAR) studies have been performed on the cannabinoid system (CB1 and CB2 receptors) [8,9,10,11,12,13], however, few structure-activity relationship studies have been performed on FAAH inhibitors, and mostly with carbamate-type structures. Dainese et al. calculated theoretical molecular descriptors in a series of naturally occurring FAAH inhibitors [14]. Käsnänen et al. reported the synthesis and 3D-QSAR studies of carbamate inhibitors [15]. Mor et al. constructed 2D-QSAR equations that could explain the inhibition activity of biphenyl-alkylcarbamates. [6]. Vacondio et al. developed structure-property relationships to explain the hydrolytic stability of carbamate inhibitors [16]. Han et al. reported a comparative molecular field analysis (CoMFA) study on a series of oleoylethanolamide structure inhibitors [17]. To date, there are no 3D-QSAR studies of irreversible inhibitors with the piperazine-carboxamides structure. This type of general structure was shown to have good physical and pharmacokinetic properties and has been reported to be capable of elevating plasma concentrations of AEA, AEP, and AEO in rats [18]. For this reason, the formulation of a QSAR model for the design and prediction of FAAH inhibitor activity based on this structural moiety is significant from a pharmacological point of view.
In the present work, three-dimensional quantitative structure-activity relationships (3D-QSAR) studies based on comparative molecular similarity indices analysis (CoMSIA) were carried out on a set of various reported urea-based FAAH inhibitors. The aim of our 3D-QSAR is to derive useful binding information in order to guide the design of future FAAH inhibitors. The importance of steric, electrostatic and hydrogen-bond characteristics can be analyzed by aligning similar analogues based on key pharmacophoric features [19]. Knowledge of binding requirements can then be used to derive predictive 3D-QSAR models that can, in turn, aid in the design of new inhibitors.

2. Results and Discussion

2.1. Statistical Results

The statistical results for CoMSIA are presented in Table 1. All possible field combinations were tested for both CoMFA and CoMSIA. In the case of CoMFA, no combination was statistically significant. The CoMSIA models with the highest q2 values were those that considered the field combinations SEDA, EDA, EHDA, and SEHDA. The SEDA and EHDA models presented a donor hydrogen-bond contribution of 0.099 and 0.093 respectively. While in the EDA model, the H-bond donor contribution was 0.111 versus 0.889 of the Electrostatic and H-bond Acceptor contributions. The imbalance in the field contribution of these models made us discard them. The final selected model SEHDA, presents a good balance between the field contributions, a high value of q2 (0.734) and r2ncv (0.937) and higher value of F (138.36). Extensive additional validation was carried out to test the predictive quality of this model.
Table 2 presents a summary of the external validation parameters of the CoMSIA model. The model has a high value of r2 (0.966), which indicates satisfactory external predictive power. Nevertheless, as described by Golbraikh and Tropsha, high values of q2 and r2 (conditions 1 and 2) are not sufficient for model validation. For a reliable predictive capacity, the line of experimental versus predicted activity should be as close as possible to the line y = x. This is observed when the conditions (3a or 3b), (4a or 4b) and (5a or 5b) shown in Table 2 are fulfilled. Condition 6 or r2m metrics, measures the agreement between the observed and predicted activity. Our derived CoMSIA model satisfied all the conditions for internal and external validation. In addition, we performed the calculation of several external validation descriptors (conditions 7 to 12) [20,21,22]. In all cases, our model passed the validation tests.
Furthermore, the Y-randomization test [23] (Table 3) was applied to assess the robustness of the model (see Table S1 of the Supplementary Material) as previously described by Lorca et al. [24]. The obtained QSAR models show low q2 and r2ncv values (Table 3).
The values of experimental activity, predicted activity, and the residual values for the best CoMSIA model are shown in Table 4. All the compounds showed low residual values and deviations of the predicted activity greater than a logarithmic unit were not observed. Figure 2A shows a plot of experimental versus predicted activity and the data distribution is close to the y = x line. The model shows a good balance in terms of predictive capacity. Forty-two compounds showed negative residual values and 48 presented positive deviations (Figure 2B). The residual range was from −0.82 to 0.89. As shown in Figure 2C the CoMSIA model shows a satisfactory predictive capability throughout the whole data set (training and test set) as well as a good predictive power for both, less active (1, 6 and 7) and most active compounds (65, 66, and 67).

2.2. Applicability Domain

In this work we used the method developed by Roy et al. [25] for determination of applicability domain (AD) as previously described by Lorca et al. [24].
The calculation was carried out using the free application available on the author’s page and all compounds were found to be within the domain of applicability, except compounds 83 and 87. These two compounds are the only ones bearing an imidazopyridine or imidazopyrimidine ring and the only difference between them is the position in which the heterocycle is connected to the urea moiety. For this reason, compounds with these heterocyclic systems connected to the urea moiety were not proposed as new molecules.
In summary, the CoMSIA model generated here presents good internal and external validation parameters (q2 = 0.734; r2 = 0.966), and meets the validation criteria of Tropsha and Roy (r2m = 0.723). All the molecules studied are within the applicability domain (except compounds 83 and 87) and the model was validated by the Y-randomization test. Therefore, reliable information can be extracted from the contour maps as discussed in the next section.

2.3. Contour Maps Analysis

The result of a 3D-QSAR study can be visualized graphically unlike a traditional 2D-QSAR equation. Contour maps represented by colored polyhedrons can be seen around the molecules. The maps obtained in our study correspond to the steric, electrostatic, hydrophobic, H-bond donor and H-bond acceptor contour maps. Regions where a molecular property is favorable or unfavorable are indicated by different colored polyhedrons. Figure 3 presents the different maps around the most active compound (66, pIC50 = 10.602; on the left) and the least active compound (6, pIC50 = 5.176; on the right).

2.3.1. Steric Contour Map

The steric contour map shows a large yellow polyhedron close to the pyridazine ring of the most active compound indicating that bulky substituents in this region should be avoided in order to favor biological activity (Figure 3A,B). Alternatively, a smaller five-membered ring could replace the pyridazine ring following the same steric requirement. This can be seen in the proposed molecules (Table 5) where the analog 2x shows a considerable increase in predicted activity when the pyridazine ring is replaced by a pyrazole ring. This relation is further supported by compounds 54 and 55 (Table 6) which bear phenyl substituents in the corresponding position and show low activity consistent with their bulkier nature.
On the other side, the yellow polyhedron that surrounds the pyrimidine ring indicates that reducing the size of this ring or replacing it with a smaller linker while maintaining the electronic properties would be beneficial for activity. Additionally, the green polyhedrons around the ortho and meta positions of the disubstituted benzene ring indicate that bulky substituents in these positions can be favorable for activity. This can also be seen in Table 5 where substitution with a methyl group considerably increases the predicted activity in all the proposed molecules. For the least active compounds (Figure 3B), the steric factor by itself does not seem to explain the lower activity values.

2.3.2. Electrostatic Contour Map

Regarding the electrostatic contour maps, the red polyhedrons around the pyridazine ring (Figure 3C) highlight the importance of nitrogenated heterocycles that can confer electronegative areas in this region. This may explain why molecules 1, 6 and 7 (Table 6) bearing benzene rings with more homogeneous charge distribution show lower activity. The blue polyhedron inside the pyridazine ring shows that an electropositive center is beneficial for activity, therefore, incorporating electro withdrawing substituents in positions 5 and 6 of the pyridazine ring could increase activity. This electron distribution with an electron rich edge and electron deficient center suggests possible pi-stacking or pi-cation interactions with the target enzyme. Similarly, the expansion of the blue contour at position 4 of the pyridazine ring indicates that electron withdrawing substituents particularly at this position are favorable for activity. In agreement with this, proposed analogues 1x, 9x and 10x (Table 5) that follow this substitution pattern display high predicted activity. The polarization of the carbon atom directly attached to the electroattractive or electronegative groups nitrile and fluorine lower the electron density right where the blue polyhedron lies.
On the other side, the red contour over the pyrimidinic nitrogen (Figure 3C) shows the importance of an electronegative atom at this position as is present in the most active molecules (65, 66, 67 and 68 from Table 6) and absent in the least active one (compound 6). Likewise, an electron rich benzene ring is favorable and thus replacing the fluorine atom for an electrodonating group would be recommended. Accordingly, the proposed molecules 9x and 10x substituted with electrodonating methylene groups show the highest predicted activities. In Figure 3D the red contour over the electron deficient nitrile carbon atom may explain the lower activity of this analog.

2.3.3. Hydrophobic Contour Map

The hydrophobic contour maps (Figure 3E,F) show a gray polyhedron over the pyridazine nitrogen atoms of the most active molecules, indicating that incorporating hydrophilic atoms in this region can favor activity. The yellow polyhedron over the pyrimidine ring shows that a hydrophobic linker region is important. Additional yellow polyhedrons surrounding the disubstituted benzene ring suggest hydrophobic substituents in the meta and ortho positions can also increase activity. Following both the hydrophobic and the previously mentioned steric requirement all proposed molecules (Table 5) possess an ortho-methyl substituent in the benzene ring.

2.3.4. Donor and Acceptor Contour Maps

In the hydrogen bond donor map (Figure 3G,H) cyan polyhedrons surrounding the pyridazine ring suggest that incorporation of hydrogen bond donor groups can favor activity. For this reason, the proposed molecule 3x (Table 5) was designed with a hydroxyl group able to form hydrogen bonds. Furthermore, cyan polyhedrons around the urea linker suggest that the urea moiety is involved in a hydrogen bond interaction with the target enzyme.
Finally, the hydrogen bond acceptor map shows red polyhedrons next to position 4 of the pyridazine ring (Figure 3I,J), position 5 of the pyrimidine ring and over the urea carbonyl. Indicating that hydrogen bond acceptor groups in these positions is unfavorable for activity. Therefore, using different linker groups without hydrogen bond acceptor groups could be advisable in order to design new inhibitors.

2.4. Design of New FAAH Inhibitors

Based on the information obtained from the contour maps, we have designed a series of compounds evaluating multiple combinations of fragments. Substituents and functional groups were proposed, taking into consideration the electronic, steric, hydrophobic and hydrogen bonding properties suggested by the contour maps. Table 5 shows the compounds that presented the best predicted inhibitory activity. All proposed molecules have better predicted activity than the most active compound in the series (66, pIC50 = 10.602). In general, the presence of 6-member rings or fused systems on the left side did not greatly increase activity (1x pIC50 = 10.889; 4x pIC50 = 10.990). However, the insertion of a pyrazole ring generated derivatives with a significant increase in the pIC50 value (best compounds: 9x, pIC50 = 12.416; 10x, pIC50 = 12.196). This is because the pyrazole system meets the electronic, hydrophilic and hydrogen bonding requirements suggested by the contour maps. On the right side, the insertion of halogens and alkyl groups slightly increased the activity.
Due to the reported potential toxicity of heterocyclic compounds similar to those presented in Table 5 [26], we conducted a predictive toxicity study using the PreAdmet online platform (https://preadmet.bmdrc.kr/adme/). The calculation showed that only compound 1x presents a high risk of toxicity mediated by affinity to the anti-target hERG.

3. Materials and Methods

3.1. Molecular Alignment

CoMSIA studies were performed with Sybyl X-1.2 software (1.2, Tripos International, St. Louis, MS, USA) [27] installed in a Windows 10 environment on a PC with an Intel core i7 CPU. Each compound was drawn in ChemDraw and then geometry optimized using MM2 molecular mechanics as implemented in ChemBio3D software (15.1.0, PerkinElmer, Waltham, MA, USA). The structures were further minimized by Tripos force field implemented in Sybyl. MMFF94 charges were assigned to each atom. Among the techniques to perform molecular alignment are: (1) atom-by-atom alignments using a common fragment, (2) rigid alignment that minimizes rms distance, (3) flexible alignments [28] and (4) receptor guided alignments. In the present study, the first two options were carried out (Figure 4). In the first case, the piperazinyl-urea nucleus was chosen as the common scaffold for alignment. In the second case, the alignment was carried out using the Distill rigid protocol, as implemented in Sybyl. The best results in terms of statistical validation (q2, N) were for the atom-by-atom alignment, which was chosen as the basis for the development of the models.

3.2. CoMSIA Field Calculation

To derive the CoMSIA descriptor fields, the aligned training set molecules were placed in a 3D cubic lattice with grid spacing of 2Å in x, y, and z directions such that the entire set was included in it. The CoMSIA analysis, the standard settings (probe with charge +1.0, radius 1Å, hydrophobicity +1.0, hydrogen-bond donating +1.0, hydrogen bond accepting +1.0 [29]) were used to calculate five different fields: steric, electrostatic, hydrophobic, donor and acceptor. Gaussian-type distance dependence was used to measure the relative attenuation of the field position of each atom in the lattice, and led to much smoother sampling of the fields around the molecules when compared to CoMFA. The default value of 0.3 was set for the attenuation factor α.

3.3. Data Set Selection and Inhibitory Activity

CoMSIA studies were performed on a set of 90 piperazinyl urea derivatives reported in the literature [18,30,31,32,33,34,35,36] (Table 6). The derivatives displayed potent fatty acid amide hydrolase (FAAH) inhibitors activity. The IC50 values were converted to pIC50 (–logIC50). The compounds were divided into training (73 compounds, 81%) and test sets (17 compounds, 19%), ensuring that both sets contained structurally diverse compounds with high, medium and low activity, and a uniform distribution to avoid possible problems during the external validation. The distribution of pIC50 values for the whole set, the training set and the test set is shown in Figure 5. In all three cases the biological activity follows a Gaussian distribution.

3.4. Internal Validation and Partial Least Squares (PLS) Analysis

PLS analysis [37] was used to construct a linear correlation between the CoMFA and CoMSIA descriptors (independent variables) and the activity values (dependent variables) as previously described by Lorca et al. [24].

3.5. External Validation

The models were subjected to external validation criteria according to the proposed test by Golbraikh and Tropsha [38,39], which considers a QSAR model predictive, if the following conditions are satisfied:
q 2 > 0.5
r 2 > 0.6
( r 2 r 0 2 ) r 2 < 0.1   or   ( r 2 r   0 2 ) r 2 < 0.1
0.85 1.15   or   0.85 k 1.15
It has been demonstrated [38] that to adequately assess the predictive ability of a QSAR model the above criteria are necessary.
Furthermore, the external predictive power of the developed 3D-QSAR models using the test set was examined by considering r2m metrics as shown below [40]:
r m 2 = r 2   ( 1 | r 2 r 0 2 | )
where r2 and r20 are squared correlation coefficients between the observed and predicted activities of the test set with and without intercept, respectively. For a significant external model validation, the value of r2m should be greater than 0.5.
Additionally, the following descriptors were calculated:
Q F 1 2 = 1   i = 1 n E X T ( y i y ^ i ) 2 i = 1 n E X T ( y i y ¯ T R ) 2
Q F 2 2 = 1   i = 1 n E X T ( y i y ^ i ) 2 i = 1 n E X T ( y i y ¯ E X T ) 2
Q F 3 2 = 1   i = 1 n E X T ( y i y ^ i ) 2 / n E X T i = 1 n E X T ( y i y ¯ T R ) 2 / n T R
C C C = 2 i = 1 n E X T ( y i y ¯ ) ( y ^ i y ^ ¯ ) i = 1 n E X T ( y i y ¯ ) 2 +   i = 1 n E X T ( y ^ i   y ^ ¯ ) 2 +   n E X T ( y ¯ y ^ ¯ ) 2
Δ r m 2 = | r m 2 r m 2 |
where, TR = training set, EXT = external prediction set, y i = experimental data values, y ^ i = predicted data values, y ¯ = average of the experimental data values, y ^ ¯ = average of the predicted data values. Finally, r m 2 is calculated using the experimental values on the ordinate axis, while r m 2 using them on the abscissa.

3.6. Applicability Domain (AD) Calculation

The AD was evaluated based on the simple standardization method reported by Roy et al. [25] and as described by Lorca et al. [24].

4. Conclusions

In this contribution, a 3D-QSAR CoMSIA study was carried out on an extensive database of 90 irreversible inhibitors of the enzyme FAAH with a pyrimidinyl-piperazine-carboxamide general structure. The best model obtained considered all the field contributions, being the electrostatic and hydrogen-bond acceptor properties the ones that contributed most to the activity (30.4% and 33.0% respectively). The model was validated internally (q2 = 0.734) and externally (r2test = 0.966) and was also submitted to Tropsha validation criteria, r2m calculation (0.723) and Y-randomization test, passing all tests. The information derived from the contour maps was used to design a series of new compounds that showed promising predicted activities (pIC50 of the most active compounds = 12.196 and 12.416). The main structure-activity relationships found in this study and summarized in Figure 6 are a useful tool to guide the future design of new FAAH inhibitors. The extensive database used in this study could motivate future work complementing the information obtained from contour maps with QSAR studies by Dragon-based descriptors [41].

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/20/10/2510/s1.

Author Contributions

Data curation, J.R.-P.; Funding acquisition, C.D.P.-M.; Methodology, M.L.; Supervision, J.M.; Validation, Y.V. and H.C.; Writing—original draft, C.D.P.-M.; Writing—review & editing, J.M.

Funding

This research was funded by CONICYT FONDECYT, grant number 1150121 and 11130701.

Acknowledgments

Marcos Lorca thanks beca Conicyt-PFCHA/Doctorado Nacional/2018-21180427.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

3D-QSARThree-dimensional Quantitative Structure-Activity Relationship
CoMSIAComparative Molecular Similarity Indices Analysis
FAAHFatty Acid Amide Hydrolase
CBCannabinoid

References

  1. Pagotto, U.; Marsicano, G.; Cota, D.; Lutz, B.; Pasquali, R. The emerging role of the endocannabinoid system in endocrine regulation and energy balance. Endocr. Rev. 2006, 27, 73–100. [Google Scholar] [CrossRef] [PubMed]
  2. Manzanares, J.; Julian, M.D.; Carrascosa, A. Role of the cannabinoid system in pain control and therapeutic implications for the management of acute and chronic pain episodes. Curr. Neuropharmacol. 2006, 4, 239–257. [Google Scholar] [CrossRef] [PubMed]
  3. Petrosino, S.; Di Marzo, V. FAAH and MAGL inhibitors: Therapeutic opportunities from regulating endocannabinoid levels. Curr. Opin. Investig. Drugs 2010, 11, 51–62. [Google Scholar] [PubMed]
  4. Long, J.Z.; Cravatt, B.F. The metabolic serine hydrolases and their functions in mammalian physiology and disease. Chem. Rev. 2011, 111, 6022–6063. [Google Scholar] [CrossRef] [PubMed]
  5. Paylor, B.; Holt, S.; Fowler, C.J. The potency of the fatty acid amide hydrolase inhibitor URB597 is dependent upon the assay pH. Pharmacol. Res. 2006, 54, 481–485. [Google Scholar] [CrossRef] [PubMed]
  6. Mor, M.; Lodola, A.; Rivara, S.; Vacondio, F.; Duranti, A.; Tontini, A.; Sanchini, S.; Piersanti, G.; Clapper, J.R.; King, A.R.; et al. Synthesis and quantitative structure-activity relationship of fatty acid amide hydrolase inhibitors: Modulation at the N-portion of biphenyl-3-yl alkylcarbamates. J. Med. Chem. 2008, 51, 3487–3498. [Google Scholar] [CrossRef] [PubMed]
  7. Saario, S.M.; Poso, A.; Juvonen, R.O.; Jarvinen, T.; Salo-Ahen, O.M. Fatty acid amide hydrolase inhibitors from virtual screening of the endocannabinoid system. J. Med. Chem. 2006, 49, 4650–4656. [Google Scholar] [CrossRef] [PubMed]
  8. Cichero, E.; Ligresti, A.; Allara, M.; di Marzo, V.; Lazzati, Z.; D’Ursi, P.; Marabotti, A.; Milanesi, L.; Spallarossa, A.; Ranise, A.; et al. Homology modeling in tandem with 3D-QSAR analyses: A computational approach to depict the agonist binding site of the human CB2 receptor. Eur. J. Med. Chem. 2011, 46, 4489–4505. [Google Scholar] [CrossRef] [PubMed]
  9. Cichero, E.; Cesarini, S.; Mosti, L.; Fossa, P. CoMFA and CoMSIA analyses on 1,2,3,4-tetrahydropyrrolo[3,4-b]indole and benzimidazole derivatives as selective CB2 receptor agonists. J. Mol. Model. 2010, 16, 1481–1498. [Google Scholar] [CrossRef] [Green Version]
  10. Cichero, E.; Cesarini, S.; Mosti, L.; Fossa, P. CoMFA and CoMSIA analyses on 4-oxo-1,4-dihydroquinoline and 4-oxo-1,4-dihydro-1,5-, -1,6- and -1,8-naphthyridine derivatives as selective CB2 receptor agonists. J. Mol. Model. 2010, 16, 677–691. [Google Scholar] [CrossRef]
  11. Romero-Parra, J.; Chung, H.; Tapia, R.A.; Faundez, M.; Morales-Verdejo, C.; Lorca, M.; Lagos, C.F.; Di Marzo, V.; David Pessoa-Mahana, C.; Mella, J. Combined CoMFA and CoMSIA 3D-QSAR study of benzimidazole and benzothiophene derivatives with selective affinity for the CB2 cannabinoid receptor. Eur. J. Pharm. Sci. 2017, 101, 1–10. [Google Scholar] [CrossRef] [PubMed]
  12. Mella-Raipan, J.A.; Lagos, C.F.; Recabarren-Gajardo, G.; Espinosa-Bustos, C.; Romero-Parra, J.; Pessoa-Mahana, H.; Iturriaga-Vasquez, P.; Pessoa-Mahana, C.D. Design, synthesis, binding and docking-based 3D-QSAR studies of 2-pyridylbenzimidazoles-a new family of high affinity CB1 cannabinoid ligands. Molecules 2013, 18, 3972–4001. [Google Scholar] [CrossRef] [PubMed]
  13. Mella-Raipan, J.; Hernandez-Pino, S.; Morales-Verdejo, C.; Pessoa-Mahana, D. 3D-QSAR/CoMFA-based structure-affinity/selectivity relationships of aminoalkylindoles in the cannabinoid CB1 and CB2 receptors. Molecules 2014, 19, 2842–2861. [Google Scholar] [CrossRef] [PubMed]
  14. Dainese, E.; Gasperi, V.; Maccarrone, M. Partial QSAR analysis of some selected natural inhibitors of FAAH suggests a working hypothesis for the development of endocannabinoid-based drugs. Curr. Drug Targets CNS Neurol. Disord. 2005, 4, 709–714. [Google Scholar] [CrossRef]
  15. Kasnanen, H.; Myllymaki, M.J.; Minkkila, A.; Kataja, A.O.; Saario, S.M.; Nevalainen, T.; Koskinen, A.M.; Poso, A. 3-Heterocycle-phenyl N-alkylcarbamates as FAAH inhibitors: Design, synthesis and 3D-QSAR studies. ChemMedChem 2010, 5, 213–231. [Google Scholar] [CrossRef] [PubMed]
  16. Vacondio, F.; Silva, C.; Lodola, A.; Fioni, A.; Rivara, S.; Duranti, A.; Tontini, A.; Sanchini, S.; Clapper, J.R.; Piomelli, D.; et al. Structure-property relationships of a class of carbamate-based fatty acid amide hydrolase (FAAH) inhibitors: Chemical and biological stability. ChemMedChem 2009, 4, 1495–1504. [Google Scholar] [CrossRef]
  17. Han, D.; Wang, B.; Jin, H.; Wang, H.; Chen, M. Design, synthesis and CoMFA studies of OEA derivatives as FAAH inhibitors. Med. Chem. Res. 2017, 26, 2951–2966. [Google Scholar] [CrossRef]
  18. Keith, J.M.; Apodaca, R.; Tichenor, M.; Xiao, W.; Jones, W.; Pierce, J.; Seierstad, M.; Palmer, J.; Webb, M.; Karbarz, M.; et al. Aryl piperazinyl ureas as inhibitors of fatty acid amide hydrolase (FAAH) in rat, dog, and primate. ACS Med. Chem. Lett. 2012, 3, 823–827. [Google Scholar] [CrossRef]
  19. Durdagi, S.; Papadopoulos, M.G.; Papahatjis, D.P.; Mavromoustakos, T. Combined 3D QSAR and molecular docking studies to reveal novel cannabinoid ligands with optimum binding activity. Bioorg. Med. Chem. Lett. 2007, 17, 6754–6763. [Google Scholar] [CrossRef]
  20. Golbraikh, A.; Shen, M.; Xiao, Z.; Xiao, Y.D.; Lee, K.H.; Tropsha, A. Rational selection of training and test sets for the development of validated QSAR models. J. Comput. Aided Mol. Des. 2003, 17, 241–253. [Google Scholar] [CrossRef]
  21. Roy, K.; Mitra, I.; Kar, S.; Ojha, P.K.; Das, R.N.; Kabir, H. Comparative studies on some metrics for external validation of QSPR models. J. Chem. Inf. Model. 2012, 52, 396–408. [Google Scholar] [CrossRef] [PubMed]
  22. Chirico, N.; Gramatica, P. Real external predictivity of QSAR models. Part 2. New intercomparable thresholds for different validation criteria and the need for scatter plot inspection. J. Chem. Inf. Model. 2012, 52, 2044–2058. [Google Scholar] [CrossRef] [PubMed]
  23. Rucker, C.; Rucker, G.; Meringer, M. y-Randomization and its variants in QSPR/QSAR. J. Chem. Inf. Model. 2007, 47, 2345–2357. [Google Scholar] [CrossRef]
  24. Lorca, M.; Morales-Verdejo, C.; Vásquez-Velásquez, D.; Andrades-Lagos, A.; Campanini-Salinas, J.; Soto-Delgado, J.; Recabarren-Gajardo, G.; Mella, J. Structure-activity relationships based on 3D-QSAR CoMFA/CoMSIA and design of aryloxypropanol-amine agonists with selectivity for the human β3-adrenergic receptor and anti-obesity and anti-diabetic profiles. Molecules 2018, 23, 1191. [Google Scholar] [CrossRef] [PubMed]
  25. Roy, K.; Kar, S.; Ambure, P. On a simple approach for determining applicability domain of QSAR models. Chemometr. Intell. Lab. Syst. 2015, 145, 22–29. [Google Scholar] [CrossRef]
  26. Aronov, A.M. Ligand structural aspects of hERG channel blockade. Curr. Top. Med. Chem. 2008, 8, 1113–1127. [Google Scholar] [CrossRef] [PubMed]
  27. Certara ENHANCES SYBYL-X Drug Design and Discovery Software Suite. Available online: https://www.certara.com/pressreleases/certara-enhances-sybyl-x-drug-design-and-discovery-software-suite/ (accessed on 9 May 2019).
  28. Xue, C.X.; Cui, S.Y.; Liu, M.C.; Hu, Z.D.; Fan, B.T. 3D QSAR studies on antimalarial alkoxylated and hydroxylated chalcones by CoMFA and CoMSIA. Eur. J. Med. Chem. 2004, 39, 745–753. [Google Scholar] [CrossRef]
  29. Klebe, G.; Abraham, U.; Mietzner, T. Molecular similarity indices in a comparative analysis (CoMSIA) of drug molecules to correlate and predict their biological activity. J. Med. Chem. 1994, 37, 4130–4146. [Google Scholar] [CrossRef]
  30. Gustin, D.J.; Ma, Z.; Min, X.; Li, Y.; Hedberg, C.; Guimaraes, C.; Porter, A.C.; Lindstrom, M.; Lester-Zeiner, D.; Xu, G.; et al. Identification of potent, noncovalent fatty acid amide hydrolase (FAAH) inhibitors. Bioorg. Med. Chem. Lett. 2011, 21, 2492–2496. [Google Scholar] [CrossRef]
  31. Kono, M.; Matsumoto, T.; Kawamura, T.; Nishimura, A.; Kiyota, Y.; Oki, H.; Miyazaki, J.; Igaki, S.; Behnke, C.A.; Shimojo, M.; et al. Synthesis, SAR study, and biological evaluation of a series of piperazine ureas as fatty acid amide hydrolase (FAAH) inhibitors. Bioorg. Med. Chem. 2013, 21, 28–41. [Google Scholar] [CrossRef]
  32. Keith, J.M.; Apodaca, R.; Xiao, W.; Seierstad, M.; Pattabiraman, K.; Wu, J.; Webb, M.; Karbarz, M.J.; Brown, S.; Wilson, S.; et al. Thiadiazolopiperazinyl ureas as inhibitors of fatty acid amide hydrolase. Bioorg. Med. Chem. Lett. 2008, 18, 4838–4843. [Google Scholar] [CrossRef]
  33. Kono, M.; Matsumoto, T.; Imaeda, T.; Kawamura, T.; Fujimoto, S.; Kosugi, Y.; Odani, T.; Shimizu, Y.; Matsui, H.; Shimojo, M.; et al. Design, synthesis, and biological evaluation of a series of piperazine ureas as fatty acid amide hydrolase inhibitors. Bioorg. Med. Chem. 2014, 22, 1468–1478. [Google Scholar] [CrossRef]
  34. Tichenor, M.S.; Keith, J.M.; Jones, W.M.; Pierce, J.M.; Merit, J.; Hawryluk, N.; Seierstad, M.; Palmer, J.A.; Webb, M.; Karbarz, M.J.; et al. Heteroaryl urea inhibitors of fatty acid amide hydrolase: Structure-mutagenicity relationships for arylamine metabolites. Bioorg. Med. Chem. Lett. 2012, 22, 7357–7362. [Google Scholar] [CrossRef]
  35. Keith, J.M.; Tichenor, M.S.; Apodaca, R.L.; Xiao, W.; Jones, W.M.; Seierstad, M.; Pierce, J.M.; Palmer, J.A.; Webb, M.; Karbarz, M.J.; et al. The SAR of brain penetration for a series of heteroaryl urea FAAH inhibitors. Bioorg. Med. Chem. Lett. 2016, 26, 3109–3114. [Google Scholar] [CrossRef]
  36. Keith, J.M.; Jones, W.M.; Pierce, J.M.; Seierstad, M.; Palmer, J.A.; Webb, M.; Karbarz, M.J.; Scott, B.P.; Wilson, S.J.; Luo, L.; et al. Heteroarylureas with spirocyclic diamine cores as inhibitors of fatty acid amide hydrolase. Bioorg. Med. Chem. Lett. 2014, 24, 737–741. [Google Scholar] [CrossRef]
  37. Clark, M.; Cramer, R.D.; Van Opdenbosch, N. Validation of the general purpose Tripos 5.2 force field. J. Comput. Chem. 1989, 10, 982–1012. [Google Scholar] [CrossRef]
  38. Golbraikh, A.; Tropsha, A. Beware of q 2! J. Mol. Graph. Model. 2002, 20, 269–276. [Google Scholar] [CrossRef]
  39. Tropsha, A. Best practices for QSAR model development, validation, and exploitation. Mol. Inform. 2010, 29, 476–488. [Google Scholar] [CrossRef] [PubMed]
  40. Roy, K.; Chakraborty, P.; Mitra, I.; Ojha, P.K.; Kar, S.; Das, R.N. Some case studies on application of “r(m)2” metrics for judging quality of quantitative structure-activity relationship predictions: Emphasis on scaling of response data. J. Comput. Chem. 2013, 34, 1071–1082. [Google Scholar] [CrossRef]
  41. Mauri, A.; Consonni, V.; Pavan, M.; Todeschini, R. Dragon software: An easy approach to molecular descriptor calculations MATCH Commun. Math. Comput. Chem. 2006, 56, 237–248. [Google Scholar]
Figure 1. Examples of reported Fatty Acid Amide Hydrolase (FAAH) inhibitors
Figure 1. Examples of reported Fatty Acid Amide Hydrolase (FAAH) inhibitors
Ijms 20 02510 g001
Figure 2. (A) Plots of experimental versus predicted pIC50 values for the training and test set compounds. (B). Residual plots between predicted and experimental values. (C) CoMSIA predictions for all molecules in the set.
Figure 2. (A) Plots of experimental versus predicted pIC50 values for the training and test set compounds. (B). Residual plots between predicted and experimental values. (C) CoMSIA predictions for all molecules in the set.
Ijms 20 02510 g002
Figure 3. CoMSIA steric (A,B), electrostatic (C,D), hydrophobic (E,F), donor (G,H) and acceptor (I,J) contour maps around compounds 66 (left) and 6 (right), the most active and least active of the series respectively. Sterically favored areas are in green and disfavored areas are in yellow. Electropositive favoured areas are in blue and electronegative favoured areas are in red. Hydrophobic favoured areas are in yellow and unfavourable areas in grey. Donor and acceptor favoured areas are in cyan and magenta respectively, and donor and acceptor unfavourable areas are in purple and red, respectively.
Figure 3. CoMSIA steric (A,B), electrostatic (C,D), hydrophobic (E,F), donor (G,H) and acceptor (I,J) contour maps around compounds 66 (left) and 6 (right), the most active and least active of the series respectively. Sterically favored areas are in green and disfavored areas are in yellow. Electropositive favoured areas are in blue and electronegative favoured areas are in red. Hydrophobic favoured areas are in yellow and unfavourable areas in grey. Donor and acceptor favoured areas are in cyan and magenta respectively, and donor and acceptor unfavourable areas are in purple and red, respectively.
Ijms 20 02510 g003
Figure 4. The superimposed structures of all compounds used in the CoMSIA model. (A) Atom fit method. (B) Distill rigid method.
Figure 4. The superimposed structures of all compounds used in the CoMSIA model. (A) Atom fit method. (B) Distill rigid method.
Ijms 20 02510 g004
Figure 5. Histogram of frequency distribution data.
Figure 5. Histogram of frequency distribution data.
Ijms 20 02510 g005
Figure 6. Main structure-activity relationships derived from this study.
Figure 6. Main structure-activity relationships derived from this study.
Ijms 20 02510 g006
Table 1. Statistical parameters and field combinations for comparative molecular similarity indices analysis (CoMSIA).
Table 1. Statistical parameters and field combinations for comparative molecular similarity indices analysis (CoMSIA).
Modelq2NSEPSEEr2ncvFField Contributions
SEHDA
CoMSIA-S0.29331.0420.9030.47020.3831
CoMSIA-E0.534100.8930.4760.86740.592 1
CoMSIA-H0.31741.0320.7880.60225.722 1
CoMSIA-D0.25381.1121.0310.3594.482 1
CoMSIA-A0.52050.8710.6380.74338.785 1
CoMSIA-SE0.519100.9070.4170.89854.7990.3140.686
CoMSIA-SEH0.53480.8790.2160.977110.7010.1830.4920.324
CoMSIA-SEHD0.62870.7790.3820.91094.4480.1590.4580.2600.123
CoMSIA-SEHA0.68870.7130.3320.933128.4860.1310.3460.177 0.347
CoMSIA-SED0.62690.7930.3820.91373.6000.2450.639 0.116
CoMSIA-SEA0.72570.6700.3500.925114.6650.1820.387 0.421
CoMSIA-SEDA0.76570.6200.3270.934132.4750.1540.357 0.0990.389
CoMSIA-SH0.31641.0330.7720.61827.5020.408 0.592
CoMSIA-SD0.364191.1280.5250.86217.4330.814 0.186
CoMSIA-SA0.57270.8360.4840.85755.4530.344 0.656
CoMSIA-SHD0.42630.9390.7990.58532.3960.219 0.4790.303
CoMSIA-SHA0.52960.8700.4830.85564.9090.201 0.303 0.496
CoMSIA-SDA0.71970.6780.4040.90083.4850.235 0.2100.555
CoMSIA-SHDA0.67370.7310.3660.918103.7440.156 0.2400.1640.440
CoMSIA-EH0.550100.8770.3910.91163.307 0.5370.427
CoMSIA-ED0.61690.8040.4070.90264.163 0.856 0.144
CoMSIA-EA0.70160.6930.4080.89695.020 0.498 0.502
CoMSIA-EHD0.64180.7710.3760.91585.695 0.5250.3430.132
CoMSIA-EHA0.69170.7100.3900.925115.138 0.3900.234 0.376
CoMSIA-EDA0.75270.6360.3660.918103.539 0.453 0.1110.437
CoMSIA-EHDA0.74280.6540.3110.942128.936 0.3410.2110.0930.355
CoMSIA-HD0.42890.9810.5280.83435.186 0.8040.196
CoMSIA-HA0.53760.8620.4930.84961.685 0.426 0.574
CoMSIA-HDA0.682100.7380.3310.93690.828 0.3560.1550.490
CoMSIA-DA0.705110.7160.4650.87639.132 0.2400.760
CoMSIA-ALL0.73470.6590.3200.937138.3600.1100.3040.1560.1000.330
q2 = the square of the LOO cross-validation (CV) coefficient; N = the optimum number of components; SEP = standard error of prediction; SEE is the standard error of estimation of non CV analysis; r2ncv is the square of the non CV coefficient; F is the F-test value; S, E, H, D and A are the steric, electrostatic, hydrophobic, hydrogen-bond donor, and hydrogen-bond acceptor contributions respectively. No statistically significant models were found for CoMFA.
Table 2. Summary of external validation parameters for CoMSIA.
Table 2. Summary of external validation parameters for CoMSIA.
ConditionParametersThreshold ValueCoMSIA
1q2>0.50.734
2r2>0.60.966
3ar02Close to value of r20.920
3br02Close to value of r20.944
4ak0.85 < k < 1.151.004
4bk0.85 < k′ < 1.150.995
5a(r2−r20)/r2<0.10.048
5b(r2−r’20)/r2<0.10.023
06|r20−r′20|<0.30.024
7 r m 2 >0.50.723
8Q2F1>0.70.944
9Q2F2>0.70.943
10Q2F3>0.70.951
11CCC>0.850.967
12∆r2m<0.20.056
q2 = the square of the LOO cross-validation (CV) coefficient; r2 is the regression coefficient for the test set exclusively; r02 and k are the correlation coefficient between the actual and predicted activities for test set and the respective slope of regression; and r02 and k’ are the correlation coefficient between the predicted and actual activities for test set and the respective slope of regression. r2m was defined in equation 5. Parameters 8–12 are defined in the methods section.
Table 3. Y-randomization test for CoMSIA model.
Table 3. Y-randomization test for CoMSIA model.
Iterationq2r2ncvIterationq2r2ncv
Random 1−0.0130.107Random 60.0060.119
Random 2−0.0300.087Random 7−0.0930.183
Random 3−0.0520.082Random 80.0850.188
Random 4−0.1980.108Random 9−0.0340.086
Random 5−0.2020.179Random 10−0.1000.073
Table 4. Experimental and predicted pIC50 and residual values for the analyzed compounds obtained with the CoMSIA model.
Table 4. Experimental and predicted pIC50 and residual values for the analyzed compounds obtained with the CoMSIA model.
MolExp. pIC50Pred. pIC50ResidualMolExp. pIC50Pred. pIC50Residual
15.3315.620−0.2946 t8.0097.7740.23
2 t6.0766.0340.04476.5756.2860.29
36.8936.3170.58487.9877.6990.29
4 t5.6075.961−0.35497.8867.915−0.03
55.5586.107−0.55507.3017.2310.07
65.1765.295−0.12515.5746.018−0.44
75.3315.438−0.1152 t7.6387.6240.01
86.4566.0320.42535.9336.164−0.23
9 t5.8156.141−0.33545.9845.8760.11
106.3106.2910.02556.4956.3810.11
117.1316.8160.31567.1557.619−0.46
127.2087.689−0.48576.4956.4360.06
137.9217.8100.1158 t7.1557.268−0.11
14 t7.3377.386−0.05597.6387.5120.13
157.8247.4970.33607.4447.467−0.02
168.8548.7180.14617.5537.4080.14
177.9218.062−0.14626.8546.7410.11
188.3197.9720.35638.6588.820−0.16
198.4328.2060.23648.8248.7400.08
20 t6.9597.170−0.216510.14310.287−0.14
218.7708.7580.016610.60210.685−0.08
228.5388.559−0.0267 t10.1439.7590.38
23 t6.9217.421−0.506810.09710.144−0.05
248.6208.2730.35698.0618.0230.04
259.0369.126−0.09708.1148.1060.01
268.6788.4150.26717.4957.3970.10
279.1749.0570.12726.9876.7870.20
288.4958.4270.0773 t7.5537.707−0.15
298.9219.236−0.32747.5537.757−0.20
30 t9.0668.8140.25758.2377.3470.89
317.4826.7430.74765.8865.7390.15
326.8186.4900.33778.1748.1530.02
335.7316.422−0.6978 t8.4698.1550.31
346.7626.6540.11797.1617.463−0.30
355.8866.710−0.82808.0008.225−0.23
366.7146.734−0.02817.6588.237−0.58
37 t6.4606.848−0.39829.3019.394−0.09
386.5906.753−0.16837.3377.2160.12
397.4997.0720.4384 t7.4827.2780.20
407.0186.9010.12857.0007.134−0.13
41 t6.8456.878−0.03867.7708.110−0.34
425.9726.412−0.44876.8546.957−0.10
436.0796.471−0.39889.0008.9260.07
447.1617.1480.01898.8868.6860.20
45 t7.5747.855−0.28909.0008.9600.04
t test set compound.
Table 5. The proposed structures of new molecules and their predicted pIC50 values using the CoMSIA model.
Table 5. The proposed structures of new molecules and their predicted pIC50 values using the CoMSIA model.
StructurePred. pIC50StructurePred. pIC50
1x Ijms 20 02510 i00110.8896x Ijms 20 02510 i00211.744
2x Ijms 20 02510 i00311.3887x Ijms 20 02510 i00411.599
3x Ijms 20 02510 i00511.4908x Ijms 20 02510 i00611.822
4x Ijms 20 02510 i00710.9909x Ijms 20 02510 i00812.416
5x Ijms 20 02510 i00911.25310x Ijms 20 02510 i01012.196
Table 6. Chemical structure, IC50 (nM) and pIC50 values of the studied molecules.
Table 6. Chemical structure, IC50 (nM) and pIC50 values of the studied molecules.
Ijms 20 02510 i011
Ijms 20 02510 i012
Ijms 20 02510 i013
Ijms 20 02510 i014
Ijms 20 02510 i015
IC50 = Half maximal inhibitory concentration; pIC50 = −logIC50; M = mol·L−1

Share and Cite

MDPI and ACS Style

Lorca, M.; Valdes, Y.; Chung, H.; Romero-Parra, J.; Pessoa-Mahana, C.D.; Mella, J. Three-Dimensional Quantitative Structure-Activity Relationships (3D-QSAR) on a Series of Piperazine-Carboxamides Fatty Acid Amide Hydrolase (FAAH) Inhibitors as a Useful Tool for the Design of New Cannabinoid Ligands. Int. J. Mol. Sci. 2019, 20, 2510. https://doi.org/10.3390/ijms20102510

AMA Style

Lorca M, Valdes Y, Chung H, Romero-Parra J, Pessoa-Mahana CD, Mella J. Three-Dimensional Quantitative Structure-Activity Relationships (3D-QSAR) on a Series of Piperazine-Carboxamides Fatty Acid Amide Hydrolase (FAAH) Inhibitors as a Useful Tool for the Design of New Cannabinoid Ligands. International Journal of Molecular Sciences. 2019; 20(10):2510. https://doi.org/10.3390/ijms20102510

Chicago/Turabian Style

Lorca, Marcos, Yudisladys Valdes, Hery Chung, Javier Romero-Parra, C. David Pessoa-Mahana, and Jaime Mella. 2019. "Three-Dimensional Quantitative Structure-Activity Relationships (3D-QSAR) on a Series of Piperazine-Carboxamides Fatty Acid Amide Hydrolase (FAAH) Inhibitors as a Useful Tool for the Design of New Cannabinoid Ligands" International Journal of Molecular Sciences 20, no. 10: 2510. https://doi.org/10.3390/ijms20102510

APA Style

Lorca, M., Valdes, Y., Chung, H., Romero-Parra, J., Pessoa-Mahana, C. D., & Mella, J. (2019). Three-Dimensional Quantitative Structure-Activity Relationships (3D-QSAR) on a Series of Piperazine-Carboxamides Fatty Acid Amide Hydrolase (FAAH) Inhibitors as a Useful Tool for the Design of New Cannabinoid Ligands. International Journal of Molecular Sciences, 20(10), 2510. https://doi.org/10.3390/ijms20102510

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop