Insight into binding of phosphodiesterase-9A selective inhibitors by crystal structures and mutagenesis

Identified Isosteric Replacements of Ligands' Glycosyl Domain by Data Mining

Biologically equivalent replacements of key moieties in molecules rationalize scaffold hopping, patent busting, or R-group enumeration. Yet, this information may depend upon the expert-defined space, and might be subjective and biased toward the chemistries they get used to. Most importantly, these practices are often informatively incomplete since they are often compromised by a try-and-error cycle, and although they depict what kind of substructures are suitable for the replacement occurrence, they fail to explain the driving forces to support such interchanges. The protein data bank (PDB) encodes a receptor-ligand interaction pattern and could be an optional source to mine structural surrogates. However, manual decoding of PDB has become almost impossible and redundant to excavate the bioisosteric know-how. Therefore, a text parsing workflow has been developed to automatically extract the local structural replacement of a specific structure from PDB by finding spatial and steric interaction overlaps between the fragments in endogenous ligands and particular ligand fragments. Taking the glycosyl domain for instance, a total of 49 520 replacements that overlap on nucleotide ribose were identified and categorized based on their SMILE codes. A predominately ring system, such as aliphatic and aromatic rings, was observed; yet, amide and sulfonamide replacements also occur. We believe these findings may enlighten medicinal chemists on the structure design and optimization of ligands using the bioisosteric replacement strategy.

8-Alkylmercaptocaffeine derivatives: antioxidant, molecular docking, and in-vitro cytotoxicity studies

Due to their unique pharmacological characteristics, methylxanthines are known as therapeutic agents in a fascinating range of medicinal scopes. In this report, we aimed to examine some biological effects of previously synthesized 8-alkylmercaptocaffeine derivatives. Cytotoxic and antioxidative activity of 8-alkylmercaptocaffeine derivatives were measured in malignant A549, MCF7, and C152 cell lines. Assessment of cGMP levels and caspase-3 activity were carried out using a colorimetric competitive ELISA kit. Computational approaches were employed to discover the inhibitory mechanism of synthesized compounds. Among the twelve synthesized derivatives, three compounds (C1, C5, and C7) bearing propyl, heptyl, and 3-methyl-butyl moieties showed higher and more desirable cytotoxic activity against all the studied cell lines (IC₅₀ < 100 µM). Furthermore, C5 synergistically enhanced cisplatin-induced cytotoxicity in MCF-7 cells (CI < 1). Both C5 and C7 significantly increased caspase-3 activity and intracellular cGMP levels at specific time intervals in all studied cell lines (P < 0.05). However, these derivatives did not elevate LDH leakage (P > 0.05) and exhibited no marked ameliorating effects on oxidative damage (P > 0.05). Computational studies showed that H-bond formation between the nitrogen atom in pyrazolo[4,3-D] pyrimidine moiety with Gln817 and creating a hydrophobic cavity result in the stability of the alkyl group in the PDE5A active site. We found that synthesized 8-alkylmercaptocaffeine derivatives induced cell death in different cancer cells through the cGMP pathway. These findings will help us to get a deeper insight into the role of methylxanthines as useful alternatives to conventional cancer therapeutics.

Shifting the paradigm in treating multi-factorial diseases: polypharmacological co-inhibitors of HDAC6

Multi-factorial diseases are illnesses that exploit multiple cellular processes, or stages within one process, and thus highly targeted therapies often succumb to the disease, losing efficacy as resistance sets in. Combination therapies have become a mainstay to battle these diseases, however these regimens are plagued with caveats. An emerging avenue to treat multi-factorial diseases is polypharmacology, wherein a single drug is rationally designed to bind multiple targets, and is widely touted to be superior to combination therapy by inherently addressing the latter's shortcomings, which include poor patient compliance, narrow therapeutic windows and spiraling healthcare costs. Through its roles in intracellular trafficking, cell motility, mitosis, protein folding and as a back-up to the proteasome pathway, HDAC6 has rapidly become an exciting new target for therapeutics, particularly in the discovery of new drugs to treat Alzheimer's disease and cancer. Herein, we describe recent efforts to marry together HDAC pharmacophores, with a particular emphasis on HDAC6 selectivity, with those of other targets towards the discovery of potent therapeutics to treat these evasive diseases. Such polypharmacological agents may supercede combination therapies through inherent synergism, permitting reduced dosing, wider therapeutic windows and improved compliance.

The long and winding road of designing phosphodiesterase inhibitors for the treatment of heart failure

Cyclic nucleotide phosphodiesterases (PDEs) are a superfamily of enzymes known to play a critical role in the indirect regulation of several intracellular metabolism pathways through the selective hydrolysis of the phosphodiester bonds of specific second messenger substrates such as cAMP (3',5'-cyclic adenosine monophosphate) and cGMP (3',5'-cyclic guanosine monophosphate), influencing the hypertrophy, contractility, apoptosis and fibroses in the cardiovascular system. The expression and/or activity of multiple PDEs is altered during heart failure (HF), which leads to changes in levels of cyclic nucleotides and function of cardiac muscle. Within the cardiovascular system, PDEs 1-5, 8 and 9 are expressed and are interesting targets for the HF treatment. In this comprehensive review we will present a briefly description of the biochemical importance of each cardiovascular related PDE to the HF, and cover almost all the "long and winding road" of designing and discovering ligands, hits, lead compounds, clinical candidates and drugs as PDE inhibitors in the last decade.

Multitarget Approach for the Treatment of Alzheimer's Disease: Inhibition of Phosphodiesterase 9 (PDE9) and Histone Deacetylases (HDACs) Covering Diverse Selectivity Profiles

Here, we present a series of dual-target phosphodiesterase 9 (PDE9) and histone deacetylase (HDAC) inhibitors devised as pharmacological tool compounds for assessing the implications of these two targets in Alzheimer's disease (AD). These novel inhibitors were designed taking into account the key pharmacophoric features of known selective PDE9 inhibitors as well as privileged chemical structures, bearing zinc binding groups (hydroxamic acids and ortho-amino anilides) that hit HDAC targets. These substituents were selected according to rational criteria and previous knowledge from our group to explore diverse HDAC selectivity profiles (pan-HDAC, HDAC6 selective, and class I selective) that were confirmed in biochemical screens. Their functional response in inducing acetylation of histone and tubulin and phosphorylation of cAMP response element binding (CREB) was measured as a requisite for further progression into complete in vitro absorption, distribution, metabolism and excretion (ADME) and in vivo brain penetration profiling. Compound 31b, a selective HDAC6 inhibitor with acceptable brain permeability, was chosen for assessing in vivo efficacy of these first-in-class inhibitors, as well as studying their mode of action (MoA).

Phosphodiesterase 9 (PDE9) regulates bovine tracheal smooth muscle relaxation

The present study was designed to clarify phosphodiesterase 9 (PDE9) expression in bovine tracheal smooth muscle tissue, and to elucidate that PDE9 may contribute to the regulation of airway relaxation. PDE9 mRNA expression was detected in bovine tracheal smooth muscle. Sodium nitroprusside (an NO donor) and BAY 73-6691 (a selective PDE9 inhibitor) reduced high K⁺- and carbachol-induced contraction. BAY 73-6691 relaxed tracheal tissue on the same level with vardenafil (a selective PDE5 inhibitor). These results support our hypothesis that PDE9 plays functional role in the tracheal smooth muscle relaxation. PDE9 inhibitors are expected to be a novel target of the add-on treatment of airway hyperresponsiveness.

Structural Asymmetry of Phosphodiesterase-9A and a Unique Pocket for Selective Binding of a Potent Enantiomeric Inhibitor

Phosphodiesterase-9 (PDE9) inhibitors have been studied as potential therapeutics for treatment of central nervous system diseases and diabetes. Here, we report the discovery of a new category of PDE9 inhibitors by rational design on the basis of the crystal structures. The best compound, (S)-6-((1-(4-chlorophenyl)ethyl)amino)-1-cyclopentyl-1,5,6,7-tetrahydro-4H-pyrazolo[3,4-day]pyrimidin-4-one [(S)-C33], has an IC50 value of 11 nM against PDE9 and the racemic C33 has bioavailability of 56.5% in the rat pharmacokinetic model. The crystal structures of PDE9 in the complex with racemic C33, (R)-C33, and (S)-C33 reveal subtle conformational asymmetry of two M-loops in the PDE9 dimer and different conformations of two C33 enantiomers. The structures also identified a small hydrophobic pocket that interacts with the tyrosyl tail of (S)-C33 but not with (R)-C33, and is thus possibly useful for improvement of selectivity of PDE9 inhibitors. The asymmetry of the M-loop and the different interactions of the C33 enantiomers imply the necessity to consider the whole PDE9 dimer in the design of inhibitors.

Discovery of a phosphodiesterase 9A inhibitor as a potential hypoglycemic agent

Phosphodiesterase 9 (PDE9) inhibitors have been studied as potential therapeutics for treatment of diabetes and Alzheimer’s disease. Here we report a potent PDE9 inhibitor 3r that has an IC₅₀ of 0.6 nM and >150-fold selectivity over other PDEs. The HepG2 cell-based assay shows that 3r inhibits the mRNA expression of phosphoenolpyruvate carboxykinase and glucose 6-phosphatase. These activities of 3r, together with the reasonable pharmacokinetic properties and no acute toxicity at 1200 mg/kg dosage, suggest its potential as a hypoglycemic agent. The crystal structure of PDE9-3r reveals significantly different conformation and hydrogen bonding pattern of 3r from those of previously published 28s. Both 3r and 28s form a hydrogen bond with Tyr424, a unique PDE9 residue (except for PDE8), but 3r shows an additional hydrogen bond with Ala452. This structure information might be useful for design of PDE9 inhibitors.

Molecular dynamics-based discovery of novel phosphodiesterase-9A inhibitors with non-pyrazolopyrimidinone scaffolds

Phosphodiesterase-9A (PDE9A) is a promising therapeutic target for the treatment of diabetes and Alzheimer's disease (AD). The Pfizer PDE9A inhibitor PF-04447943 has completed Phase II clinical trials in subjects with mild to moderate AD in 2013. However, most of the reported PDE9A inhibitors share the same scaffold as pyrazolopyrimidinone, which lacks structural diversity and is unfavorable for the development of novel PDE9A inhibitors. In the present study, a combinatorial method including pharmacophores, molecular docking, molecular dynamics simulations, binding free energy calculations, and bioassay was used to discover novel PDE9A inhibitors with new scaffolds rather than pyrazolopyrimidinones from the SPECS database containing about 200,000 compounds. As a result, 15 hits out of 29 molecules (a hit rate of 52%) with five novel scaffolds were identified to be PDE9A inhibitors with inhibitory affinities no more than 50 μM to enrich the structural diversity, different from the pyrazolopyrimidinone-derived family. The high hit ratio of 52% for this virtual screening method indicated that the combinatorial method is a good compromise between computational cost and accuracy. Binding pattern analyses indicate that those hits with non-pyrazolopyrimidinone scaffolds can bind the same active site pocket of PDE9A as classical PDE9A inhibitors. In addition, structural modification of compound AG-690/40135604 (IC50=8.0 μM) led to a new one, 16, with an improved inhibitory affinity of 2.1 μM as expected. The five novel scaffolds discovered in the present study can be used for the rational design of PDE9A inhibitors with higher affinities.

Structure-based discovery of highly selective phosphodiesterase-9A inhibitors and implications for inhibitor design

A new series of phosphodiesterase-9 (PDE9) inhibitors that contain a scaffold of 6-amino-pyrazolopyrimidinone have been discovered by a combination of structure-based design and computational docking. This procedure significantly saved the load of chemical synthesis and is an effective method for the discovery of inhibitors. The best compound 28 has an IC(50) of 21 nM and 3.3 μM, respectively, for PDE9 and PDE5 and about 3 orders of magnitude of selectivity against other PDE families. The crystal structure of the PDE9 catalytic domain in complex with 28 has been determined and shows a hydrogen bond between 28 and Tyr424. This hydrogen bond may account for the 860-fold selectivity of 28 against PDE1B, in comparison with about 30-fold selectivity of BAY73-6691. Thus, our studies suggest that Tyr424, a unique residue of PDE8 and PDE9, is a potential target for improvement of selectivity of PDE9 inhibitors.

Structural asymmetry of phosphodiesterase-9, potential protonation of a glutamic acid, and role of the invariant glutamine

PDE9 inhibitors show potential for treatment of diseases such as diabetes. To help with discovery of PDE9 inhibitors, we performed mutagenesis, kinetic, crystallographic, and molecular dynamics analyses on the active site residues of Gln453 and its stabilizing partner Glu406. The crystal structures of the PDE9 Q453E mutant (PDE9Q453E) in complex with inhibitors IBMX and (S)-BAY73-6691 showed asymmetric binding of the inhibitors in two subunits of the PDE9Q453E dimer and also the significant positional change of the M-loop at the active site. The kinetic analysis of the Q453E and E406A mutants suggested that the invariant glutamine is critical for binding of substrates and inhibitors, but is unlikely to play a key role in the differentiation between substrates of cGMP and cAMP. The molecular dynamics simulations suggest that residue Glu406 may be protonated and may thus explain the hydrogen bond distance between two side chain oxygens of Glu453 and Glu406 in the crystal structure of the PDE9Q453E mutant. The information from these studies may be useful for design of PDE9 inhibitors.