Inherent efficacies of pyrazole-based derivatives for cancer therapy: the interface between experiment and in silico.

There has been an increasing trend in the design of novel pyrazole derivatives for desired biological applications. For a cost-effective strategy, scientists have implemented various computational drug design tools to go hand in hand with experiments for the design and discovery of potentially effective pyrazole-based therapeutics. This review highlights the milestones of pyrazole-containing inhibitors and the use of molecular modeling techniques in conjunction with experimental studies to provide a view of the binding mechanism of these compounds. The review focuses on the established targets that play a key role in cancer therapy, including proteins involved in tubulin polymerization, carbonic anhydrase and tyrosine kinase. Overall, using both experimental and computational methods in drug design represents a promising approach to cancer therapy.

Piecing the Fragments Together: Dynamical Insights into the Enhancement of BRD4-BD1 (BET Protein) Druggability in Cancer Chemotherapy Using Novel 8-methyl-pyrrolo[1,2-a]pyrazin-1(2H)-one Derivatives.

BACKGROUND: Fragment-based drug discovery in recent times has been explored in the design of highly potent therapeutics. METHODS: In this study, we explored the inhibitory dynamics of Compound 38 (Cpd38), a newly synthesized Bromodomain-containing protein 4 bromodomain 1 (BRD4-BD1) protein inhibitor derived from the synthetic coupling of Fragment 47 (Fgt47) into ABBV-075 scaffold. Using dynamic simulation methods, we unraveled the augmentative effects of chemical fragmentation on improved BRD4- BD1 inhibition. RESULTS: Findings from this study revealed that although Fgt47 exhibited a considerable ΔGbind, its incorporation into the difluoro-phenoxy pyridine scaffold (Cpd38) notably enhanced the binding affinity. Time-based analyses of interaction dynamics further revealed that the bulkiness of Cpd38 favored its interaction at the BRD4-BD1 active site relative to the fragment. Strikingly, compared to Fgt47, Cpd38 demonstrated high mobility, which could have enabled it to bind optimally and complementarily with key residues of the active site such as Ile146, Asn140, Cys136, Tyr98, Leu94, Val87, Phe83, and Trp81. DISCUSSION: On the contrary, the majority of these interactions were gradually lost in Fgt47, which could further indicate the essence of coupling it with the difluoro-phenoxy pyridine scaffold. Furthermore, Cpd38 had a more altering effect on BRD4-BDI relative to Fgt47, which could also be a result of its higher inhibitory activity. CONCLUSION: Conclusively, the design of highly potent therapeutics could be facilitated by the incorporation of pharmacologically active small molecule fragments into the scaffold of existing drugs.

Probing Protein-protein Interactions and Druggable Site Identification: Mechanistic Binding Events Between Ubiquitin and Zinc Finger with UFM1-specific Peptidase Domain Protein (ZUFSP).

BACKGROUND: Deubiquitinating enzymes (DUBs) protein family have been implicated in some deregulated pathways involved in carcinogeneses, such as cell cycle, gene expression, and DNA damage response (DDR). Zinc finger with UFM1-specific peptidase domain protein (ZUFSP) is one of the recently discovered members of the DUBs. OBJECTIVES: To identify and cross-validate the ZUFSP binding site using the bioinformatic tools, including SiteMap&Metapocket, respectively. To understand the molecular basis of complementary ZUFSP-Ub interaction and associated structural events using MD Simulation. METHODS: In this study, four binding pockets were predicted, characterized, and cross-validated based on physiochemical features such as site score, druggability score, site volume, and site size. Also, a molecular dynamics simulation technique was employed to determine the impact of ubiquitin-binding on ZUFSP. RESULTS: Site 1 with a site score 1.065, Size 102, D scores 1.00, and size volume 261 was predicted to be the most druggable site. Structural studies revealed that upon ubiquitin-binding, the motional movement of ZUFSP was reduced when compared to the unbound ZUFSP. Also, the ZUFSP helical arm (ZHA) domain orient in such a way that it moves closer to the Ub; this orientation enables the formation of a UBD which is very peculiar to ZUFSP. CONCLUSION: The impact of ubiquitin on ZUFSP movement and the characterization of its predicted druggable site can be targeted in the development of therapeutics.

Prospecting the therapeutic edge of a novel compound (B12) over berberine in the selective targeting of Retinoid X Receptor in colon cancer.

The Retinoid X Receptor (RXR) is an attractive target in the treatment of colon cancer. Different therapeutic binders with high potency have been used to specifically target RXR. Among these compounds is a novel analogue of berberine, B12. We provided structural and molecular insights into the therapeutic activity properties of B12 relative to its parent compound, berberine, using force field estimations and thermodynamic calculations. Upon binding of B12 to RXR, the high instability elicited by RXR was markedly reduced; similar observation was seen in the berberine-bound RXR. However, our analysis revealed that B12 could have a more stabilizing effect on RXR when compared to berberine. Interestingly, the mechanistic behaviour of B12 in the active site of RXR opposed its impact on RXR protein. This disparity could be due to the bond formation and breaking elicited between B12/berberine and the active site residues. We observed that B12 and berberine could induce a disparate conformational change in regions Gly250-Asp258 located on the His-RXRα/LBD domain. Comparatively, the high agonistic and activation potential reported for B12 compared to berberine might be due to its superior binding affinity as evidenced in the thermodynamic estimations. The total affinity for B12 (-25.76 kcal/mol) was contributed by electrostatic interactions from Glu243 and Glu239. Also, Arg371, which plays a crucial role in the activity of RXR, formed a strong hydrogen bond with B12; however, a weak interaction was elicited between Arg371 and berberine. Taken together, our study has shown the RXRα activating potential of B12, and findings from this study could provide a framework in the future design of RXRα binders specifically tailored in the selective treatment of colon cancer.

East to West not North-West: Structure-Based Mechanistic Resolution of 8-Hydroxyl Replacement and Resulting Effects on the Activities of Imidazole-Based Heme Oxygenase-1 Inhibitors.

Upregulation of Heme Oxygenase-1 (HO-1) has been widely implicated in cancer growth and chemoresistance. This explains the numerous drug discovery efforts aimed at mitigating its pro-carcinogenic roles till date. In a recent study, two selective azole-based HO-1 inhibitors (Cpd1 and Cpd2) were synthesized, which exhibited differential inhibitory potencies of ~200µm. Interestingly, variations in the affinities of these compounds were determined by their positioning across specific regions of the HO-1 binding domain, pin-pointing a pharmacological interrelationship that remains unresolved. Therefore, in this study, using molecular dynamics simulations and binding free energy calculations, we investigate how dynamical orientations of these compounds influence their binding affinities at the active HO-1 domain. Findings revealed favorable binding for the bromobenzene and imidazole substituents of Cpd1 at the western and eastern regions of the HO-1 active domain. The constituent hydroxyl group was coordinated by residues Asp140 and Arg136 over the simulation period. On the contrary, stable binding of the bromobenzene and imidazole substituents were negated by the optimal orientations of the benzyl substituent, which extended into the northeastern region. These were supported by the displacement of Asp140 and Arg136, crucial for hydrogen bond formation in Cpd1. Also, we observed that Cpd2 exhibited high deviations indicative of an unstable binding relative to Cpd1. This further supports the presumption that Cpd2 was systematically oriented away from the active HO-1 region, a phenomenon that was due to the optimal motions of the benzyl group at the northeastern regions. The highlight of our findings is that the benzyl substituent in Cpd2 elicited negative effects on HO-1, vis a vis, instability, displacement of crucial residues, and low binding energy when compared to Cpd1. Findings pave the way for future drug discovery efforts related to HO-1 inhibition in cancer therapy.

Functional Analysis of Single Nucleotide Polymorphism in ZUFSP Protein and Implication in Pathogenesis.

Researches have revealed that functional non-synonymous Single Nucleotide Polymorphism (nsSNPs) present in the Zinc-finger with UFM1-Specific Peptidase domain protein (ZUFSP) may be involved in genetic instability and carcinogenesis. For the first time, we employed in-silico approach using predictive tools to identify and validate potential nsSNPs that could be pathogenic. Our result revealed that 8 nsSNPs (rs 112738382, rs 140094037, rs 201652589, rs 201847265, rs 202076827, rs 373634906, rs 375114528, rs 772591104) are pathogenic after being subjected to rigorous filtering process. The structural impact of the nsSNPs on ZUFSP structure indicated that the nsSNPs affect the stability of the protein by lowering ZUFSP protein stability. Furthermore, conservation analysis showed that rs 201652589, rs 140094037, rs 201847265, and rs 772591104 were highly conserved. Interestingly, the protein-protein affinity between ZUFSP and Ubiquitin was altered rs 201652589, rs 140094037, rs 201847265, and rs 772591104 had a binding affinity of - 0.46, - 0.83, - 1.62, and - 1.12 kcal/mol respectively. Our study has been able to identify potential nsSNPs that could be used as genetic biomarkers for some diseases arising as a result of aberration in the ZUFSP structure, however, being a predictive study, the identified nsSNPs need to be experimentally investigated.

Update and Potential Opportunities in CBP [Cyclic Adenosine Monophosphate (cAMP) Response Element-Binding Protein (CREB)-Binding Protein] Research Using Computational Techniques.

CBP [cyclic adenosine monophosphate (cAMP) response element-binding protein (CREB)-binding protein] is one of the most researched proteins for its therapeutic function. Several studies have identified its vast functions and interactions with other transcription factors to initiate cellular signals of survival. In cancer and other diseases such as Alzheimer's, Rubinstein-taybi syndrome, and inflammatory diseases, CBP has been implicated and hence an attractive target in drug design and development. In this review, we explore the various computational techniques that have been used in CBP research, furthermore we identified computational gaps that could be explored to facilitate the development of highly therapeutic CBP inhibitors.

Leveraging on Active Site Similarities; Identification of Potential Inhibitors of Zinc-Finger and UFSP domain Protein (ZUFSP).

BACKGROUND: ZUFSP (Zinc-finger and UFSP domain protein) is a novel representative member of the recently characterized seventh class of deubiquitinating enzymes (DUBs). Due to the roles DUBs play in genetic instability, they have become a major drug target in cancer and neurodegenerative diseases. ZUFSP, being a DUB enzyme has also been implicated in genetic stability. However, no lead compound has been developed to target ZUFSP. OBJECTIVE/METHODS: Therefore, in this study, we used a combined drug repurposing, virtual screening and per-Residue Energy Decomposition (PRED) to identify ZUFSP inhibitors with therapeutic potential. 3-bromo-6-{[4-hydroxy-1-3(3-phenylbutanoyl)piperidin-4-yl]methyl}-4H,5H,6H,7H-thieno[2,3- C]pyridine-7-one (BHPTP) which is an inhibitor of USP7 was repurposed to target ZUFSP. The rationale behind this is based on the similarity of the active between USP7 and ZUFSP. RESULTS: PRED of the binding between BHPTP and ZUFSP revealed Cys223, Arg408, Met410, Asn460, and Tyr465 as the crucial residues responsible for this interaction. The pharmacophoric moieties of BHPTP responsible for this binding along with other physiochemical properties were used as a filter to retrieve potential ligands. 799 compounds were retrieved, ZINC083241427, ZINC063648749, and ZINC063648753 were selected due to the binding energy they exhibited. Cheminformatics analysis revealed that the compounds possess high membrane permeability, however, BHPTP had a low membrane permeability. Furthermore, the compounds are drug like, having obeyed Lipinski's rule of five. CONCLUSION: Taken together, findings from this study put ZINC083241427, ZINC063648749, and ZINC063648753 as potential ZUFSP inhibitor, however, more experimental validation is required to unravel the mechanism of actions of these compounds.

Therapeutic Path to Double Knockout: Investigating the Selective Dual-Inhibitory Mechanisms of Adenosine Receptors A1 and A2 by a Novel Methoxy-Substituted Benzofuran Derivative in the Treatment of Parkinson's Disease.

The dual inhibition of adenosine receptors A1 (A1 AR) and A2 (A2A AR) has been considered as an efficient strategy in the treatment of Parkinson's disease (PD). This led to the recent development of a series of methoxy-substituted benzofuran derivatives among which compound 3j exhibited dual-inhibitory potencies in the micromolar range. Therefore, in this study, we seek to resolve the mechanisms by which this novel compound elicits its selective dual targeting against A1 AR and A2A AR. Unique to the binding of 3j in both proteins, from our findings, is the ring-ring interaction elicited by A1Phe275 (→ A2Phe170) with the benzofuran ring of the compound. As observed, this π-stacking interaction contributes notably to the stability of 3j at the active sites of A1 and A2A AR. Besides, conserved active site residues in the proteins such as A1Ala170 (→ A2Ala65), A1Ile173 (→ A2Ile68), A1Val191 (→ A2Val86), A1Leu192 (→ A2Leu87), A1Ala195 (→ A2Ala90), A1Met284 (→ A2Met179), A1Tyr375 (→ A2Tyr369), A1Ile378 (→ A2Ile372), and A1His382 (→ A2His376) were commonly involved with other ring substituents which further complement the dual binding and stability of 3j. This reflects a similar interaction mechanism that involved aromatic (π) interactions. Consequentially, vdW energies contributed immensely to the dual binding of the compound, which culminated in high ΔGbinds that were homogenous in both proteins. Furthermore, 3j commonly disrupted the stable and compact conformation of A1 and A2A AR, coupled with their active sites where Cα deviations were relatively high. Ligand mobility analysis also revealed that both compounds exhibited a similar motion pattern at the active site of the proteins relative to their optimal dual binding. We believe that findings from this study with significantly aid the structure-based design of highly selective dual-inhibitors of A1 and A2A AR.

Investigating the Mechanistic Inhibitory Discrepancies of Novel Halogen and Alkyl Di-Substituted Oxadiazole-Based Dibenzo-Azepine-Dione Derivatives on Poly (ADP-Ribose) Polymerase-1.

Numerous studies have established the involvement of Poly (ADP-ribose) Polymerase-1 (PARP-1) in cancer presenting it as an important therapeutic target over recent years. Although homology among the PARP protein family makes selective targeting difficult, two compounds [d11 (0.939 µM) and d21 (0.047 µM)] with disparate inhibitory potencies against PARP-1 were recently identified. In this study, free energy calculations and molecular simulations were used to decipher underlying mechanisms of differential PARP-1 inhibition exhibited by the two compounds. The thermodynamics calculation revealed that compound d21 had a relatively higher ΔGbind than d11. High involvement of van der Waal and electrostatic effects potentiated the affinity of d21 at PARP-1 active site. More so, incorporated methyl moiety in d11 accounted for steric hindrance which, in turn, prevented complementary interactions of key site residues such as TYR889, MET890, TYR896, TYR907. Conformational studies also revealed that d21 is more stabilized for interactions in the active site compared to d11. We believe that findings from this study would provide an important avenue for the development of selective PARP-1 inhibitors.

Identification and classification of differentially expressed genes reveal potential molecular signature associated with SARS-CoV-2 infection in lung adenocarcinomal cells.

Genomic techniques such as next-generation sequencing and microarrays have facilitated the identification and classification of molecular signatures inherent in cells upon viral infection, for possible therapeutic targets. Therefore, in this study, we performed a differential gene expression analysis, pathway enrichment analysis, and gene ontology on RNAseq data obtained from SARS-CoV-2 infected A549 cells. Differential expression analysis revealed that 753 genes were up-regulated while 746 down-regulated. SNORA81, OAS2, SYCP2, LOC100506985, and SNORD35B are the top 5 upregulated genes upon SARS-Cov-2 infection. Expectedly, these genes have been implicated in the immune response to viral assaults. In the Ontology of protein classification, a high percentage of the genes are classified as Gene-specific transcriptional regulator, metabolite interconversion enzyme, and Protein modifying enzymes. Twenty pathways with P-value lower than 0.05 were enriched in the up-regulated genes while 18 pathways are enriched in the down-regulated DEGs. The toll-like receptor signalling pathway is one of the major pathways enriched. This pathway plays an important role in the innate immune system by identifying the pathogen-associated molecular signature emanating from various microorganisms. Taken together, our results present a novel understanding of genes and corresponding pathways upon SARS-Cov-2 infection, and could facilitate the identification of novel therapeutic targets and biomarkers in the treatment of COVID-19.

Delving into the Characteristic Features of "Menace" Mycobacterium tuberculosis Homologs: A Structural Dynamics and Proteomics Perspectives.

The global increase in the morbidity/mortality rate of Mycobacterial infections, predominantly renascent tuberculosis, leprosy, and Buruli ulcers have become worrisome over the years. More challenging is the incidence of resistance mediated by mutant Mycobacterium strains against front-line antitubercular drugs. Homologous to all Mycobacteria species is the GlcNAc-6-phosphate deacetylase (NagA) which catalyzes essential amino sugars synthesis required for cell wall architecture, hence, metamorphosing into an important pharmacological target for curtailing virulence and drug-resistance. This study used integrated bioinformatics methods, MD simulations, and DynaMut and PolyPhen2 to; explore unique features, monitor dynamics, and analyze the functional impact of non-synonymous single-nucleotide polymorphisms of the six NagA of most ruinous Mycobacterium species; tuberculosis (Mtb), smegmatis (MS), marinum (MM), ulcerans, africanum, and microti respectively. This approach is essential for multi-targeting and could result in the identification of potential polypharmacological antitubercular compounds. Comparative sequential analyses revealed ≤ 50% of the overall structure, including the catalytic Asp267 and reactive Cys131, remained conserved. Interestingly, MS-NagA and MM-NagA possess unique hydrophobic isoleucine (Ile) residues at their active sites in contrast to leucine (Leu) found in other variants. More so, unique to the active sites of the NagA is a 'subunit loop' that covers the active site; probably crucial in binding (entry and exit) mechanisms of targeted NagA inhibitors. Relatively, nsSNP mutations exerted a destabilizing effect on the native NagA conformation. Structural and dynamical insights provided, basically pin-pointed the "Achilles' heel" explorable for the rational drug design of target-specific 'NagA' inhibitors potent against a wide range of mycobacterial diseases.

Exploring the Role of Asp1116 in Selective Drug Targeting of CREBcAMP- Responsive Element-binding Protein Implicated in Prostate Cancer.

BACKGROUND: The selective targeting of CREB-cAMP-responsive element-binding protein (CBP) has recently evolved as a vital therapeutic approach for curtailing its aberrant upregulation associated with the development of prostate cancer. Inhibition of CBP has been discovered to be an important therapeutic option in androgen receptor signalling pathway mediated prostate cancer. Y08197, a novel selective inhibitor of CBP, has shown promising therapeutic outcome in prostate carcinogenesis over non-selective analogues such as CPI-637. METHODS/RESULTS: Herein, we used molecular dynamics simulation to gain insights into the mechanistic and selective targeting of Y08197 at the bromodomain active site. Molecular Mechanics/ Poisson-Boltzmann Surface Area (MM/PBSA) analysis revealed a similar inhibitory effect between Y08197 and CPI-637. Furthermore, in exploring the selective affinity of Y08197 towards CBP in combination with Bromodomain and PHD finger-containing protein 1(BRPF1), our findings highlighted Asp1116 as the 'culprit' residue responsible for this selective targeting. Upon binding, Asp1116 assumed a conformation that altered the architecture of the bromodomain active site, thereby orienting the helices around the active site in a more compacted position. In addition to some specific structural perturbations mediated by Asp1116 on the dynamics of CBP, our study revealed that the strong hydrogen bond interaction (N-H...O) elicited in CBP-Y08197 sequestered Y08197 tightly into the CBP bromodomain active site. CONCLUSION: Conclusively, the inhibition and selective pattern of Y08197 can be replicated in future structure-based CBP inhibitors and other bromodomain implicated in carcinogenesis.

From the Explored to the Unexplored: Computer-Tailored Drug Design Attempts in the Discovery of Selective Caspase Inhibitors.

The pathophysiological roles of caspases have made them attractive targets in the treatment and amelioration of neurologic diseases. In normal conditions, the expression of caspases is regulated in the brain, while at the onset of neurodegeneration, such as in Alzheimer's disease, they are typically overexpressed. Till date, several therapeutic efforts that include the use of small endogenous binders have been put forward to curtail dysfunctionalities that drive aberrant death in neuronal cells. Caspases are highly homologous, both in structure and in sequence, which leaves us with the question: is it possible to specifically and individually target caspases, while multiple therapeutic attempts to achieve selective targeting have failed! Based on antecedent events, the use of Computer-Aided Drug Design (CADD) methods has significantly contributed to the design of small molecule inhibitors, especially with selective target ability and reduced off-target therapeutic effects. Interestingly, we found out that there still exists an enormous room for the integration of structure/ligand-based drug design techniques towards the development of highly specific reversible and irreversible caspase inhibitors. Therefore, in this review, we highlight drug discovery approaches that have been directed towards caspase inhibition in addition to an insightful focus on applicable CADD techniques for achieving selective targeting in caspase research.

Drug promiscuity: Exploring the polypharmacology potential of 1, 3, 6-trisubstituted 1, 4-diazepane-7-ones as an inhibitor of the 'god father' of immune checkpoint.

High production cost, instability, low tumor penetration are some of the shortcomings that have characterized and undermined the use of antibodies as a target for Cytotoxic T-lymphocytes associated protein 4 (CTLA-4). Design and discovery of small molecule inhibitors have therefore become a sine qua non in targeting immune proteins implicated in immune disorders. In this study, we utilized a drug repositioning approach to explore the characteristic feature of unrelated proteins to have similar binding sites and the promiscuity of drugs to repurpose an existing drug to target CTLA-4. CTLA-4 and Kallikrein-7 were found to have similar binding sites, we therefore used 1, 3, 6-trisubstituted 1, 4-diazepane-7-ones (TDSO) which is an inhibitor of Kallikrein-7 as our lead compound. High throughput screening using TDSO as a lead compound resulted in 9 hits with ZINC04515726 and ZINC08985213 having the highest binding score. We went ahead to investigate the interaction of these compounds with CTLA-4 by conducting a molecular dynamic simulation. Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) estimations revealed that TDSO had the highest binding energy value of -28.51Kcal/mol, with ZINC04515726 and ZINC08985213 having -23.76Kcal/mol and -21.03Kcal/mol respectively. The per-residue decomposition highlighted Tyr24, Ala25, Gly28, Ala30, Tyr53 and Asn72 as having significantly high electrostatic energy contributions and the main contributing residues to the binding of TDSO, ZINC04515726 and ZINC08985213 to Cytotoxic T lymphocytes CTLA-4. Summarily, from the results gathered, we proposed that TDSO can be an effective immune check point small molecule inhibitor against the suppression of T-cell activation, proliferation, and tumor cell eradication.