Search

The lack of publicly available, large, and unbiased datasets is a key bottleneck for the application of machine learning (ML) methods in synthetic chemistry. Data from electronic laboratory notebooks (ELNs) could provide less biased, large datasets, but no such datasets have been made publicly available. The first real-world dataset from the ELNs of a large pharmaceutical company is disclosed and its relationship to high-throughput experimentation (HTE) datasets is described. For chemical yield predictions, a key task in chemical synthesis, an attributed graph neural network (AGNN) performs as good or better than the best previous models on two HTE datasets for the Suzuki and Buchwald-Hartwig reactions. However, training of the AGNN on the ELN dataset does not lead to a predictive model. The implications of using ELN data for training ML-based models are discussed in the context of yield predictions.

Supplementary Data from EP300 Selectively Controls the Enhancer Landscape of MYCN-Amplified Neuroblastoma

Gene expression is regulated by promoters and enhancers marked by histone H3 lysine 27 acetylation (H3K27ac), which is established by the paralogous histone acetyltransferases (HAT) EP300 and CBP. These enzymes display overlapping regulatory roles in untransformed cells, but less characterized roles in cancer cells. We demonstrate that the majority of high-risk pediatric neuroblastoma (NB) depends on EP300, whereas CBP has a limited role. EP300 controls enhancer acetylation by interacting with TFAP2β, a transcription factor member of the lineage-defining transcriptional core regulatory circuitry (CRC) in NB. To disrupt EP300, we developed a proteolysis-targeting chimera (PROTAC) compound termed “JQAD1” that selectively targets EP300 for degradation. JQAD1 treatment causes loss of H3K27ac at CRC enhancers and rapid NB apoptosis, with limited toxicity to untransformed cells where CBP may compensate. Furthermore, JQAD1 activity is critically determined by cereblon (CRBN) expression across NB cells.Significance:EP300, but not CBP, controls oncogenic CRC-driven transcription in high-risk NB by binding TFAP2β. We developed JQAD1, a CRBN-dependent PROTAC degrader with preferential activity against EP300 and demonstrated its activity in NB. JQAD1 has limited toxicity to untransformed cells and is effective in vivo in a CRBN-dependent manner.This article is highlighted in the In This Issue feature, p. 587

The noninflamed microenvironment in prostate cancer represents a barrier to immunotherapy. Genetic alterations underlying cancer cell–intrinsic oncogenic signaling are increasingly appreciated for their role in shaping the immune landscape. Recently, we identified Pygopus 2 ( PYGO2 ) as the driver oncogene for the amplicon at 1q21.3 in prostate cancer. Here, using transgenic mouse models of metastatic prostate adenocarcinoma, we found that Pygo2 deletion decelerated tumor progression, diminished metastases, and extended survival. Pygo2 loss augmented the activation and infiltration of cytotoxic T lymphocytes (CTLs) and sensitized tumor cells to T cell killing. Mechanistically, Pygo2 orchestrated a p53/Sp1/Kit/Ido1 signaling network to foster a microenvironment hostile to CTLs. Genetic or pharmacological inhibition of Pygo2 enhanced the antitumor efficacy of immunotherapies using immune checkpoint blockade (ICB), adoptive cell transfer, or agents inhibiting myeloid-derived suppressor cells. In human prostate cancer samples, Pygo2 expression was inversely correlated with the infiltration of CD8 + T cells. Analysis of the ICB clinical data showed association between high PYGO2 level and worse outcome. Together, our results highlight a potential path to improve immunotherapy using Pygo2-targeted therapy for advanced prostate cancer.

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is the causative agent of the coronavirus disease of 2019 (COVID-19). Its genome encodes two open reading frames for two large proteins, PP1a and PP1ab. Within the two polypeptide stretches, there are two proteases that process the large proteins into 15 discrete proteins essential for the assembly of the virion during its replication. We describe herein the cloning of the genes for these discrete proteins optimized for expression in Escherichia coli, production of the proteins, and their purification to homogeneity. These included all but six: NSP6, which possesses eight transmembrane regions, and five that are small proteins/peptides (E, ORF3b, ORF6, ORF7b, and ORF10). These proteins are intended for experimental validation of small-molecule binders as molecular template hits. The proof of concept was established with the ADP-ribosylhydrolase (ARH) domain of NSP3 in discovery of small-molecule templates that could serve as the basis for further optimization. The hit molecules include one submicromolar and a few low-micromolar binders to the ARH domain. Availability of these proteins in soluble forms opens up the opportunity for discoveries of novel templates with the potential for anti-COVID-19 pharmaceuticals.

Ferrocene derivatives have a wide range of applications, including as ligands in asymmetric catalysis, due to their chemical stability, rigid backbone, steric bulk, and ability to encode stereochemical information via planar chirality. Unfortunately, few of the available molecular mechanics force fields incorporate parameters for the accurate study of this important building block. Here, we present a MM3* force field for ferrocenyl ligands, which was generated using the quantum-guided molecular mechanics (Q2MM) method. Detailed validation by comparison to DFT calculations and crystal structures demonstrates the accuracy of the parameters and uncovers the physical origin of deviations through excess energy analysis. Combining the ferrocene force field with a force field for Pd-allyl complexes and comparing the crystal structures shows the compatibility with previously developed MM3* force fields. Finally, the ferrocene force field was combined with a previously published transition-state force field to predict the stereochemical outcomes of the aminations of Pd-allyl complexes with different amines and different chiral ferrocenyl ligands, with an

A phosphorus(III)-mediated dearomatization of ortho-substituted dianiline squaraine dyes results in an unusual skeletal rearrangement to provide exotic, highly conjugated benzofuranone and oxindole scaffolds bearing a C3 side chain comprised of a linear conflagration of an enol, a phosphorus ylide, and 2,4-disubstituted aniline. Employing experimental and computational analysis, a mechanistic evaluation revealed a striking dependence on the acidity of the aniline ortho substituent. Notably, the rearrangement adducts underwent rapid and complete reversion to the parent squaraine in the presence of a Brønsted acid.

Understanding the mechanisms of enzymatic catalysis requires a detailed understanding of the complex interplay of structure and dynamics of large systems that is a challenge for both experimental and computational approaches. More importantly, the computational demands of QM/MM simulations mean that the dynamics of the reaction can only be considered on a timescale of nanoseconds even though the conformational changes needed to reach the catalytically active state happen on a much slower timescale. Here we demonstrate an alternative approach that uses transition state force fields (TSFFs) derived by the quantum-guided molecular mechanics (Q2MM) method that provides a consistent treatment of the entire system at the classical molecular mechanics level and allows simulations at the microsecond timescale. Application of this approach to the second hydride transfer transition state of HMG-CoA reductase from Pseudomonas mevalonii (PmHMGR) identified three remote residues, R396, E399 and L407, (15–27 Å away from the active site) that have a remote dynamic effect on enzyme activity. The predictions were subsequently validated experimentally via site-directed mutagenesis. These results show that microsecond timescale MD simulations of transition states are possible and can predict rather than just rationalize remote allosteric residues., Transition state force fields enable MD simulations at the transition state of HMGCoA reductase that sample the transition state ensemble on the μs timescale to identify remote residues that affect the reaction rate.

The noninflamed microenvironment in prostate cancer represents a barrier to immunotherapy. Genetic alterations underlying cancer cell-intrinsic oncogenic signaling have emerging roles in shaping the immune landscape. Recently, we identified Pygopus 2 (PYGO2) as the driver oncogene for the amplicon at 1q21.3. Here, using transgenic mouse models of metastatic prostate adenocarcinoma, we found that Pygo2 deletion decelerated tumor progression, diminished metastases, and extended survival. Pygo2 loss augmented the infiltration of cytotoxic T lymphocytes (CTLs) and sensitized tumor cells to T cell killing. Mechanistically, Pygo2 orchestrated a p53/Sp1/Kit/Ido1 signaling network to foster a microenvironment hostile to CTLs. Genetic or pharmacological inhibition of Pygo2 enhanced the anti-tumor efficacy of immunotherapies using immune checkpoint blockade, adoptive cell transfer, or myeloid-derived suppressor cell inhibitors. In human prostate cancer samples, Pygo2 expression was inversely correlated with CD8+ T cells. Our results highlight a promising path to improving immunotherapy with targeted therapy for lethal prostate cancer.

The nucleophilic addition of alkyl, benzyl, and allyl organozinc reagents to protected pyridinium riboses proceeds under mild conditions and with yields of >90% in several cases and complete regioselectivity for the 4-position.

Deregulation of dual-specificity tyrosine phosphorylation-regulated kinase 1A (DYRK1A) plays a significant role in developmental brain defects, early-onset neurodegeneration, neuronal cell loss, dementia, and several types of cancer. Herein, we report the discovery of three new classes of N-heterocyclic DYRK1A inhibitors based on the potent, yet toxic kinase inhibitors, harmine and harmol. An initial in vitro evaluation of the small molecule library assembled revealed that the core heterocyclic motifs benzofuranones, oxindoles, and pyrrolones, showed statistically significant DYRK1A inhibition. Further, the utilization of a low cost, high-throughput functional genomic in vivo model system to identify small molecule inhibitors that normalize DYRK1A overexpression phenotypes is described. This in vivo assay substantiated the in vitro results, and the resulting correspondence validates generated classes as architectural motifs that serve as potential DYRK1A inhibitors. Further expansion and analysis of these core compound structures will allow discovery of safe, more effective chemical inhibitors of DYRK1A to ameliorate phenotypes caused by DYRK1A overexpression.

Identity determining transcription factors (TFs), or core regulatory (CR) TFs, are governed by cell-type specific super enhancers (SEs). Drugs to selectively inhibit CR circuitry are of high interest for cancer treatment. In alveolar rhabdomyosarcoma, PAX3-FOXO1 activates SEs to induce the expression of other CR TFs, providing a model system for studying cancer cell addiction to CR transcription. Using chemical genetics, the systematic screening of chemical matter for a biological outcome, here we report on a screen for epigenetic chemical probes able to distinguish between SE-driven transcription and constitutive transcription. We find that chemical probes along the acetylation-axis, and not the methylation-axis, selectively disrupt CR transcription. Additionally, we find that histone deacetylases (HDACs) are essential for CR TF transcription. We further dissect the contribution of HDAC isoforms using selective inhibitors, including the newly developed selective HDAC3 inhibitor LW3. We show HDAC1/2/3 are the co-essential isoforms that when co-inhibited halt CR transcription, making CR TF sites hyper-accessible and disrupting chromatin looping., Core regulatory transcription factors are usually regulated by cell-type specific super enhancers (SEs). Here, the authors screen for chemical probes able to distinguish between SE-driven and promoter-driven transcription and find that histone deacetylases are selectively required for core regulatory transcription.

Chemical reactions are a complex process, as they involve interaction between several molecular compounds. As a result, predicting the success of a reaction is a non-trivial task, which often requires running several experiments in the lab. This process is is expensive, time consuming, and inefficient. As a result, in recent years, researchers have explored the use of machine learning algorithms to predict reaction success. These methods mainly rely on chemical properties of the molecules involved in the reactions. Despite their promising success, none of existing methods explored the use of structural properties of molecules in predicting reaction success. In this work, we develop an Attributed Graph Neural Network model that integrates both structural properties as well as chemical properties of molecules for predicting reaction success. Our model shows remarkable performance on two hand-crafted datasets obtained from high-throughput experiments, as well as one real-world dataset.

Deregulation of dual-specificity tyrosine phosphorylation-regulated kinase 1A (DYRK1A) plays a significant role in developmental brain defects, early-onset neurodegeneration, neuronal cell loss, and dementia. Herein, we report the discovery of three new classes of N-heterocyclic DYRK1A inhibitors based on the potent, yet toxic kinase inhibitors, harmine and harmol. An initial in vitro evaluation of the small molecule collection assembled revealed that the core heterocyclic motifs benzofuranones, oxindoles, and pyrrolones, showed statistically significant DYRK1A inhibition. Further, the utilization of a low cost, high-throughput functional genomic in vivo model system to identify small molecule inhibitors that normalize DYRK1A overexpression phenotypes is described. This in vivo assay confirmed the in vitro results, and the resulting correspondence validates generated classes as architectural motifs that serve as potential DYRK1A inhibitors. Further expansion and analysis of these core compound structures will allow discovery of safe, more effective chemical inhibitors of DYRK1A to ameliorate phenotypes caused by DYRK1A overexpression.

Deregulation of dual-specificity tyrosine phosphorylation-regulated kinase 1A (DYRK1A) plays a significant role in developmental brain defects, early-onset neurodegeneration, neuronal cell loss, and dementia. Herein, we report the discovery of three new classes of N-heterocyclic DYRK1A inhibitors based on the potent, yet toxic kinase inhibitors, harmine and harmol. An initial in vitro evaluation of the small molecule collection assembled revealed that the core heterocyclic motifs benzofuranones, oxindoles, and pyrrolones, showed statistically significant DYRK1A inhibition. Further, the utilization of a low cost, high-throughput functional genomic in vivo model system to identify small molecule inhibitors that normalize DYRK1A overexpression phenotypes is described. This in vivo assay confirmed the in vitro results, and the resulting correspondence validates generated classes as architectural motifs that serve as potential DYRK1A inhibitors. Further expansion and analysis of these core compound structures will allow discovery of safe, more effective chemical inhibitors of DYRK1A to ameliorate phenotypes caused by DYRK1A overexpression.

The conjugate addition of aryl boronic acids to enones is a powerful synthetic tool to introduce quaternary chiral centers, but the experimentally observed stereoselectivities vary widely, and the identification of suitable substrate–ligand combinations requires significant effort. We describe the development and application of a transition-state force field (TSFF) by the quantum-guided molecular mechanics (Q2MM) method that is validated using an automated screen of 9 ligands, 38 aryl boronic acids, and 22 enones, leading to a MUE of 1.8 kJ/mol and a R(2) value of 0.877 over 82 examples. A detailed error analysis identified the structural origin for the deviations in the small group of outliers. The TSFF was then used to predict the stereoselectivity for 27 ligands and 59 enones. The vast majority of the virtual screening results are in line with the expected results. Selected results for 6-substituted pyrox ligands, which were not part of the training set, were followed up by density functional theory and experimental studies.

Gene expression is regulated by promoters and enhancers marked by histone H3 lysine 27 acetylation (H3K27ac), which is established by the paralogous histone acetyltransferases (HAT) EP300 and CBP. These enzymes display overlapping regulatory roles in untransformed cells, but less characterized roles in cancer cells. We demonstrate that the majority of high-risk pediatric neuroblastoma (NB) depends on EP300, whereas CBP has a limited role. EP300 controls enhancer acetylation by interacting with TFAP2β, a transcription factor member of the lineage-defining transcriptional core regulatory circuitry (CRC) in NB. To disrupt EP300, we developed a proteolysis-targeting chimera (PROTAC) compound termed “JQAD1” that selectively targets EP300 for degradation. JQAD1 treatment causes loss of H3K27ac at CRC enhancers and rapid NB apoptosis, with limited toxicity to untransformed cells where CBP may compensate. Furthermore, JQAD1 activity is critically determined by cereblon (CRBN) expression across NB cells. Significance: EP300, but not CBP, controls oncogenic CRC-driven transcription in high-risk NB by binding TFAP2β. We developed JQAD1, a CRBN-dependent PROTAC degrader with preferential activity against EP300 and demonstrated its activity in NB. JQAD1 has limited toxicity to untransformed cells and is effective in vivo in a CRBN-dependent manner. This article is highlighted in the In This Issue feature, p. 587

Niemann-Pick disease type C1 (NPC1) is a rare genetic cholesterol storage disorder caused by mutations in the NPC1 gene. Mutations in this transmembrane late endosome protein lead to loss of normal cholesterol efflux from late endosomes and lysosomes. It has been shown that broad spectrum histone deacetylase inhibitors (HDACi's) such as Vorinostat correct the cholesterol accumulation phenotype in the majority of NPC1 mutants tested in cultured cells. In order to determine the optimal specificity for HDACi correction of the mutant NPC1s, we screened 76 HDACi's of varying specificity. We tested the ability of these HDACi's to correct the excess accumulation of cholesterol in patient fibroblast cells that homozygously express NPC1 I1061T , the most common mutation. We determined that inhibition of HDACs 1, 2, and 3 is important for correcting the defect, and combined inhibition of all three is needed to achieve the greatest effect, suggesting a need for multiple effects of the HDACi treatments. Identifying the specific HDACs involved in the process of regulating cholesterol trafficking in NPC1 will help to focus the search for more specific druggable targets.

The parameterization of a transition state force field for the title reaction is described. Validation for 82 literature examples leads to a MUE of 1.8 kJ/mol and an R2 of 0.877 between computed and experimental stereoselectivities. The use if the TSFF is demonstrated for a virtual library of 27 ligands and 59 enones.

The development of new chemical tools with improved properties is essential to chemical and cell biology. Of particular interest is the development of mimics of small molecules with important cellular function that allow the direct observation of their trafficking in a cell. To this end, a novel 15-azasterol has been designed and synthesized as a luminescent cholesterol mimic for the monitoring of cholesterol trafficking. The brightness of this probe, which is ∼32-times greater than the widely used dehydroergosterol probe, is combined with resistance to photobleaching in solution and in human fibroblasts and an exceptionally large Stokes-like shift of ∼150−200 nm. The photophysical properties of the probe have been studied experimentally and computationally, suggesting an intersystem crossing to the triplet excited state with subsequent phosphorescent decay. Molecular dynamics simulations show a similar binding mode of cholesterol and the azasterol probe to NPC proteins, demonstrating the structural similarity of the probe to cholesterol.

The generation of surrogate potential energy functions (PEF) that are orders of magnitude faster to compute but as accurate as the underlying training data from high-level electronic structure methods is one of the most promising applications of fitting procedures in chemistry. In previous work, we have shown that transition state force fields (TSFFs), fitted to the functional form of MM3* force fields using the quantum guided molecular mechanics (Q2MM) method, provide an accurate description of transition states that can be used for stereoselectivity predictions of small molecule reactions. Here, we demonstrate the applicability of the method for fit TSFFs to the well-established Amber force field, which could be used for molecular dynamics studies of enzyme reaction. As a case study, the fitting of a TSFF to the second hydride transfer in Pseudomonas mevalonii 3-hydroxy-3-methylglutaryl coenzyme A reductase (PmHMGR) is used. The differences and similarities to fitting of small molecule TSFFs are discussed.

HMG-CoA reductase (Pseudomonas mevalonii) utilizes mevalonate, coenzyme A (CoA) and the cofactor NAD in a complex mechanism involving two hydride transfers with cofactor exchange, accompanied by large conformational changes by a 50 residue subdomain, to generate HMG-CoA. Details about this mechanism such as the conformational changes that allow intermediate formation, cofactor exchange and product release remain unknown. The formation of the proposed intermediates has also not been observed in structural studies with natural substrates. Having been shown to be an essential enzyme for the survival of gram-positive antibiotic resistant pathogenic bacteria, studying its mechanism in detail will be beneficial in developing novel antibacterials. The enzyme has been shown to be catalytically active inside the crystal with dithio-HMG-CoA and NADH but curiously is found to be inactive in the reverse direction in the structure bound to mevalonate, CoA and NAD.To understand the factors limiting activity in the HMGR crystal with mevalonate, CoA and NAD, we studied the effect of crystallization components and pH on enzymatic activity. We observed a strong inhibition in the crystallization buffer and an increase in activity with increasing pH. We attribute this inhibitive effect to the presence of ammonium ions present in the crystal since inhibition is also observed with several other ammonium salt buffers. Additionally, the lack of inhibition was observed in the absence of ammonium. The effect of each ligand (mevalonate, CoA and NAD) on the rate of the enzymatic reaction in the crystallization environment was further investigated by measuring their Kmin the crystallization buffer. The Kmmeasurements indicate that the hydride transfer step between NAD and mevalonate is inhibited in the crystallization environment. To test this further, we solved a crystal structure of pmHMGR bound to the post-hydride transfer intermediate (mevaldehyde) and cofactor Coenzyme A. The resulting turnover with the formation of a thiohemiacetal indicated that the crystallization environment inhibited the oxidative acylation of mevalonate and the reaction intermediate mevaldyl-CoA.

The nucleophilic addition of alkyl, benzyl, and allyl organozinc reagents to protected pyridinium riboses proceeds under mild conditions and with yields of >90% in several cases and complete regioselectivity for the 4-position.

A transition state force field (TSFF) was developed using the quantum-guided molecular mechanics (Q2MM) method to describe the stereodetermining migratory insertion step of the enantioselective redox-relay Heck reaction for a range of multisubstituted alkenes. We show that the TSFF is highly predictive through an external validation of the TSFF against 151 experimentally determined stereoselectivities resulting in an R(2) of 0.89 and MUE of 1.8 kJ/mol. In addition, limitations in the underlying force field were identified by comparison of the TSFF results to DFT level calculations. A novel application of the TSFF was demonstrated for 31 cases where the enantiomer predicted by the TSFF differed from the originally published values. Experimental determination of the absolute configuration demonstrated that the computational predictions were accurate, suggesting that TSFFs can be used for the rapid prediction of the absolute stereochemistry for a class of reactions. Finally, a virtual ligand screen was conducted utilizing both the TSFF and a simple molecular correlation method. Both methods were similarly predictive, but the TSFF was able to show greater utility through transferability, speed, and interpretability.