A chemogenomic analysis of ionization constants--implications for drug discovery

Kinetically-Controlled Ni-Catalyzed Direct Carboxylation of Unactivated Secondary Alkyl Bromides without Chain Walking

Herein, we report the direct carboxylation of unactivated secondary alkyl bromides enabled by the merger of photoredox and nickel catalysis, a previously inaccessible endeavor in the carboxylation arena. Site-selectivity is dictated by a kinetically controlled insertion of CO2 at the initial C(sp3)-Br site by the rapid formation of Ni(I)-alkyl species, thus avoiding undesired β-hydride elimination and chain-walking processes. Preliminary mechanistic experiments reveal the subtleties of stereoelectronic effects for guiding the reactivity and site-selectivity.

Protolytic equilibria of ACE inhibitors in micellar solution of nonionic surfactant Brij 35

The acid-base equilibria of six ACE inhibitors (ACEIs), captopril, cilazapril, enalapril, lisinopril, quinapril, and ramipril, were investigated in the presence of micelles of nonionic surfactant Brij 35. The pKa values were potentiometrically determined at 25 °C and at a constant ionic strength (0.1 M NaCl). The obtained potentiometric data were evaluated in the computer program Hyperquad. On the basis of the shift in the pKa values (ΔpKa) determined in micellar media in relation to the pKa values previously determined in "pure" water, the effect of Brij 35 micelles on ACEIs ionization was estimated. The presence of nonionic Brij 35 micelles caused a shift in the pKa values of all ionizable groups of the investigated ACEIs (ΔpKa from - 3.44 to + 1.9) while shifting the protolytic equilibria of both acidic and basic groups toward the molecular form. The Brij 35 micelles expressed the most pronounced effect on the ionization of captopril among the investigated ACEIs and stronger effect on the ionization of amino than on the ionization of carboxyl groups. The obtained results suggest that ionizable functional groups of ACEIs are involved in interactions with palisade layer of nonionic Brij 35 micelles, which potentially can be considered in physiological conditions. Distribution diagrams of the investigated ACEIs equilibrium forms as a function of pH indicate that the change in distribution is most strongly expressed in pH range 4-8, which includes biopharmaceutically important pH values.

Graphical abstract

Natural Lipid Extracts as an Artificial Membrane for Drug Permeability Assay: In Vitro and In Silico Characterization

In vitro non-cellular permeability models such as the parallel artificial membrane permeability assay (PAMPA) are widely applied tools for early-phase drug candidate screening. In addition to the commonly used porcine brain polar lipid extract for modeling the blood–brain barrier’s permeability, the total and polar fractions of bovine heart and liver lipid extracts were investigated in the PAMPA model by measuring the permeability of 32 diverse drugs. The zeta potential of the lipid extracts and the net charge of their glycerophospholipid components were also determined. Physicochemical parameters of the 32 compounds were calculated using three independent forms of software (Marvin Sketch, RDKit, and ACD/Percepta). The relationship between the lipid-specific permeabilities and the physicochemical descriptors of the compounds was investigated using linear correlation, Spearman correlation, and PCA analysis. While the results showed only subtle differences between total and polar lipids, permeability through liver lipids highly differed from that of the heart or brain lipid-based models. Correlations between the in silico descriptors (e.g., number of amide bonds, heteroatoms, and aromatic heterocycles, accessible surface area, and H-bond acceptor–donor balance) of drug molecules and permeability values were also found, which provides support for understanding tissue-specific permeability.

A New Unnatural Amino Acid Derived from the Modification of 4′-(p-tolyl)-2,2′:6′,2″-terpyridine and Its Mixed-Ligand Complexes with Ruthenium: Synthesis, Characterization, and Photophysical Properties

The modification of the methyl group of 4′-(p-tolyl)-2,2′:6′,2″-terpyridine produced the novel unnatural amino acid 3-(4-([2,2′:6′,2″-terpyridin]-4′-yl)phenyl)-2-aminopropanoic acid (phet). Mononuclear heteroleptic ruthenium complexes of the general formulae [Ru(L1)(L2)](PF6)2 (L1 = 2-acetylamino-2-(4-[2,2′:6′,2″]terpyridine-4′-yl-benzyl)-malonic acid diethyl ester, (phem), 3-(4-([2,2′:6′,2″-terpyridin]-4′-yl)phenyl)-2-aminopropanoic acid, (phet), and L2 = 2,2′:6′,2″-terpyridine (tpy), 4′-phenyl-2,2′:6′,2″-terpyridine (ptpy), 4′-(p-tolyl)-2,2′:6′,2″-terpyridine (mptpy)), as well as the homoleptic [Ru(phem)2](PF6)2 and [Ru(phet)2](PF6)2, were synthesized and characterized by means of NMR spectroscopic techniques, elemental analysis, and high-resolution mass spectrometry. The photophysical properties of the synthesized complexes were also studied.

Trends in small molecule drug properties: A developability molecule assessment perspective

Developability molecule assessment is a key interfacial capability across the biopharmaceutical industry, screening and staging molecules discovered by medicinal chemists for successful chemistry manufacturing controls (CMC) development and launch. The breadth of responsibility and expertise such teams possess puts them in a unique position to understand the impact of the physicochemical properties of a drug during its initial discovery and subsequent development. However, most of the publications describing trends in physicochemical properties are written from a medicinal chemistry perspective with the aim to identify molecules with better ADMET profiles that are either lead-like or drug-like, failing to describe the impact these properties have on CMC development. To systematically uncover knowledge obtained from recent trends in physicochemical properties and the corresponding impact on CMC development, a comprehensive analysis was conducted on molecules in the drug repurposing hub dataset. The only physicochemical property that seems to have been preserved in FDA-approved oral molecules over the decades (1900-2020) is a constant H-bond donor count, highlighting the importance this property has on cell permeability and lattice energy. Pharmaceutical attrition analysis suggests that partition-distribution coefficient, H-bond acceptors, polar surface area and the fraction of sp³ carbons are properties that are associated with compound attrition. Looking at pharmaceutical attrition asynchronously with the temporal analysis of FDA-approved oral molecules highlights the opposing trends, risks and diminishing effects some of these physiochemical properties (cLogP, cLogD and Fsp3) have on describing compound attrition during the past decade. Trellising the dataset by target class suggests that certain formulation and drug delivery strategies can be anticipated or put into place based on target class of a molecule. For example, molecules binding to nuclear hormone receptors are amenable to lipid-based drug delivery systems with proven commercial success. Although the poor solubility of kinase inhibitors is a combination of hydrophobicity (due to aromaticity) required to bind to its target and high lattice energy (melting point), they are a challenging target class to formulate. The influence of drug targets on physicochemical properties and the temporal nature of these properties is highlighted when comparing molecules in the drug repurposing dataset to those developed at Amgen. An improved understanding of the impact of molecular properties on performance attributes can accelerate decisions and facilitate risk assessments during candidate selection and development.

Applications of Isosteres of Piperazine in the Design of Biologically Active Compounds: Part 2

Applications of piperazine and homopiperazine in drug design are well-established, and these heterocycles have found use as both scaffolding and terminal elements and also as a means of introducing a water-solubilizing element into a molecule. In the accompanying review (10.1021/acs.jafc.2c00726), we summarized applications of piperazine and homopiperazine and their fused ring homologues in bioactive compound design along with illustrations of the use of 4-substituted piperidines and a sulfoximine-based mimetic. In this review, we discuss applications of pyrrolidine- and fused-pyrrolidine-based mimetics of piperazine and homopiperazine and illustrate derivatives of azetidine that include stretched and spirocyclic motifs, along with applications of a series of diaminocycloalkanes.

Direct Aerobic α-Hydroxylation of Arylacetates for the Synthesis of Mandelates

Aerobic α-hydroxylation of α-methylene esters has proven challenging due to overoxidation and hydrolysis of the materials. In this article, KO^tBu-promoted TBAB-catalyzed α-hydroxylation of α-methylene aryl esters using O₂ as the oxygen source has been developed. Both low reaction temperature and catalyst TBAB are keys to success. This reaction provides an environmentally friendly and low-cost approach to mandelates, which are valuable building blocks and widely present in pharmaceuticals.

Controlled meta-Selective C-H Mono- and Di-Olefination of Mandelic Acid Derivatives

Mandelic acids represent a key structural motif present in many drug molecules. Herein, we report the controlled meta-selective mono- and diolefination of mandelic acids by the careful design of the substrate and oxidant. Furthermore, free meta-functionalized mandelic acid was generated by selectively removing the template under mild basic conditions. The synthesis of functionalized homatropine and cyclandelate drug derivatives was demonstrated. Kinetic isotope effects revealed C-H activation as the rate-limiting step.

The Acid/Base Characterization of Molecules with Epigenetic Activity

The acidic and basic functional groups in a molecule strongly influence its physicochemical properties, affinity for a macromolecule, pharmacokinetics, and toxicity. For instance, basicity has been correlated with molecular promiscuity, hERG blockade, and phospholipidosis. Nonetheless, no systematic characterization of the acid/base profile of epigenomic databases has been reported. This study describes an analysis of the acidic ionization constant distribution of a library of 7820 compounds with reported activity against epigenetic targets. Furthermore, the epigenomics database's acid/base profile was compared to the reference libraries of food chemicals, natural products, and approved drugs. It was found that the acid/base profile of histone lysine demethylase ligands is more similar to previously approved drugs, and histone acetyltransferase ligands have acidic and basic functional groups largely found in food chemicals and natural products; this support the potential of these libraries for finding new epigenetic inhibitors.

Into the Fray! A Beginner's Guide to Medicinal Chemistry

Modern medicinal chemistry is a complex, multidimensional discipline that operates at the interface of the chemical and biological sciences. The medicinal chemistry contribution to drug discovery is typically described in the context of the well-recited linear progression of the drug discovery pipeline. However, compound optimization is idiosyncratic to each project, and clear definitions of hit and lead molecules and the subsequent progress along the pipeline becomes easily blurred. In addition, this description lacks insight into the entangled relationship between chemical and pharmacological properties, and thus provides limited guidance on how innovative medicinal chemistry strategies can be applied to solve optimization problems, regardless of the stage in the pipeline. Through discussion and illustrative examples, this article seeks to provide insights into the finesse of medicinal chemistry and the subtlety of balancing chemical properties pharmacology. In so doing, it aims to serve as an accessible and simple-to-digest guide for anyone who wishes to learn about the underlying principles of medicinal chemistry, in a context that has been decoupled from the pipeline description.

Overview of the SAMPL6 pKa challenge: evaluating small molecule microscopic and macroscopic pKa predictions

The prediction of acid dissociation constants (pKa) is a prerequisite for predicting many other properties of a small molecule, such as its protein-ligand binding affinity, distribution coefficient (log D), membrane permeability, and solubility. The prediction of each of these properties requires knowledge of the relevant protonation states and solution free energy penalties of each state. The SAMPL6 pKa Challenge was the first time that a separate challenge was conducted for evaluating pKa predictions as part of the Statistical Assessment of Modeling of Proteins and Ligands (SAMPL) exercises. This challenge was motivated by significant inaccuracies observed in prior physical property prediction challenges, such as the SAMPL5 log D Challenge, caused by protonation state and pKa prediction issues. The goal of the pKa challenge was to assess the performance of contemporary pKa prediction methods for drug-like molecules. The challenge set was composed of 24 small molecules that resembled fragments of kinase inhibitors, a number of which were multiprotic. Eleven research groups contributed blind predictions for a total of 37 pKa distinct prediction methods. In addition to blinded submissions, four widely used pKa prediction methods were included in the analysis as reference methods. Collecting both microscopic and macroscopic pKa predictions allowed in-depth evaluation of pKa prediction performance. This article highlights deficiencies of typical pKa prediction evaluation approaches when the distinction between microscopic and macroscopic pKas is ignored; in particular, we suggest more stringent evaluation criteria for microscopic and macroscopic pKa predictions guided by the available experimental data. Top-performing submissions for macroscopic pKa predictions achieved RMSE of 0.7-1.0 pKa units and included both quantum chemical and empirical approaches, where the total number of extra or missing macroscopic pKas predicted by these submissions were fewer than 8 for 24 molecules. A large number of submissions had RMSE spanning 1-3 pKa units. Molecules with sulfur-containing heterocycles or iodo and bromo groups were less accurately predicted on average considering all methods evaluated. For a subset of molecules, we utilized experimentally-determined microstates based on NMR to evaluate the dominant tautomer predictions for each macroscopic state. Prediction of dominant tautomers was a major source of error for microscopic pKa predictions, especially errors in charged tautomers. The degree of inaccuracy in pKa predictions observed in this challenge is detrimental to the protein-ligand binding affinity predictions due to errors in dominant protonation state predictions and the calculation of free energy corrections for multiple protonation states. Underestimation of ligand pKa by 1 unit can lead to errors in binding free energy errors up to 1.2 kcal/mol. The SAMPL6 pKa Challenge demonstrated the need for improving pKa prediction methods for drug-like molecules, especially for challenging moieties and multiprotic molecules.

In vitro pH dependent passive transport of ketoprofen and metformin

The kinetics of passive transport of ketoprofen and metformin, as model substances for high and low permeability, respectively, across the artificial membrane under the influence of the pH of donor solution was investigated. There was an upward trend in the apparent permeation coefficient (P _app) of ketoprofen with the decrease in pH to a value close to pKa. At the pH value below pKa the permeation coefficient had lower value, due to the higher retention of ketoprofen in the artificial membrane. Metformin is a low permeable compound, and the highest permeation values were recorded at pH 7.4. Two dissociation constants determine that metformin at physiological pH exists as a hydrophilic cationic molecule, i.e. predominantly in ionized form. At pH values below 2.8, metformin mainly exists in diprotonated form, and it was, thus, very poorly permeable. The highest retention, i.e. affinity of both ketoprofen and metformin to the membrane, was at the lowest pH values, which is explained by different mechanisms. At higher pH values of donor compartment the substances showed significantly less affinity to the membrane. The obtained values of apparent permeation coefficients at studied pH values showed good correlation with the obtained experimental values by other in vitro methods.

Binding interactions of halo-benzoic acids, halo-benzenesulfonic acids and halo-phenylboronic acids with human transthyretin

The anionic form-dependent binding interaction of halo-phenolic substances with human transthyretin (hTTR) has been observed previously. This indicates that ionizable compounds should be the primary focus in screening potential hTTR disruptors. Here, the potential binding potency of halo-benzoic acids, halo-benzenesulfonic acids/sulfates and halo-phenylboronic acids with hTTR was determined and analyzed by competitive fluorescence displacement assay integrated with computational methods. The laboratorial results indicated that the three test groups of model compounds exhibited a distinct binding affinity to hTTR. All the tested halo-phenylboronic acids, some of the tested halo-benzoic acids and halo-benzenesulfonic acids/sulfates were shown to be inactive with hTTR. Other halo-benzoic acids and halo-benzenesulfonic acids/sulfates were moderate and/or weak hTTR binders. The binding affinity of halo-benzoic acids and halo-benzenesulfonic acids/sulfates with hTTR was similar. The low distribution ability of the model compounds from water to hTTR may be the reason why they exhibited the binding potency observed with hTTR. By introducing other highly hydrophobic compounds, we observed that the binding affinity between compounds and hTTR increased with increasing molecular hydrophobicity. Those results indicated that the highly hydrophobic halo-benzoic acids and halo-benzenesulfonic acids/sulfates may be high-priority hTTR disruptors. Finally, a binary classification model was constructed employing three predictive variables. The sensitivity (S_n), specificity (S_p), predictive accuracy (Q) values of the training set and validation set were >0.83, indicating that the model had good classification performance. Thus, the binary classification model developed here could be used to distinguish whether a given ionizable compound is a potential hTTR binder or not.

Analysis of the Acid/Base Profile of Natural Products from Different Sources

For many years drug discovery and other areas in chemistry have successfully relied on natural products. Recent advances in computational methods have made possible to study the chemical space of natural products from different sources. Ionizable acidic and basic functional groups heavily influence physicochemical properties and thus a molecule's absorption, distribution, metabolism, excretion, and toxicity characteristics as well as their affinity for biological targets. This work reports the generation and critical comparison of the acid/base profiles of ten chemical databases including seven natural products sets from different origins, a set of semisynthetic compounds, a collection of approved drugs, and a compendium of food chemicals. Similarities were found in the proportion of the main charge state categories among the natural products databases with few differences in their pK_a distributions. Clear differences were observed between natural products and the approved drugs and semi-synthetic natural products databases, whereas natural products share some trends with the food chemical database. We noted that the natural products collections comprise around 45 % of neutral compounds. The proportion of single acids was approximately twice that found for FDA drugs, and they demonstrated a similar distribution of pK_a values. In contrast to drugs, only 5 % of compounds among the natural products sets had a single basic group. Likewise, simple ampholytes were less prevalent in the natural products databases relative to drugs.

Dynamic Protonation Dramatically Affects the Membrane Permeability of Drug-like Molecules

Permeability (P_m) across biological membranes is of fundamental importance and a key factor in drug absorption, distribution, and development. Although the majority of drugs will be charged at some point during oral delivery, our understanding of membrane permeation by charged species is limited. The canonical model assumes that only neutral molecules partition into and passively permeate across membranes, but there is mounting evidence that these processes are also facile for certain charged species. However, it is unknown whether such ionizable permeants dynamically neutralize at the membrane surface or permeate in their charged form. To probe protonation-coupled permeation in atomic detail, we herein apply continuous constant-pH molecular dynamics along with free energy sampling to study the permeation of a weak base propranolol (PPL), and evaluate the impact of including dynamic protonation on P_m. The simulations reveal that PPL dynamically neutralizes at the lipid-tail interface, which dramatically influences the permeation free energy landscape and explains why the conventional model overestimates the assigned intrinsic permeability. We demonstrate how fixed-charge-state simulations can account for this effect, and propose a revised model that better describes pH-coupled partitioning and permeation. Our results demonstrate how dynamic changes in protonation state may play a critical role in the permeation of ionizable molecules, including pharmaceuticals and drug-like molecules, thus requiring a revision of the standard picture.

Molecular Dynamics Simulations of Membrane Permeability

This Review illustrates the evaluation of permeability of lipid membranes from molecular dynamics (MD) simulation primarily using water and oxygen as examples. Membrane entrance, translocation, and exit of these simple permeants (one hydrophilic and one hydrophobic) can be simulated by conventional MD, and permeabilities can be evaluated directly by Fick's First Law, transition rates, and a global Bayesian analysis of the inhomogeneous solubility-diffusion model. The assorted results, many of which are applicable to simulations of nonbiological membranes, highlight the limitations of the homogeneous solubility diffusion model; support the utility of inhomogeneous solubility diffusion and compartmental models; underscore the need for comparison with experiment for both simple solvent systems (such as water/hexadecane) and well-characterized membranes; and demonstrate the need for microsecond simulations for even simple permeants like water and oxygen. Undulations, subdiffusion, fractional viscosity dependence, periodic boundary conditions, and recent developments in the field are also discussed. Last, while enhanced sampling methods and increasingly sophisticated treatments of diffusion add substantially to the repertoire of simulation-based approaches, they do not address directly the critical need for force fields with polarizability and multipoles, and constant pH methods.

The influence and manipulation of acid/base properties in drug discovery

There is a growing awareness of the importance of acid/base properties in medicinal chemistry research. In many drug classes, ionisable groups are present that make critical interactions with the receptor and are essential for potency. Yet the presence of these groups may cause problems with oral bioavailability, pharmacokinetics, or toxicity. Manipulating pK_a values during drug development or applying pro-drug techniques are strategies that can overcome potential deficits in a variety of these areas. Knowledge of drug ionisation states coupled with a consideration of pH-specific cellular environments can be used advantageously to target chemoresistance. As modern drug research ventures into drug candidates that exceed the rule of 5, such exploration requires an understanding of drug acid/base properties and how these factors affect ADMET characteristics.

А theoretical study on ionization of sartans in aqueous media and on interactions with surfactant micelles

The ionization order of sartans in aqueous media and possible way of interactions between their equilibrium forms and surfactant micelles have been theoretically investigated. The examined sartans are ampholytes (irbesartan and losartan) and a diacid (valsartan) with the close values of ionization constants. In order to get a better insight in the overlapped protolytic equilibria of sartans and to predict an affinity of the equilibrium forms interacting with micelles as biomembrane mimetic systems, the theoretical study was performed. Energy calculation of the optimized structures of the equilibrium forms was performed at the B3LYP/6-31G (d,p) level of the Density Functional Theory (DFT). The results of the theoretical study helped to assign the experimentally determined pK_a values to the corresponding ionizable centers and confirmed that in all examined compounds, the higher pK_a values can be attributed to ionization of tetrazole. The molecular descriptor values showed that sartans interact predominantly with the micelle surfaces. The equilibrium forms of ampholytes demonstrate higher affinity to the micelles, as compared to the forms of the diprotic acid. Additionally, it was shown that the uncharged molecular forms of ampholytes are more lipophylic then their zwitterionic forms.

Synthesis of unnatural α-amino acid derivatives via selective o-C–H functionalization

Selective o-C–H functionalization of aryl based amino acids including arylation, alkylation, alkynylation, halogenation, alkoxylation, and acyloxylation were developed.

Proton dissociation properties of arylphosphonates: Determination of accurate Hammett equation parameters

Determination of the proton dissociation constants of several arylphosphonic acid derivatives was carried out to investigate the accuracy of the Hammett equations available for this family of compounds. For the measurement of the pK_a values modern, accurate methods, such as the differential potentiometric titration and NMR-pH titration were used. We found our results significantly different from the pK_a values reported before (pK_a1: MAE = 0.16 pK_a2: MAE=0.59). Based on our recently measured pK_a values, refined Hammett equations were determined that might be used for predicting highly accurate ionization constants of newly synthesized compounds (pK_a1=1.70-0.894σ, pK_a2=6.92-0.934σ).

The acid/base profile of agrochemicals

Drug-likeness has long been studied in the pursuit of finding new medicines. Similarly, in the agrochemical field there is value in exploring the properties of the chemicals involved. Patterns that emerge can potentially influence future discovery campaigns to improve the probability of commercial success. In this study we investigate the acid/base properties and physicochemical characteristics of three classes of compounds: fungicides, herbicides and insecticides. In comparison with FDA-approved drugs, it was noted that the pesticides were generally smaller, possessed a neutral charge state and were more lipophilic. The results are discussed in the light of their intended targets.

High throughput methods to measure the propensity of compounds to form intramolecular hydrogen bonding

Implementation of IMHB considerations in drug discovery needs robust and validated descriptors to experimentally verify the propensity of compounds to exhibit IMHBs. The first part of the paper presents an overview of the most common techniques to measure the propensity of compounds to form IMHBs. Then we review and discuss recently proposed high throughput (HT) physicochemical descriptors (i.e. Δlog Poct-tol, EPSA and log k'80 PLRP-S) which provide the same information. Analysis of the available data enabled us to extract guidelines for the application of these descriptors in drug discovery programs.

Lipid bilayer permeation of aliphatic amine and carboxylic acid drugs: rates of insertion, translocation and dissociation from MD simulations

The overabundance of drugs containing aliphatic amine and carboxylic acid groups is rationalized in terms of their membrane permeability.

Monoselective o-C-H Functionalizations of Mandelic Acid and α-Phenylglycine

Pd-catalyzed C-H functionalization of mandelic acid and α-phenylglycine is reported. We have developed different protocols for the arylation, iodination, acetoxylation, and olefination of these substrates based on two different (Pd(II)/Pd(IV) and Pd(II)/Pd(0)) catalytic cycles. Four crucial features of these protocols are advantageous for practical applications. First, the α-hydroxyl and amino groups are protected with simple protecting groups such as acetates (Ac, Piv) and carbamates (Boc, Fmoc), respectively. Second, these protocols do not involve installation and removal of a directing group. Third, monoselectivity is accomplished. Fourth, no epimerization occurs at the vulnerable α-chiral centers.

Exploring experimental and computational markers of cyclic peptides: Charting islands of permeability

An increasing number of macrocyclic peptides that cross biological membranes are being reported, suggesting that it might be possible to develop peptides into orally bioavailable therapeutics; however, current understanding of what makes macrocyclic peptides cell permeable is still limited. Here, we synthesized 62 cyclic hexapeptides and characterized their permeability using in vitro assays commonly used to predict in vivo absorption rates, i.e. the Caco-2 and PAMPA assays. We correlated permeability with experimentally measured parameters of peptide conformation obtained using rapid methods based on chromatography and nuclear magnetic resonance spectroscopy. Based on these correlations, we propose a model describing the interplay between peptide permeability, lipophilicity and hydrogen bonding potential. Specifically, peptides with very high permeability have high lipophilicity and few solvent hydrogen bond interactions, whereas peptides with very low permeability have low lipophilicity or many solvent interactions. Our model is supported by molecular dynamics simulations of the cyclic peptides calculated in explicit solvent, providing a structural basis for the observed correlations. This prospective exploration into biomarkers of peptide permeability has the potential to unlock wider opportunities for development of peptides into drugs.

Target-based analysis of ionization states of bioactive compounds

Ionization states within a chemical neighborhood. Shown are an acidic and a basic analog of a neutral compound.

The CARLSBAD database: a confederated database of chemical bioactivities

Many bioactivity databases offer information regarding the biological activity of small molecules on protein targets. Information in these databases is often hard to resolve with certainty because of subsetting different data in a variety of formats; use of different bioactivity metrics; use of different identifiers for chemicals and proteins; and having to access different query interfaces, respectively. Given the multitude of data sources, interfaces and standards, it is challenging to gather relevant facts and make appropriate connections and decisions regarding chemical-protein associations. The CARLSBAD database has been developed as an integrated resource, focused on high-quality subsets from several bioactivity databases, which are aggregated and presented in a uniform manner, suitable for the study of the relationships between small molecules and targets. In contrast to data collection resources, CARLSBAD provides a single normalized activity value of a given type for each unique chemical-protein target pair. Two types of scaffold perception methods have been implemented and are available for datamining: HierS (hierarchical scaffolds) and MCES (maximum common edge subgraph). The 2012 release of CARLSBAD contains 439 985 unique chemical structures, mapped onto 1,420 889 unique bioactivities, and annotated with 277 140 HierS scaffolds and 54 135 MCES chemical patterns, respectively. Of the 890 323 unique structure-target pairs curated in CARLSBAD, 13.95% are aggregated from multiple structure-target values: 94 975 are aggregated from two bioactivities, 14 544 from three, 7 930 from four and 2214 have five bioactivities, respectively. CARLSBAD captures bioactivities and tags for 1435 unique chemical structures of active pharmaceutical ingredients (i.e. 'drugs'). CARLSBAD processing resulted in a net 17.3% data reduction for chemicals, 34.3% reduction for bioactivities, 23% reduction for HierS and 25% reduction for MCES, respectively. The CARLSBAD database supports a knowledge mining system that provides non-specialists with novel integrative ways of exploring chemical biology space to facilitate knowledge mining in drug discovery and repurposing. Database URL: http://carlsbad.health.unm.edu/carlsbad/.