Identification of Novel Phosphodiesterase-4D Inhibitors Prescreened by Molecular Dynamics-Augmented Modeling and Validated by Bioassay

Hepatotoxicity assessment investigations on PFASs targeting L-FABP using binding affinity data and machine learning-based QSAR model

Per- and polyfluoroalkyl substances (PFASs) are persistent organic pollutants that have been detected in various environmental media and human serum, but their safety assessment remains challenging. PFASs may accumulate in liver tissues and cause hepatotoxicity by binding to liver fatty acid binding protein (L-FABP). Therefore, evaluating the binding affinity of PFASs to L-FABP is crucial in assessing the potential hepatotoxic effects. In this study, two binding sites of L-FABP were evaluated, results suggested that the outer site possessed high affinity to polyfluoroalkyl sulfates and the inner site preferred perfluoroalkyl sulfonamides, overall, the inner site of L-FABP was more sensitive to PFASs. The binding affinity data of PFASs to L-FABP were used as training set to develop a machine learning model-based quantitative structure-activity relationship (QSAR) for efficient prediction of potentially hazardous PFASs. Further Bayesian Kernel Machine Regression (BKMR) model disclosed flexibility as the determinant molecular property on PFASs-induced hepatotoxicity. It can influence affinity of PFASs to target protein through affecting binding conformations directly (individual effect) as well as integrating with other molecular properties (joint effect). Our present work provided more understanding on hepatotoxicity of PFASs, which could be significative in hepatotoxicity gradation, administration guidance, and safer alternatives development of PFASs.

A label-free LC/MS-based enzymatic activity assay for the detection of PDE5A inhibitors

Phosphodiesterase type 5 (PDE5), a cyclic nucleotide phosphodiesterase, controls the duration of the cyclic guanosine monophosphate (cGMP) signal by hydrolyzing cGMP to GMP. Inhibiting the activity of PDE5A has proven to be an effective strategy for treating pulmonary arterial hypertension and erectile dysfunction. Current enzymatic activity assay methods for PDE5A mainly use fluorescent or isotope-labeled substrates, which are expensive and inconvenient. Here, we developed an LC/MS-based enzymatic activity assay for PDE5A without labeling, which detects the enzymatic activity of PDE5A by quantifying the substrate cGMP and product GMP at a concentration of 100 nM. The accuracy of this method was verified by a fluorescently labeled substrate. Moreover, a new inhibitor of PDE5A was identified by this method and virtual screening. It inhibited PDE5A with an IC₅₀ value of 870 nM. Overall, the proposed strategy provides a new method for screening PDE5A inhibitors.

Free energy perturbation-based large-scale virtual screening for effective drug discovery against COVID-19

As a theoretically rigorous and accurate method, FEP-ABFE (Free Energy Perturbation-Absolute Binding Free Energy) calculations showed great potential in drug discovery, but its practical application was difficult due to high computational cost. To rapidly discover antiviral drugs targeting SARS-CoV-2 Mpro and TMPRSS2, we performed FEP-ABFE-based virtual screening for ∼12,000 protein-ligand binding systems on a new generation of Tianhe supercomputer. A task management tool was specifically developed for automating the whole process involving more than 500,000 MD tasks. In further experimental validation, 50 out of 98 tested compounds showed significant inhibitory activity towards Mpro, and one representative inhibitor, dipyridamole, showed remarkable outcomes in subsequent clinical trials. This work not only demonstrates the potential of FEP-ABFE in drug discovery but also provides an excellent starting point for further development of anti-SARS-CoV-2 drugs. Besides, ∼500 TB of data generated in this work will also accelerate the further development of FEP-related methods.

Influence of heat processing on the anti-inflammatory activity of fresh Smilax glabra based on PDE4 inhibition

•

Heat processing plays a key role in chemical profiles and health benefits of fresh SG.

•

Fresh SG exhibits significant anti‐inflammatory effect based on PDE4 inhibition.

•

The heat-labile quality and safety aspects of four astilbin isomers are compared.

Smilax glabra Roxb. (SG) is widely used as functional food with various beneficial effects. Fresh SG without processing has been eaten directly for anti-inflammation from ancient China, while the underlying mechanism remains underexplored. This study aims to investigate the anti-inflammatory activity of fresh SG by using metabolites profiles, affinity ultrafiltration mass spectrometry, PDE4 enzyme inhibition assay, and in silico analysis. Encouragingly, fresh SG showed promising anti-inflammatory effect with IC₅₀ value (0.009 μg/μL) on PDE4 was about 12 times higher than that of processed SG (0.110 μg/μL). Astilbin was identified as the main bioactive compound of fresh SG responsible for PDE4 inhibitory activity. We found that heat processing strongly affected astilbin isomerization, leading to significant changes in contents and PDE4 inhibitory activities of four astilbin isomers, resulting in decreased anti-inflammatory activity of fresh SG. This finding will provide theoretical basis for systematic research and food/nutraceutical applications of fresh Smilax glabra in the future.

Free energy perturbation (FEP)-guided scaffold hopping

Scaffold hopping refers to computer-aided screening for active compounds with different structures against the same receptor to enrich privileged scaffolds, which is a topic of high interest in organic and medicinal chemistry. However, most approaches cannot efficiently predict the potency level of candidates after scaffold hopping. Herein, we identified potent PDE5 inhibitors with a novel scaffold via a free energy perturbation (FEP)-guided scaffold-hopping strategy, and FEP shows great advantages to precisely predict the theoretical binding potencies ΔG FEP between ligands and their target, which were more consistent with the experimental binding potencies ΔG EXP (the mean absolute deviations| Δ G FEP - Δ G EXP |  < 2 kcal/mol) than those ΔG MM-PBSA or ΔG MM-GBSA predicted by the MM-PBSA or MM-GBSA method. Lead L12 had an IC50 of 8.7 nmol/L and exhibited a different binding pattern in its crystal structure with PDE5 from the famous starting drug tadalafil. Our work provides the first report via the FEP-guided scaffold hopping strategy for potent inhibitor discovery with a novel scaffold, implying that it will have a variety of future applications in rational molecular design and drug discovery.

A novel inhibitor of N 6-methyladenosine demethylase FTO induces mRNA methylation and shows anti-cancer activities

N⁶-methyladenosine (m⁶A) modification is critical for mRNA splicing, nuclear export, stability and translation. Fat mass and obesity-associated protein (FTO), the first identified m⁶A demethylase, is critical for cancer progression. Herein, we developed small-molecule inhibitors of FTO by virtual screening, structural optimization, and bioassay. As a result, two FTO inhibitors namely 18077 and 18097 were identified, which can selectively inhibit demethylase activity of FTO. Specifically, 18097 bound to the active site of FTO and then inhibited cell cycle process and migration of cancer cells. In addition, 18097 reprogrammed the epi-transcriptome of breast cancer cells, particularly for genes related to P53 pathway. 18097 increased the abundance of m⁶A modification of suppressor of cytokine signaling 1 (SOCS1) mRNA, which recruited IGF2BP1 to increase mRNA stability of SOCS1 and subsequently activated the P53 signaling pathway. Further, 18097 suppressed cellular lipogenesis via downregulation of peroxisome proliferator-activated receptor gamma (PPARγ), CCAAT/enhancer-binding protein alpha (C/EBPα), and C/EBPβ. Animal studies confirmed that 18097 can significantly suppress in vivo growth and lung colonization of breast cancer cells. Collectively, we identified that FTO can work as a potential drug target and the small-molecule inhibitor 18097 can serve as a potential agent against breast cancer.

Probing Allosteric Regulation Mechanism of W7.35 on Agonist-Induced Activity for μOR by Mutation Simulation

The residue located at 15 positions before the most conserved residue in TM7 (7.35 of Ballesteros-Weinstein number) plays important roles in ligand binding and the receptor activity for class A GPCRs. Nevertheless, its regulation mechanism has not been clearly clarified in experiments, and some controversies also exist for its impact on μ-opioid receptors (μOR) bound by agonists. Thus, we chose the μ-opioid receptor (μOR) of class A GPCRs as a representative and utilized a microsecond accelerated molecular dynamics simulation (aMD) coupled with a protein structure network (PSN) to explore the effect of W318^7.35 on its functional activity induced by the agonist endomorphin2 mainly by a comparison of the wild system and its W7.35A mutant. When endomorphin2 binds to the wild-type μOR, TM6 in μOR moves outward to form an open intracellular conformation that is beneficial to accommodating the β-arrestin transducer, rather than the G-protein transducer due to the clash with the α5 helix of G-protein, thus acting as a β-arrestin biased agonist. However, the W318A mutation induces the intracellular part of μOR to form a closed state, which disfavors coupling with either G-protein or β-arrestin. The allosteric pathway analysis further reveals that the binding of endomorphin2 to the wild-type μOR transmits more activation signals to the β-arrestin binding site while the W318A mutation induces more structural signals to transmit to common binding residues of the G protein and β-arrestin. More interestingly, the residue at the 7.35 position regulates the shortest allosteric pathway in indirect ways by influencing the interactions between other ligand-binding residues and endomorphin2. W293^6.48 and F289^6.44 are important for regulating the different activities of μOR induced either by the agonist or by the mutation. Y336^7.53, F343^8.50, and D340^8.47 play crucial roles in activating the β-arrestin biased signal induced by the agonist endomorphin2, while L158^3.43 and V286^6.41 devote important contributions to the change in the activity of endomorphin2 from the β-arrestin biased agonist to the antagonist upon the W318A mutation.

Rational Design of 2-Chloroadenine Derivatives as Highly Selective Phosphodiesterase 8A Inhibitors

To validate the hypothesis that Tyr748 is a crucial residue to aid the discovery of highly selective phosphodiesterase 8A (PDE8A) inhibitors, we identified a series of 2-chloroadenine derivatives based on the hit clofarabine. Structure-based design targeting Tyr748 in PDE8 resulted in the lead compound 3a (IC₅₀ = 0.010 μM) with high selectivity with a reasonable druglike profile. In the X-ray crystal structure, 3a bound to PDE8A with a different mode from 3-isobutyl-1-methylxanthine (a pan-PDE inhibitor) and gave a H-bond of 2.7 Å with Tyr748, which possibly interprets the 220-fold selectivity of 3a against PDE2A. Additionally, oral administration of compound 3a achieved remarkable therapeutic effects against vascular dementia (VaD), indicating that PDE8 inhibitors could serve as potential anti-VaD agents.

Discovery and Optimization of Chromone Derivatives as Novel Selective Phosphodiesterase 10 Inhibitors

Phosphodiesterase 10 (PDE10) inhibitors have received much attention as promising therapeutic agents for central nervous system (CNS) disorders such as schizophrenia and Huntington's disease. Recently, a hit compound 1 with a novel chromone scaffold has shown moderate inhibitory activity against PDE10A (IC₅₀ = 500 nM). Hit-to-lead optimization has resulted in compound 3e with an improved inhibitory activity (IC₅₀ = 6.5 nM), remarkable selectivity (>95-fold over other PDEs), and good metabolic stability (RLM t_1/2 = 105 min) by using an integrated strategy (molecular modeling, chemical synthesis, bioassay, and cocrystal structure). The cocrystal structural information provides insights into the binding pattern of 3e in the PDE10A catalytic domain to highlight the key role of the halogen and hydrogen bonds toward Tyr524 and Tyr693, respectively, thereby resulting in high selectivity against other PDEs. These new observations are of benefit for the rational design of the next generation PDE10 inhibitors for CNS disorders.

Inverse Molecular Docking as a Novel Approach to Study Anticarcinogenic and Anti-Neuroinflammatory Effects of Curcumin

Research efforts are placing an ever increasing emphasis on identifying signal transduction pathways related to the chemopreventive activity of curcumin. Its anticarcinogenic effects are presumably mediated by the regulation of signaling cascades, including nuclear factor κB (NF-κB), activator protein 1 (AP-1), and mitogen-activated protein kinases (MAPK). By modulating signal transduction pathways, curcumin induces apoptosis in malignant cells, thus inhibiting cancer development and progression. Due to the lack of mechanistic insight in the scientific literature, we developed a novel inverse molecular docking protocol based on the CANDOCK algorithm. For the first time, we performed inverse molecular docking of curcumin into a collection of 13,553 available human protein structures from the Protein Data Bank resulting in prioritized target proteins of curcumin. Our predictions were in agreement with the scientific literature and confirmed that curcumin binds to folate receptor β, DNA (cytosine-5)-methyltransferase 3A, metalloproteinase-2, mitogen-activated protein kinase 9, epidermal growth factor receptor and apoptosis-inducing factor 1. We also identified new potential protein targets of curcumin, namely deoxycytidine kinase, NAD-dependent protein deacetylase sirtuin-1 and -2, ecto-5'-nucleotidase, core histone macro-H2A.1, tyrosine-protein phosphatase non-receptor type 11, macrophage colony-stimulating factor 1 receptor, GTPase HRas, aflatoxin B1 aldehyde reductase member 3, aldo-keto reductase family 1 member C3, amiloride-sensitive amine oxidase, death-associated protein kinase 2 and tryptophan-tRNA ligase, that may all play a crucial role in its observed anticancer effects. Moreover, our inverse docking results showed that curcumin potentially binds also to the proteins cAMP-specific 3',5'-cyclic phosphodiesterase 4D and 17-β-hydroxysteroid dehydrogenase type 10, which provides a new explanation for its efficiency in the treatment of Alzheimer's disease. We firmly believe that our computational results will complement and direct future experimental studies on curcumin's anticancer activity as well as on its therapeutic effects against Alzheimer's disease.

Structural Features and Ligand Selectivity for 10 Intermediates in the Activation Process of β2-Adrenergic Receptor

It has already been suggested by researchers that there should be multiple intermediate states in the activation process for G-protein-coupled receptors (GPCRs). However, the intermediate states are very short-lived and hardly captured by the experiments, leading to very limited understanding of their structural features and drug efficacies. In this work, a novel joint strategy of targeted molecular dynamics simulation, conventional molecular dynamics simulation, and virtual screening is developed to address the problems. The results from 10 intermediate conformations obtained from the work reveal that the ligand pocket is very unstable and fluctuates between the inactive state and the active one in the case of ligand-free, in particular for ECL2 as a gate-keeper of the ligand-binding. The ligand-binding site could be stable in the active state with a small volume and a completely closed ECL2, only when the G-protein-binding region is fully activated. In addition, the activations of the ligand-binding pocket and G-protein-binding site are relatively independent and exhibit a loose allosteric coupling, which contributes to the existence of multiple intermediate conformations. Interestingly, the screening performance of the agonists does not increase on increasing the overall activity of the intermediate state, but is dependent on the activated extent of the ligand pocket. The receptor is prone to bind the agonist when closing ECL2 and reducing the ligand-binding pocket volume, whereas it is more favorable for binding the antagonist when opening ECL2 and increasing the pocket volume. These observations added to previous studies could help us better understand the activation mechanism of GPCRs and provide valuable information for drug design.

Phosphodiesterase Inhibitors as a Therapeutic Approach to Neuroprotection and Repair

A wide diversity of perturbations of the central nervous system (CNS) result in structural damage to the neuroarchitecture and cellular defects, which in turn are accompanied by neurological dysfunction and abortive endogenous neurorepair. Altering intracellular signaling pathways involved in inflammation and immune regulation, neural cell death, axon plasticity and remyelination has shown therapeutic benefit in experimental models of neurological disease and trauma. The second messengers, cyclic adenosine monophosphate (cyclic AMP) and cyclic guanosine monophosphate (cyclic GMP), are two such intracellular signaling targets, the elevation of which has produced beneficial cellular effects within a range of CNS pathologies. The only known negative regulators of cyclic nucleotides are a family of enzymes called phosphodiesterases (PDEs) that hydrolyze cyclic nucleotides into adenosine monophosphate (AMP) or guanylate monophosphate (GMP). Herein, we discuss the structure and physiological function as well as the roles PDEs play in pathological processes of the diseased or injured CNS. Further we review the approaches that have been employed therapeutically in experimental paradigms to block PDE expression or activity and in turn elevate cyclic nucleotide levels to mediate neuroprotection or neurorepair as well as discuss both the translational pathway and current limitations in moving new PDE-targeted therapies to the clinic.

Exploring the origin of phosphodiesterase inhibition via proteochemometric modeling

A proteochemometric study of a set of phosphodiesterase 4B and 4D inhibitors sheds light on the origin of their inhibition and selectivities.

Novel anti-cancer agents based on germacrone: design, synthesis, biological activity, docking studies and MD simulations

Germacrone is a major activity component found in Curcuma zedoaria oil product, which is extracted from Curcuma zedoaria.

Further Insights in the Binding Mode of Selective Inhibitors to Human PDE4D Enzyme Combining Docking and Molecular Dynamics

Alzheimer's disease has recently emerged as a possible field of application for PDE4D inhibitors (PDE4DIs). The great structure similarity among the various PDE4 isoforms and, furthermore, the lack of the full length crystal structure of the enzyme, impaired the rational design of new selective PDE4DIs. In this paper, with the aim of exploring new insights into the PDE4D binding, we tackled the problem by performing a computational study based on docking simulations combined with molecular dynamics (D-MD). Our work uniquely identified the binding mode and the key residues involved in the interaction with a number of in-house catechol iminoether derivatives, acting as PDE4DIs. Moreover, the new binding mode was tested using a series of analogues previously reported by us and it was used to confirm their key structural features to allow PDE4D inhibition. The binding model disclosed within the current computational study may prove to be useful to further advance the design and synthesis of novel, more potent and selective, PDE4D inhibitors.

New insights into selective PDE4D inhibitors: 3-(Cyclopentyloxy)-4-methoxybenzaldehyde O-(2-(2,6-dimethylmorpholino)-2-oxoethyl) oxime (GEBR-7b) structural development and promising activities to restore memory impairment

Phosphodiesterase type 4D (PDE4D) has been indicated as a promising target for treating neurodegenerative pathologies such as Alzheimer's Disease (AD). By preventing cAMP hydrolysis, PDE4 inhibitors (PDE4Is) increase the cAMP response element-binding protein (CREB) phosphorylation, synaptic plasticity and long-term memory formation. Pharmacological and behavioral studies on our hit GEBR-7b demonstrated that selective PDE4DIs could improve memory without causing emesis and sedation. The hit development led to new molecule series, herein reported, characterized by a catechol structure bonded to five member heterocycles. Molecular modeling studies highlighted the pivotal role of a polar alkyl chain in conferring selective enzyme interaction. Compound 8a showed PDE4D3 selective inhibition and was able to increase intracellular cAMP levels in neuronal cells, as well as in the hippocampus of freely moving rats. Furthermore, 8a was able to readily cross the blood-brain barrier and enhanced memory performance in mice without causing any emetic-like behavior. These data support the view that PDE4D is an adequate molecular target to restore memory deficits in different neuropathologies, including AD, and also indicate compound 8a as a promising candidate for further preclinical development.

Discover natural compounds as potential phosphodiesterase-4B inhibitors via computational approaches

cAMP, intracellular cyclic adenosine monophosphate, is a ubiquitous second messenger that plays a key role in many physiological processes. PDE4B which can reduce the cAMP level by hydrolyzing cAMP to 5'-AMP has become a therapeutic target for the treatment of human diseases such as respiratory disorders, inflammation diseases, neurological and psychiatric disorders. However, the use of currently available PDE4B inhibitors is restricted due to serious side effects caused by targeting PDE4D. Hence, we are attempting to find out subfamily-selective PDE4B inhibitors from natural products, using computer-aided approaches such as virtual screening, docking, and molecular dynamics simulation. Finally, four potential PDE4B-selective inhibitors (ZINC67912770, ZINC67912780, ZINC72320169, and ZINC28882432) were found. Compared to the reference drug (roflumilast), they scored better during the virtual screening process. Binding free energy for them was -317.51, -239.44, -215.52, and -165.77 kJ/mol, better than -129.05 kJ/mol of roflumilast. The pharmacophore model of the four candidate inhibitors comprised six features, including one hydrogen bond donor, four hydrogen bond acceptors, and one aromatic ring feature. It is expected that our study will pave the way for the design of potent PDE4B-selective inhibitors of new drugs to treat a wide variety of diseases such as asthma, COPD, psoriasis, depression, etc.

Discovery and modelling studies of natural ingredients from Gaultheria yunnanensis (FRANCH.) against phosphodiesterase-4

Phosphodiesterase-4 (PDE4) is an anti-inflammatory target for treatment of asthma and chronic obstructive pulmonary disease (COPD). Here, we report the isolation and characterization of 13 compounds (G1-G13) by bioassay-guided fractionation of the ethyl acetate extraction of Gaultheria yunnanensis (FRANCH.), one of which pentacyclic triterpene (G1) has never been reported. Four of them (G1, G2, G4, and G5) inhibit PDE4 with the IC50 values < 20 μM and G1 is the most potent ingredient with an IC50 of 245 nM and moderate selectivity over other PDE families. Molecular dynamics simulations suggest that G1 forms a hydrogen bond with Asn362, in addition to the hydrogen bond with Gln369 and π-π interactions with Phe372, which are commonly observed in the binding of most PDE4 inhibitors. The calculated binding free energies for the interactions of PDE4-G1 and PDE4-G2 are -19.4 and -18.8 kcal/mol, in consistence with the bioassay that G1 and G2 have IC50 of 245 nM and 542 nM, respectively. The modelling results of these active compounds may aid the rational design of novel PDE4 inhibitors as anti-inflammatory agents.

Understanding the effects on constitutive activation and drug binding of a D130N mutation in the β2 adrenergic receptor via molecular dynamics simulation

G-protein-coupled receptors (GPCRs) are currently one of the largest families of drug targets. The constitutive activation induced by mutation of key GPCR residues is associated closely with various diseases. However, the structural basis underlying such activation and its role in drug binding has remained unclear. Herein, we used all-atom molecular dynamics simulations and free energy calculations to study the effects of a D130N mutation on the structure of β2 adrenergic receptor (β2AR) and its binding of the agonist salbutamol. The results indicate that the mutation caused significant changes in some key helices. In particular, the mutation leads to the departure of transmembrane 3 (TM3) from transmembrane 6 (TM6) and marked changes in the NPxxY region as well as the complete disruption of a key ionic lock, all of which contribute to the observed constitutive activation. In addition, the D130N mutation weakens some important H-bonds, leading to structural changes in these regions. Binding free energy calculations indicate that van der Waals and electrostatic interactions are the main driving forces in binding salbutamol; however, binding strength in the mutant β2AR is significantly enhanced mainly through modifying electrostatic interactions. Further analysis revealed that the increase in binding energy upon mutation stems mainly from the H-bonds formed between the hydroxyl group of salbutamol and the serine residues of TM5. This observation suggests that modifications of the H-bond groups of this drug could significantly influence drug efficacy in the treatment of diseases associated with this mutation.

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.

Prenylated coumarins: natural phosphodiesterase-4 inhibitors from Toddalia asiatica

Bioassay-guided fractionation of the ethanolic extract of the roots of Toddalia asiatica led to the isolation of seven new prenylated coumarins (1-7) and 14 known analogues (8-21). The structures of 1-7 were elucidated by spectroscopic analysis, and their absolute configurations were determined by combined chemical methods and chiral separation analysis. Compounds 1-5, named toddalin A, 3‴-O-demethyltoddalin A, and toddalins B-D, represent an unusual group of phenylpropenoic acid-coupled prenylated coumarins. Compounds 1-21 and four modified analogues, 10a, 11a, 13a, and 17a, were screened by using tritium-labeled adenosine 3',5'-cyclic monophosphate ([3H]-cAMP) as substrate for their inhibitory activity against phosphodiesterase-4 (PDE4), which is a drug target for the treatment of asthma and chronic obstructive pulmonary disease. Compounds 3, 8, 10, 10a, 11, 11a, 12, 13, 17, and 21 exhibited inhibition with IC50 values less than 10 μM. Toddacoumalone (8), the most active compound (IC50=0.14 μM), was more active than the positive control, rolipram (IC50=0.59 μM). In addition, the structure-activity relationship and possible inhibitory mechanism of the active compounds are also discussed.

Natural phosphodiesterase-4 (PDE4) inhibitors from Crotalaria ferruginea

Bioassay-guided fractionation of the ethanol extract of the Chinese folk medicine Crotalaria ferruginea led to the isolation of a new isoflavonoid, 4'-hydroxy-2'-methylalpinum-isoflavone (1), and eight known analogs (2-9). Their structures were elucidated by spectroscopic analysis. Compounds 1, 2, 5, and 8 showed inhibitory activities against phosphodiesterase-4 (PDE4), a therapeutic target of asthma, with IC50 values ranging from 2.57 to 8.94 μM. The possible action mechanism and the structure-activity relationship (SAR) of the active isoflavonoids were explored by using molecular docking and molecular dynamics (MD) simulation methods. Our study herein may explain the anti-inflammatory function of this plant in Chinese folk medicine.

Six new tetraprenylated alkaloids from the South China Sea gorgonian Echinogorgia pseudossapo

Six new tetraprenylated alkaloids, designated as malonganenones L–Q (1–6), were isolated from the gorgonian Echinogorgia pseudossapo, collected in Daya Bay of Guangdong Province, China. The structures of 1–6 featuring a methyl group at N-3 and a tetraprenyl chain at N-7 in the hypoxanthine core were established by extensive spectroscopic analyses. Compounds 1–6 were tested for their inhibitory activity against the phosphodiesterases (PDEs)-4D, 5A, and 9A, and compounds 1 and 6 exhibited moderate inhibitory activity against PDE4D with IC₅₀ values of 8.5 and 20.3 µM, respectively.

The molecular basis for the inhibition of phosphodiesterase-4D by three natural resveratrol analogs. Isolation, molecular docking, molecular dynamics simulations, binding free energy, and bioassay

The phosphodiesterase-4 (PDE4) enzyme is a promising therapeutic target for several diseases. Our previous studies found resveratrol and moracin M to be natural PDE4 inhibitors. In the present study, three natural resveratrol analogs [pterostilbene, (E)-2',3,5',5-tetrahydroxystilbene (THSB), and oxyresveratrol] are structurally related to resveratrol and moracin M, but their inhibition and mechanism against PDE4 are still unclear. A combined method consisting of molecular docking, molecular dynamics (MD) simulations, binding free energy, and bioassay was performed to better understand their inhibitory mechanism. The binding pattern of pterostilbene demonstrates that it involves hydrophobic/aromatic interactions with Phe340 and Phe372, and forms hydrogen bond(s) with His160 and Gln369 in the active site pocket. The present work also reveals that oxyresveratrol and THSB can bind to PDE4D and exhibits less negative predicted binding free energies than pterostilbene, which was qualitatively validated by bioassay (IC50=96.6, 36.1, and 27.0μM, respectively). Additionally, a linear correlation (R(2)=0.953) is achieved for five PDE4D/ligand complexes between the predicted binding free energies and the experimental counterparts approximately estimated from their IC50 values (≈RT ln IC50). Our results imply that hydrophobic/aromatic forces are the primary factors in explaining the mechanism of inhibition by the three products. Results of the study help to understand the inhibitory mechanism of the three natural products, and thus help the discovery of novel PDE4 inhibitors from resveratrol, moracin M, and other natural products.