2D/3D-QSAR Model Development Based on a Quinoline Pharmacophoric Core for the Inhibition of Plasmodium falciparum: An In Silico Approach with Experimental Validation.

Malaria is an infectious disease caused by Plasmodium spp. parasites, with widespread drug resistance to most antimalarial drugs. We report the development of two 3D-QSAR models based on comparative molecular field analysis (CoMFA), comparative molecular similarity index analysis (CoMSIA), and a 2D-QSAR model, using a database of 349 compounds with activity against the P. falciparum 3D7 strain. The models were validated internally and externally, complying with all metrics (q2 > 0.5, r2test > 0.6, r2m > 0.5, etc.). The final models have shown the following statistical values: r2test CoMFA = 0.878, r2test CoMSIA = 0.876, and r2test 2D-QSAR = 0.845. The models were experimentally tested through the synthesis and biological evaluation of ten quinoline derivatives against P. falciparum 3D7. The CoMSIA and 2D-QSAR models outperformed CoMFA in terms of better predictive capacity (MAE = 0.7006, 0.4849, and 1.2803, respectively). The physicochemical and pharmacokinetic properties of three selected quinoline derivatives were similar to chloroquine. Finally, the compounds showed low cytotoxicity (IC50 > 100 µM) on human HepG2 cells. These results suggest that the QSAR models accurately predict the toxicological profile, correlating well with experimental in vivo data.

Py-CoMFA, docking, and molecular dynamics simulations of Leishmania (L.) amazonensis arginase inhibitors.

Leishmaniasis is a disease caused by a protozoan of the genus Leishmania, affecting millions of people, mainly in tropical countries, due to poor social conditions and low economic development. First-line chemotherapeutic agents involve highly toxic pentavalent antimonials, while treatment failure is mainly due to the emergence of drug-resistant strains. Leishmania arginase (ARG) enzyme is vital in pathogenicity and contributes to a higher infection rate, thus representing a potential drug target. This study helps in designing ARG inhibitors for the treatment of leishmaniasis. Py-CoMFA (3D-QSAR) models were constructed using 34 inhibitors from different chemical classes against ARG from L. (L.) amazonensis (LaARG). The 3D-QSAR predictions showed an excellent correlation between experimental and calculated pIC50 values. The molecular docking study identified the favorable hydrophobicity contribution of phenyl and cyclohexyl groups as substituents in the enzyme allosteric site. Molecular dynamics simulations of selected protein-ligand complexes were conducted to understand derivatives' interaction modes and affinity in both active and allosteric sites. Two cinnamide compounds, 7g and 7k, were identified, with similar structures to the reference 4h allosteric site inhibitor. These compounds can guide the development of more effective arginase inhibitors as potential antileishmanial drugs.

Tetrahydrocurcumin Derivatives Enhanced the Anti-Inflammatory Activity of Curcumin: Synthesis, Biological Evaluation, and Structure-Activity Relationship Analysis.

Tetrahydrocurcumin, the most abundant curcumin transformation product in biological systems, can potentially be a new alternative therapeutic agent with improved anti-inflammatory activity and higher bioavailability than curcumin. In this article, we describe the synthesis and evaluation of the anti-inflammatory activities of tetrahydrocurcumin derivatives. Eleven tetrahydrocurcumin derivatives were synthesized via Steglich esterification on both sides of the phenolic rings of tetrahydrocurcumin with the aim of improving the anti-inflammatory activity of this compound. We showed that tetrahydrocurcumin (2) inhibited TNF-α and IL-6 production but not PGE2 production. Three tetrahydrocurcumin derivatives inhibited TNF-α production, five inhibited IL-6 production, and three inhibited PGE2 production. The structure-activity relationship analysis suggested that two factors could contribute to the biological activities of these compounds: the presence or absence of planarity and their structural differences. Among the tetrahydrocurcumin derivatives, cyclic compound 13 was the most active in terms of TNF-α production, showing even better activity than tetrahydrocurcumin. Acyclic compound 11 was the most effective in terms of IL-6 production and retained the same effect as tetrahydrocurcumin. Moreover, acyclic compound 12 was the most active in terms of PGE2 production, displaying better inhibition than tetrahydrocurcumin. A 3D-QSAR analysis suggested that the anti-inflammatory activities of tetrahydrocurcumin derivatives could be increased by adding bulky groups at the ends of compounds 2, 11, and 12.

QSAR Studies, Synthesis, and Biological Evaluation of New Pyrimido-Isoquinolin-Quinone Derivatives against Methicillin-Resistant Staphylococcus aureus.

According to the WHO, antimicrobial resistance is among the top 10 threats to global health. Due to increased resistance rates, an increase in the mortality and morbidity of patients has been observed, with projections of more than 10 million deaths associated with infections caused by antibacterial resistant microorganisms. Our research group has developed a new family of pyrimido-isoquinolin-quinones showing antibacterial activities against multidrug-resistant Staphylococcus aureus. We have developed 3D-QSAR CoMFA and CoMSIA studies (r2 = 0.938; 0.895), from which 13 new derivatives were designed and synthesized. The compounds were tested in antibacterial assays against methicillin-resistant Staphylococcus aureus and other bacterial pathogens. There were 12 synthesized compounds active against Gram-positive pathogens in concentrations ranging from 2 to 32 µg/mL. The antibacterial activity of the derivatives is explained by the steric, electronic, and hydrogen-bond acceptor properties of the compounds.

In silico approaches to develop new phenyl-pyrimidines as glycogen synthase kinase 3 (GSK-3) inhibitors with halogen-bonding capabilities: 3D-QSAR CoMFA/CoMSIA, molecular docking and molecular dynamics studies.

Glycogen synthase kinase 3 (GSK-3) is involved in different diseases, such as manic-depressive illness, Alzheimer's disease and cancer. Studies have shown that insulin inhibits GSK-3 to keep glycogen synthase active. Inhibiting GSK-3 may have an indirect pro-insulin effect by favouring glycogen synthesis. Therefore, the development of GSK-3 inhibitors can be a useful alternative for the treatment of type II diabetes. Aminopyrimidine derivatives already proved to be interesting GSK-3 inhibitors. In the current study, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) have been performed on a series of 122 aminopyrimidine derivatives in order to generate a robust model for the rational design of new compounds with promising antidiabetic activity. The q2 values obtained for the best CoMFA and CoMSIA models have been 0.563 and 0.598, respectively. In addition, the r2 values have been 0.823 and 0.925 for CoMFA and CoMSIA, respectively. The models were statistically validated, and from the contour maps analysis, a proposal of 10 new compounds has been generated, with predicted pIC50 higher than 9. The final contribution of our work is that: (a) we provide an extensive structure-activity relationship for GSK-3 inhibitory pyrimidines; and (b) these models may speed up the discovery of GSK-3 inhibitors based on the aminopyrimidine scaffold. Finally, we carried out docking and molecular dynamics studies of the two best candidates, which were shown to establish halogen-bond interactions with the enzyme.Communicated by Ramaswamy H. Sarma.

Synthesis and Phytotoxic Evaluation of Isatin Derivatives Supported by 3D-QSAR Study.

Concerned about weed infestation, a major threat to food production and herbicide resistance that interferes in the mechanism of action of the main herbicides, we have synthesized eight isatin derivatives using the "Click Chemistry" approach through copper-catalyzed azide-alkyne cycloadditions (CuAAC). Sixteen isatin derivatives were evaluated for their phytotoxic activity against the seed culture of the model plants, Lactuca sativa and Allium cepa. Six of them showed phytotoxic activity similar to the positive control, trifluralin. Hypocotyl length measurement analysis in L. sativa revealed that triazole derivative 8 is more active than trifluralin. For A. cepa, root length measurement analyses revealed that 3, 10, 14, 16, and 17 were similar to the positive control trifluralin. Three-dimensional quantitative structure-activity relationship (3D-QSAR) comparative molecular field analysis (CoMFA) model construction using the acetolactate synthase (ALS) crystallographic structure displayed pki values of predicted inhibitory activity and contour maps revealing sterically bulky groups for 11, the CF3 group in ortho, and for 17, Br in ortho, favoring the inhibitory ALS activity.

An in silico molecular modeling approach of halolactone derivatives as potential inhibitors for human immunodeficiency virus type-1 reverse transcriptase enzyme.

Acquired Immune Deficiency Syndrome (AIDS) is an infectious disease caused by Human Immunodeficiency Virus (HIV) infection and its replication requires the Reverse Transcriptase (RT) enzyme. RT plays a key role in the HIV life cycle, making it one of the most important targets for designing new drugs. Thus, in order to increase therapeutic options against AIDS, halolactone derivatives (D-halolactone) that have been showed as potential non-nucleoside inhibitors of the RT enzyme were studied. In the present work, a series of D-halolactone were investigated by molecular modeling studies, combining Three-dimensional Quantitative Structure-Activity Relationship (3 D-QSAR), molecular docking and Molecular Dynamics (MD) techniques, to understand the molecular characteristics that promote biological activity. The internal and external validation parameters indicated that the 3 D-QSAR model has good predictive capacity and statistical significance. Contour maps provided useful information on the structural characteristics of compounds for anti-HIV-1 activity. The docking results showed that D-halolactone present good complementarity by the RT allosteric site. In MD simulations it was observed that the formation of enzyme-ligand complexes were favorable, and from the free energy decomposition it was found that Leu100, Val106, Tyr181, Try188, and Trp229 are key residues for stabilization in the enzymatic site. Thus, the results showed that the proposed models can be used to design promising HIV-1 RT inhibitors. Communicated by Ramaswamy H. Sarma.

3D-QSARpy: Combining variable selection strategies and machine learning techniques to build QSAR models

Abstract Quantitative Structure-Activity Relationship (QSAR) is a computer-aided technology in the field of medicinal chemistry that seeks to clarify the relationships between molecular structures and their biological activities. Such technologies allow for the acceleration of the development of new compounds by reducing the costs of drug design. This work presents 3D-QSARpy, a flexible, user-friendly and robust tool, freely available without registration, to support the generation of QSAR 3D models in an automated way. The user only needs to provide aligned molecular structures and the respective dependent variable. The current version was developed using Python with packages such as scikit-learn and includes various techniques of machine learning for regression. The diverse techniques employed by the tool is a differential compared to known methodologies, such as CoMFA and CoMSIA, because it expands the search space of possible solutions, and in this way increases the chances of obtaining relevant models. Additionally, approaches for select variables (dimension reduction) were implemented in the tool. To evaluate its potentials, experiments were carried out to compare results obtained from the proposed 3D-QSARpy tool with the results from already published works. The results demonstrated that 3D-QSARpy is extremely useful in the field due to its expressive results.

Development of 3D-QSAR and pharmacophoric models to design new anti-Trypanosoma cruzi agents based on 2-aryloxynaphthoquinone scaffold.

In this work we have collected a set of 30 trypanosomicidal naphthoquinones and developed pharmacophoric and 3D-QSAR models as tools for the design of new potential anti-Chagasic compounds. Firstly, qualitative information was obtained from SAR and pharmacophoric models identifying some fragments around the 2-aryloxynaphthoquinone scaffold important for the antiparasitic activity. Then, 3D-QSAR CoMFA and CoMSIA models were developed. The models showed adequate statistical parameters where the steric, electrostatic, and hydrophobic features explain the trypanosomicidal effect. Therefore, to validate our models, we carried out the design, synthesis, and biological evaluation on T. cruzi epimastigotes of five new compounds (33a-e). According to CoMFA model, three out of five compounds showed pIC50 values within one logarithmic unit of deviation. The two compounds that did not fit the predictions were those with high lipophilicity, which agreed with the SAR and pharmacophore models. Docking and molecular dynamic studies were performed on T. cruzi trypanothione reductase, in a proposed binding site for this type of naphthoquinone. Interestingly, 33a-e showed the same interaction pattern as a naphthoquinone inhibitor (2). Finally, predicted drug-likeness properties indicated that 33a-e have optimal oral bioavailability. Thus, this study provides new in silico models for obtaining novel trypanosomicidal compounds.

Design, Synthesis, In Silico Studies and Inhibitory Activity towards Bcr-Abl, BTK and FLT3-ITD of New 2,6,9-Trisubstituted Purine Derivatives as Potential Agents for the Treatment of Leukaemia.

We report 31 new compounds designed, synthesized and evaluated on Bcr-Abl, BTK and FLT3-ITD as part of our program to develop 2,6,9-trisubstituted purine derivatives as inhibitors of oncogenic kinases. The design was inspired by the chemical structures of well-known kinase inhibitors and our previously developed purine derivatives. The synthesis of these purines was simple and used a microwave reactor for the final step. Kinase assays showed three inhibitors with high selectivity for each protein that were identified: 4f (IC50 = 70 nM for Bcr-Abl), 5j (IC50 = 0.41 µM for BTK) and 5b (IC50 = 0.38 µM for FLT-ITD). The 3D-QSAR analysis and molecular docking studies suggested that two fragments are potent and selective inhibitors of these three kinases: a substitution at the 6-phenylamino ring and the length and volume of the alkyl group at N-9. The N-7 and the N-methyl-piperazine moiety linked to the aminophenyl ring at C-2 are also requirements for obtaining the activity. Furthermore, most of these purine derivatives were shown to have a significant inhibitory effect in vitro on the proliferation of leukaemia and lymphoma cells (HL60, MV4-11, CEM, K562 and Ramos) at low concentrations. Finally, we show that the selected purines (4i, 5b and 5j) inhibit the downstream signalling of the respective kinases in cell models. Thus, this study provides new evidence regarding how certain chemical modifications of purine ring substituents provide novel inhibitors of target kinases as potential anti-leukaemia drugs.

QSAR-driven synthesis of antiproliferative chalcones against SH-SY5Y cancer cells: Design, biological evaluation, and redesign.

Neuroblastoma is one of the most frequent types of cancer found in infants, and traditional chemotherapy has limited efficacy against this pathology. Thus, the development of new compounds with higher activity and selectivity than traditional drugs is a current challenge in medicinal chemistry research. In this study, we report the synthesis of 21 chalcones with antiproliferative activity and selectivity against the neuroblastoma cell line SH-SY5Y. Then, we developed three-dimensional quantitative structure-activity relationship models (comparative molecular field analysis and comparative molecular similarity index analysis) with high-quality statistical values (q2 > 0.7; r2 > 0.8; r2 pred > 0.7), using IC50 and selectivity index (SI) data as dependent variables. With the information derived from these theoretical models, we designed and synthesized 16 new molecules to prove their consistency, finding good antiproliferative activity against SH-SY5Y cells on these derivatives, with three of them showing higher SI than the referential drugs 5-fluorouracil and cisplatin, displaying also a proapoptotic effect comparable to these drugs, as proven by measuring their effects on executor caspases 3/7 activity induction, Bcl-2/Bax messenger RNA levels alteration, and DNA fragmentation promotion.

Is conformation relevant for QSAR purposes? 2D Chemical representation in a 3D-QSAR perspective.

Conformation has a key role in the mechanism of interaction between small molecules and biological receptors. However, encoding this type of information in molecular descriptors for the construction of robust quantitative structure-activity relationships (QSAR) models is not an easy task and, so far, the dependence of these models on such feature has not been thoroughly investigated. In the present study, the authors explore the effects of conformational information on a 3D-QSAR technique by comparing models built with descriptors that encode fully described tridimensional aspects (structures docked inside a biological target), with descriptors in which this information is suppressed (flat structures) or not fully described (structures with quantum-chemically optimized geometries). As a result, the validation parameters indicate that the robustness of the models seems to be more related to the alignment aspect of the structures than to how well their tridimensional features are described.

Coumarin-Resveratrol-Inspired Hybrids as Monoamine Oxidase B Inhibitors: 3-Phenylcoumarin versus trans-6-Styrylcoumarin.

Monoamine oxidases (MAOs) are attractive targets in drug design. The inhibition of one of the isoforms (A or B) is responsible for modulating the levels of different neurotransmitters in the central nervous system, as well as the production of reactive oxygen species. Molecules that act selectively on one of the MAO isoforms have been studied deeply, and coumarin has been described as a promising scaffold. In the current manuscript we describe a comparative study between 3-phenylcoumarin (endo coumarin-resveratrol-inspired hybrid) and trans-6-styrylcoumarin (exo coumarin-resveratrol-inspired hybrid). Crystallographic structures of both compounds were obtained and analyzed. 3D-QSAR models, in particular CoMFA and CoMSIA, docking simulations and molecular dynamics simulations have been performed to support and better understand the interaction of these molecules with both MAO isoforms. Both molecules proved to inhibit MAO-B, with trans-6-styrylcoumarin being 107 times more active than 3-phenylcoumarin, and 267 times more active than trans-resveratrol.

Understanding the Molecular Basis of 5-HT4 Receptor Partial Agonists through 3D-QSAR Studies.

Alzheimer's disease (AD) is a neurodegenerative disorder whose prevalence has an incidence in senior citizens. Unfortunately, current pharmacotherapy only offers symptom relief for patients with side effects such as bradycardia, nausea, and vomiting. Therefore, there is a present need to provide other therapeutic alternatives for treatments for these disorders. The 5-HT4 receptor is an attractive therapeutic target since it has a potential role in central and peripheral nervous system disorders such as AD, irritable bowel syndrome, and gastroparesis. Quantitative structure-activity relationship analysis of a series of 62 active compounds in the 5-HT4 receptor was carried out in the present work. The structure-activity relationship was estimated using three-dimensional quantitative structure-activity relationship (3D-QSAR) techniques based on these structures' field molecular (force and Gaussian field). The best force-field QSAR models achieve a value for the coefficient of determination of the training set of R2training = 0.821, and for the test set R2test = 0.667, while for Gaussian-field QSAR the training and the test were R2training = 0.898 and R2test = 0.695, respectively. The obtained results were validated using a coefficient of correlation of the leave-one-out cross-validation of Q2LOO = 0.804 and Q2LOO = 0.886 for force- and Gaussian-field QSAR, respectively. Based on these results, novel 5-HT4 partial agonists with potential biological activity (pEC50 8.209-9.417 for force-field QSAR and 9.111-9.856 for Gaussian-field QSAR) were designed. In addition, for the new analogues, their absorption, distribution, metabolism, excretion, and toxicity properties were also analyzed. The results show that these new derivatives also have reasonable pharmacokinetics and drug-like properties. Our findings suggest novel routes for the design and development of new 5-HT4 partial agonists.

Behavior of Chemokine Receptor 6 (CXCR6) in Complex with CXCL16 Soluble form Chemokine by Molecular Dynamic Simulations: General Protein‒Ligand Interaction Model and 3D-QSAR Studies of Synthetic Antagonists.

The CXCR6‒CXCL16 axis is involved in several pathological processes, and its overexpression has been detected in different types of cancer, such as prostate, breast, ovary, and lung cancer, along with schwannomas, in which it promotes invasion and metastasis. Moreover, this axis is involved in atherosclerosis, type 1 diabetes, primary immune thrombocytopenia, vitiligo, and other autoimmune diseases, in which it is responsible for the infiltration of different immune system cells. The 3D structure of CXCR6 and CXCL16 has not been experimentally resolved; therefore, homology modeling and molecular dynamics simulations could be useful for the study of this signaling axis. In this work, a homology model of CXCR6 and a soluble form of CXCL16 (CXCR6‒CXCL16s) are reported to study the interactions between CXCR6 and CXCL16s through coarse-grained molecular dynamics (CG-MD) simulations. CG-MD simulations showed the two activation steps of CXCR6 through a decrease in the distance between the chemokine and the transmembrane region (TM) of CXCR6 and transmembrane rotational changes and polar interactions between transmembrane segments. The polar interactions between TM3, TM5, and TM6 are fundamental to functional conformation and the meta-active state of CXCR6. The interactions between D77-R280 and T243-TM7 could be related to the functional conformation of CXCR6; alternatively, the interaction between Q195-Q244 and N248 could be related to an inactive state due to the loss of this interaction, and an arginine cage broken in the presence of CXCL16s allows the meta-active state of CXCR6. A general protein‒ligand interaction supports the relevance of TM3‒TM5‒TM6 interactions, presenting three relevant pharmacophoric features: HAc (H-bond acceptor), HDn (H-bond donator), and Hph (hydrophobic), distributed around the space between extracellular loops (ECLs) and TMs. The HDn feature is close to TM3 and TM6; likewise, the HAc and Hph features are close to ECL1 and ECL2 and could block the rotation and interactions between TM3‒TM6 and the interactions of CXCL16s with the ECLs. Tridimensional quantitative structure-activity relationships (3D-QSAR) models show that the positive steric (VdW) and electrostatic fields coincide with the steric and positive electrostatic region of the exo-azabicyclo[3.3.1]nonane scaffold in the best pIC50 ligands. This substructure is close to the E274 residue and therefore relevant to the activity of CXCR6. These data could help with the design of new molecules that inhibit chemokine binding or antagonize the receptor based on the activation mechanism of CXCR6 and provoke a decrease in chemotaxis caused by the CXCR6‒CXCL16 axis.

Identification of the naphthoquinone derivative inhibitors binding site in heat shock protein 90: an induced-fit docking, molecular dynamics and 3D-QSAR study.

The combination of molecular modeling methods to identify the putative binding site of inhibitors constitutes an important tool in drug discovery. In this work, we used these analyses to understand the potent inhibitory effect of naphthoquinone derivatives on heat shock protein 90 (Hsp90), one of the proteins involved in many types of cancer. Molecular docking results indicated that some favorable interactions of key amino acid residues at the binding site of Hsp90 with these quinones would be responsible for the inhibition of Hsp90 activity. Molecular docking and molecular dynamics simulation were carried out to further understand the binding modes and the interactions between the protein and these inhibitors. The main residues of the internal cavity were Val136, Phe138, Tyr139, Val150, Trp162 and Val186. The high concordance between the docking results and 3D-QSAR contour maps gives us helpful information about the environment of the binding site. Our results provide the bases for a rational modification of new molecules based in quinone scaffold, in order to design more potent Hsp90 inhibitors, which would exhibit highly potent antitumor activity.Communicated by Ramaswamy H. Sarma.

3-Arylcoumarins as highly potent and selective monoamine oxidase B inhibitors: Which chemical features matter?

Monoamine oxidase B inhibitory activity is closely regulated by the interaction of the small molecules with the enzyme. It is therefore desirable to use theoretical approaches to design rational methods to develop new molecules to modulate specific interactions with the protein. Here, we report such methods, and we illustrate their successful implementation by studying six synthetized 3-arylcoumarins (71-76) based on them. Monoamine oxidase B inhibition is essential to maintain the balance of dopamine, and one of its major functions is to combat dopamine degradation, a phenomenon linked to Parkinson's disease. In this work, we study small-molecule inhibitors based on the 3-arylcoumarin scaffold and their monoamine oxidase B selective inhibition. We show that 3D-QSAR models, in particular CoMFA and CoMSIA, and molecular docking approaches, enhance the probability to find new interesting inhibitors, avoiding very costly and time-consuming synthesis and biological evaluations.

Promising 2,6,9-Trisubstituted Purine Derivatives for Anticancer Compounds: Synthesis, 3D-QSAR, and Preliminary Biological Assays.

We designed, synthesized, and evaluated novel 2,6,9-trisubstituted purine derivatives for their prospective role as antitumor compounds. Using simple and efficient methodologies, 31 compounds were obtained. We tested these compounds in vitro to draw conclusions about their cell toxicity on seven cancer cells lines and one non-neoplastic cell line. Structural requirements for antitumor activity on two different cancer cell lines were analyzed with SAR and 3D-QSAR. The 3D-QSAR models showed that steric properties could better explain the cytotoxicity of compounds than electronic properties (70% and 30% of contribution, respectively). From this analysis, we concluded that an arylpiperazinyl system connected at position 6 of the purine ring is beneficial for cytotoxic activity, while the use of bulky systems at position C-2 of the purine is not favorable. Compound 7h was found to be an effective potential agent when compared with a currently marketed drug, cisplatin, in four out of the seven cancer cell lines tested. Compound 7h showed the highest potency, unprecedented selectivity, and complied with all the Lipinski rules. Finally, it was demonstrated that 7h induced apoptosis and caused cell cycle arrest at the S-phase on HL-60 cells. Our study suggests that substitution in the purine core by arylpiperidine moiety is essential to obtain derivatives with potential anticancer activity.

Structural Requirements of N-alpha-Mercaptoacetyl Dipeptide (NAMdP) Inhibitors of Pseudomonas Aeruginosa Virulence Factor LasB: 3D-QSAR, Molecular Docking, and Interaction Fingerprint Studies.

The zinc metallopeptidase Pseudomonas elastase (LasB) is a virulence factor of Pseudomonas aeruginosa (P. aeruginosa), a pathogenic bacterium that can cause nosocomial infections. The present study relates the structural analysis of 118 N-alpha-mercaptoacetyl dipeptides (NAMdPs) as LasB inhibitors. Field-based 3D-QSAR and molecular docking methods were employed to describe the essential interactions between NAMdPs and LasB binding sites, and the chemical features that determine their differential activities. We report a predictive 3D-QSAR model that was developed according to the internal and external validation tests. The best model, including steric, electrostatic, hydrogen bond donor, hydrogen bond acceptor, and hydrophobic fields, was found to depict a three-dimensional map with the local positive and negative effects of these chemotypes on the LasB inhibitory activities. Furthermore, molecular docking experiments yielded bioactive conformations of NAMdPs inside the LasB binding site. The series of NAMdPs adopted a similar orientation with respect to phosphoramidon within the LasB binding site (crystallographic reference), where the backbone atoms of NAMdPs are hydrogen-bonded to the LasB residues N112, A113, and R198, similarly to phosphoramidon. Our study also included a deep description of the residues involved in the protein-ligand interaction patterns for the whole set of NAMdPs, through the use of interaction fingerprints (IFPs).

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.

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).