The mechanism of cyclic nucleotide hydrolysis in the phosphodiesterase catalytic site

Combined ligand-based and structure-based design of PDE 9A inhibitors against Alzheimer's disease

PDE9 enzyme hydrolyzes cGMP, which is involved in the regulation of synaptic plasticity through the NMDA pathway (a well-known excitotoxic target for AD) via activation of calcium/calmodulin-dependent neuronal NO synthases in the postsynaptic neurons. The inhibition of PDE9 leads to elevated cGMP levels, causing enhanced NMDA signaling and thus contributing to an increase in synaptic plasticity and stabilization. Therefore, it could be considered a pertinent target for AD drug discovery. PF-04447943 and BI-409306 targeting PDE9 are undergoing clinical trials (Phase II). The present study encompasses a pharmacophoric approach to identify potent PDE9 inhibitors using various computational methods. Pharmacophores generated from the PDB 6A3N yielded 37,554 virtual hits, which underwent drug likeliness and PAINS filtering to arrive at a few virtual leads. The leads were further subjected to extra precision docking, ADMET predictions, and molecular dynamics. The final hits, ZINC000001305675 and ZINC000000377099, showed superior docking scores of - 10.90 and - 10.30 kcal/mol and satisfactory predicted ADMET scores. The hits were subjected to molecular dynamics (MD) studies, wherein they formed stable complexes with PDE9 protein and had ligand RMSDs within acceptable limits. The processes involved in the combined ligand and structure-based strategies.

The stereospecific interaction sites and target specificity of cGMP analogs in mouse cortex

cGMP interactors play a role in several pathologies and may be targets for cGMP analog-based drugs, but the success of targeting depends on the biochemical stereospecificity between the cGMP-analog and the interactor. The stereospecificity between general cGMP analogs-or such that are selectivity-modified to obtain, for example, inhibitory actions on a specific target, like the cGMP-dependent protein kinase-have previously been investigated. However, the importance of stereospecificity for cGMP-analog binding to interactors is not known. We, therefore, applied affinity chromatography on mouse cortex proteins utilizing analogs with cyclic phosphate (8-AET-cGMP, 2-AH-cGMP, 2'-AHC-cGMP) and selectivity-modified analogs with sulfur-containing cyclic phosphorothioates (Rp/Sp-8-AET-cGMPS, Rp/Sp-2'-AHC-cGMPS) immobilized to agaroses. The results illustrate the cGMP analogs' stereospecific binding for PKG, PKA regulatory subunits and PKA catalytic subunits, PDEs, and EPAC2 and the involvement of these in various KEGG pathways. For the seven agaroses, PKG, PKA regulatory subunits, and PKA catalytic subunits were more prone to be enriched by 2-AH-, 8-AET-, Rp-8-AET-, and Sp-8-AET-cGMP, whereas PDEs and EPAC2 were more likely to be enriched by 2-AH-, Rp-2'-AHC-, and Rp-8-AET-cGMP. Our findings help elucidate the stereospecific-binding sites essential for the interaction between individual cGMP analogs and cGMP-binding proteins, as well as the cGMP analogs' target specificity, which are two crucial parameters in drug design.

Photoreceptor Phosphodiesterase (PDE6): Structure, Regulatory Mechanisms, and Implications for Treatment of Retinal Diseases

The photoreceptor phosphodiesterase (PDE6) is a member of large family of Class I phosphodiesterases responsible for hydrolyzing the second messengers cAMP and cGMP. PDE6 consists of two catalytic subunits and two inhibitory subunits that form a tetrameric protein. PDE6 is a peripheral membrane protein that is localized to the signal-transducing compartment of rod and cone photoreceptors. As the central effector enzyme of the G-protein coupled visual transduction pathway, activation of PDE6 catalysis causes a rapid decrease in cGMP levels that results in closure of cGMP-gated ion channels in the photoreceptor plasma membrane. Because of its importance in the phototransduction pathway, mutations in PDE6 genes result in various retinal diseases that currently lack therapeutic treatment strategies due to inadequate knowledge of the structure, function, and regulation of this enzyme. This review focuses on recent progress in understanding the structure of the regulatory and catalytic domains of the PDE6 holoenzyme, the central role of the multi-functional inhibitory γ-subunit, the mechanism of activation by the heterotrimeric G protein, transducin, and future directions for pharmacological interventions to treat retinal degenerative diseases arising from mutations in the PDE6 genes.

PDE9 Inhibitor PF-04447943 Attenuates DSS-Induced Colitis by Suppressing Oxidative Stress, Inflammation, and Regulating T-Cell Polarization

Ulcerative colitis (UC) is a form of inflammatory bowel disease, which manifests as irritation or swelling and sores in the large intestine in a relapsing and remitting manner. In a dextran sulfate sodium sulfate (DSS)-induced UC model in female mice, we found that the levels of cyclic guanosine monophosphate (cGMP) are reduced, while the expression of phosphodiesterase 9A (PDE9A) is highest among all phosphodiesterase (PDEs). Since PDE9 has the highest affinity toward cGMP, we evaluated the selective PDE9 inhibitor PF-04447943 (PF) as a potential candidate for UC treatment. PF has been extensively studies in cognitive function and in sickle cell disease, but not in models for inflammatory bowel disease (IBD). Therefore, we used female C57BL/6 mice treated with 3% DSS alone or co-treated with PF or sulfasalazine (SASP) to study the body weight, colon length, histopathology, and measure superoxide dismutase (SOD), malondialdehyde (MDA), and cGMP level, as well as cytokines such as tumor necrosis factor-alpha (TNF-α), interleukin-6 (IL-6), interleukin-17 (IL-17), interleukin-12/23 (IL-12/23), interleukin-10 (IL-10), and pathways including nuclear factor kappa B (NF-κB), signal transducer and activator of transcription 3 (STAT3), and inflammasome activation. In addition, the number of dendritic cells (DC) and regulatory T cells (Treg cell) was assessed in the spleen, lymph node, and colon using flow cytometry. DSS reduced the number of goblet cells, decreased colon lengths and body weights, all of them were attenuated by PF treatment. It also suppressed the elevated level of inflammatory cytokines and increased level the anti-inflammatory cytokine, IL-10. PF treatment also reduced the DSS-induced inflammation by suppressing oxidative stress, NF-κB, STAT3, and inflammasome activation, by upregulating nuclear factor erythroid 2-related factor 2 (Nrf-2) and its downstream proteins via extracellular signal-regulated kinase (ERK) phosphorylation. Importantly, PF reversed imbalance in Treg/T helper 17 cells (Th17) cells ratio, possibly by regulating dendritic cells and Treg developmental process. In summary, this study shows the protective effect of a PDE9A inhibitor in ulcerative colitis by suppressing oxidative stress and inflammation as well as reversing the Treg/Th17 cells imbalance.

Crystal structures of SAMHD1 inhibitor complexes reveal the mechanism of water-mediated dNTP hydrolysis

SAMHD1 regulates cellular 2'-deoxynucleoside-5'-triphosphate (dNTP) homeostasis by catalysing the hydrolysis of dNTPs into 2'-deoxynucleosides and triphosphate. In CD4⁺ myeloid lineage and resting T-cells, SAMHD1 blocks HIV-1 and other viral infections by depletion of the dNTP pool to a level that cannot support replication. SAMHD1 mutations are associated with the autoimmune disease Aicardi-Goutières syndrome and hypermutated cancers. Furthermore, SAMHD1 sensitises cancer cells to nucleoside-analogue anti-cancer therapies and is linked with DNA repair and suppression of the interferon response to cytosolic nucleic acids. Nevertheless, despite its requirement in these processes, the fundamental mechanism of SAMHD1-catalysed dNTP hydrolysis remained unknown. Here, we present structural and enzymological data showing that SAMHD1 utilises an active site, bi-metallic iron-magnesium centre that positions a hydroxide nucleophile in-line with the P^α-O^5' bond to catalyse phosphoester bond hydrolysis. This precise molecular mechanism for SAMHD1 catalysis, reveals how SAMHD1 down-regulates cellular dNTP and modulates the efficacy of nucleoside-based anti-cancer and anti-viral therapies.

Synthesis and degradation of cAMP in Giardia lamblia: possible role and characterization of a nucleotidyl cyclase with a single cyclase homology domain

Despite its importance in the regulation of growth and differentiation processes of a variety of organisms, the mechanism of synthesis and degradation of cAMP (cyclic AMP) has not yet been described in Giardia lamblia In this work, we measured significant quantities of cAMP in trophozoites of G. lamblia incubated in vitro and later detected how it increases during the first hours of encystation, and how it then returns to basal levels at 24 h. Through an analysis of the genome of G. lamblia, we found sequences of three putative enzymes - one phosphodiesterase (gPDE) and two nucleotidyl cyclases (gNC1 and gNC2) - that should be responsible for the regulation of cAMP in G. lamblia Later, an RT-PCR assay confirmed that these three genes are expressed in trophozoites. The bioinformatic analysis indicated that gPDE is a transmembrane protein of 154 kDa, with a single catalytic domain in the C-terminal end; gNC1 is predicted to be a transmembrane protein of 74 kDa, with only one class III cyclase homology domain (CHD) at the C-terminal end; and gNC2 should be a transmembrane protein of 246 kDa, with two class III CHDs. Finally, we cloned and enriched the catalytic domain of gNC1 (gNC1cd) from bacteria. After that, we confirmed that gNC1cd has adenylyl cyclase (AC) activity. This enzymatic activity depends on the presence of Mn²⁺ and Ca²⁺, but no significant activity was displayed in the presence of Mg²⁺ Additionally, the AC activity of gNC1cd is competitively inhibited with GTP, so it is highly possible that gNC1 has guanylyl cyclase activity as well.

Assessing protein-ligand binding modes with computational tools: the case of PDE4B

In a first step in the discovery of novel potent inhibitor structures for the PDE4B family with limited side effects, we present a protocol to rank newly designed molecules through the estimation of their IC[Formula: see text] values. Our protocol is based on reproducing the linear relationship between the logarithm of experimental IC[Formula: see text] values [[Formula: see text](IC[Formula: see text])] and their calculated binding free energies ([Formula: see text]). From 13 known PDE4B inhibitors, we show here that (1) binding free energies obtained after a docking process by AutoDock are not accurate enough to reproduce this linear relationship; (2) MM-GB/SA post-processing of molecular dynamics (MD) trajectories of the top ranked AutoDock pose improves the linear relationship; (3) by taking into account all representative structures obtained by AutoDock and by averaging MM-GB/SA computations on a series of 40 independent MD trajectories, a linear relationship between [Formula: see text](IC[Formula: see text]) and the lowest [Formula: see text] is achieved with [Formula: see text].

Suppression of β-catenin/TCF transcriptional activity and colon tumor cell growth by dual inhibition of PDE5 and 10

Previous studies suggest the anti-inflammatory drug, sulindac inhibits tumorigenesis by a COX independent mechanism involving cGMP PDE inhibition. Here we report that the cGMP PDE isozymes, PDE5 and 10, are elevated in colon tumor cells compared with normal colonocytes, and that inhibitors and siRNAs can selectively suppress colon tumor cell growth. Combined treatment with inhibitors or dual knockdown suppresses tumor cell growth to a greater extent than inhibition from either isozyme alone. A novel sulindac derivative, ADT-094 was designed to lack COX-1/-2 inhibitory activity but have improved potency to inhibit PDE5 and 10. ADT-094 displayed >500 fold higher potency to inhibit colon tumor cell growth compared with sulindac by activating cGMP/PKG signaling to suppress proliferation and induce apoptosis. Combined inhibition of PDE5 and 10 by treatment with ADT-094, PDE isozyme-selective inhibitors, or by siRNA knockdown also suppresses β-catenin, TCF transcriptional activity, and the levels of downstream targets, cyclin D1 and survivin. These results suggest that dual inhibition of PDE5 and 10 represents novel strategy for developing potent and selective anticancer drugs.

Analysis of proton wires in the enzyme active site suggests a mechanism of c-di-GMP hydrolysis by the EAL domain phosphodiesterases

We report for the first time a hydrolysis mechanism of the cyclic dimeric guanosine monophosphate (c-di-GMP) by the EAL domain phosphodiesterases as revealed by molecular simulations. A model system for the enzyme-substrate complex was prepared on the base of the crystal structure of the EAL domain from the BlrP1 protein complexed with c-di-GMP. The nucleophilic hydroxide generated from the bridging water molecule appeared in a favorable position for attack on the phosphorus atom of c-di-GMP. The most difficult task was to find a pathway for a proton transfer to the O3' atom of c-di-GMP to promote the O3'P bond cleavage. We show that the hydrogen bond network extended over the chain of water molecules in the enzyme active site and the Glu359 and Asp303 side chains provides the relevant proton wires. The suggested mechanism is consistent with the structural, mutagenesis, and kinetic experimental studies on the EAL domain phosphodiesterases. Proteins 2016; 84:1670-1680. © 2016 Wiley Periodicals, Inc.

Cyclic di-GMP: the first 25 years of a universal bacterial second messenger

Twenty-five years have passed since the discovery of cyclic dimeric (3'→5') GMP (cyclic di-GMP or c-di-GMP). From the relative obscurity of an allosteric activator of a bacterial cellulose synthase, c-di-GMP has emerged as one of the most common and important bacterial second messengers. Cyclic di-GMP has been shown to regulate biofilm formation, motility, virulence, the cell cycle, differentiation, and other processes. Most c-di-GMP-dependent signaling pathways control the ability of bacteria to interact with abiotic surfaces or with other bacterial and eukaryotic cells. Cyclic di-GMP plays key roles in lifestyle changes of many bacteria, including transition from the motile to the sessile state, which aids in the establishment of multicellular biofilm communities, and from the virulent state in acute infections to the less virulent but more resilient state characteristic of chronic infectious diseases. From a practical standpoint, modulating c-di-GMP signaling pathways in bacteria could represent a new way of controlling formation and dispersal of biofilms in medical and industrial settings. Cyclic di-GMP participates in interkingdom signaling. It is recognized by mammalian immune systems as a uniquely bacterial molecule and therefore is considered a promising vaccine adjuvant. The purpose of this review is not to overview the whole body of data in the burgeoning field of c-di-GMP-dependent signaling. Instead, we provide a historic perspective on the development of the field, emphasize common trends, and illustrate them with the best available examples. We also identify unresolved questions and highlight new directions in c-di-GMP research that will give us a deeper understanding of this truly universal bacterial second messenger.

C-di-GMP hydrolysis by Pseudomonas aeruginosa HD-GYP phosphodiesterases: analysis of the reaction mechanism and novel roles for pGpG

In biofilms, the bacterial community optimizes the strategies to sense the environment and to communicate from cell to cell. A key player in the development of a bacterial biofilm is the second messenger c-di-GMP, whose intracellular levels are modulated by the opposite activity of diguanylate cyclases and phosphodiesterases. Given the huge impact of bacterial biofilms on human health, understanding the molecular details of c-di-GMP metabolism represents a critical step in the development of novel therapeutic approaches against biofilms. In this study, we present a detailed biochemical characterization of two c-di-GMP phosphodiesterases of the HD-GYP subtype from the human pathogen Pseudomonas aeruginosa, namely PA4781 and PA4108. Upstream of the catalytic HD-GYP domain, PA4781 contains a REC domain typical of two-component systems, while PA4108 contains an uncharacterized domain of unknown function. Our findings shed light on the activity and catalytic mechanism of these phosphodiesterases. We show that both enzymes hydrolyse c-di-GMP in a two-step reaction via the linear intermediate pGpG and that they produce GMP in vitro at a surprisingly low rate. In addition, our data indicate that the non-phosphorylated REC domain of PA4781 prevents accessibility of c-di-GMP to the active site. Both PA4108 and phosphorylated PA4781 are also capable to use pGpG as an alternative substrate and to hydrolyse it into GMP; the affinity of PA4781 for pGpG is one order of magnitude higher than that for c-di-GMP. These results suggest that these enzymes may not work (primarily) as genuine phosphodiesterases. Moreover, the unexpected affinity of PA4781 for pGpG may indicate that pGpG could also act as a signal molecule in its own right, thus further widening the c-di-GMP-related signalling scenario.

Synthesis, docking study and relaxant effect of 2-alkyl and 2-naphthylchromones on rat aorta and guinea-pig trachea through phosphodiesterase inhibition

Chromone (4), which form the base structure of various flavonoids isolated as natural products, is capable of relaxing smooth muscle. This is relevant to the treatment of high blood pressure, asthma and chronic obstructive pulmonary disease. The former disorder involves the contraction of vascular smooth muscle (VSM), and the latter two bronchoconstriction of airway smooth muscle (ASM). One of the principal mechanisms by which flavonoids relax muscle tissue is the inhibition of phosphodiesterases (PDEs), present in both VSM and ASM. Therefore, a study was designed to analyze the structure-activity relationship of chromone derivatives in vaso- and bronchorelaxation through the inhibition of PDE. Docking studies showed that these chromones bind at the catalytic site of PDEs. Consequently, we synthesized analogs of chromones substituted at position C-2 with alkyl and naphthyl groups. These compounds were synthesized from 2-hydroxyacetophenone and acyl chlorides in the presence of DBU and pyridine, modifying the methodology reported for the synthesis of 3-acylchromones by changing the reaction temperature from 80 to 30°C and using methylene chloride as solvent, yielding the corresponding phenolic esters 10a-10h. These compounds were cyclized with an equivalent of DBU, pyridine as solvent, and heated at reflux temperature, yielding the chromones 11a-11h. Evaluation of the vasorelaxant effect of 4, 11a-11h on rat aorta demonstrated that potency decreases with branched alkyl groups. Whereas the EC50 of compound 11d (substituted by an n-hexyl group) was 8.64±0.39 μM, that of 11f (substituted by an isobutyl group) was 14.58±0.64 μM. Contrarily, the effectiveness of the compound is directly proportional to the length of the alkyl chain, as evidenced by the increase in maximal effect of compound 11c versus 11d (66% versus 100%) and 11e versus 11f (60% versus 96%). With an aromatic group like naphthyl as the C-2 substituent, the effectiveness was only 43%. All compounds tested on guinea pig trachea showed less than 55% effectiveness. Compounds 4, 11a-11h were evaluated as PDE inhibitors in vitro, with 11d showing the greatest effect (73%), corroborating the importance of a long alkyl chain, which inhibits the decomposition of cGMP. Docking studies showed that the compound 11d was selective for the inhibition of PDE-5.

An update view on the substrate recognition mechanism of phosphodiesterases: a computational study of PDE10 and PDE4 bound with cyclic nucleotides

Early studies strongly implied that the specificity of cyclic nucleotide phosphodiesterases (PDEs) toward its endogenous substrates can be uniquely determined by the amido orientation of the invariant glutamine locating in the binding pocket of the enzyme. However, recently solved crystal structures of PDE4 (cAMP specific) and PDE10 (dual specific) in the presence of endogenous substrates have revealed that their invariant glutamine orientations are very similar despite exhibiting different substrate specificities proven physiologically. To understand this subtle specificity issue in the PDE family, here several experimentally inaccessible PDE-substrate complex models have been studied computationally, and the results are juxtaposed and compared in detail. Modeling results show that PDE10 in fact favors cAMP energetically but still can bind to cGMP owing to the robust hydrogen-bond network in the vicinity of the invariant glutamine side chain. PDE4 fails to accommodate cGMP is correlated to the weakening of this same hydrogen-bond network but not owing to any steric strain in the binding pocket. An Asn residue in the binding pocket of PDE4 has enhanced the specificity of the binding to cAMP sideway as observed in our computer simulation. Further to the previously studied syn- versus anti-conformational specificity of cAMP in PDE10, the unexpected substrate-binding mode in PDE10 versus PDE4 as reported here strongly suggested that there are remaining uncertainties in the substrate orientation and recognition mechanism in the PDE families. The molecular details of the binding pocket observed in this study provide hints for more optimal PDE4 and PDE10 inhibitor design.

THEORETICAL STUDIES ON THE MECHANISM OF CYCLIC NUCLEOTIDE MONOPHOSPHATE HYDROLYSIS WITHIN PHOSPHODIESTERASES

As the only metabolizing enzyme for the degradation of second messenger cAMP and cGMP, phosphodiesterase (PDE) has been the clinical target of various human diseases. But the hydrolysis procedure of PDE is still unclear. To investigate the mechanism of PDE catalysis, three types of PDE (PDE4d, PDE5a and PDE10a) were selected and studied by using molecular dynamics (MD) simulation and quantum mechanics (QM) calculation methods. MD Simulation results indicate that different PDEs share a similar hydrolysis area in the active sites, and the phosphate parts of cyclic nucleotides take the same orientation and are partly surrounded by water molecules. Based on the statistical data of MD simulation, the QM calculation models were built. The calculation results indicate that in aqueous solution, the nucleophile hydroxide ion that attacks the phosphor atom of the cyclic nucleotide in the hydrolysis may migrate between the two metal ions in the active site. To help the ring-open reaction, it is the water molecule that provides proton to the O3′ atom of cyclic nucleotide, and generates another hydroxide ion complexed with the metal ion.

Fundamental reaction pathway and free energy profile for hydrolysis of intracellular second messenger adenosine 3',5'-cyclic monophosphate (cAMP) catalyzed by phosphodiesterase-4

As important drug targets for a variety of human diseases, cyclic nucleotide phosphodiesterases (PDEs) are a superfamily of enzymes sharing a similar catalytic site. We have performed pseudobond first-principles quantum mechanical/molecular mechanical-free energy perturbation (QM/MM-FE) and QM/MM-Poisson-Boltzmann surface area (PBSA) calculations to uncover the detailed reaction mechanism for PDE4-catalyzed hydrolysis of adenosine 3',5'-cyclic monophosphate (cAMP). This is the first report on QM/MM reaction-coordinate calculations including the protein environment of any PDE-catalyzed reaction system, demonstrating a unique catalytic reaction mechanism. The QM/MM-FE and QM/MM-PBSA calculations revealed that the PDE4-catalyzed hydrolysis of cAMP consists of two reaction stages: cAMP hydrolysis (stage 1) and bridging hydroxide ion regeneration (stage 2). The stage 1 includes the binding of cAMP in the active site, nucleophilic attack of the bridging hydroxide ion on the phosphorus atom of cAMP, cleavage of O3'-P phosphoesteric bond of cAMP, protonation of the departing O3' atom, and dissociation of hydrolysis product (AMP). The stage 2 includes the binding of solvent water molecules with the metal ions in the active site and regeneration of the bridging hydroxide ion. The dissociation of the hydrolysis product is found to be rate-determining for the enzymatic reaction process. The calculated activation Gibbs free energy of ≥16.0 and reaction free energy of -11.1 kcal/mol are in good agreement with the experimentally derived activation free energy of 16.6 kcal/mol and reaction free energy of -11.5 kcal/mol, suggesting that the catalytic mechanism obtained from this study is reliable and provides a solid base for future rational drug design.

Characterization of the deoxynucleotide triphosphate triphosphohydrolase (dNTPase) activity of the EF1143 protein from Enterococcus faecalis and crystal structure of the activator-substrate complex

The EF1143 protein from Enterococcus faecalis is a distant homolog of deoxynucleotide triphosphate triphosphohydrolases (dNTPases) from Escherichia coli and Thermus thermophilus. These dNTPases are important components in the regulation of the dNTP pool in bacteria. Biochemical assays of the EF1143 dNTPase activity demonstrated nonspecific hydrolysis of all canonical dNTPs in the presence of Mn(2+). In contrast, with Mg(2+) hydrolysis required the presence of dGTP as an effector, activating the degradation of dATP and dCTP with dGTP also being consumed in the reaction with dATP. The crystal structure of EF1143 and dynamic light scattering measurements in solution revealed a tetrameric oligomer as the most probable biologically active unit. The tetramer contains four dGTP specific allosteric regulatory sites and four active sites. Examination of the active site with the dATP substrate suggests an in-line nucleophilic attack on the α-phosphate center as a possible mechanism of the hydrolysis and two highly conserved residues, His-129 and Glu-122, as an acid-base catalytic dyad. Structural differences between EF1143 apo and holo forms revealed mobility of the α3 helix that can regulate the size of the active site binding pocket and could be stabilized in the open conformation upon formation of the tetramer and dGTP effector binding.

Insight into the phosphodiesterase mechanism from combined QM/MM free energy simulations

Molecular dynamics simulations employing a combined quantum mechanical and molecular mechanical potential have been carried out to elucidate the reaction mechanism of the hydrolysis of a cyclic nucleotide cAMP substrate by phosphodiesterase 4B (PDE4B). PDE4B is a member of the PDE superfamily of enzymes that play crucial roles in cellular signal transduction. We have determined a two-dimensional potential of mean force (PMF) for the coupled phosphoryl bond cleavage and proton transfer through a general acid catalysis mechanism in PDE4B. The results indicate that the ring-opening process takes place through an S(N)2 reaction mechanism, followed by a proton transfer to stabilize the leaving group. The computed free energy of activation for the PDE4B-catalyzed cAMP hydrolysis is about 13 kcal·mol(-1) and an overall reaction free energy is about -17 kcal·mol(-1), both in accord with experimental results. In comparison with the uncatalyzed reaction in water, the enzyme PDE4B provides a strong stabilization of the transition state, lowering the free energy barrier by 14 kcal·mol(-1). We found that the proton transfer from the general acid residue His234 to the O3' oxyanion of the ribosyl leaving group lags behind the nucleophilic attack, resulting in a shallow minimum on the free energy surface. A key contributing factor to transition state stabilization is the elongation of the distance between the divalent metal ions Zn(2+) and Mg(2+) in the active site as the reaction proceeds from the Michaelis complex to the transition state.

Mammalian Cyclic Nucleotide Phosphodiesterases: Molecular Mechanisms and Physiological Functions

The superfamily of cyclic nucleotide (cN) phosphodiesterases (PDEs) is comprised of 11 families of enzymes. PDEs break down cAMP and/or cGMP and are major determinants of cellular cN levels and, consequently, the actions of cN-signaling pathways. PDEs exhibit a range of catalytic efficiencies for breakdown of cAMP and/or cGMP and are regulated by myriad processes including phosphorylation, cN binding to allosteric GAF domains, changes in expression levels, interaction with regulatory or anchoring proteins, and reversible translocation among subcellular compartments. Selective PDE inhibitors are currently in clinical use for treatment of erectile dysfunction, pulmonary hypertension, intermittent claudication, and chronic pulmonary obstructive disease; many new inhibitors are being developed for treatment of these and other maladies. Recently reported x-ray crystallographic structures have defined features that provide for specificity for cAMP or cGMP in PDE catalytic sites or their GAF domains, as well as mechanisms involved in catalysis, oligomerization, autoinhibition, and interactions with inhibitors. In addition, major advances have been made in understanding the physiological impact and the biochemical basis for selective localization and/or recruitment of specific PDE isoenzymes to particular subcellular compartments. The many recent advances in understanding PDE structures, functions, and physiological actions are discussed in this review.

Computational determination of binding structures and free energies of phosphodiesterase-2 with benzo[1,4]diazepin-2-one derivatives

Phosphodiesterase-2 (PDE2) is a key enzyme catalyzing hydrolysis of both cyclic adenosine monophosphate (cAMP) and cyclic guanosine monophosphate (cGMP) that serve as intracellular second messengers. PDE2 has been recognized as an attractive drug target, and selective inhibitors of PDE2 are expected to be promising candidates for the memory enhancer, antidepressant, and anxiolytic agent. In the present study, we examined the detailed binding structures and free energies for PDE2 interacting with a promising series of inhibitors, i.e., benzo[1,4]diazepin-2-one derivatives, by carrying out molecular docking, molecular dynamics (MD) simulations, binding free energy calculations, and binding energy decompositions. The computational results provide valuable insights into the detailed enzyme-inhibitor binding modes including important intermolecular interactions, e.g., the π-π stacking interactions with the common benzo[1,4]diazepin-2-one scaffold of the inhibitors, hydrogen bonding and hydrophobic interactions with the substituents on the benzo[1,4]diazepin-2-one scaffold. Future rational design of new, more potent inhibitors of PDE2 should carefully account for all of these favorable intermolecular interactions. By use of the MD-simulated binding structures, the calculated binding free energies are in good agreement with the experimental activity data for all of the examined benzo[1,4]diazepin-2-one derivatives. The enzyme-inhibitor binding modes determined and the agreement between the calculated and experimental results are expected to be valuable for future rational design of more potent inhibitors of PDE2.

Structural insight into the mechanism of c-di-GMP hydrolysis by EAL domain phosphodiesterases

Cyclic diguanylate (or bis-(3'-5') cyclic dimeric guanosine monophosphate; c-di-GMP) is a ubiquitous second messenger that regulates diverse cellular functions, including motility, biofilm formation, cell cycle progression, and virulence in bacteria. In the cell, degradation of c-di-GMP is catalyzed by highly specific EAL domain phosphodiesterases whose catalytic mechanism is still unclear. Here, we purified 13 EAL domain proteins from various organisms and demonstrated that their catalytic activity is associated with the presence of 10 conserved EAL domain residues. The crystal structure of the TBD1265 EAL domain was determined in free state (1.8 Å) and in complex with c-di-GMP (2.35 A), and unveiled the role of conserved residues in substrate binding and catalysis. The structure revealed the presence of two metal ions directly coordinated by six conserved residues, two oxygens of c-di-GMP phosphate, and potential catalytic water molecule. Our results support a two-metal-ion catalytic mechanism of c-di-GMP hydrolysis by EAL domain phosphodiesterases.

A substrate selectivity and inhibitor design lesson from the PDE10-cAMP crystal structure: a computational study

Phosphodiesterases (PDEs) catalyze the hydrolysis of second messengers cAMP and cGMP in regulating many important cellular signals and have been recognized as important drug targets. Experimentally, a range of specificity/selectivity toward cAMP and cGMP is well-known for the individual PDE families. The study reported here reveals that PDEs might also exhibit selectivity toward conformations of the endogenous substrates cAMP and cGMP. Molecular dynamics simulations and free energy study have been applied to study the binding of the cAMP torsional conformers about the glycosyl bond in PDE10A2. The computational results elucidated that PDE10A2 is energetically more favorable in complex with the syn cAMP conformer (as reported in the crystal structure) and the binding of anti cAMP to PDE10A2 would lead to either a nonreactive configuration or significant perturbation on the catalytic pocket of the enzyme. This experimentally inaccessible information provides important molecular insights for the development of effective PDE10 ligands.

Structure and mechanism of a bacterial light-regulated cyclic nucleotide phosphodiesterase

The ability to respond to light is crucial for most organisms. BLUF is a recently identified photoreceptor protein domain that senses blue light using a FAD chromophore. BLUF domains are present in various proteins from the Bacteria, Euglenozoa and Fungi. Although structures of single-domain BLUF proteins have been determined, none are available for a BLUF protein containing a functional output domain; the mechanism of light activation in this new class of photoreceptors has thus remained poorly understood. Here we report the biochemical, structural and mechanistic characterization of a full-length, active photoreceptor, BlrP1 (also known as KPN_01598), from Klebsiella pneumoniae. BlrP1 consists of a BLUF sensor domain and a phosphodiesterase EAL output domain which hydrolyses cyclic dimeric GMP (c-di-GMP). This ubiquitous second messenger controls motility, biofilm formation, virulence and antibiotic resistance in the Bacteria. Crystal structures of BlrP1 complexed with its substrate and metal ions involved in catalysis or in enzyme inhibition provide a detailed understanding of the mechanism of the EAL-domain c-di-GMP phosphodiesterases. These structures also sketch out a path of light activation of the phosphodiesterase output activity. Photon absorption by the BLUF domain of one subunit of the antiparallel BlrP1 homodimer activates the EAL domain of the second subunit through allosteric communication transmitted through conserved domain-domain interfaces.

Structural basis for the catalytic mechanism of human phosphodiesterase 9

The phosphodiesterases (PDEs) are metal ion-dependent enzymes that regulate cellular signaling by metabolic inactivation of the ubiquitous second messengers cAMP and cGMP. In this role, the PDEs are involved in many biological and metabolic processes and are proven targets of successful drugs for the treatments of a wide range of diseases. However, because of the rapidity of the hydrolysis reaction, an experimental knowledge of the enzymatic mechanisms of the PDEs at the atomic level is still lacking. Here, we report the structures of reaction intermediates accumulated at the reaction steady state in PDE9/crystal and preserved by freeze-trapping. These structures reveal the catalytic process of a PDE and explain the substrate specificity of PDE9 in an actual reaction and the cation requirements of PDEs in general.