Comparative molecular field analysis (CoMFA) of phthalazine derivatives as phosphodiesterase IV inhibitors

The pyridazine heterocycle in molecular recognition and drug discovery

The pyridazine ring is endowed with unique physicochemical properties, characterized by weak basicity, a high dipole moment that subtends π-π stacking interactions and robust, dual hydrogen-bonding capacity that can be of importance in drug-target interactions. These properties contribute to unique applications in molecular recognition while the inherent polarity, low cytochrome P450 inhibitory effects and potential to reduce interaction of a molecule with the cardiac hERG potassium channel add additional value in drug discovery and development. The recent approvals of the gonadotropin-releasing hormone receptor antagonist relugolix (24) and the allosteric tyrosine kinase 2 inhibitor deucravacitinib (25) represent the first examples of FDA-approved drugs that incorporate a pyridazine ring. In this review, the properties of the pyridazine ring are summarized in comparison to the other azines and its potential in drug discovery is illustrated through vignettes that explore applications that take advantage of the inherent physicochemical properties as an approach to solving challenges associated with candidate optimization.

Graphical Abstract

α-Aminophosphonates as novel anti-leishmanial chemotypes: synthesis, biological evaluation, and CoMFA studies

α-Aminophosphonates have been identified as novel anti-leishmanial chemotypes against theL. donovanipromastigote with low μM range activity.

Synthesis and biological evaluation of 4-arylphthalazones bearing benzenesulfonamide as anti-inflammatory and anti-cancer agents

Nine 4-arylphthalazones bearing benzenesulfonamide (2a-i) were synthesized by the condensation of the appropriate 2-aroylbenzoic acid (1a-i) and 4-hydrazinobenzenesulfonamide in ethanol. The structures of these compounds were elucidated by elemental analysis, IR, ¹H NMR, ¹³C NMR, and MS spectroscopy. Two compounds, 2b and 2i, showed significant anti-inflammatory activity comparable to that of the standard drug celecoxib in the carrageenan-induced rat paw edema model. These compounds (2b and 2i) had selective inhibitory activity towards the COX-2 enzyme. Compound 2b had a better selectivity ratio (COX-1/COX-2) compared to that of celecoxib and can be used as a novel template for the design of selective COX-2 inhibitors. Compounds 2d and 2i were screened for their antiproliferative activity toward 60 human cancer cell lines by the National Cancer Institute (USA). The compounds 2d and 2i displayed mild activity toward the renal cancer cell line UO-31.

Novel thieno[2,3-d]pyrimidines: their design, synthesis, crystal structure analysis and pharmacological evaluation

Novel thieno[2,3-d]pyrimidines containing a cyclohexane ring fused with a six- or five-membered heterocyclic moiety along with a benzylic nitrile were designed as potential inhibitors of PDE4. Expeditious synthesis of these compounds was carried out via a multi-step sequence consisting of a few key steps such as Gewald reaction, Dieckmann type cyclisation and Krapcho decarboxylation. This newly developed strategy involved construction of the thienopyrimidine ring followed by the cyclohexanone moiety and subsequently the fused heterocyclic ring. A number of thieno[2,3-d]pyrimidine based derivatives were synthesized using this method some of which showed promising PDE4B inhibitory properties. One of them was tested for PDE4D inhibition in vitro and dose dependent inhibition of TNF-α. A few selected molecules were docked into the PE4B protein the results of which showed good overall correlations to their observed PDE4B inhibitory properties in vitro. The crystal structure analysis of representative compounds along with hydrogen bonding patterns and molecular arrangement present within the molecule is described.

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) of some benzimidazole derivatives with trichomonicidal activity

Trichomonosis is a common sexually transmitted infectious disease linked to reproductive health complications. Recently, the benzimidazole nucleus has emerged as a promising scaffold to develop new trichomonicidal agents. Despite the fact that large amounts of experimental data have been accumulated over the past eight years, no quantitative studies have yet been reported on this class of compounds. In our effort to develop new antiparasitic benzimidazole derivatives, we report in this paper CoMFA and CoMSIA studies with an initial set of 70 benzimidazole derivatives with trichomonicidal activity. Four CoMFA models and eight CoMSIA models were generated; ten of these models had values of r(2) > 0.6 and q(2) > 0.5. The best CoMFA model had r(2) = 0.936 and q(2) = 0.634, and the best CoMSIA model had r(2) = 0.858 and q(2) = 0.642. These models were generated by using two conformer selection methodologies (minimum energy conformations and 3D similarity), and three charge types (Mulliken, Gasteiger-Hükel and electrostatic potential atomic charges). The putative active tautomers of 1H-benzimidazole derivatives were selected using 3D-QSAR calculations. All models were validated via an external test set with 13 molecules. The best models satisfied additional validation criteria. The contour maps generated show the most important features that a benzimidazole derivative should have for trichomonicidal activity; they also, suggest that substituents at the 2- and 6-positions are important in the generation of derivatives with strong activity.

Development of improved models for phosphodiesterase-4 inhibitors with a multi-conformational structure-based QSAR method

Phosphodiesterase-4 (PDE-4) is an important drug target for several diseases, including COPD (chronic obstructive pulmonary disorder) and neurodegenerative diseases. In this paper, we describe the development of improved QSAR (quantitative structure-activity relationship) models using a novel multi-conformational structure-based pharmacophore key (MC-SBPPK) method. Similar to our previous work, this method calculates molecular descriptors based on the matching of a molecule's pharmacophore features with those of the target binding pocket. Therefore, these descriptors are PDE4-specific, and most relevant to the problem under study. Furthermore, this work expands our previous SBPPK QSAR method by explicitly including multiple conformations of the PDE-4 inhibitors in the regression analysis, and thus addresses the issue of molecular flexibility. The nonlinear regression problem resulted from including multiple conformations has been transformed into a linear equation and solved by an iterative partial least square (iPLS) procedure, according to the Lukacova-Balaz scheme. 35 PDE-4 inhibitors have been analyzed with this new method, and predictive models have been developed. Based on the prediction statistics for both the training set and the test set, these new models are more robust and predictive than those obtained by traditional ligand-based QSAR techniques as well as that obtained with the SBPPK method reported in our previous work. As a result, multiple predictive models have been added to the collection of QSAR models for PDE4 inhibitors. Collectively, these models will be useful for the discovery of new drug candidates targeting the PDE-4 enzyme.

A new structure-based QSAR method affords both descriptive and predictive models for phosphodiesterase-4 inhibitors

We describe the application of a new QSAR (quantitative structure-activity relationship) formalism to the analysis and modeling of PDE-4 inhibitors. This new method takes advantage of the X-ray structural information of the PDE-4 enzyme to characterize the small molecule inhibitors. It calculates molecular descriptors based on the matching of their pharmacophore feature pairs with those (the reference) of the target binding pocket. Since the reference is derived from the X-ray crystal structures of the target under study, these descriptors are target-specific and easy to interpret. We have analyzed 35 indole derivative-based PDE-4 inhibitors where Partial Least Square (PLS) analysis has been employed to obtain the predictive models. Compared to traditional QSAR methods such as CoMFA and CoMSIA, our models are more robust and predictive measured by statistics for both the training and test sets of molecules. Our method can also identify critical pharmacophore features that are responsible for the inhibitory potency of the small molecules. Thus, this structure-based QSAR method affords both descriptive and predictive models for phosphodiesterase-4 inhibitors. The success of this study has also laid a solid foundation for systematic QSAR modeling of the PDE family of enzymes, which will ultimately contribute to chemical genomics research and drug discovery targeting the PDE enzymes.

Phthalazine PDE IV inhibitors: Conformational study of some 6-methoxy-1,4-disubstituted derivatives

This report describes the detailed conformational analysis and synthesis of a series of phthalazine phosphodiesterase-type (IV) (PDE IV) inhibitors bearing either mono- or dichloro 4-methylenepyridine at the 'down' position of the phthalazine nucleus or different heterocycles at the 'top' position 4 of the phthalazine moiety. Both the mono- and dichloro 4-methylenepyridine units linked at carbon C1 of the phthalazine nucleus show identical conformational behaviour with the substituent preferentially oriented towards the external part of the molecule and the pyridine plane almost orthogonal to that of phthalazine ring. The heterocyclic five-membered rings linked at carbon C4 of phthalazine show quite different conformational behaviour. The 1,3-thiazole ring exists in a well-defined conformation almost coplanar with respect to the phthalazine nucleus while the 1,2,4-triazole and the 1,3-diazole heterocycles show a great conformational freedom with large torsion angles. Compound 3 bearing the thiazole ring at C4 displays the major biological activity, thus suggesting that a planar and rather rigid conformation of the pentacycle should favour the PDE IV inhibition capacity of this class of compounds.

Three-dimensional quantitative structure (3-D QSAR) activity relationship studies on imidazolyl and N-pyrrolyl heptenoates as 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR) inhibitors by comparative molecular similarity indices analysis (CoMSIA)

A comparative molecular similarity indices analysis (CoMSIA) of a set of 29 imidazolyl and N-pyrrolyl heptenoates have been performed to find out the structural requirements for 3-hydroxy-3-methylglutaryl-CoA reductase (HMGR) inhibitory activity. The HMG like side chain, a common moiety of statins, was used to align the molecules. The results guide to design new chemical entities with high potency.