In silicoinvestigation of PARP-1 catalytic domains inholoandapostates for the design of high-affinity PARP-1 inhibitors

Mathematical modeling of cardio-oncology: Modeling the systemic effects of cancer therapeutics on the cardiovascular system

Cardiotoxicity is a common side-effect of many cancer therapeutics; however, to-date there has been very little push to understand the mechanisms underlying this group of pathologies. This has led to the emergence of cardio-oncology, a field of medicine focused on understanding the effects of cancer and its treatment on the human heart. Here, we describe how mechanistic modeling approaches have been applied to study open questions in the cardiovascular system and how these approaches are being increasingly applied to advance knowledge of the underlying effects of cancer treatments on the human heart. A variety of mechanistic, mathematical modeling techniques have been applied to explore the link between common cancer treatments, such as chemotherapy, radiation, targeted therapy, and immunotherapy, and cardiotoxicity, nevertheless there is limited coverage in the different types of cardiac dysfunction that may be associated with these treatments. Moreover, cardiac modeling has a rich heritage of mathematical modeling and is well suited for the further development of novel approaches for understanding the cardiotoxicities associated with cancer therapeutics. There are many opportunities to combine mechanistic, bottom-up approaches with data-driven, top-down approaches to improve personalized, precision oncology to better understand, and ultimately mitigate, cardiac dysfunction in cancer patients.

In-silico method for elucidation of prodigiosin as PARP-1 inhibitor a prime target of Triple-negative breast cancer

Triple-Negative Breast Cancer (TNBC) is found to be one of the life-threatening cancer. Poly (ADP-Ribose) Polymerase-1 (PARP-1) is overexpressed by those tumour cells, which become resistant to chemotherapies. Inhibition of PARP-1 has a considerable effect on treating TNBC. Prodigiosin is a valuable pharmaceutical compound that exhibits anticancer properties. The present study aims to virtually evaluate prodigiosin as a potent PARP-1 inhibitor using Molecular docking and Molecular Dynamics (MD) simulation studies. The PASS (Prediction of Activity Spectra for Substances) prediction tool evaluated the biological properties of prodigiosin. Then the drug-likeness and pharmacokinetic properties of prodigiosin were determined using Swiss-ADME software. It was suggested that prodigiosin obeyed Lipinski's rule of five and thus could act as a drug with good pharmacokinetic properties. Moreover, molecular docking was done with AutoDock 4.2 to identify the critical amino acids of the protein-ligand complex. It was indicated that prodigiosin has a docking score of -8.08 kcal/mol, which showed its effective interaction with crucial amino acid, His201A of PARP-1 protein. Further, MD simulation was performed using Gromacs software to validate the stability of the prodigiosin-PARP-1 complex. Prodigiosin was found to have good structural stability and affinity at the active site of PARP-1 protein. Additionally, PCA and MM-PBSA were calculated for the prodigiosin-PARP-1 complex, which revealed that prodigiosin has an excellent binding affinity towards PARP-1 protein. Prodigiosin can possibly be used as oral drug due to its PARP-1 inhibition through high binding affinity, structural stability, and receptor flexibility towards crucial amino acid residue His201A of PARP-1 protein. In-addition, in-vitro cytotoxicity, and apoptosis analysis of prodigiosin-treated TNBC cell line-MDA-MB-231 revealed that prodigiosin exhibited significant anticancer activity in 101.1 µg/mL concentration, when compared to commercially available synthetic drug cisplatin. Thus, prodigiosin could act as a potential candidate for treatment of TNBC than the commercially available synthetic drugs.

How ligands regulate the binding of PARP1 with DNA: Deciphering the mechanism at the molecular level

The catalytic (CAT) domain is a key region of poly (ADP-ribose) polymerase 1 (PARP1), which has crucial interactions with inhibitors, DNA, and other domains of PARP1. To facilitate the development of potential inhibitors of PARP1, it is of great significance to clarify the differences in structural dynamics and key residues between CAT/inhibitors and DNA/PARP1/inhibitors through structure-based computational design. In this paper, conformational changes in PAPR1 and differences in key residue interactions induced by inhibitors were revealed at the molecular level by comparative molecular dynamics (MD) simulations and energy decomposition. On one hand, PARP1 inhibitors indirectly change some residues of the CAT domain which interact with DNA and other domains. Furthermore, the interaction between ligands and catalytic binding sites can be transferred to the DNA recognition domain of PARP1 by a strong negative correlation movement among multi-domains of PARP1. On the other hand, it is not reliable to use the binding energy of CAT/ligand as a measure of ligand activity, because it may in some cases differs greatly from the that of PARP1/DNA/ligand. For PARP1/DNA/ligand, the stronger the binding stability between the ligand and PARP1, the stronger the binding stability between PARP1 and DNA. The findings of this work can guide further novel inhibitor design and the structural modification of PARP1 through structure-based computational design.

Unveiling the mechanistic roles of chlorine substituted phthalazinone-based compounds containing chlorophenyl moiety towards the differential inhibition of poly (ADP-ribose) polymerase-1 in the treatment of lung cancer

PARP-1 has become an attractive target in cancer treatment owing to its significant role in breast and ovarian cancers. The design of highly selective and effective poly (ADP ribose) polymerase-1 inhibitors has significant therapeutic advantages and has remained the core of several PARP-1-based drug discovery research. The pharmacophoric relevance of a chlorine substituent in a recent study led to the design of compounds 11c (meta-chlorophenyl) and 11d (para-chlorophenyl). In this study, we resolved the mechanistic effects of the changes in chlorine positional orientation, which underlie the inhibitory potencies and selectivity exhibited disparately by 11c and 11d. Compared to 11d, among other multiple higher-affinity complementary interactions with key site residues, the meta-Cl positioning in 11c facilitated its optimal motion and orientation towards conserved residues Arg878 and Asp766 with consistent pi-cation and pi-anion interactions, respectively, thereby favoring the stability of the ligand towards PARP-1. These could account for the higher inhibitory potency exhibited by 11c relative to 11d against PARP-1. The thermodynamics calculation revealed that 11c had a relatively higher total binding energy (ΔGbind) than 11d. We also observed that 11d displayed high deviations, compared to 11c, indicative of its unstable binding orientation. Furthermore, we reported in this study that the high involvement of electrostatic and van der Waal effects potentiated the binding affinity and strength of 11c (ΔE_vdW = -50.58 and ΔE_ele = -27.20) relative to 11d (ΔE_vdW = -49.46 and ΔE_ele = -19.96) at PARP-1 binding pocket. We believe the findings in this current study would provide valuable insights into the design of selective PARP-1 inhibitors containing chlorine substituent for cancer treatment, including lung cancer.Communicated by Ramaswamy H. Sarma.

Deciphering the functional mechanism of zinc ions of PARP1 binding with single strand breaks and double strand breaks

Poly(ADP-ribose)polymerase 1 (PARP1) is a key target for the treatment of cancer-related diseases, and plays an important role in biological processes such as DNA repair, regulating a variety of metabolic and signal transduction processes. Understanding the dynamic binding mechanisms between each domain of PARP1 and DNA is of great significance to deepen the understanding on the function of PARP1 and to facilitate the design of inhibitors. Herein, strategies such as classical molecular dynamics simulation, conformational analysis, binding free energy calculation and energy decomposition were used to shed light on the binding mechanisms of different DNA binding domains (DBDs, including ZnF1, ZnF2 and ZnF3) in PARP1 with DNA and on the influences of zinc ions on the binding process. On one hand, during binding with DNA, ZnF2 tends to expand its space to identify the DNA damage sites and ZnF1/ZnF2 recognizes the interfaces on both sides of DNA damage rather than one side during the process of DNA repair. More importantly, the stable secondary structure of L₂ of ZnF2 (PRO146 to MET153) is the key conformational change for ZnF1 and ZnF2 to recognize DNA damage. Meanwhile, ZnF3 has little effect on the binding mechanisms of PARP1. On the other hand, for the structural differences of DBD domains, zinc ions in ZnF1 and ZnF2 (Zn1 and Zn2) have an impact not only on the conformational changes of PARP1, but also on the conformational changes brought by the interaction of double strand breaks (DSB) and single strand breaks (SSB). And meanwhile, Zn3 also has little effect on ZnF3 for the system of ZnF3/DSB. The findings presented in this work deepen the understanding on the functional mechanism of PARP1 and provide a theoretical basis for further study on the interaction between different inhibitors and DBD domains to design more potential inhibitors.

Poly(ADP-ribose)polymerase 1 (PARP1) is a key target for treatment of cancer-related diseases. Detailed structural changes DBD in PARP1 during the binding process with DNA were investigated and the dynamic conformational differences of DBD caused by zinc ions were revealed.

Biochemical, Structural Analysis, and Docking Studies of Spiropyrazoline Derivatives

In this study, we evaluated the antiproliferative potential, DNA damage, crystal structures, and docking calculation of two spiropyrazoline derivatives. The main focus of the research was to evaluate the antiproliferative potential of synthesized compounds towards eight cancer cell lines. Compound I demonstrated promising antiproliferative properties, especially toward the HL60 cell line, for which IC₅₀ was equal to 9.4 µM/L. The analysis of DNA damage by the comet assay showed that compound II caused DNA damage to tumor lineage cells to a greater extent than compound I. The level of damage to tumor cells of the HEC-1-A lineage was 23%. The determination of apoptotic and necrotic cell fractions by fluorescence microscopy indicated that cells treated with spiropyrazoline-based analogues were entering the early phase of programmed cell death. Compounds I and II depolarized the mitochondrial membranes of cancer cells. Furthermore, we performed simple docking calculations, which indicated that the obtained compounds are able to bind to the PARP1 active site, at least theoretically (the free energy of binding values for compound I and II were −9.7 and 8.7 kcal mol⁻¹, respectively). In silico studies of the influence of the studied compounds on PARP1 were confirmed in vitro with the use of eight cancer cell lines. The degradation of the PARP1 enzyme was observed, with compound I characterized by a higher protein degradation activity.

An Insight of RuBisCO Evolution through a Multilevel Approach

RuBisCO is the most abundant enzyme on earth; it regulates the organic carbon cycle in the biosphere. Studying its structural evolution will help to develop new strategies of genetic improvement in order to increase food production and mitigate CO₂ emissions. In the present work, we evaluate how the evolution of sequence and structure among isoforms I, II and III of RuBisCO defines their intrinsic flexibility and residue-residue interactions. To do this, we used a multilevel approach based on phylogenetic inferences, multiple sequence alignment, normal mode analysis, and molecular dynamics. Our results show that the three isoforms exhibit greater fluctuation in the loop between αB and βC, and also present a positive correlation with loop 6, an important region for enzymatic activity because it regulates RuBisCO conformational states. Likewise, an increase in the flexibility of the loop structure between αB and βC, as well as Lys330 (form II) and Lys322 (form III) of loop 6, is important to increase photosynthetic efficiency. Thus, the cross-correlation dynamics analysis showed changes in the direction of movement of the secondary structures in the three isoforms. Finally, key amino acid residues related to the flexibility of the RuBisCO structure were indicated, providing important information for its enzymatic engineering.

Discovery of New Protein Targets of BPA Analogs and Derivatives Associated with Noncommunicable Diseases: A Virtual High-Throughput Screening

BACKGROUND

Bisphenol A analogs and derivatives (BPs) have emerged as new contaminants with little or no information about their toxicity. These have been found in numerous everyday products, from thermal paper receipts to plastic containers, and measured in human samples.

OBJECTIVES

The objectives of this research were to identify in silico new protein targets of BPs associated with seven noncommunicable diseases (NCDs), and to study their protein-ligand interactions using computer-aided tools.

METHODS

Fifty BPs were identified by a literature search and submitted to a virtual high-throughput screening (vHTS) with 328 proteins associated with NCDs. Protein-protein interactions between predicted targets were examined using STRING, and the protocol was validated in terms of binding site recognition and correlation between in silico affinities and in vitro data.

RESULTS

According to the vHTS, several BPs may target proteins associated with NCDs, some of them with stronger affinities than bisphenol A (BPA). The best affinity score (the highest in silico affinity absolute value) was obtained after docking 4,4'-bis(N-carbamoyl-4-methylbenzensulfonamide)diphenylmethane (BTUM) on estradiol 17-beta-dehydrogenase 1 (-13.7 kcal/mol). However, other molecules, such as bisphenol A bis(diphenyl phosphate) (BDP), bisphenol PH (BPPH), and Pergafast 201 also exhibited great affinities (top 10 affinity scores for each disease) with proteins related to NCDs.

DISCUSSION

Molecules such as BTUM, BDP, BPPH, and Pergafast 201 could be targeting key signaling pathways related to NCDs. These BPs should be prioritized for in vitro and in vivo toxicity testing and to further assess their possible role in the development of these diseases. https://doi.org/10.1289/EHP7466.

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 ΔG_bind 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.

Identifying new piperazine-based PARP1 inhibitors using text mining and integrated molecular modeling approaches

One of the important molecular targets for antitumor drug discovery is the polyadenosine diphosphate-ribose polymerase-1 (PARP1) enzyme. It is linked with various biological functions including DNA repair and apoptosis. It is primarily a nuclear enzyme linked to chromatin, which is activated by DNA damage. Improved expression of PARP1 in melanomas, breast cancer, lung cancer and other neoplastic diseases is often observed. A tremendous PARP research concerning cancer and ischemia is progressing very rapidly. There are currently four PARP1 inhibitors approved by the FDA on the market, namely Olaparib, Rucaparib, Niraparib and Talazoparib. All of these molecules are non-selective inhibitors of PARP1. Currently there is an urgent need for novel and selective PARP1 inhibitors. In this work, asmall molecule database (Specs SC) were used to identify the new selective lead inhibitors of PARP1. Piperazine scaffold is an important fragment that is used in many currently used FDA approved drugs in different diseases including PARP1 inhibitor Olaparib. Thus, based on text mining studies, 4674 compounds thatinclude piperazine fragments were identified and virtually screened at the binding pocket of target protein PARP1. Compounds that have high docking scores were used in molecular dynamics (MD) simulations. Free energy calculations were also performed to compare the predicted binding energies with known PARP1 inhibitors. The critical amino acid interactions of these newly identified hits in the binding pocket were also investigated in detail for better understanding of the structural features required for next generation PARP1 inhibitors. Thus, here together with combination of text-mining and integrated molecular modeling approaches, we identified novel piperazine-based hits against PARP1 enzyme.Communicated by Ramaswamy H. Sarma.

Insight into the binding of a synthetic nitro-flavone derivative with human poly (ADP-ribose) polymerase 1

Flavones are important bioactive compounds, many of which are effective in cancer therapy for their ability to target enzymes related to DNA repair and cell proliferation. In this report, the interaction of a synthetic nitroflavone, 2,4-nitrophenylchromen-4-one (4NCO) with human poly (ADP-ribose) polymerase 1 (hPARP1) was investigated to explore its inhibitory action. Its interaction with hPARP1 was compared with that of other inhibitors through molecular docking studies. Further insight into the 4NCO-hPARP1 interaction was obtained from competitive docking and molecular dynamic simulation studies. In silico mutagenesis studies and per-residue interaction energy calculations were carried out. Quantitative Structure Activity Relationship analysis was also performed to calculate its predictive percent inhibitory activity. Our results indicated that 4NCO exhibited competitive mode of binding to hPARP1. It formed a stable interaction with the protein thereby hindering any further molecular interaction to render it inactive with a predictive inhibition of 96%. It also had good ADMET properties and showed best Autodock binding free energy values compared to other known inhibitors. 4NCO showed good hPARP1 inhibitory properties with higher bioavailability and lower probability of getting effluxed. Development of inhibitors against hPARP1 is important for cell proliferative disorders, where 4NCO can be predicted as a potential new drug.

Target-based drug discovery through inversion of quantitative structure-drug-property relationships and molecular simulation: CA IX-sulphonamide complexes

In this work, a target-based drug screening method is proposed exploiting the synergy effect of ligand-based and structure-based computer-assisted drug design. The new method provides great flexibility in drug design and drug candidates with considerably lower risk in an efficient manner. As a model system, 45 sulphonamides (33 training, 12 testing ligands) in complex with carbonic anhydrase IX were used for development of quantitative structure-activity-lipophilicity (property)-relationships (QSPRs). For each ligand, nearly 5,000 molecular descriptors were calculated, while lipophilicity (logk_w) and inhibitory activity (logK_i) were used as drug properties. Genetic algorithm-partial least squares (GA-PLS) provided a QSPR model with high prediction capability employing only seven molecular descriptors. As a proof-of-concept, optimal drug structure was obtained by inverting the model with respect to reference drug properties. 3509 ligands were ranked accordingly. Top 10 ligands were further validated through molecular docking. Large-scale MD simulations were performed to test the stability of structures of selected ligands obtained through docking complemented with biophysical experiments.

Examination of Diazaspiro Cores as Piperazine Bioisosteres in the Olaparib Framework Shows Reduced DNA Damage and Cytotoxicity

Development of poly(ADP-ribose) polymerase inhibitors (PARPi's) continues to be an attractive area of research due to synthetic lethality in DNA repair deficient cancers; however, PARPi's also have potential as therapeutics to prevent harmful inflammation. We investigated the pharmacological impact of incorporating spirodiamine motifs into the phthalazine architecture of FDA approved PARPi olaparib. Synthesized analogues were screened for PARP-1 affinity, enzyme specificity, catalytic inhibition, DNA damage, and cytotoxicity. This work led to the identification of 10e (12.6 ± 1.1 nM), which did not induce DNA damage at similar drug concentrations as olaparib. Interestingly, several worst in class compounds with low PARP-1 affinity, including 15b (4397 ± 1.1 nM), induced DNA damage at micromolar concentrations, which can explain the cytotoxicity observed in vitro. This work provides further evidence that high affinity PARPi's can be developed without DNA damaging properties offering potential new drugs for treating inflammatory related diseases.

Identification of potential inhibitors of PARP-1, a regulator of caspase-independent cell death pathway, from Withania somnifera phytochemicals for combating neurotoxicity: A structure-based in-silico study

Poly (ADP-ribose) polymerase-1 (PARP-1) reverses DNA damage by repairing DNA nicks and breaks in the normal cellular environment. However, during abnormal conditions like stroke and other neurological disorders, overactivation of PARP-1 leads to neuronal cell death via a caspase-independent programmed cell death pathway. Strategies involving inhibition or knockout of PARP-1 have proved beneficial in combating neuro-cytotoxicity. In this study, we performed in-silico analysis of 27 phytochemicals of Withania somnifera (Ashwagandha), to investigate their inhibition efficiency against PARP-1. Out of 27 phytochemicals, we report 12 phytochemicals binding to the catalytic domain of PARP-1 with an affinity higher than FR257517, PJ34 and Talazoparib (highly potent inhibitors of the enzyme). Among these 12 compounds, five phytochemicals namely Stigmasterol, Withacnistin, Withaferin A, Withanolide G and Withanolide B show an exceptionally high binding affinity for the catalytic domain of PARP-1 and bind to the enzyme with similar hydrogen bond formation and hydrophobic interaction pattern as their inhibitors. All of these phytochemicals are BBB permeable so that they can be further developed into potential future neuro-therapeutic drugs against neurodegenerative disorders involving neuronal cell death.

Synthesis, biological activity and multiscale molecular modeling studies for coumaryl-carboxamide derivatives as selective carbonic anhydrase IX inhibitors

New coumaryl-carboxamide derivatives with the thiourea moiety as a linker between the alkyl chains and/or the heterocycle nucleus were synthesized and their inhibitory activity against the human carbonic anhydrase (hCA) isoforms hCA I, II, VII and IX were evaluated. While the hCA I, II and VII isoforms were not inhibited by the investigated compounds, the tumour-associated isoform hCA IX was inhibited in the high nanomolar range. 2-Oxo-N-((2-(pyrrolidin-1-yl)ethyl)carbamothioyl)-2H-chromene-3-carboxamide (e11) exhibited a selective inhibitory action against hCA IX with the K_i of 107.9 nM. In order to better understand the inhibitory profiles of studied molecules, multiscale molecular modeling approaches were used. Different molecular docking algorithms were used to investigate binding poses and predicted binding energies of studied compounds at the active sites of the CA I, II, VII and IX isoforms.