Retrospective assessment of rat liver microsomal stability at NCATS: data and QSAR models.

Hepatic metabolic stability is a key pharmacokinetic parameter in drug discovery. Metabolic stability is usually assessed in microsomal fractions and only the best compounds progress in the drug discovery process. A high-throughput single time point substrate depletion assay in rat liver microsomes (RLM) is employed at the National Center for Advancing Translational Sciences. Between 2012 and 2020, RLM stability data was generated for ~ 24,000 compounds from more than 250 projects that cover a wide range of pharmacological targets and cellular pathways. Although a crucial endpoint, little or no data exists in the public domain. In this study, computational models were developed for predicting RLM stability using different machine learning methods. In addition, a retrospective time-split validation was performed, and local models were built for projects that performed poorly with global models. Further analysis revealed inherent medicinal chemistry knowledge potentially useful to chemists in the pursuit of synthesizing metabolically stable compounds. In addition, we deposited experimental data for ~ 2500 compounds in the PubChem bioassay database (AID: 1508591). The global prediction models are made publicly accessible ( https://opendata.ncats.nih.gov/adme ). This is to the best of our knowledge, the first publicly available RLM prediction model built using high-quality data generated at a single laboratory.

The Critical Role of Passive Permeability in Designing Successful Drugs.

Passive permeability is a key property in drug disposition and delivery. It is critical for gastrointestinal absorption, brain penetration, renal reabsorption, defining clearance mechanisms and drug-drug interactions. Passive diffusion rate is translatable across tissues and animal species, while the extent of absorption is dependent on drug properties, as well as in vivo physiology/pathophysiology. Design principles have been developed to guide medicinal chemistry to enhance absorption, which combine the balance of aqueous solubility, permeability and the sometimes unfavorable compound characteristic demanded by the target. Permeability assays have been implemented that enable rapid development of structure-permeability relationships for absorption improvement. Future advances in assay development to reduce nonspecific binding and improve mass balance will enable more accurately measurement of passive permeability. Design principles that integrate potency, selectivity, passive permeability and other ADMET properties facilitate rapid advancement of successful drug candidates to patients.

Quantification of cyclocreatine in mouse and rat plasma using hydrophilic-interaction ultra-performance liquid chromatography-tandem mass spectrometry.

An accurate, rapid and selective method was developed to quantify cyclocreatine in mouse and rat plasma using hydrophilic interaction (HILIC) ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). The plasma samples were prepared by protein precipitation with acetonitrile:methanol (70:30). Chromatographic separation was performed on a HILIC BEH amide column (2.1mm×50mm, 1.7µm) with a 3min gradient elution at a flow rate of 0.5mL/min. For mass spectrometric detection, selected reaction monitoring (SRM) was used; the SRM transitions were m/z 144→98 and m/z 144→56 for cyclocreatine and m/z 148→102 for the internal standard (D4-cyclocreatine) in the positive ionization mode. No endogenous components interfered with the analysis of cyclocreatine and the internal standard in mouse and rat plasma. Plasma calibration curves were constructed in the range of 0.01-25µM. The correlation coefficient of the calibration curves was greater than 0.99. The mean intraday assay accuracy for all quality control (QC) replicates was between 93 and 105%. The mean intraday assay precision (CV%) was 1.9-11% for all QC levels. The HILIC-UPLC-MS/MS method was successfully applied in pharmacokinetic (PK) studies of cyclocreatine in mice and rats for the first time. After a single 30mg/kg oral administration in mice and rats, the AUC0-∞ (area under the curve) was 84.1µgh/mL and 91.7±18.0µgh/mL, respectively.

Highly predictive and interpretable models for PAMPA permeability.

Cell membrane permeability is an important determinant for oral absorption and bioavailability of a drug molecule. An in silico model predicting drug permeability is described, which is built based on a large permeability dataset of 7488 compound entries or 5435 structurally unique molecules measured by the same lab using parallel artificial membrane permeability assay (PAMPA). On the basis of customized molecular descriptors, the support vector regression (SVR) model trained with 4071 compounds with quantitative data is able to predict the remaining 1364 compounds with the qualitative data with an area under the curve of receiver operating characteristic (AUC-ROC) of 0.90. The support vector classification (SVC) model trained with half of the whole dataset comprised of both the quantitative and the qualitative data produced accurate predictions to the remaining data with the AUC-ROC of 0.88. The results suggest that the developed SVR model is highly predictive and provides medicinal chemists a useful in silico tool to facilitate design and synthesis of novel compounds with optimal drug-like properties, and thus accelerate the lead optimization in drug discovery.

An Automated High-Throughput Metabolic Stability Assay Using an Integrated High-Resolution Accurate Mass Method and Automated Data Analysis Software.

Advancement of in silico tools would be enabled by the availability of data for metabolic reaction rates and intrinsic clearance (CLint) of a diverse compound structure data set by specific metabolic enzymes. Our goal is to measure CLint for a large set of compounds with each major human cytochrome P450 (P450) isozyme. To achieve our goal, it is of utmost importance to develop an automated, robust, sensitive, high-throughput metabolic stability assay that can efficiently handle a large volume of compound sets. The substrate depletion method [in vitro half-life (t1/2) method] was chosen to determine CLint The assay (384-well format) consisted of three parts: 1) a robotic system for incubation and sample cleanup; 2) two different integrated, ultraperformance liquid chromatography/mass spectrometry (UPLC/MS) platforms to determine the percent remaining of parent compound, and 3) an automated data analysis system. The CYP3A4 assay was evaluated using two long t1/2 compounds, carbamazepine and antipyrine (t1/2 > 30 minutes); one moderate t1/2 compound, ketoconazole (10 < t1/2 < 30 minutes); and two short t1/2 compounds, loperamide and buspirone (t½ < 10 minutes). Interday and intraday precision and accuracy of the assay were within acceptable range (∼12%) for the linear range observed. Using this assay, CYP3A4 CLint and t1/2 values for more than 3000 compounds were measured. This high-throughput, automated, and robust assay allows for rapid metabolic stability screening of large compound sets and enables advanced computational modeling for individual human P450 isozymes.

Synthesis and structure-activity relationship studies of N-benzyl-2-phenylpyrimidin-4-amine derivatives as potent USP1/UAF1 deubiquitinase inhibitors with anticancer activity against nonsmall cell lung cancer.

Deregulation of ubiquitin conjugation or deconjugation has been implicated in the pathogenesis of many human diseases including cancer. The deubiquitinating enzyme USP1 (ubiquitin-specific protease 1), in association with UAF1 (USP1-associated factor 1), is a known regulator of DNA damage response and has been shown as a promising anticancer target. To further evaluate USP1/UAF1 as a therapeutic target, we conducted a quantitative high throughput screen of >400000 compounds and subsequent medicinal chemistry optimization of small molecules that inhibit the deubiquitinating activity of USP1/UAF1. Ultimately, these efforts led to the identification of ML323 (70) and related N-benzyl-2-phenylpyrimidin-4-amine derivatives, which possess nanomolar USP1/UAF1 inhibitory potency. Moreover, we demonstrate a strong correlation between compound IC50 values for USP1/UAF1 inhibition and activity in nonsmall cell lung cancer cells, specifically increased monoubiquitinated PCNA (Ub-PCNA) levels and decreased cell survival. Our results establish the druggability of the USP1/UAF1 deubiquitinase complex and its potential as a molecular target for anticancer therapies.

Passive lipoidal diffusion and carrier-mediated cell uptake are both important mechanisms of membrane permeation in drug disposition.

Recently, it has been proposed that drug permeation is essentially carrier-mediated only and that passive lipoidal diffusion is negligible. This opposes the prevailing hypothesis of drug permeation through biological membranes, which integrates the contribution of multiple permeation mechanisms, including both carrier-mediated and passive lipoidal diffusion, depending on the compound's properties, membrane properties, and solution properties. The prevailing hypothesis of drug permeation continues to be successful for application and prediction in drug development. Proponents of the carrier-mediated only concept argue against passive lipoidal diffusion. However, the arguments are not supported by broad pharmaceutics literature. The carrier-mediated only concept lacks substantial supporting evidence and successful applications in drug development.

Synthesis and structure-activity relationship studies of 4-((2-hydroxy-3-methoxybenzyl)amino)benzenesulfonamide derivatives as potent and selective inhibitors of 12-lipoxygenase.

Human lipoxygenases (LOXs) are a family of iron-containing enzymes which catalyze the oxidation of polyunsaturated fatty acids to provide the corresponding bioactive hydroxyeicosatetraenoic acid (HETE) metabolites. These eicosanoid signaling molecules are involved in a number of physiologic responses such as platelet aggregation, inflammation, and cell proliferation. Our group has taken a particular interest in platelet-type 12-(S)-LOX (12-LOX) because of its demonstrated role in skin diseases, diabetes, platelet hemostasis, thrombosis, and cancer. Herein, we report the identification and medicinal chemistry optimization of a 4-((2-hydroxy-3-methoxybenzyl)amino)benzenesulfonamide-based scaffold. Top compounds, exemplified by 35 and 36, display nM potency against 12-LOX, excellent selectivity over related lipoxygenases and cyclooxygenases, and possess favorable ADME properties. In addition, both compounds inhibit PAR-4 induced aggregation and calcium mobilization in human platelets and reduce 12-HETE in ß-cells.

Evidence-based approach to assess passive diffusion and carrier-mediated drug transport.

Evidence supporting the action of passive diffusion and carrier-mediated (CM) transport in drug bioavailability and disposition is discussed to refute the recently proposed theory that drug transport is CM-only and that new transporters will be discovered that possess transport characteristics ascribed to passive diffusion. Misconceptions and faulty speculations are addressed to provide reliable guidance on choosing appropriate tools for drug design and optimization.

Identification of 3-sulfonylindazole derivatives as potent and selective 5-HT(6) antagonists.

As part of our efforts to develop agents for cognitive enhancement, we have been focused on the 5-HT(6) receptor in order to identify potent and selective ligands for this purpose. Herein we report the identification of a novel series of 3-sulfonylindazole derivatives with acyclic amino side chains as potent and selective 5-HT(6) antagonists. The synthesis and detailed SAR of this class of compounds are reported.

The effect of plasma protein binding on in vivo efficacy: misconceptions in drug discovery.

Data from in vitro plasma protein binding experiments that determine the fraction of protein-bound drug are frequently used in drug discovery to guide structure design and to prioritize compounds for in vivo studies. However, we consider that these practices are usually misleading, because in vivo efficacy is determined by the free (unbound) drug concentration surrounding the therapeutic target, not by the free drug fraction. These practices yield no enhancement of the in vivo free drug concentration. So, decisions based on free drug fraction could result in the wrong compounds being advanced through drug discovery programmes. This Perspective provides guidance on the application of plasma protein binding information in drug discovery.

Insights for predicting blood-brain barrier penetration of CNS targeted molecules using QSPR approaches.

Due to the high attrition rate of central nervous system drug candidates during clinical trials, the assessment of blood-brain barrier (BBB) penetration in early research is particularly important. A genetic approximation (GA)-based regression model was developed for predicting in vivo blood-brain partitioning data, expressed as logBB (log[brain]/[blood]). The model was built using an in-house data set of 193 compounds assembled from 22 different therapeutic projects. The final model (cross-validated r(2) = 0.72) with five molecular descriptors was selected based on validation using several large internal and external test sets. We demonstrate the potential utility of the model by applying it to a set of literature reported secretase inhibitors. In addition, we describe a rule-based approach for rapid assessment of brain penetration with several simple molecular descriptors.

Simultaneous determination of 2-arachidonoylglycerol, 1-arachidonoylglycerol and arachidonic acid in mouse brain tissue using liquid chromatography/tandem mass spectrometry.

Endocannabinoids (ECs), such as anandamide (AEA) and 2-arachidonoylglycerol (2-AG), modulate a number of physiological processes, including pain, appetite and emotional state. Levels of ECs are tightly controlled by enzymatic biosynthesis and degradation in vivo. However, there is limited knowledge about the enzymes that terminate signaling of the major brain EC, 2-AG. Identification and quantification of 2-AG, 1-AG and arachidonic acid (AA) is important for studying the enzymatic hydrolysis of 2-AG. We have developed a sensitive and specific quantification method for simultaneous determination of 2-AG, 1-AG and AA from mouse brain and adipose tissues by liquid chromatography/tandem mass spectrometry (LC/MS/MS) using a simple brain sample preparation method. The separations were carried out based on reversed phase chromatography. Optimization of electrospray ionization conditions established the limits of detection (S/N = 3) at 50, 25 and 65 fmol for 2-AG, 1-AG and AA, respectively. The methods were selective, precise (%R.S.D. < 10%) and sensitive over a range of 0.02-20, 0.01-10 and 0.05-50 ng/mg tissue for 2-AG, 1-AG and AA, respectively. The quantification method was validated with consideration of the matrix effects and the mass spectrometry (MS) responses of the analytes and the deuterium labeled internal standard (IS). The developed methods were applied to study the hydrolysis of 2-AG from mouse brain extracts containing membrane bound monoacylglycerol lipase (MAGL), and to measure the basal levels of 2-AG, 1-AG and AA in mouse brain and adipose tissues.

Stability challenges in drug discovery.

Stability is one of the most important properties of drug candidates. Instable compounds can lead to false positive high-throughput screening (HTS) hits, incorrect bioassay results, erroneous structure-activity relationships (SAR), low oral bioavailability, drug withdrawal, toxic reactions from degradation products, and difficult formulation development. Screening of stability has been implemented early in drug discovery to identify labile chemotypes and guide structural modification. The most commonly applied stability studies in drug discovery are stability-pH profile, stability in gastrointestinal fluids, stability in bioassay media, excipient compatibility, and prodrug screening. The strategy enhances the quality of drug development candidates and reduces the risks.

Drug-like property concepts in pharmaceutical design.

The pharmaceutical industry is facing an ever increasing challenge to deliver safer and more effective medicines. Traditionally, drug discovery programs were driven solely by potency, regardless of the properties. As a result, the development of non-drug-like molecules was costly, had high risk and low success rate. To meet the challenges, the bar has been rising higher for drug candidates. They not only need to be active, but also drug-like to be advanced to clinical development. Drug-like properties, such as solubility, permeability, metabolic stability and transporter effects are of critical importance for the success of drug candidates. They affect oral bioavailability, metabolism, clearance, toxicity, as well as in vitro pharmacology. Insoluble and impermeable compounds can result in erroneous biological data and unreliable SAR in enzyme and cell-based assays. Rapid metabolism by enzymes and high efflux by transporters can lead to high clearance, short half-life, low systemic exposure and inadequate efficacy. Early property information helps teams make informed decisions and avoids wasting precious resources. Structure-property relationships are essential to guide structural modification to improve properties. High throughput ADME/TOX assays have been implemented and are being widely used to drive drug discovery projects in parallel with activity screening. Property design has become an integrated and inseparable part of the modern drug discovery paradigm. The approach has been proven to be a winning strategy.

Identification of a novel series of 3-piperidinyl-5-sulfonylindazoles as potent 5-HT6 ligands.

Cognitive dysfunction is a characteristic of various forms of dementia such as Alzheimer's disease (AD) and a core feature of schizophrenia. As part of our continuing efforts to develop agents for cognitive enhancement, we have been focused on the 5-HT(6) receptor-one of the emerging therapeutic targets in this area. Herein, we report the identification of a novel series of 3-piperidinyl-5-sulfonylindazole derivatives as potent 5-HT(6) antagonists. The synthesis and SAR of this class of compounds are reported.

1-Sulfonylindazoles as potent and selective 5-HT6 ligands.

As part of our continuing efforts to identify therapeutics for CNS diseases, such as schizophrenia and Alzheimer's disease (AD), we have been focused on the 5-HT(6) receptor in an attempt to identify ligands as a potential treatment for cognitive dysfunction. Herein we report the identification of a novel series of 1-sulfonylindazole derivatives as potent and selective 5-HT(6) antagonists. The synthesis and SAR of this class of compounds are reported. Several potent compounds in both binding and cyclase functional assays also display good selectivity, microsomal stability, solubility, and brain penetration as well as low cytochrome P450 inhibition. One compound exemplified in this series showed 24% oral bioavailability and in vivo efficacy in a NOR cognition model at 10mg/kg following an oral administration in rats.

Identification of a series of benzoxazoles as potent 5-HT6 ligands.

As part of our continuing efforts to identify therapeutics for CNS diseases such as schizophrenia and Alzheimer's disease (AD), we have been focused on the 5-HT(6) receptor in order to identify potent and selective ligands as a potential treatment for cognitive dysfunction. Herein we report the identification of a novel series of benzoxazole derivatives as potent 5-HT(6) ligands. The synthesis and detailed SAR of this class of compounds are reported. The compounds have been shown to be full antagonists in a cyclic AMP functional assay.

Discovery of (2S,4R)-1-(2-aminoacetyl)-4-benzamidopyrrolidine-2-carboxylic acid hydrochloride (GAP-134)13, an orally active small molecule gap-junction modifier for the treatment of atrial fibrillation.

Rotigaptide (3) is an antiarrhythmic peptide that improves cardiac conduction by modifying gap-junction communication. Small molecule gap-junction modifiers with improved physical properties were identified from a Zealand Pharma peptide library using pharmaceutical profiling, established SAR around 3, and a putative pharmacophore model for rotigaptide. Activity of the compounds was confirmed in a mouse cardiac conduction block model of arrhythmia. Dipeptide 9f (GAP-134) was identified as a potent, orally active gap-junction modifier for clinical development.