Rapid Affinity and Microsomal Stability Ranking of Crude Mixture Libraries of Histone Deacetylase Inhibitors.

The science of drug discovery involves multiparameter optimization of molecular structures through iterative design-make-test cycles. For medicinal chemistry library synthesis, traditional workflows involve the isolation of each individual compound, gravimetric quantitation, and preparation of a standard concentration solution for biological assays. In this work, we explore ways to expedite this process by testing unpurified library mixtures using a combination of mass spectrometry-based assays for affinity selection and microsomal metabolic stability. Utilizing this approach, microgram quantities of crude library mixtures can be used to identify high affinity, metabolically stable library members for isolation and full characterization. This streamlined approach was demonstrated for the synthesis and evaluation of two libraries of histone deacetylase inhibitors and was shown to generate decision-making data in line with traditional workflows. The advantages of this paradigm include greatly reduced cycle time, reduced material requirements, and concentration of resources on the most promising compounds.

Strategy for improving circular dichroism spectra deconvolution accuracy for macrocyclic peptides in drug discovery.

Peptide therapeutics have emerged as an appealing modality in the pharmaceutical industry. Understanding peptide conformation in solution remains one of the most critical areas for peptide drug development. Circular dichroism (CD) spectroscopy is a useful technique to study the secondary structure of proteins and peptides, but the current approaches are limited to protein-focused models to predict high-order structures of peptides, and the models were built based on X-ray crystallography instead of solution-based technique, as a result, such models may have poor predictions for peptides. In this study, we present a novel CD deconvolution model to determine peptide conformation in solution. To quantitatively obtain secondary structure information using CD, a calibration model is needed beforehand to establish the relationship between each secondary structure feature and the corresponding CD response. A reference set containing the majority of cyclic peptides with known structures from solution-state NMR spectroscopy was used to build the calibration model for CD deconvolution. Improved prediction accuracy on the secondary structure determination for cyclic peptides was achieved by this model compared to the commercial standard model using commercially available platforms. This new CD deconvolution method is crucial for peptide conformational analysis in solution, and has the potential to greatly accelerate peptide drug candidate optimization in the pharmaceutical drug discovery field.

Cerastecin Inhibition of the Lipooligosaccharide Transporter MsbA to Combat Acinetobacter baumannii: From Screening Impurity to In Vivo Efficacy.

Acinetobacter baumannii, a commonly multidrug-resistant Gram-negative bacterium responsible for large numbers of bloodstream and lung infections worldwide, is increasingly difficult to treat and constitutes a growing threat to human health. Structurally novel antibacterial chemical matter that can evade existing resistance mechanisms is essential for addressing this critical medical need. Herein, we describe our efforts to inhibit the essential A. baumannii lipooligosaccharide (LOS) ATP-binding cassette (ABC) transporter MsbA. An unexpected impurity from a phenotypic screening was optimized as a series of dimeric compounds, culminating with 1 (cerastecin D), which exhibited antibacterial activity in the presence of human serum and a pharmacokinetic profile sufficient to achieve efficacy against A. baumannii in murine septicemia and lung infection models.

CASE-DFT Structure Elucidation of Proton-Deficient Chlorodepsidones from the Indonesian Lichen Teloschistes flavicans and Structure Revision of Flavicansone.

Biodiscovery efforts in Indonesia have aimed to explore the understudied chemical diversity of its rich lichen flora, seeking to find new products endowed with significant biological properties. The chemical screening of a Teloschistes flavicans extract led to selection of this species for further investigation. LC/MS and 1H NMR-based dereplication pinpointed six chlorodepsidones from the thallus of a sample of this lichen. This led to the streamlined isolation and the subsequent structure elucidation of the three new compounds norflavicansone 1, flavicansone 2, and isocaloploicin 3, along with the known chlorodepsidones 4-6, stictic acid 7, aurantiamide acetate 8, and parietin 9. The challenging structure elucidation of these proton-deficient metabolites benefited from a state-of-the-art workflow involving a synergistic combination of Computer-Assisted Structure Elucidation (CASE) and Density Functional Theory (DFT) calculations of the top-ranked candidates. This investigation also led to the revision of flavicansone's structure, previously described from this species. The three new molecules that are being reported here are remarkable in that they represent hybrid depsidones in which one of the aromatic rings is derived from orsellinic acid and the other is derived from ß-orcinol, a rare structural feature for lichen depsidones.

DFT investigation of coupling constant anomalies in substituted ß-lactams.

ß-lactams are a chemically diverse group of molecules with a wide range of biological activities. Having recently observed curious trends in 2JHH coupling values in studies on this structural class, we sought to obtain a more comprehensive understanding of these diagnostic NMR parameters, specifically interrogating 1JCH, 2JCH, and 2JHH, to differentiate 3- and 4-monosubstituted ß-lactams. Further investigation using computational chemistry methods was employed to explore the geometric and electronic origins for the observed and calculated differences between the two substitution patterns.

Unifying the configuration of historical alkaloids from Borreria capitata through an extensive spectroscopic reinvestigation.

Reinvestigation of the structure of borrecapine and borreline through extensive spectroscopic analysis of their authentic samples led to the assignment of their absolute configurations. Newly acquired spectroscopic data determined that the previously assigned relative configuration for borrecapine was incorrect and that the claimed absolute configuration of borreline should be revised to its enantiomer.

Enhancing Efficiency of Natural Product Structure Revision: Leveraging CASE and DFT over Total Synthesis.

Natural products remain one of the major sources of coveted, biologically active compounds. Each isolated compound undergoes biological testing, and its structure is usually established using a set of spectroscopic techniques (NMR, MS, UV-IR, ECD, VCD, etc.). However, the number of erroneously determined structures remains noticeable. Structure revisions are very costly, as they usually require extensive use of spectroscopic data, computational chemistry, and total synthesis. The cost is particularly high when a biologically active compound is resynthesized and the product is inactive because its structure is wrong and remains unknown. In this paper, we propose using Computer-Assisted Structure Elucidation (CASE) and Density Functional Theory (DFT) methods as tools for preventive verification of the originally proposed structure, and elucidation of the correct structure if the original structure is deemed to be incorrect. We examined twelve real cases in which structure revisions of natural products were performed using total synthesis, and we showed that in each of these cases, time-consuming total synthesis could have been avoided if CASE and DFT had been applied. In all described cases, the correct structures were established within minutes of using the originally published NMR and MS data, which were sometimes incomplete or had typos.

DELTA50: A Highly Accurate Database of Experimental 1H and 13C NMR Chemical Shifts Applied to DFT Benchmarking.

Density functional theory (DFT) benchmark studies of 1H and 13C NMR chemical shifts often yield differing conclusions, likely due to non-optimal test molecules and non-standardized data acquisition. To address this issue, we carefully selected and measured 1H and 13C NMR chemical shifts for 50 structurally diverse small organic molecules containing atoms from only the first two rows of the periodic table. Our NMR dataset, DELTA50, was used to calculate linear scaling factors and to evaluate the accuracy of 73 density functionals, 40 basis sets, 3 solvent models, and 3 gauge-referencing schemes. The best performing DFT methodologies for 1H and 13C NMR chemical shift predictions were WP04/6-311++G(2d,p) and ωB97X-D/def2-SVP, respectively, when combined with the polarizable continuum solvent model (PCM) and gauge-independent atomic orbital (GIAO) method. Geometries should be optimized at the B3LYP-D3/6-311G(d,p) level including the PCM solvent model for the best accuracy. Predictions of 20 organic compounds and natural products from a separate probe set had root-mean-square deviations (RMSD) of 0.07 to 0.19 for 1H and 0.5 to 2.9 for 13C. Maximum deviations were less than 0.5 and 6.5 ppm for 1H and 13C, respectively.

Ominoxanthone-The First Xanthone Linearly Fused to a γ-Lactone from Cortinarius ominosus Bidaud Basidiomata. CASE- and DFT-Based Structure Elucidation.

The UHPLC-HRMS analysis of Cortinarius ominosus basidiomata extract revealed that this mushroom accumulated elevated yields of an unreported specialized metabolite. The molecular formula of this unknown compound, C17H10O8, indicated that a challenging structure elucidation lay ahead, owing to its critically low H/C atom ratio. The structure of this new isolate, namely ominoxanthone (1), could not be solved from the interpretation of the usual set of 1D/2D NMR data that conveyed too limited information to afford a single, unambiguous structure. To remedy this, a Computer-Assisted Structure Elucidation (CASE) workflow was used to rank the different possible structure candidates consistent with our scarce spectroscopic data. DFT-based chemical shift calculations on a limited set of top-ranked structures further ascertained the determined structure for ominoxanthone. Although the determined scaffold of ominoxanthone is unprecedented as a natural product, a plausible biosynthetic scenario involving a precursor known from cortinariaceous sources and classical biogenetic reactions could be proposed.

Rapid Evolution of a Fragment-like Molecule to Pan-Metallo-Beta-Lactamase Inhibitors: Initial Leads toward Clinical Candidates.

With the emergence and rapid spreading of NDM-1 and existence of clinically relevant VIM-1 and IMP-1, discovery of pan inhibitors targeting metallo-beta-lactamases (MBLs) became critical in our battle against bacterial infection. Concurrent with our fragment and high-throughput screenings, we performed a knowledge-based search of known metallo-beta-lactamase inhibitors (MBLIs) to identify starting points for early engagement of medicinal chemistry. A class of compounds exemplified by 11, discovered earlier as B. fragilis metallo-beta-lactamase inhibitors, was selected for in silico virtual screening. From these efforts, compound 12 was identified with activity against NDM-1 only. Initial exploration on metal binding design followed by structure-guided optimization led to the discovery of a series of compounds represented by 23 with a pan MBL inhibition profile. In in vivo studies, compound 23 in combination with imipenem (IPM) robustly lowered the bacterial burden in a murine infection model and became the lead for the invention of MBLI clinical candidates.

Lead Optimization to Advance Protease-Activated Receptor-1 Antagonists in Early Discovery.

Vorapaxar is an approved drug for the reduction of thrombotic cardiovascular events in patients with a history of myocardial infarction or with peripheral arterial disease. Subsequent to the discovery of Vorapaxar, medicinal chemistry efforts were continued to identify structurally differentiated leads. Toward this goal, extensive structure-activity relationship studies using a C-ring-truncated version of Vorapaxar culminated in the discovery of three leads, represented as 13, 14, and 23. Among these leads, compound 14 possessed favorable pharmacokinetic properties and an off-target profile, which supported additional profiling in an exploratory rat toxicology study.

Driving to a Better Understanding of Acyl Glucuronide Transformations Using NMR and Molecular Modeling.

Acyl glucuronide (AG) metabolites of carboxylic acid-containing drugs and products of their transformations have long been implicated in drug-induced liver injury (DILI). To inform on the DILI risk arising from AG reactive intermediates, a comprehensive mechanistic study of enzyme-independent AG rearrangements using nuclear magnetic resonance (NMR) and density functional theory (DFT) was undertaken. NMR spectroscopy was utilized for structure elucidation and kinetics measurements of nine rearrangement and hydrolysis products of 1ß-O-acyl glucuronide of ibufenac. To extract rate constants of rearrangement, mutarotation, and hydrolysis from kinetic data, 11 different kinetic models were examined. Model selection and estimated rate constant verification were supported by measurements of H/D kinetic isotope effects. DFT calculations of ground and transition states supported the proposed kinetic mechanisms and helped to explain the unusually fast intramolecular transacylation rates found for some of the intermediates. The findings of the current study reinforce the notion that the short half-life of parent AG and slow hydrolysis rates of AG rearrangement products are the two key factors that can influence the in vivo toxicity of AGs.

19 F-detected dual-optimized inverted 1 JCC 1,n-ADEQUATE.

Modification of the recently reported 19 F-detected 1,1-ADEQUATE experiment that incorporates dual-optimization to selectively invert a wide range of 1 JCC correlations in a 1,n-ADEQUATE experiment is reported. Parameters for the dual-optimization segment of the pulse sequence were modified to accommodate the increased size of 1 JCC homonuclear coupling constants of poly- and perfluorinated molecules relative to protonated molecules to allow broadband inversion of the 1 JCC correlations. The observation and utility of isotope shifts are reported for the first time for 1,1- and 1,n-ADEQUATE correlations.

Development of ProTx-II Analogues as Highly Selective Peptide Blockers of Nav1.7 for the Treatment of Pain.

Inhibitor cystine knot peptides, derived from venom, have evolved to block ion channel function but are often toxic when dosed at pharmacologically relevant levels in vivo. The article describes the design of analogues of ProTx-II that safely display systemic in vivo blocking of Nav1.7, resulting in a latency of response to thermal stimuli in rodents. The new designs achieve a better in vivo profile by improving ion channel selectivity and limiting the ability of the peptides to cause mast cell degranulation. The design rationale, structural modeling, in vitro profiles, and rat tail flick outcomes are disclosed and discussed.

Synthesis of HDAC Inhibitor Libraries via Microscale Workflow.

An integrated workflow has been established that enables the synthesis, purification, and subsequent biological testing of compound libraries on a microgram scale. This approach utilizes mass directed preparative HPLC in conjunction with charged aerosol detection (CAD) to generate solutions of investigational compounds at high purity and standardized concentrations, facilitating high fidelity biological testing. This new workflow successfully delivered libraries of histone deacetylase (HDAC) inhibitors that afforded biological data consistent with that obtained from standard scale parallel medicinal chemistry techniques. The advantages of this new approach to library synthesis include greatly reduced material requirements and amenability to high-throughput experimentation.

Development of 19 F-detected 1,1-ADEQUATE for the characterization of polyfluorinated and perfluorinated compounds.

Polyfluorinated and perfluorinated compounds in the environment are a growing health concern. 19 F-detected variants of commonly employed heteronuclear shift correlation experiments such as heteronuclear single quantum correlation (HSQC) and heteronuclear multiple bond correlation (HMBC) are available; 19 F-detected experiments that employ carbon-carbon homonuclear coupling, in contrast, have never been reported. Herein, we report the measurement of the 1 JCC and n JCC coupling constants of a simple perfluorinated phthalonitrile and the first demonstration of a 19 F-detected 1,1-ADEQUATE experiment.

An orally available non-nucleotide STING agonist with antitumor activity.

Pharmacological activation of the STING (stimulator of interferon genes)-controlled innate immune pathway is a promising therapeutic strategy for cancer. Here we report the identification of MSA-2, an orally available non-nucleotide human STING agonist. In syngeneic mouse tumor models, subcutaneous and oral MSA-2 regimens were well tolerated and stimulated interferon-ß secretion in tumors, induced tumor regression with durable antitumor immunity, and synergized with anti-PD-1 therapy. Experimental and theoretical analyses showed that MSA-2 exists as interconverting monomers and dimers in solution, but only dimers bind and activate STING. This model was validated by using synthetic covalent MSA-2 dimers, which were potent agonists. Cellular potency of MSA-2 increased upon extracellular acidification, which mimics the tumor microenvironment. These properties appear to underpin the favorable activity and tolerability profiles of effective systemic administration of MSA-2.

Molecular packing of pharmaceuticals analyzed with paramagnetic relaxation enhancement and ultrafast magic angle pinning NMR.

Understanding the relationship between the structure and the physicochemical attributes of crystalline pharmaceuticals requires high-resolution molecular details. Solid-state nuclear magnetic resonance (ssNMR) spectroscopy is an indispensable tool for analyzing molecular structures, but often experiences challenges of low spectral resolution and sensitivity, particularly in the characterization of unlabeled pharmaceutical materials. Besides, the relatively long spin-lattice relaxation times in pharmaceutical crystals result in time-consuming data collections. In this study, we utilize ultrafast magic angle spinning (UF-MAS) of the sample at 60 and 110 kHz to enable proton and fluorine spectroscopies for probing the structural details of crystalline posaconazole. Paramagnetic relaxation enhancement (PRE), obtained by doping Cu(ii) ions into the crystalline lattice and coating on particle surface, is implemented to shorten the spin-lattice relaxation time for speeding up the ssNMR acquisition. Our results demonstrate a remarkably improved 1H and 19F resolution and sensitivity, which enables multi-dimensional 1H-1H and heteronuclear 1H-19F correlations. In combination with density functional theory (DFT) calculations of chemical shifts, molecular details of posaconazole are established in terms of 1H and 19F networks for identifying "head-to-tail" and "head-to-head" intermolecular packings, with presumably critical contacts that stabilize the crystalline structure. The PRE and UF-MAS techniques enable the high-resolution structure characterization of fluorinated drug molecules in pharmaceutical formulations at natural abundance.

Enhancing computer-assisted structure elucidation with DFT analysis of J-couplings.

Computer-assisted structure elucidation (CASE) is the class of expert systems that derives molecular structures primarily from one-dimensional and two-dimensional nuclear magnetic resonance data. Contemporary CASE systems, including Advanced Chemistry Development/Structure Elucidator (ACD/SE), consider cross-peaks in heteronuclear multiple bond coherence (HMBC) and correlation spectroscopy (COSY) spectra as two- or three-bond correlations by default. However, four and more bond correlations (nonstandard correlations [NSCs]) could be present in these spectra too. The indiscriminate addition of NSCs to the CASE computations is prohibitively expensive. To address this problem, the ACD/SE program performs a logical analysis of observed correlations and determines the minimum number of NSCs. Guided by this information, a more efficient fuzzy structure generation (FSG) algorithm is subsequently applied. Until now, the FSG algorithm was utilized without any verification of the reliability of found NSCs. Here, we report a verification method for NSCs based on the relationship between NSCs and J-couplings computed with high accuracy density functional theory (DFT) methods. We used the example of strychnine to show that 41 (32%) of 8-Hz HMBC cross-peaks were NSCs and were consistent with 4-6 JCH couplings greater than 0.3 Hz. This cutoff value was largely confirmed by the analysis of NSCs in 11 real-world natural products elucidated by ACD/SE. Additionally, utilizing the example of the CASE study of cleospinol A, we showed that the DFT-computed J-couplings of NSCs can distinctively differentiate the correct structure among six proposed isomers. The proposed approach of NSC verification should further improve the robustness of CASE analysis and can help reveal potential problems with reported experimental data.

Bioactivation of GPR40 Agonist MK-8666: Formation of Protein Adducts in Vitro from Reactive Acyl Glucuronide and Acyl CoA Thioester.

MK-8666, a selective GPR40 agonist developed for the treatment of type 2 diabetes mellitus, was discontinued in phase I clinical trials due to liver safety concerns. To address whether chemically reactive metabolites played a causative role in the observed drug induced liver injury (DILI), we characterized the metabolism, covalent binding to proteins, and amino acid targets of MK-8666 in rat and human hepatocytes or cofactor-fortified liver microsomes. MK-8666 was primarily metabolized to an acyl glucuronide in hepatocytes of both species and a taurine conjugate in rat hepatocytes. Similar levels of covalent binding to proteins were observed in rat and human hepatocytes following incubation with [3H]MK-8666. After protease digestion of hepatocyte pellets, amino acid adducts A1, A2, and A3 were identified as transacylated products with lysine, serine, and cysteine residues, respectively. Amino acid adducts A4a-c were identified as glycation adducts resulting from rearrangement of MK-8666-1-O-ß-acyl glucuronide to ring-opened aldehydes which further condensed with lysine residues of proteins into imine adducts. Adducts A1-A3 and A4a-c were detected in rat and human liver microsomes fortified with UDPGA. Adducts A1-A3 were detected in rat and human liver microsomes fortified with CoA and ATP. Additionally, a trace amount of CoA thioester metabolite of MK-8666 and its transacylated GSH adduct were detected in human liver microsomes fortified with CoA, ATP, and GSH. Higher levels of covalent binding to protein were observed when [3H]MK-8666 was incubated in liver microsomes supplemented with CoA and ATP compared to UDPGA. Addition of GSH attenuated levels of CoA thioester-mediated covalent binding by 41-45%. Collectively, these studies indicated that metabolism of the -COOH moiety of MK-8666 can form a reactive acyl glucuronide and an acyl CoA thioester, which covalently modifies proteins and may represent one causative mechanism of the observed DILI.