In silico drug repositioning against human NRP1 to block SARS-CoV-2 host entry

Despite COVID-19 turned into a pandemic, no approved drug for the treatment or globally available vaccine is out yet. In such a global emergency, drug repurposing approach that bypasses a costly and long-time demanding drug discovery process is an effective way in search of finding drugs for the COVID-19 treatment. Recent studies showed that SARS-CoV-2 uses neuropilin-1 (NRP1) for host entry. Here we took advantage of structural information of the NRP1 in complex with C-terminal of spike (S) protein of SARS-CoV-2 to identify drugs that may inhibit NRP1 and S protein interaction. U.S. Food and Drug Administration (FDA) approved drugs were screened using docking simulations. Among top drugs, well-tolerated drugs were selected for further analysis. Molecular dynamics (MD) simulations of drugs-NRP1 complexes were run for 100 ns to assess the persistency of binding. MM/GBSA calculations from MD simulations showed that eltrombopag, glimepiride, sitagliptin, dutasteride, and ergotamine stably and strongly bind to NRP1. In silico Alanine scanning analysis revealed that Tyr²⁹⁷, Trp³⁰¹, and Tyr³⁵³ amino acids of NRP1 are critical for drug binding. Validating the effect of drugs analyzed in this paper by experimental studies and clinical trials will expedite the drug discovery process for COVID-19.     

Computational modeling of human coreceptor CCR5 antagonist as a HIV-1 entry inhibitor: using an integrated homology modeling,docking, and membrane molecular dynamics simulation analysis approach

Chemokine receptor 5 (CCR5) is an integral membrane protein that is utilized during human immunodeficiency virus type-1 entry into host cells. CCR5 is a G-protein coupled receptor that contains seven transmembrane (TM) helices. However, the crystal structure of CCR5 has not been reported. A homology model of CCR5 was developed based on the recently reported CXCR4 structure as template. Automated docking of the most potent (14), medium potent (37), and least potent (25) CCR5 antagonists was performed using the CCR5 model. To characterize the mechanism responsible for the interactions between ligands (14, 25, and 37) and CCR5, membrane molecular dynamic (MD) simulations were performed. The position and orientation of ligands (14, 25, and 37) were found to be changed after MD simulations, which demonstrated the ability of this technique to identify binding modes. Furthermore, at the end of simulation, it was found that residues identified by docking were changed and some new residues were introduced in the proximity of ligands. Our results are in line with the majority of previous mutational reports. These results show that hydrophobicity is the determining factor of CCR5 antagonism. In addition, salt bridging and hydrogen bond contacts between ligands (14, 25, and 37) and CCR5 are also crucial for inhibitory activity. The residues newly identified by MD simulation are Ser160, Phe166, Ser180, His181, and Trp190, and so far no site-directed mutagenesis studies have been reported. To determine the contributions made by these residues, additional mutational studies are suggested. We propose a general binding mode for these derivatives based on the MD simulation results of higher (14), medium (37), and lower (25) potent inhibitors. Interestingly, we found some trend for these inhibitors such as, salt bridge interaction between basic nitrogen of ligand and acidic Glu283 seemed necessary for inhibitory activity. Also, two aromatic pockets (pocket I – TM1-3 and pocket II – TM3-6) were linked by the central polar region in TM7, and the simulated inhibitors show important interactions with the Trp86, Tyr89, Tyr108, Phe112, Ile198, Tyr251, Leu255, and Gln280 and Glu283 residues. These results shed light on the usage of MD simulation to identify more stable, optimal binding modes of the inhibitors.     

In silico characterization of binding mode of CCR8 inhibitor: homology modeling,docking and membrane based MD simulation study

Human CC-chemokine receptor 8 (CCR8) is a crucial drug target in asthma that belongs to G-protein-coupled receptor superfamily, which is characterized by seven transmembrane helices. To date, there is no X-ray crystal structure available for CCR8; this hampers active research on the target. Molecular basis of interaction mechanism of antagonist with CCR8 remains unclear. In order to provide binding site information and stable binding mode, we performed modeling, docking and molecular dynamics (MD) simulation of CCR8. Docking study of biaryl-ether-piperidine derivative (13C) was performed inside predefined CCR8 binding site to get the representative conformation of 13C. Further, MD simulations of receptor and complex (13C-CCR8) inside dipalmitoylphosphatidylcholine lipid bilayers were performed to explore the effect of lipids. Results analyses showed that the Gln91, Tyr94, Cys106, Val109, Tyr113, Cys183, Tyr184, Ser185, Lys195, Thr198, Asn199, Met202, Phe254, and Glu286 were conserved in both docking and MD simulations. This indicated possible role of these residues in CCR8 antagonism. However, experimental mutational studies on these identified residues could be effective to confirm their importance in CCR8 antagonism. Furthermore, calculated Coulombic interactions represented the crucial roles of Glu286, Lys195, and Tyr113 in CCR8 antagonism. Important residues identified in this study overlap with the previous non-peptide agonist (LMD-009) binding site. Though, the non-peptide agonist and currently studied inhibitor (13C) share common substructure, but they differ in their effects on CCR8. So, to get more insight into their agonist and antagonist effects, further side-by-side experimental studies on both agonist (LMD-009) and antagonist (13C) are suggested.     

Binding of histamine to the H1 receptor—a molecular dynamics study

Binding of histamine to the G-protein coupled histamine H₁ receptor plays an important role in the context of allergic reactions; however, no crystal structure of the resulting complex is available yet. To deduce the histamine binding site, we performed unbiased molecular dynamics (MD) simulations on a microsecond time scale, which allowed to monitor one binding event, in which particularly the residues of the extracellular loop 2 were involved in the initial recognition process. The final histamine binding pose in the orthosteric pocket is characterized by interactions with Asp107^3.32, Tyr108^3.33, Thr194^5.43, Asn198^5.46, Trp428^6.48, Tyr431^6.51, Phe432^6.52, and Phe435^6.55, which is in agreement with existing mutational data. The conformational stability of the obtained complex structure was subsequently confirmed in 2 μs equilibrium MD simulations, and a metadynamics simulation proved that the detected binding site represents an energy minimum. A complementary investigation of a D107A mutant, which has experimentally been shown to abolish ligand binding, revealed that this exchange results in a significantly weaker interaction and enhanced ligand dynamics. This finding underlines the importance of the electrostatic interaction between the histamine ammonium group and the side chain of Asp107^3.32 for histamine binding.     

Theoretical calculations of the catalytic triad in short-chain alcohol dehydrogenases/reductases

Three highly conserved active site residues (Ser, Tyr, and Lys) of the family of short-chain alcohol dehydrogenases/reductases (SDRs) were demonstrated to be essential for catalytic activity and have been denoted the catalytic triad of SDRs. In this study computational methods were adopted to study the ionization properties of these amino acids in SDRs from Drosophila melanogaster and Drosophila lebanonensis. Three enzyme models, with different ionization scenarios of the catalytic triad that might be possible when inhibitors bind to the enzyme cofactor complex, were constructed. The binding of the two alcohol competitive inhibitors were studied using automatic docking by the Internal Coordinate Mechanics program, molecular dynamic (MD) simulations with the AMBER program package, calculation of the free energy of ligand binding by the linear interaction energy method, and the hydropathic interactions force field. The calculations indicated that deprotonated Tyr acts as a strong base in the binary enzyme-NAD⁺ complex. Molecular dynamic simulations for 5 ns confirmed that deprotonated Tyr is essential for anchoring and orientating the inhibitors at the active site, which might be a general trend for the family of SDRs. The findings here have implications for the development of therapeutically important SDR inhibitors.     

Regulatory phosphorylation of cyclin-dependent kinase 2: insights from molecular dynamics simulations

The structures of fully active cyclin-dependent kinase-2 (CDK2) complexed with ATP and peptide substrate, CDK2 after the catalytic reaction, and CDK2 inhibited by phosphorylation at Thr14/Tyr15 were studied using molecular dynamics (MD) simulations. The structural details of the CDK2 catalytic site and CDK2 substrate binding box were described. Comparison of MD simulations of inhibited complexes of CDK2 was used to help understand the role of inhibitory phosphorylation at Thr14/Tyr15. Phosphorylation at Thr14/Tyr15 causes ATP misalignment for the phosphate-group transfer, changes in the Mg²⁺ coordination sphere, and changes in the H-bond network formed by CDK2 catalytic residues (Asp127, Lys129, Asn132). The inhibitory phosphorylation causes the G-loop to shift from the ATP binding site, which leads to opening of the CDK2 substrate binding box, thus probably weakening substrate binding. All these effects explain the decrease in kinase activity observed after inhibitory phosphorylation at Thr14/Tyr15 in the G-loop. Interaction of the peptide substrate, and the phosphorylated peptide product, with CDK2 was also studied and compared. These results broaden hypotheses drawn from our previous MD studies as to why a basic residue (Arg/Lys) is preferred at the P₊₂ substrate position. Figure View of the substrate binding site of the fully active cyclin-dependent kinase-2 (CDK2) (pT160-CDK2/cyclin A/ATP). The pThr160 activation site is located in the T-loop (yellow secondary structure). The G-loop, which partly forms the ATP binding site, is shown in blue. The Thr14 and Tyr15 inhibitory phosphorylation sites located in the G-loop are shown in licorice representation     

Synthesis and structure–activity relationship of non-peptidic antagonists of neuropilin-1 receptor

Neuropilins (NRPs) are VEGF-A₁₆₅ co-receptors over-expressed in tumor cells, and considered as targets in angiogenic-related pathologies. We previously identified compound 1, the first non-peptidic antagonist of the VEGF-A₁₆₅/NRP binding, which exhibits in vivo anti-angiogenic and anti-tumor activities. We report here the synthesis and biological evaluations of new antagonists structurally-related to compound 1. Among these molecules, 4a, 4c and 4d show cytotoxic effects on HUVEC and MDA-MB-31 cells, and antagonize VEGF-A₁₆₅/NRP-1 binding. This study confirmed our key structure–activity relationships hypothesis and paved the way to compound 1 ‘hit to lead’ optimization.     

Identification of the first small-molecule inhibitor of the REV7 DNA repair protein interaction

DNA interstrand crosslink (ICL) repair (ICLR) has been implicated in the resistance of cancer cells to ICL-inducing chemotherapeutic agents. Despite the clinical significance of ICL-inducing chemotherapy, few studies have focused on developing small-molecule inhibitors for ICLR. The mammalian DNA polymerase ζ, which comprises the catalytic subunit REV3L and the non-catalytic subunit REV7, is essential for ICLR. To identify small-molecule compounds that are mechanistically capable of inhibiting ICLR by targeting REV7, high-throughput screening and structure–activity relationship (SAR) analysis were performed. Compound 1 was identified as an inhibitor of the interaction of REV7 with the REV7-binding sequence of REV3L. Compound 7 (an optimized analog of compound 1) bound directly to REV7 in nuclear magnetic resonance analyses, and inhibited the reactivation of a reporter plasmid containing an ICL in between the promoter and reporter regions. The normalized clonogenic survival of HeLa cells treated with cisplatin and compound 7 was lower than that for cells treated with cisplatin only. These findings indicate that a small-molecule inhibitor of the REV7/REV3L interaction can chemosensitize cells by inhibiting ICLR.     

Study on the binding mode of a pyrrolotriazin derivative with JAK2 by docking and MD simulation

The pyrrolotriazin derivative 2-(4-(4-((7-(3-(N-methylmethylsulfonamido)phenyl)pyrrolo [2,1-f][1,2,4]triazin-2-yl)amino)phenyl)piperidin-1-yl)acetamide (PPA) is a potential Janus kinase 2 (JAK2) inhibitor. The binding mode between PPA and JAK2 was investigated by using a combined method of docking, molecular dynamics (MD) simulation and binding free-energy calculation. The docking calculations preliminarily indicated that there were two possible binding modes 1 and 2; MD simulations and binding free-energy calculations identified that binding mode 1 was more stable and favourable, with the lower MM-PBSA binding free energy of ?34.00?±?0.17?kcal/mol. Moreover, some valuable binding information is revealed as follows: the inhibitor PPA is suitably located at the ATP-binding site of JAK2 and the hydrophobic interaction plays an essential role. PPA not only interacts with residues Leu855, Val863, Ala880, Tyr931, Leu932 and Leu983 via hydrophobic interaction but also interacts with Ser936 and Asp994 by hydrogen bonds. These two factors are advantageous for PPA to strongly bind to JAK2. These results help to understand the action mechanisms and designing new compounds with a higher affinity to JAK2.     

VEGF binding to NRP1 is essential for VEGF stimulation of endothelial cell migration, complex formation between NRP1 and VEGFR2, and signaling via FAK Tyr407 phosphorylation

In endothelial cells, neuropilin-1 (NRP1) binds vascular endothelial growth factor (VEGF)-A and is thought to act as a coreceptor for kinase insert domain-containing receptor (KDR) by associating with KDR and enhancing VEGF signaling. Here we report mutations in the NRP1 b1 domain (Y297A and D320A), which result in complete loss of VEGF binding. Overexpression of Y297A and D320A NRP1 in human umbilical vein endothelial cells reduced high-affinity VEGF binding and migration toward a VEGF gradient, and markedly inhibited VEGF-induced angiogenesis in a coculture cell model. The Y297A NRP1 mutant also disrupted complexation between NRP1 and KDR and decreased VEGF-dependent phosphorylation of focal adhesion kinase at Tyr407, but had little effect on other signaling pathways. Y297A NRP1, however, heterodimerized with wild-type NRP1 and NRP2 indicating that nonbinding NRP1 mutants can act in a dominant-negative manner through formation of NRP1 dimers with reduced binding affinity for VEGF. These findings indicate that VEGF binding to NRP1 has specific effects on endothelial cell signaling and is important for endothelial cell migration and angiogenesis mediated via complex formation between NRP1 and KDR and increased signaling to focal adhesions. Identification of key residues essential for VEGF binding and biological functions provides the basis for a rational design of antagonists of VEGF binding to NRP1.     

Agonism activities of lyso-phosphatidylcholines (LPC) Ligands binding to peroxisome proliferator-activated receptor gamma (PPARγ)

Abstract

PPARγ is an isoform of peroxisome proliferator-activated receptor (PPAR) belonging to a super family of nuclear receptors and is a primary target of the effective drug to treat the type II diabetes. The experiments found that Lyso-phosphatidylcholines (LPC) could bind to PPARγ, but the binding modes remain unknown. We used the Molecular Docking and Molecular Dynamic (MD) simulations to study the binding of four LPC ligands (LPC16:0, LPC18:0, LPC18:1-1 and LPC18:1-2) to PPARγ. The two-step MD simulations were employed to determine the final binding modes. The 20?ns MD simulations for four final LPC-PPARγ complexes were performed to analyze their structures, the binding key residues, and agonism activities. The results reveal that three LPC ligands (LPC16:0, LPC18:0 and LPC18:1-1) bind to Arm II and III regions of the Ligand Binding Domain (LBD) pocket, whereas they do not interact with Tyr473 of Helix 12 (H12). In contrast, LPC18:1-2 can form the hydrogen bonds with Tyr473 and bind into Arm I and II regions. Comparing with the paradigm systems of the full agonist (Rosiglitazone–PPARγ) and the partial agonist (MRL24–PPARγ), our results indicate that LPC16:0, LPC18:0 and LPC18:1-1 could be the potential partial agonists and LPC18:1-2 could be a full agonist. The in-depth analysis of the residue fluctuations and structure alignment confirm the present prediction of the LPC agonism activities.

Communicated by Ramaswamy H. Sarma     

Optimization and bioevaluation of Cdc37-derived peptides: An insight into Hsp90-Cdc37 protein-protein interaction modulators

Targeting Hsp90-Cdc37 protein-protein interaction (PPI) is becoming an alternative approach for future anti-cancer drug development. We previously reported the discovery of an eleven-residue peptide (Pep-1) with micromolar activity for the disruption of Hsp90-Cdc37 PPI. Efforts to improve upon the Pep-1 led to the discovery of more potent modulators for Hsp90-Cdc37 PPI. Through the analysis of peptides binding patterns, more peptides were designed for further verification which resulted in Pep-5, the shortest peptide targeting Hsp90-Cdc37, exerting the optimal structure and the most efficient binding mode. Subsequent MD simulation analysis also confirmed that Pep-5 could perform more stable binding ability and better ligand properties than Pep-1. Under the premise of retentive binding capacity, Pep-5 exhibited lower molecular weight and higher ligand efficiency with a K_d value of 5.99 μM (Pep-1 K_d = 6.90 μM) in both direct binding determination and biological evaluation. The optimal and shortest Pep-5 might provide a breakthrough and a better model for the future design of small molecule inhibitors targeting Hsp90-Cdc37 PPI.     

In silico study of carvone derivatives as potential neuraminidase inhibitors

Recent outbreaks of highly pathogenic influenza strains have highlighted the need to develop new anti-influenza drugs. Here, we report an in silico study of carvone derivatives to analyze their binding modes with neuraminidase (NA) active sites. Two proposed carvone analogues, CV(A) and CV(B), with 36 designed ligands were predicted to inhibit NA (PDB ID: 3TI6) using molecular docking. The design is based on structural resemblance with the commercial inhibitor, oseltamivir (OTV), ligand polarity, and amino acid residues in the NA active sites. Docking simulations revealed that ligand A18 has the lowest energy binding (?G_bind) value of ?8.30 kcal mol^-1, comparable to OTV with ?G_bind of ?8.72 kcal mol^-1. A18 formed seven hydrogen bonds (H-bonds) at residues Arg292, Arg371, Asp151, Trp178, Glu227, and Tyr406, while eight H-bonds were formed by OTV with amino acids Arg118, Arg292, Arg371, Glu119, Asp151, and Arg152. Molecular dynamics (MD) simulation was conducted to compare the stability between ligand A18 and OTV with NA. Our simulation study showed that the A18-NA complex is as stable as the OTV-NA complex during the MD simulation of 50 ns through the analysis of RMSD, RMSF, total energy, hydrogen bonding, and MM/PBSA free energy calculations.     

Identification of a novel selective small-molecule inhibitor of protein arginine methyltransferase 5 (PRMT5) by virtual screening,resynthesis and biological evaluations

As one of the most promising anticancer target in protein arginine methyltransferase (PRMT) family, PRMT5 has been drawing more and more attentions, and many efforts have been devoted to develop its inhibitors. In this study, three PRMT5 inhibitors (9, 16, and 23) with novel scaffolds were identified by performing pharmacophore- and docking-based virtual screening combined with in vitro radiometric-based scintillation proximity assay (SPA). Substructure search based on the scaffold of the most active 9 afforded 26 additional analogues, and SPA results indicated that two analogues (9–1 and 9–2) showed increased PRMT5 inhibitory activity compared with the parental compound. Resynthesis of 9, 9–1, and 9–2 confirmed their PRMT5 enzymatic inhibition activity. In addition, compound 9–1 displayed selectivity against PRMT5 over other key homological members (PRMT1 and CARM1 (PRMT4)). While the structure–activity relationship (SAR) of this series of compounds was discussed to provide clues for further structure optimization, the probable binding modes of active compounds were also probed by molecular docking and molecular dynamics simulations. Finally, the antiproliferative effect of 9–1 on MV4-11 leukemia cell line was confirmed and its impact on regulating the target gene of PRMT5 was also validated. The hit compounds identified in this work have provided more novel scaffolds for future hit-to-lead optimization of small-molecule PRMT5 inhibitors.     

Interaction of a calix[4]arene derivative with a DOPC bilayer: Biomolecular simulations towards chloride transport

The ability of a calix[4]arene derivative (CX-1), bearing four protonated NH₃⁺ groups located in the upper rim and aliphatic tails in the lower rim, to interact with a 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) model bilayer and promote transmembrane chloride transport was investigated by molecular dynamics (MD) simulations. Unconstrained MD simulations show that the interaction of CX-1 with DOPC occurs via the NH₃⁺ groups, which are able to establish electrostatic interactions and multiple hydrogen bonds with the DOPC phosphate groups, while the aliphatic tails point towards the water phase (when CX-1 starts from the water phase) or to the membrane (when CX-1 is initially positioned within the bilayer). The interaction does not induce any relevant perturbation on the biophysical properties of the bilayer system (area per lipid, thickness, and hydration) apart from a systematic increase in the order parameter of the C2 carbon atom of the sn-1 lipid tail, meaning that the bilayer conserves its integrity. Since total internalization of CX-1 was not observed in the unconstrained MD time-scale, constant velocity steered molecular dynamics (SMD) simulations were performed in order to simulate the CX-1 permeation across the bilayer. At pulling velocities lower than 0.0075 nm ps^− 1, chloride transport was observed. The Potential of Mean Force (PMF), calculated with the weighted histogram analysis method, indicates a barrier of ca. 58 kJ mol^− 1 for this mobile carrier to cross the membrane.     

Tyr74 is essential for the formation, stability and function of Plasmodium falciparum triosephosphate isomerase dimer

Plasmodium falciparum triosephosphate isomerase (PfTIM) is known to be functional only as a homodimer. Although many studies have shown that the interface Cys13 plays a major role in the stability of the dimer, a few reports have demonstrated that structurally conserved Tyr74 may be essential for the stability of PfTIM dimer. To understand the role of Tyr74, we have performed molecular dynamics (MD) simulations of monomeric and dimeric PfTIM mutated to glycine and cysteine at position 74. Simulations of the monomer revealed that mutant Tyr74Gly does not produce changes in folding and stability of the monomer. Interestingly, comparison of the flexibility of Tyr74 in the monomer and dimer revealed that this residue possesses an intrinsic restricted mobility, indicating that Tyr74 is an anchor residue required for homodimerization. Tyr74 also appears to play an important role in binding by facilitating the disorder-to-order transitions of loops 1 and 3, which allows Cys13 to form favorable interactions with loop 3 and Lys12 to be locked in a favorable position for catalysis. High-temperature MD simulations of the wild-type and Tyr74Gly PfTIM dimers showed that the aromatic moiety of Tyr74 is necessary to preserve the geometry and native contacts between loops 1 and 3 at the interface of the dimer. Disulfide cross-linking between mutant Tyr74Cys and Cys13 further revealed that Tyr74 stabilizes the geometry of loop 1 (which contains the catalytic residue Lys12) and the interactions between loops 1 and 3 via aromatic-aromatic interactions with residues Phe69, Tyr101, and Phe102. Principal component analysis showed that Tyr74 is also necessary to preserve the collective motions in the dimer that contribute to the catalytic efficiency of PfTIM dimer. We conclude that Tyr74 not only plays a role in the stability of the dimer, but also participates in the dimerization process and collective motions via coupled disorder-to-order transitions of intrinsically disordered regions, necessary for efficiency in the catalytic function of PfTIM.     

Variations in the Binding Pocket of an Inhibitor of the Bacterial Division Protein FtsZ across Genotypes and Species

The recent increase in antibiotic resistance in pathogenic bacteria calls for new approaches to drug-target selection and drug development. Targeting the mechanisms of action of proteins involved in bacterial cell division bypasses problems associated with increasingly ineffective variants of older antibiotics; to this end, the essential bacterial cytoskeletal protein FtsZ is a promising target. Recent work on its allosteric inhibitor, PC190723, revealed in vitro activity on Staphylococcus aureus FtsZ and in vivo antimicrobial activities. However, the mechanism of drug action and its effect on FtsZ in other bacterial species are unclear. Here, we examine the structural environment of the PC190723 binding pocket using PocketFEATURE, a statistical method that scores the similarity between pairs of small-molecule binding sites based on 3D structure information about the local microenvironment, and molecular dynamics (MD) simulations. We observed that species and nucleotide-binding state have significant impacts on the structural properties of the binding site, with substantially disparate microenvironments for bacterial species not from the Staphylococcus genus. Based on PocketFEATURE analysis of MD simulations of S. aureus FtsZ bound to GTP or with mutations that are known to confer PC190723 resistance, we predict that PC190723 strongly prefers to bind Staphylococcus FtsZ in the nucleotide-bound state. Furthermore, MD simulations of an FtsZ dimer indicated that polymerization may enhance PC190723 binding. Taken together, our results demonstrate that a drug-binding pocket can vary significantly across species, genetic perturbations, and in different polymerization states, yielding important information for the further development of FtsZ inhibitors.     

Discovery of a novel allosteric inhibitor scaffold for polyadenosine-diphosphate-ribose polymerase 14 (PARP14) macrodomain 2

The polyadenosine-diphosphate-ribose polymerase 14 (PARP14) has been implicated in DNA damage response pathways for homologous recombination. PARP14 contains three (ADP ribose binding) macrodomains (MD) whose exact contribution to overall PARP14 function in pathology remains unclear. A medium throughput screen led to the identification of N-(2(-9H-carbazol-1-yl)phenyl)acetamide (GeA-69, 1) as a novel allosteric PARP14 MD2 (second MD of PARP14) inhibitor. We herein report medicinal chemistry around this novel chemotype to afford a sub-micromolar PARP14 MD2 inhibitor. This chemical series provides a novel starting point for further development of PARP14 chemical probes.