Computational methods applied to drug discovery: the rational design of dual inhibitors of FAAH and COX

The search for effective and safe drugs in pain-relief treatment represents a great challenge for medicinal chemists. Lipid derived mediators, such as endocannabinoids, may have different roles as agonists of cannabinoid receptors, relieving pain, or as substrates of cyclooxygenase (COX), generating the pro-inflammatory prostamides. Moreover, the tissue-protective endocannabinoid anandamide is metabolised by fatty acid amide hydrolase (FAAH). Therefore, a new challenging approach in pain-relief might be the development of dual action FAAH/COX inhibitors. The purpose of this thesis is to apply computational methods in drug discovery to assist medicinal chemistry studies targeting the rational design of novel FAAH/COX inhibitors, and to exploit structural studies relative to two side projects on other biological targets. The wider project of this thesis explores the mechanism of action and the rational design of novel FAAH/COX dual inhibitors. The reversible mixed type inhibitors Flu-AM1 and Ibu-AM5, derivatives of flurbiprofen and ibuprofen, respectively, retain similar COX inhibitory properties and are more potent FAAH inhibitors than the parent compounds. Applying a combination of molecular docking, MD simulations and free energy evaluation of the ligand-receptor complex, the binding mode of the enantiomer forms of Flu-AM1 and Ibu-AM5 has been found in the substrate access channel of FAAH and has been supported by studies of site-directed mutagenesis. The substitution of the isobutyl group of Ibu-AM5 with 4-(2-(trifluoromethyl)pyridin-4-yl)amino group led to the design of TPA5 derivative, which showed an inhibitory activity (IC50 = 0.59 μM) similar to the lead compound (Ibu-AM5, IC50 = 0.52 μM). Kinetic studies of TPA5 revealed that it is a pure competitive inhibitor of rat FAAH and molecular modeling studies supported a binding mode that overlap the anandamide analog MAFP. Among TPA5 derivatives, compound TPA27 exhibited a 10-fold enhancement in the inhibitory profile against FAAH (IC50 = 0.058 μM). Thermodynamic Integration calculations performed to complete the transformation of TPA5 in TPA27 yielded a free energy difference of 0.3 kcal/mol, which indicates a slight lower affinity of TPA27 with respect to TPA5, in the competitive binding site. Kinetics studies showed that TPA27 could be considered the first non-competitive reversible FAAH inhibitor reported so far, and that it more likely binds to an allosteric site. Differences in the inhibitory potency against rat and mouse FAAH for all compounds studied suggested different aminoacid composition of both competitive and non-competitive binding sites. This information was used as criteria of selection for a putative allosteric site found between the cytosolic port and the interface of the FAAH monomers. Computational studies in the allosteric site allowed the definition of the binding mode of Ibu-AM5 and TPA27. Nevertheless, a series of derivatives of Ibu-AM5 and Flu-AM1 were designed in order to get more information on the structure-activity relationships, leading to the identification of novel derivatives with improved activity against FAAH (Ibu-AM56, IC50 = 0.08 μM; Ibu-AM57, IC50 = 0.1 μM; Flu-AM3, IC50 = 0.02 μM. Finally, the thesis also reports the results of two other projects: i) the design of potential anticancer peptides that interfere in the formation of the tetrameric complex hUbA1/UbcH10/Ub2, key intermediate of the ubiquitination cascade.; ii) structural studies on the hybridization of PNA of different length with miR-509-3p, involved in regulating the expression of the CFTR gene, as a way to validate a potential new strategy for the treatment of Cystic Fibrosis

Design, synthesis and <i>in vitro</i> and <i>in vivo</i> biological evaluation of flurbiprofen amides as new fatty acid amide hydrolase/cyclooxygenase-2 dual inhibitory potential analgesic agents

Compounds combining dual inhibitory action against FAAH and cyclooxygenase (COX) may be potentially useful analgesics. Here, we describe a novel flurbiprofen analogue, N-(3-bromopyridin-2-yl)-2-(2-fluoro-(1,1'-biphenyl)-4-yl)propanamide (Flu-AM4). The compound is a competitive, reversible inhibitor of FAAH with a Ki value of 13 nM and which inhibits COX activity in a substrate-selective manner. Molecular modelling suggested that Flu-AM4 optimally fits a hydrophobic pocket in the ACB region of FAAH, and binds to COX-2 similarly to flurbiprofen. In vivo studies indicated that at a dose of 10 mg/kg, Flu-AM4 was active in models of prolonged (formalin) and neuropathic (chronic constriction injury) pain and reduced the spinal expression of iNOS, COX-2, and NFκB in the neuropathic model. Thus, the present study identifies Flu-AM4 as a dual-action FAAH/substrate-selective COX inhibitor with anti-inflammatory and analgesic activity in animal pain models. These findings underscore the potential usefulness of such dual-action compounds

Peptide Nucleic Acids as miRNA Target Protectors for the Treatment of Cystic Fibrosis

Cystic Fibrosis (CF) is one of the most common life shortening conditions in Caucasians. CF is caused by mutations in the CF Transmembrane Conductance Regulator (CFTR) gene which result in reduced or altered CFTR functionality. Several microRNAs (miRNAs) downregulate the expression of CFTR, thus causing or exacerbating the symptoms of CF. In this context, the design of anti-miRNA agents represents a valid functional tool, but its translation to the clinic might lead to unpredictable side effects because of the interference with the expression of other genes regulated by the same miRNAs. Herein, for the first time, is proposed the use of peptide nucleic acids (PNAs) to protect specific sequences in the 3'UTR (untranslated region) of the CFTR messenger RNA (mRNA) by action of miRNAs. Two PNAs (7 and 13 bases long) carrying the tetrapeptide Gly-SerP-SerP-Gly at their C-end, fully complementary to the 3'UTR sequence recognized by miR-509-3p, have been synthesized and the structural features of target PNA/RNA heteroduplexes have been investigated by spectroscopic and molecular dynamics studies. The co-transfection of the pLuc-CFTR-3´UTR vector with different combinations of PNAs, miR-509-3p, and controls in A549 cells demonstrated the ability of the longer PNA to rescue the luciferase activity by up to 70% of the control, thus supporting the use of suitable PNAs to counteract the reduction in the CFTR expression

Molecular Basis for Non-Covalent, Non-Competitive FAAH Inhibition

Fatty acid amide hydrolase (FAAH) plays a key role in the control of cannabinoid signaling and it represents a promising therapeutic strategy for the treatment of a wide range of diseases, including neuropathic pain and chronic inflammation. Starting from kinetics experiments carried out in our previous work for the most potent inhibitor 2-amino-3-chloropyridine amide (TPA14), we have investigated its non-competitive mechanism of action using molecular dynamics, thermodynamic integration and QM-MM/GBSA calculations. The computational studies highlighted the impact of mutations on the receptor binding pockets and elucidated the molecular basis of the non-competitive inhibition mechanism of TPA14, which prevents the endocannabinoid anandamide (AEA) from reaching its pro-active conformation. Our study provides a rationale for the design of non-competitive potent FAAH inhibitors for the treatment of neuropathic pain and chronic inflammation

Structural Model of the hUbA1-UbcH10 Quaternary Complex: In Silico and Experimental Analysis of the Protein-Protein Interactions between E1, E2 and Ubiquitin

UbcH10 is a component of the Ubiquitin Conjugation Enzymes (Ubc; E2) involved in the ubiquitination cascade controlling the cell cycle progression, whereby ubiquitin, activated by E1, is transferred through E2 to the target protein with the involvement of E3 enzymes. In this work we propose the first three dimensional model of the tetrameric complex formed by the human UbA1 (E1), two ubiquitin molecules and UbcH10 (E2), leading to the transthiolation reaction. The 3D model was built up by using an experimentally guided incremental docking strategy that combined homology modeling, protein-protein docking and refinement by means of molecular dynamics simulations. The structural features of the in silico model allowed us to identify the regions that mediate the recognition between the interacting proteins, revealing the active role of the ubiquitin crosslinked to E1 in the complex formation. Finally, the role of these regions involved in the E1–E2 binding was validated by designing short peptides that specifically interfere with the binding of UbcH10, thus supporting the reliability of the proposed model and representing valuable scaffolds for the design of peptidomimetic compounds that can bind selectively to Ubcs and inhibit the ubiquitylation process in pathological disorders

Structural model of the hUbA1-UbcH10 quaternary complex: In silico and experimental analysis of the protein-protein interactions between E1, E2 and ubiquitin

UbcH10 is a component of the Ubiquitin Conjugation Enzymes (Ubc; E2) involved in the ubiquitination cascade controlling the cell cycle progression, whereby ubiquitin, activated by E1, is transferred through E2 to the target protein with the involvement of E3 enzymes. In this work we propose the first three dimensional model of the tetrameric complex formed by the human UbA1 (E1), two ubiquitin molecules and UbcH10 (E2), leading to the transthiolation reaction. The 3D model was built up by using an experimentally guided incremental docking strategy that combined homology modeling, protein-protein docking and refinement by means of molecular dynamics simulations. The structural features of the in silico model allowed us to identify the regions that mediate the recognition between the interacting proteins, revealing the active role of the ubiquitin crosslinked to E1 in the complex formation. Finally, the role of these regions involved in the E1-E2 binding was validated by designing short peptides that specifically interfere with the binding of UbcH10, thus supporting the reliability of the proposed model and representing valuable scaffolds for the design of peptidomimetic compounds that can bind selectively to Ubcs and inhibit the ubiquitylation process in pathological disorders

Interaction of the N-(3-Methylpyridin-2-yl)amide Derivatives of Flurbiprofen and Ibuprofen with FAAH: Enantiomeric Selectivity and Binding Mode

Background Combined fatty acid amide hydrolase (FAAH) and cyclooxygenase (COX) inhibition is a promising approach for pain-relief. The Flu-AM1 and Ibu-AM5 derivatives of flurbiprofen and ibuprofen retain similar COX-inhibitory properties and are more potent inhibitors of FAAH than the parent compounds. However, little is known as to the nature of their interaction with FAAH, or to the importance of their chirality. This has been explored here. Methodology/Principal Findings FAAH inhibitory activity was measured in rat brain homogenates and in lysates expressing either wild-type or FAAHT488A-mutated enzyme. Molecular modelling was undertaken using both docking and molecular dynamics. The (R)- and (S)-enantiomers of Flu-AM1 inhibited rat FAAH with similar potencies (IC50 values of 0.74 and 0.99 μM, respectively), whereas the (S)-enantiomer of Ibu-AM5 (IC50 0.59 μM) was more potent than the (R)-enantiomer (IC50 5.7 μM). Multiple inhibition experiments indicated that both (R)-Flu-AM1 and (S)-Ibu-AM5 inhibited FAAH in a manner mutually exclusive to carprofen. Computational studies indicated that the binding site for the Flu-AM1 and Ibu-AM5 enantiomers was located between the acyl chain binding channel and the membrane access channel, in a site overlapping the carprofen binding site, and showed a binding mode in line with that proposed for carprofen and other non-covalent ligands. The potency of (R)-Flu-AM1 was lower towards lysates expressing FAAH mutated at the proposed carprofen binding area than in lysates expressing wild-type FAAH. Conclusions/Significance The study provides kinetic and structural evidence that the enantiomers of Flu-AM1 and Ibu-AM5 bind in the substrate channel of FAAH. This information will be useful in aiding the design of novel dual-action FAAH: COX inhibitors

Structural model of the hUbA1-UbcH10 quaternary complex: In silico and experimental analysis of the protein-protein interactions between E1, E2 and ubiquitin

UbcH10 is a component of the Ubiquitin Conjugation Enzymes (Ubc; E2) involved in the ubiquitination cascade controlling the cell cycle progression, whereby ubiquitin, activated by E1, is transferred through E2 to the target protein with the involvement of E3 enzymes. In this work we propose the first three dimensional model of the tetrameric complex formed by the human UbA1 (E1), two ubiquitin molecules and UbcH10 (E2), leading to the transthiolation reaction. The 3D model was built up by using an experimentally guided incremental docking strategy that combined homology modeling, protein-protein docking and refinement by means of molecular dynamics simulations. The structural features of the in silico model allowed us to identify the regions that mediate the recognition between the interacting proteins, revealing the active role of the ubiquitin crosslinked to E1 in the complex formation. Finally, the role of these regions involved in the E1-E2 binding was validated by designing short peptides that specifically interfere with the binding of UbcH10, thus supporting the reliability of the proposed model and representing valuable scaffolds for the design of peptidomimetic compounds that can bind selectively to Ubcs and inhibit the ubiquitylation process in pathological disorders

Structural model of the hUbA1-UbcH10 quaternary complex: In silico and experimental analysis of the protein-protein interactions between E1, E2 and ubiquitin

UbcH10 is a component of the Ubiquitin Conjugation Enzymes (Ubc; E2) involved in the ubiquitination cascade controlling the cell cycle progression, whereby ubiquitin, activated by E1, is transferred through E2 to the target protein with the involvement of E3 enzymes. In this work we propose the first three dimensional model of the tetrameric complex formed by the human UbA1 (E1), two ubiquitin molecules and UbcH10 (E2), leading to the transthiolation reaction. The 3D model was built up by using an experimentally guided incremental docking strategy that combined homology modeling, protein-protein docking and refinement by means of molecular dynamics simulations. The structural features of the in silico model allowed us to identify the regions that mediate the recognition between the interacting proteins, revealing the active role of the ubiquitin crosslinked to E1 in the complex formation. Finally, the role of these regions involved in the E1-E2 binding was validated by designing short peptides that specifically interfere with the binding of UbcH10, thus supporting the reliability of the proposed model and representing valuable scaffolds for the design of peptidomimetic compounds that can bind selectively to Ubcs and inhibit the ubiquitylation process in pathological disorders